From 79841a78421cf52ceea86eb5660563a46967b49d Mon Sep 17 00:00:00 2001
From: "Eric T. Johnson" <yut23@users.noreply.github.com>
Date: Mon, 8 Jul 2024 08:19:17 -0400
Subject: [PATCH] Check for a singular matrix when using a templated network
 (#2903)

This checks the result of RHS::dgefa() for templated networks.
---
 Source/sdc/sdc_newton_solve.H | 3 +--
 1 file changed, 1 insertion(+), 2 deletions(-)

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