From c0eee83fc6f2dade69924e138adebb5d6a76f9a2 Mon Sep 17 00:00:00 2001 From: ParamThakkar123 Date: Wed, 21 Aug 2024 00:26:03 +0530 Subject: [PATCH] init added to reducer --- lib/OrdinaryDiffEqRosenbrock/src/rosenbrock_perform_step.jl | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/lib/OrdinaryDiffEqRosenbrock/src/rosenbrock_perform_step.jl b/lib/OrdinaryDiffEqRosenbrock/src/rosenbrock_perform_step.jl index a869a02fa3..af173ef9e8 100644 --- a/lib/OrdinaryDiffEqRosenbrock/src/rosenbrock_perform_step.jl +++ b/lib/OrdinaryDiffEqRosenbrock/src/rosenbrock_perform_step.jl @@ -1239,16 +1239,16 @@ end integrator.stats.nf += 1 if mass_matrix === I - linsolve_tmp = du + dtd[i] * dT + sum(dtC[i][j] * k[j] for j in 1:i-1) + linsolve_tmp = du + dtd[i] * dT + sum(dtC[i][j] * k[j] for j in 1:i-1; init=0) else - linsolve_tmp = du + dtd[i] * dT + mass_matrix * sum(dtC[i][j] * k[j] for j in 1:i-1) + linsolve_tmp = du + dtd[i] * dT + mass_matrix * sum(dtC[i][j] * k[j] for j in 1:i-1; init=0) end k[i] = _reshape(W \ -_vec(linsolve_tmp), axes(uprev)) integrator.stats.nsolve += 1 if i < 6 - u = uprev + sum(a[i+1][j] * k[j] for j in 1:i) + u = uprev + sum(a[i+1][j] * k[j] for j in 1:i; init=0) end end