fix the loop order over Jac in BackwardEuler (#1640)

the MathArray is Fortran-ordered
AMReX-Astro · Aug 5, 2024 · 115aec8 · 115aec8
1 parent 14b8b0e
commit 115aec8
Showing 1 changed file with 2 additions and 2 deletions.
diff --git 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);