fix identity tmp sparse matrix for ida solver

OpenModelica · Jan 12, 2017 · bbc489c · bbc489c
1 parent 4d49ee5
commit bbc489c
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/SimulationRuntime/c/simulation/solver/ida_solver.c b/SimulationRuntime/c/simulation/solver/ida_solver.c
@@ -423,7 +423,7 @@ ida_solver_initial(DATA* data, threadData_t *threadData, SOLVER_INFO* solverInfo
       idaData->NNZ = data->simulationInfo->analyticJacobians[data->callback->INDEX_JAC_A].sparsePattern.numberOfNoneZeros;
       flag = IDAKLU(idaData->ida_mem, idaData->N, idaData->NNZ);
       /* to add a cj identety matrix */
-      idaData->tmpJac = NewSparseMat(idaData->N, idaData->N, idaData->NNZ);
+      idaData->tmpJac = NewSparseMat(idaData->N, idaData->N, idaData->N);
     }
     if (checkIDAflag(flag)){
       throwStreamPrint(threadData, "##IDA## Setup of linear solver KLU failed!");