@@ -291,11 +291,11 @@ end
291
291
function KrylovKitJL (args... ;
292
292
KrylovAlg = KrylovKit. GMRES, gmres_restart = 0 ,
293
293
kwargs... )
294
- return KrylovJL (KrylovAlg, gmres_restart, args, kwargs)
294
+ return KrylovKitJL (KrylovAlg, gmres_restart, args, kwargs)
295
295
end
296
296
297
297
function KrylovKitJL_CG (args... ; kwargs... )
298
- KrylovKitJL (args... ; KrylovAlg = KrylovKit. CG, kwargs... )
298
+ KrylovKitJL (args... ; KrylovAlg = KrylovKit. CG, kwargs... , isposdef = true )
299
299
end
300
300
function KrylovKitJL_GMRES (args... ; kwargs... )
301
301
KrylovKitJL (args... ; KrylovAlg = KrylovKit. GMRES, kwargs... )
@@ -306,12 +306,12 @@ function SciMLBase.solve(cache::LinearCache, alg::KrylovKitJL, kwargs...)
306
306
rtol = float (cache. reltol)
307
307
maxiter = cache. maxiters
308
308
verbosity = cache. verbose ? 1 : 0
309
- krylovdim = (alg. gmres_restart == 0 ) ? min (20 , size (A, 1 )) : alg. gmres_restart
309
+ krylovdim = (alg. gmres_restart == 0 ) ? min (20 , size (cache . A, 1 )) : alg. gmres_restart
310
310
311
311
kwargs = (atol = atol, rtol = rtol, maxiter = maxiter, verbosity = verbosity,
312
312
krylovdim = krylovdim, alg. kwargs... )
313
313
314
- x, info = KrylovKit. linsolve (cache. A, cache. b, cache. u, alg . KrylovAlg )
314
+ x, info = KrylovKit. linsolve (cache. A, cache. b, cache. u; kwargs ... )
315
315
316
316
copy! (cache. u, x)
317
317
resid = info. normres
0 commit comments