diff --git a/Source/sdc/sdc_newton_solve.H b/Source/sdc/sdc_newton_solve.H
index 7a796e270c..da0c6feaad 100644
--- 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)