From 4c6aae536cb187668f67a98c2650f4f9c1b39149 Mon Sep 17 00:00:00 2001
From: Michael Zingale <michael.zingale@stonybrook.edu>
Date: Thu, 11 Jul 2024 19:44:09 -0400
Subject: [PATCH] fix linpack pivot templating for true SDC

---
 Source/sdc/sdc_newton_solve.H | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

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)