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

Insert correct conversion sections for KLU matrices

This commit is contained in:
Alessio Cacciatori 2024-09-27 10:51:44 +02:00 committed by dwarning
parent b2802e5eb9
commit ac4b7efbe3
1 changed files with 48 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

62
src/spicelib/analysis/span.c
View File

@ -662,6 +662,20 @@ SPan(CKTcircuit* ckt, int restart)
return(error); return(error);
} }
} }
#ifdef KLU
if (ckt->CKTmatrix->CKTkluMODE)
{
/* Conversion from Real Matrix to Complex Matrix */
if (!ckt->CKTmatrix->SMPkluMatrix->KLUmatrixIsComplex)
{
for (i = 0; i < DEVmaxnum; i++)
if (DEVices[i] && DEVices[i]->DEVbindCSCComplex && ckt->CKThead[i])
DEVices[i]->DEVbindCSCComplex(ckt->CKThead[i], ckt);
ckt->CKTmatrix->SMPkluMatrix->KLUmatrixIsComplex = KLUMatrixComplex;
}
}
#endif
ckt->CKTactivePort = 0; ckt->CKTactivePort = 0;
/* main loop through all scheduled frequencies */ /* main loop through all scheduled frequencies */
@ -678,6 +692,19 @@ SPan(CKTcircuit* ckt, int restart)
/* Update opertating point, if variable 'hertz' is given */ /* Update opertating point, if variable 'hertz' is given */
if (ckt->CKTvarHertz) { if (ckt->CKTvarHertz) {
#ifdef KLU
if (ckt->CKTmatrix->CKTkluMODE)
{
/* Conversion from Complex Matrix to Real Matrix */
for (i = 0; i < DEVmaxnum; i++)
if (DEVices[i] && DEVices[i]->DEVbindCSCComplexToReal && ckt->CKThead[i])
DEVices[i]->DEVbindCSCComplexToReal(ckt->CKThead[i], ckt);
ckt->CKTmatrix->SMPkluMatrix->KLUmatrixIsComplex = KLUmatrixReal;
}
#endif
#ifdef XSPICE #ifdef XSPICE
/* Call EVTop if event-driven instances exist */ /* Call EVTop if event-driven instances exist */
@ -709,7 +736,17 @@ SPan(CKTcircuit* ckt, int restart)
tfree(data); return(error); tfree(data); return(error);
} }
} }
#ifdef KLU
if (ckt->CKTmatrix->CKTkluMODE)
{
/* Conversion from Real Matrix to Complex Matrix */
for (i = 0; i < DEVmaxnum; i++)
if (DEVices[i] && DEVices[i]->DEVbindCSCComplex && ckt->CKThead[i])
DEVices[i]->DEVbindCSCComplex(ckt->CKThead[i], ckt);
ckt->CKTmatrix->SMPkluMatrix->KLUmatrixIsComplex = KLUMatrixComplex;
}
#endif
// Store previous rhs // Store previous rhs
if (rhswoPorts == NULL) if (rhswoPorts == NULL)
rhswoPorts = (double*)TMALLOC(double, ckt->CKTmaxEqNum); rhswoPorts = (double*)TMALLOC(double, ckt->CKTmaxEqNum);
@ -760,20 +797,6 @@ SPan(CKTcircuit* ckt, int restart)
else else
vsrcRoot = ckt->CKTVSRCid; vsrcRoot = ckt->CKTVSRCid;
#ifdef KLU
if (ckt->CKTmatrix->CKTkluMODE)
{
/* Conversion from Real Matrix to Complex Matrix */
if (!ckt->CKTmatrix->SMPkluMatrix->KLUmatrixIsComplex)
{
for (i = 0 ; i < DEVmaxnum ; i++)
if (DEVices [i] && DEVices [i]->DEVbindCSCComplex && ckt->CKThead [i])
DEVices [i]->DEVbindCSCComplex (ckt->CKThead [i], ckt) ;
ckt->CKTmatrix->SMPkluMatrix->KLUmatrixIsComplex = KLUMatrixComplex ;
}
}
#endif
// Pre-load everything but RF Ports (these will be updated in the next cycle). // Pre-load everything but RF Ports (these will be updated in the next cycle).
error = NIspPreload(ckt); error = NIspPreload(ckt);
@ -977,6 +1000,17 @@ endsweep:
tfree(irhswoPorts); tfree(irhswoPorts);
deleteSPmatrix(ckt); deleteSPmatrix(ckt);
tfree(data); tfree(data);
#ifdef KLU
if (ckt->CKTmatrix->CKTkluMODE)
{
/* Conversion from Complex Matrix to Real Matrix */
for (i = 0; i < DEVmaxnum; i++)
if (DEVices[i] && DEVices[i]->DEVbindCSCComplexToReal && ckt->CKThead[i])
DEVices[i]->DEVbindCSCComplexToReal(ckt->CKThead[i], ckt);
ckt->CKTmatrix->SMPkluMatrix->KLUmatrixIsComplex = KLUmatrixReal;
}
#endif
return(0); return(0);
} }