@@ -634,13 +634,9 @@ function alg_cache(alg::Rodas3P, u, rate_prototype, ::Type{uEltypeNoUnits},
634
634
tf = TimeGradientWrapper (f, uprev, p)
635
635
uf = UJacobianWrapper (f, t, p)
636
636
linsolve_tmp = zero (rate_prototype)
637
- linprob = LinearProblem (W, _vec (linsolve_tmp); u0 = _vec (tmp))
638
- Pl, Pr = wrapprecs (
639
- alg. precs (W, nothing , u, p, t, nothing , nothing , nothing ,
640
- nothing )... , weight, tmp)
641
- linsolve = init (linprob, alg. linsolve, alias_A = true , alias_b = true ,
642
- Pl = Pl, Pr = Pr,
643
- assumptions = LinearSolve. OperatorAssumptions (true ))
637
+ linprob = LinearProblem (W, _vec (linsolve_tmp), (nothing ,u,p,t); u0 = _vec (tmp))
638
+ linsolve = init (linprob, wrapprecs (alg. linsolve, W, weight), alias_A = true , alias_b = true ,
639
+ assumptions = LinearSolve. OperatorAssumptions (true ))
644
640
grad_config = build_grad_config (alg, f, tf, du1, t)
645
641
jac_config = build_jac_config (alg, f, uf, du1, uprev, u, tmp, du2)
646
642
Rodas3PCache (u, uprev, dense1, dense2, dense3, du, du1, du2, k1, k2, k3, k4, k5,
0 commit comments