luke/ngspice
1
0
Fork 0
mirror of https://git.code.sf.net/p/ngspice/ngspice synced 2026-08-22 13:57:32 +02:00
Code Issues Packages Projects Releases Wiki Activity

Maintain the KLU related changes under '#ifdef KLU'

Browse Source
This commit is contained in:
Francesco Lannutti
2020-10-28 23:49:48 +01:00
parent d2df13a63f
commit f623e62ce4
1 changed files with 16 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files
src/spicelib/analysis/cktic.c
+16
View File
@@ -25,6 +25,8 @@ CKTic(CKTcircuit *ckt)
for(node = ckt->CKTnodes;node != NULL; node = node->next) { for(node = ckt->CKTnodes;node != NULL; node = node->next) {
if(node->nsGiven) { if(node->nsGiven) {
#ifdef KLU
if (ckt->CKTkluMODE) { if (ckt->CKTkluMODE) {
node->ptr = (double *) SMPfindElt (ckt->CKTmatrix, node->number, node->number, 0) ; node->ptr = (double *) SMPfindElt (ckt->CKTmatrix, node->number, node->number, 0) ;
if (node->ptr == NULL) { if (node->ptr == NULL) {
@@ -32,14 +34,22 @@ CKTic(CKTcircuit *ckt)
printf ("Please specify an existing element for .nodeset\n") ; printf ("Please specify an existing element for .nodeset\n") ;
} }
} else { } else {
#endif
node->ptr = SMPmakeElt(ckt->CKTmatrix,node->number,node->number); node->ptr = SMPmakeElt(ckt->CKTmatrix,node->number,node->number);
#ifdef KLU
} }
#endif
if(node->ptr == NULL) return(E_NOMEM); if(node->ptr == NULL) return(E_NOMEM);
ckt->CKThadNodeset = 1; ckt->CKThadNodeset = 1;
ckt->CKTrhsOld[node->number] = ckt->CKTrhs[node->number] = node->nodeset; ckt->CKTrhsOld[node->number] = ckt->CKTrhs[node->number] = node->nodeset;
} }
if(node->icGiven) { if(node->icGiven) {
if(! ( node->ptr)) { if(! ( node->ptr)) {
#ifdef KLU
if (ckt->CKTkluMODE) { if (ckt->CKTkluMODE) {
node->ptr = (double *) SMPfindElt (ckt->CKTmatrix, node->number, node->number, 0) ; node->ptr = (double *) SMPfindElt (ckt->CKTmatrix, node->number, node->number, 0) ;
if (node->ptr == NULL) { if (node->ptr == NULL) {
@@ -47,8 +57,14 @@ CKTic(CKTcircuit *ckt)
printf ("Please specify an existing element for .ic\n") ; printf ("Please specify an existing element for .ic\n") ;
} }
} else { } else {
#endif
node->ptr = SMPmakeElt(ckt->CKTmatrix,node->number, node->number); node->ptr = SMPmakeElt(ckt->CKTmatrix,node->number, node->number);
#ifdef KLU
} }
#endif
if(node->ptr == NULL) return(E_NOMEM); if(node->ptr == NULL) return(E_NOMEM);
} }
ckt->CKTrhsOld[node->number] = ckt->CKTrhs[node->number] = node->ic; ckt->CKTrhsOld[node->number] = ckt->CKTrhs[node->number] = node->ic;