@@ -634,13 +634,9 @@ function alg_cache(alg::Rodas3P, u, rate_prototype, ::Type{uEltypeNoUnits},
634634 tf = TimeGradientWrapper (f, uprev, p)
635635 uf = UJacobianWrapper (f, t, p)
636636 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 ))
644640 grad_config = build_grad_config (alg, f, tf, du1, t)
645641 jac_config = build_jac_config (alg, f, uf, du1, uprev, u, tmp, du2)
646642 Rodas3PCache (u, uprev, dense1, dense2, dense3, du, du1, du2, k1, k2, k3, k4, k5,
0 commit comments