Check for a singular matrix when using a templated network (#2903)

This checks the result of RHS::dgefa() for templated networks.
AMReX-Astro · Jul 8, 2024 · 79841a7 · 79841a7
1 parent 0538b05
commit 79841a7
Showing 1 changed file with 1 addition and 2 deletions.
diff --git a/Source/sdc/sdc_newton_solve.H b/Source/sdc/sdc_newton_solve.H
@@ -159,8 +159,7 @@ sdc_newton_solve(const Real dt_m,
 
         // solve the linear system: Jac dU = -f
 #ifdef NEW_NETWORK_IMPLEMENTATION
-        RHS::dgefa(Jac);
-        info = 0;
+        info = RHS::dgefa(Jac);
 #else
         IArray1D ipvt;
         dgefa<NumSpec+1>(Jac, ipvt, info);