Update to Arpack v0.5 (#286)

* Update to Arpack v0.5 * Update to arpack v0.5.1
qojulia · Jan 4, 2021 · f3f0c65 · f3f0c65
1 parent 50e29cb
commit f3f0c65
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 3 deletions.
diff --git a/Project.toml b/Project.toml
@@ -17,7 +17,7 @@ StochasticDiffEq = "789caeaf-c7a9-5a7d-9973-96adeb23e2a0"
 WignerSymbols = "9f57e263-0b3d-5e2e-b1be-24f2bb48858b"
 
 [compat]
-Arpack = "0.4"
+Arpack = "0.5.1"
 DiffEqCallbacks = "2"
 FFTW = "1"
 OrdinaryDiffEq = "5"

diff --git a/src/steadystate.jl b/src/steadystate.jl
@@ -68,7 +68,7 @@ function liouvillianspectrum(L::SparseSuperOpType; nev::Int = min(10, length(L.b
     d, v, nconv, niter, nmult, resid = try
         eigs(L.data; nev = nev, which = which, kwargs...)
     catch err
-        if isa(err, SingularException) || isa(err, Arpack.ARPACKException)
+        if isa(err, SingularException) || isa(err, Arpack.XYEUPD_Exception) || isa(err, Arpack.XYAUPD_Exception)
             error("Arpack's eigs() algorithm failed; try using DenseOperators or change nev.")
         else
             rethrow(err)

diff --git a/test/test_steadystate.jl b/test/test_steadystate.jl
@@ -52,7 +52,7 @@ tss, ρss = steadystate.master(Hdense, Jdense; tol=1e-4)
 ρss = steadystate.eigenvector(Hdense, Jdense, nev = 1)
 @test tracedistance(ρss, ρt[end]) < 1e-6
 
-ρss = steadystate.eigenvector(H, sqrt(2).*J; rates=0.5.*ones(length(J)))
+ρss = steadystate.eigenvector(H, sqrt(2).*J; rates=0.5.*ones(length(J)), tol=1e-8)
 @test tracedistance(ρss, ρt[end]) < 1e-3
 
 ρss = steadystate.eigenvector(H, sqrt(2).*J; rates=0.5.*ones(length(J)), nev = 1)