luke
/
ngspice
mirror of https://git.code.sf.net/p/ngspice/ngspice
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Fix SMPfindElement for KLU

This commit is contained in:
Francesco Lannutti 2020-04-05 19:56:28 +02:00 committed by Holger Vogt
parent 11fb209ee6
commit fdd815183a
1 changed files with 8 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
src/maths/KLU/klusmp.c
View File

@ -1462,8 +1462,14 @@ SMPfindElt (SMPmatrix *eMatrix, int Row, int Col, int CreateIfMissing)
Row = Matrix->ExtToIntRowMap [Row] ; Row = Matrix->ExtToIntRowMap [Row] ;
Col = Matrix->ExtToIntColMap [Col] ; Col = Matrix->ExtToIntColMap [Col] ;
for (i = eMatrix->CKTkluAp [Col - 1] ; i < eMatrix->CKTkluAp [Col] ; i++) { Row = Row - 1 ;
if (eMatrix->CKTkluAi [i] == Row - 1) { Col = Col - 1 ;
if ((Row < 0) || (Col < 0)) {
printf ("Error: Cannot find an element with row '%d' and column '%d' in the KLU matrix\n", Row, Col) ;
return NULL ;
}
for (i = eMatrix->CKTkluAp [Col] ; i < eMatrix->CKTkluAp [Col + 1] ; i++) {
if (eMatrix->CKTkluAi [i] == Row) {
if (eMatrix->CKTkluMatrixIsComplex == CKTkluMatrixReal) { if (eMatrix->CKTkluMatrixIsComplex == CKTkluMatrixReal) {
return (SMPelement *) &(eMatrix->CKTkluAx [i]) ; return (SMPelement *) &(eMatrix->CKTkluAx [i]) ;
} else if (eMatrix->CKTkluMatrixIsComplex == CKTkluMatrixComplex) { } else if (eMatrix->CKTkluMatrixIsComplex == CKTkluMatrixComplex) {