From 115aec8846801361d790190c74e1c8c2b5addcaa Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Mon, 5 Aug 2024 18:19:27 -0400 Subject: [PATCH] fix the loop order over Jac in BackwardEuler (#1640) the MathArray is Fortran-ordered --- integration/BackwardEuler/be_integrator.H | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/integration/BackwardEuler/be_integrator.H b/integration/BackwardEuler/be_integrator.H index 1a8e5c954..f82cab471 100644 --- a/integration/BackwardEuler/be_integrator.H +++ b/integration/BackwardEuler/be_integrator.H @@ -76,8 +76,8 @@ int single_step (BurnT& state, BeT& be, const amrex::Real dt) // construct the matrix for the linear system // (I - dt J) dy^{n+1} = rhs - for (int m = 1; m <= int_neqs; m++) { - for (int n = 1; n <= int_neqs; n++) { + for (int n = 1; n <= int_neqs; n++) { + for (int m = 1; m <= int_neqs; m++) { be.jac(m, n) *= -dt; if (m == n) { be.jac(m, n) = 1.0_rt + be.jac(m, n);