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

Fixed the SMPcSolve routine for KLU

This commit is contained in:
Francesco Lannutti 2016-06-12 21:47:31 +02:00 committed by Holger Vogt
parent df9439ebe4
commit f8ad3e6f70
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -339,8 +339,8 @@ SMPcSolve (SMPmatrix *Matrix, double RHS[], double iRHS[], double Spare[], doubl
pExtOrder = &Matrix->SPmatrix->IntToExtRowMap [Matrix->CKTkluN] ; pExtOrder = &Matrix->SPmatrix->IntToExtRowMap [Matrix->CKTkluN] ;
for (i = 2 * Matrix->CKTkluN - 1 ; i > 0 ; i -= 2) for (i = 2 * Matrix->CKTkluN - 1 ; i > 0 ; i -= 2)
{ {
Matrix->CKTkluIntermediate_Complex [i] = RHS [*(pExtOrder)] ; Matrix->CKTkluIntermediate_Complex [i] = iRHS [*(pExtOrder)] ;
Matrix->CKTkluIntermediate_Complex [i - 1] = iRHS [*(pExtOrder--)] ; Matrix->CKTkluIntermediate_Complex [i - 1] = RHS [*(pExtOrder--)] ;
} }
ret = klu_z_solve (Matrix->CKTkluSymbolic, Matrix->CKTkluNumeric, Matrix->CKTkluN, 1, Matrix->CKTkluIntermediate_Complex, Matrix->CKTkluCommon) ; ret = klu_z_solve (Matrix->CKTkluSymbolic, Matrix->CKTkluNumeric, Matrix->CKTkluN, 1, Matrix->CKTkluIntermediate_Complex, Matrix->CKTkluCommon) ;
@ -348,8 +348,8 @@ SMPcSolve (SMPmatrix *Matrix, double RHS[], double iRHS[], double Spare[], doubl
pExtOrder = &Matrix->SPmatrix->IntToExtColMap [Matrix->CKTkluN] ; pExtOrder = &Matrix->SPmatrix->IntToExtColMap [Matrix->CKTkluN] ;
for (i = 2 * Matrix->CKTkluN - 1 ; i > 0 ; i -= 2) for (i = 2 * Matrix->CKTkluN - 1 ; i > 0 ; i -= 2)
{ {
RHS [*(pExtOrder)] = Matrix->CKTkluIntermediate_Complex [i] ; iRHS [*(pExtOrder)] = Matrix->CKTkluIntermediate_Complex [i] ;
iRHS [*(pExtOrder--)] = Matrix->CKTkluIntermediate_Complex [i - 1] ; RHS [*(pExtOrder--)] = Matrix->CKTkluIntermediate_Complex [i - 1] ;
} }
} else { } else {