Use O(N^2) matrix-vector multiply in out-of-place IRKC

SciML · Mar 18, 2019 · ae63437 · ae63437
1 parent 5d62d61
commit ae63437
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/perform_step/rkc_perform_step.jl b/src/perform_step/rkc_perform_step.jl
@@ -602,7 +602,7 @@ function perform_step!(integrator,cache::IRKCConstantCache,repeat_step=false)
   if isnewton(nlsolver) && integrator.opts.adaptive
     update_W!(integrator, cache, dt, false)
     LU_W = get_W(nlsolver)
-    tmp = LU_W.L * LU_W.U * (@. dt*(0.5*(cache.du₂ - du₂) + (0.5 - μs₁)*(cache.du₁ - du₁)))
+    tmp = LU_W.L * (LU_W.U * (@. dt*(0.5*(cache.du₂ - du₂) + (0.5 - μs₁)*(cache.du₁ - du₁))))
     atmp = calculate_residuals(tmp, uprev, u, integrator.opts.abstol, integrator.opts.reltol, integrator.opts.internalnorm, t)
     integrator.EEst = integrator.opts.internalnorm(atmp,t)
   end