Fix segmentation fault when using KINSOL with LAPACK as linear solver.

Allocate memory for dense Jacobian when LAPACK is selected.
OpenModelica · Sep 3, 2021 · 0f49ff7 · 0f49ff7
1 parent 5eb2c51
commit 0f49ff7
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/OMCompiler/SimulationRuntime/c/simulation/solver/kinsolSolver.c b/OMCompiler/SimulationRuntime/c/simulation/solver/kinsolSolver.c
@@ -165,7 +165,7 @@ static void resetKinsolMemory(NLS_KINSOL_DATA *kinsolData,
   checkReturnFlag_SUNDIALS(flag, SUNDIALS_KIN_FLAG, "KINInit");
 
   /* Create matrix object */
-  if (kinsolData->linearSolverMethod == NLS_LS_DEFAULT) {
+  if (kinsolData->linearSolverMethod == NLS_LS_DEFAULT || kinsolData->linearSolverMethod == NLS_LS_LAPACK) {
     kinsolData->J = SUNDenseMatrix(size, size);
   } else if (kinsolData->linearSolverMethod == NLS_LS_KLU) {
     kinsolData->nnz = nlsData->sparsePattern->numberOfNoneZeros;