Skip to content
This repository was archived by the owner on May 18, 2019. It is now read-only.

Commit eae1de1

Browse files
Willi BraunOpenModelica-Hudson
authored andcommitted
use sparsePattern.nnz instead of the extended
1 parent 860bb91 commit eae1de1

File tree

1 file changed

+3
-3
lines changed

1 file changed

+3
-3
lines changed

SimulationRuntime/c/simulation/solver/ida_solver.c

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -1346,7 +1346,7 @@ static void setJacElementKluSparse(int row, int col, double value, int nth, SlsM
13461346
}
13471347

13481348
/* finish sparse matrix, by fixing colprts */
1349-
static void finishSparseColPtr(SlsMat mat)
1349+
static void finishSparseColPtr(SlsMat mat, int nnz)
13501350
{
13511351
int i;
13521352
for(i=1; i<mat->N+1; ++i){
@@ -1355,7 +1355,7 @@ static void finishSparseColPtr(SlsMat mat)
13551355
}
13561356
}
13571357
/* finish matrix colptrs */
1358-
mat->colptrs[mat->N] = mat->NNZ;
1358+
mat->colptrs[mat->N] = nnz;
13591359
}
13601360

13611361
/*
@@ -1458,7 +1458,7 @@ int jacobianSparseNumIDA(double tt, N_Vector yy, N_Vector yp, N_Vector rr, SlsMa
14581458
}
14591459
}
14601460
}
1461-
finishSparseColPtr(Jac);
1461+
finishSparseColPtr(Jac, sparsePattern->numberOfNoneZeros);
14621462

14631463

14641464
unsetContext(data);

0 commit comments

Comments
 (0)