Change default solver detection method

Author: albertomercurio

src/qobj/eigsolve.jl

@@ -378,9 +378,16 @@ function eigsolve(
         !haskey(kwargs2, :assumptions) && (kwargs2 = merge(kwargs2, (assumptions = OperatorAssumptions(true),)))
 
         prob = LinearProblem{true}(Aₛ, v0)
-        linsolve =
-            typeof(A) <: SparseMatrixCSC ? init(prob, UMFPACKFactorization(); kwargs2...) :
+
+        linsolve = if isnothing(solver)
+            if typeof(A) <: SparseMatrixCSC && isprimitivetype(T)
+                init(prob, UMFPACKFactorization(); kwargs2...)
+            else
+                init(prob; kwargs2...)
+            end
+        else
             init(prob, solver; kwargs2...)
+        end
 
         Amap = EigsolveInverseMap(T, size(A), linsolve)
 

test/ext-test/cpu/arbitrary_precision/Project.toml

@@ -1,3 +1,4 @@
 [deps]
 GenericSchur = "c145ed77-6b09-5dd9-b285-bf645a82121e"
 QuantumToolbox = "6c2fb7c5-b903-41d2-bc5e-5a7c320b9fab"
+Sparspak = "e56a9233-b9d6-4f03-8d0f-1825330902ac"

test/runtests.jl

@@ -89,7 +89,9 @@ if (GROUP == "Arbitrary-Precision")
     Pkg.update()
 
     using QuantumToolbox
+    using LinearAlgebra
     using SparseArrays
+    using Sparspak
     using GenericSchur
 
     QuantumToolbox.about()