fix

Jutho · Jun 7, 2024 · 45e7d7c · 45e7d7c
1 parent 18a042a
commit 45e7d7c
Show file tree

Hide file tree

Showing 2 changed files with 12 additions and 3 deletions.
diff --git a/src/eigsolve/eigsolve.jl b/src/eigsolve/eigsolve.jl
@@ -206,7 +206,16 @@ function eigsolve(f, x₀, howmany::Int=1, which::Selector=:LM; kwargs...)
             `λ` and `conj(λ)` equally: work in complex arithmetic by providing a complex starting vector `x₀`")
         end
     end
-    return eigsolve(f, x₀, howmany, which, alg)
+    if haskey(kwargs, :alg_rrule)
+        alg_rrule = kwargs[:alg_rrule]
+    else
+        alg_rrule = Arnoldi(; tol=alg.tol,
+                            krylovdim=alg.krylovdim,
+                            maxiter=alg.maxiter,
+                            eager=alg.eager,
+                            orth=alg.orth)
+    end
+    return eigsolve(f, x₀, howmany, which, alg; alg_rrule=alg_rrule)
 end
 
 function eigselector(f,

diff --git a/src/linsolve/linsolve.jl b/src/linsolve/linsolve.jl
@@ -111,8 +111,8 @@ function linsolve(f, b, x₀, a₀::Number=0, a₁::Number=1; kwargs...)
     T = promote_type(Core.Compiler.return_type(inner, Tuple{Tb,Tfx}), typeof(a₀),
                      typeof(a₁))
     alg = linselector(f, b, T; kwargs...)
-    if haskey(kwargs, alg_rrule)
-        alg_rrule = kwargs[alg_rrule]
+    if haskey(kwargs, :alg_rrule)
+        alg_rrule = kwargs[:alg_rrule]
     else
         alg_rrule = alg
     end