Fix the error with Hermtian matrices

exanauts · Dec 8, 2023 · 206bbdd · 206bbdd
1 parent 006e180
commit 206bbdd
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/src/generic.jl b/src/generic.jl
@@ -12,6 +12,7 @@ function LinearAlgebra.ldlt(A::CuSparseMatrixCSR{T,Cint}) where T <: BlasFloat
   n = LinearAlgebra.checksquare(A)
   structure = T <: Real ? 'S' : 'H'
   solver = CudssSolver(A, structure, 'F')
+  (T <: Complex) && cudss_set(solver, "pivot_type", 'N')
   x = CudssMatrix(T, n)
   b = CudssMatrix(T, n)
   cudss("analysis", solver, x, b)