diff --git a/Source/sdc/sdc_newton_solve.H b/Source/sdc/sdc_newton_solve.H index 580b3e6c43..bdc06fc113 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; + RHS::dgefa(Jac, info); #else IArray1D ipvt; dgefa(Jac, ipvt, info);