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

Revert "enable removal of dangling nodes also when MODETRANOP is set." 

This reverts commit 4fb403a590.
This commit is contained in:
Holger Vogt 2026-07-31 14:11:47 +02:00
parent 47c77e421d
commit ee38134f63
1 changed files with 13 additions and 18 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

31
src/spicelib/devices/cap/capload.c
View File

@ -49,25 +49,20 @@ CAPload(GENmodel *inModel, CKTcircuit *ckt)
vcap = *(ckt->CKTrhsOld+here->CAPposNode) - vcap = *(ckt->CKTrhsOld+here->CAPposNode) -
*(ckt->CKTrhsOld+here->CAPnegNode) ; *(ckt->CKTrhsOld+here->CAPnegNode) ;
} }
if (here->CAPdangling) {
/* Topology reduction: this cap hangs on a
* floating (degree-1) node. Remove it from the system:
* pin the floating node(s) with a unit diagonal so the
* matrix stays nonsingular, and contribute no charge or
* current. This eliminates the spurious LTE pressure
* that otherwise drives "Timestep too small" at the
* dangling node (set in CKTtopologyReduce()). */
if (here->CAPdangling & 1)
*(here->CAPposPosPtr) += 1.0;
if (here->CAPdangling & 2)
*(here->CAPnegNegPtr) += 1.0;
*(ckt->CKTstate0+here->CAPqcap) = 0.0;
continue;
}
if(ckt->CKTmode & (MODETRAN | MODEAC)) { if(ckt->CKTmode & (MODETRAN | MODEAC)) {
if (here->CAPdangling) {
/* Topology reduction: this cap hangs on a
* floating (degree-1) node. Remove it from the system:
* pin the floating node(s) with a unit diagonal so the
* matrix stays nonsingular, and contribute no charge or
* current. This eliminates the spurious LTE pressure
* that otherwise drives "Timestep too small" at the
* dangling node (set in CKTtopologyReduce()). */
if (here->CAPdangling & 1) *(here->CAPposPosPtr) += 1.0;
if (here->CAPdangling & 2) *(here->CAPnegNegPtr) += 1.0;
*(ckt->CKTstate0+here->CAPqcap) = 0.0;
continue;
}
#ifndef PREDICTOR #ifndef PREDICTOR
if(ckt->CKTmode & MODEINITPRED) { if(ckt->CKTmode & MODEINITPRED) {
*(ckt->CKTstate0+here->CAPqcap) = *(ckt->CKTstate0+here->CAPqcap) =