fix linpack pivot templating for true SDC (#2912)

AMReX-Astro · Jul 22, 2024 · 56f51f6 · 56f51f6
1 parent ab2b6d8
commit 56f51f6
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/Source/sdc/sdc_newton_solve.H b/Source/sdc/sdc_newton_solve.H
@@ -162,7 +162,8 @@ sdc_newton_solve(const Real dt_m,
         info = RHS::dgefa(Jac);
 #else
         IArray1D ipvt;
-        dgefa<NumSpec+1>(Jac, ipvt, info);
+        constexpr bool allow_pivot{true};
+        dgefa<NumSpec+1, allow_pivot>(Jac, ipvt, info);
 #endif
         if (info != 0) {
             ierr = newton::SINGULAR_MATRIX;
@@ -172,7 +173,7 @@ sdc_newton_solve(const Real dt_m,
 #ifdef NEW_NETWORK_IMPLEMENTATION
         RHS::dgesl(Jac, f);
 #else
-        dgesl<NumSpec+1>(Jac, ipvt, f);
+        dgesl<NumSpec+1, allow_pivot>(Jac, ipvt, f);
 #endif
 
         // on output, f is the solution (dU)