From bbc489cd20a4376c757699be92bad4f7b710720e Mon Sep 17 00:00:00 2001 From: Willi Braun Date: Thu, 12 Jan 2017 11:40:49 +0100 Subject: [PATCH] fix identity tmp sparse matrix for ida solver --- SimulationRuntime/c/simulation/solver/ida_solver.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/SimulationRuntime/c/simulation/solver/ida_solver.c b/SimulationRuntime/c/simulation/solver/ida_solver.c index 0d8dccebd81..6c4b3ac6934 100644 --- 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!");