[numerics] fix bug with log file for ma57
@@ -6330,6 +6330,8 @@ int NM_LDLT_solve(NumericsMatrix* Ao, double *b, unsigned int nrhs)
case NSM_HSL:
{
LBL_Data * lbl = (LBL_Data *)p->linear_solver_data;
+ FILE * logfile = fopen("lbl.log", "a");
+ lbl->ma57->logfile = logfile;
// Solve.
for (int irhs=0; irhs <nrhs ; irhs++)
{
@@ -6339,6 +6341,7 @@ int NM_LDLT_solve(NumericsMatrix* Ao, double *b, unsigned int nrhs)
fprintf(stderr, "NM_LDLT_solve. LBL_Solve error return from Solve: %d\n", info);
}
}
+ fclose(logfile);
break;
}
#endif
