/********** Copyright 1990 Regents of the University of California. All rights reserved. Author: 1998 Samuel Fung, Dennis Sinitsky and Stephen Tang File: b3soipdacld.c 98/5/01 Modified by Pin Su 99/4/30 Modified by Pin Su 99/9/27 Modified by Paolo Nenzi 2002 **********/ /* * Revision 2.2.3 02/3/5 Pin Su * BSIMPD2.2.3 release */ #include "ngspice/ngspice.h" #include "ngspice/cktdefs.h" #include "b3soipddef.h" #include "ngspice/sperror.h" #include "ngspice/suffix.h" int B3SOIPDacLoad(GENmodel *inModel, CKTcircuit *ckt) { B3SOIPDmodel *model = (B3SOIPDmodel*)inModel; B3SOIPDinstance *here; int selfheat; double xcggb, xcgdb, xcgsb, xcgeb, xcgT; double xcdgb, xcddb, xcdsb, xcdeb, xcdT; double xcsgb, xcsdb, xcssb, xcseb, xcsT; double xcbgb, xcbdb, xcbsb, xcbeb, xcbT; double xcegb, xceeb, xceT; double gdpr, gspr, gds; double cggb, cgdb, cgsb, cgT; double cdgb, cddb, cdsb, cdeb, cdT; double cbgb, cbdb, cbsb, cbeb, cbT; double ceeb, ceT; double GSoverlapCap, GDoverlapCap, GEoverlapCap, FwdSum, RevSum, Gm, Gmbs, GmT; double omega; double dxpart, sxpart; double gbbg, gbbdp, gbbb, gbbp, gbbsp, gbbT; double gddpg, gddpdp, gddpsp, gddpb, gddpT; double gsspg, gsspdp, gsspsp, gsspb, gsspT; double gppdp, gppb, gppp, gppT; double xcTt, cTt, gTtt, gTtg, gTtb, gTtdp, gTtsp; double EDextrinsicCap, ESextrinsicCap; double xcedb, xcesb; double m; omega = ckt->CKTomega; for (; model != NULL; model = model->B3SOIPDnextModel) { for (here = model->B3SOIPDinstances; here!= NULL; here = here->B3SOIPDnextInstance) { selfheat = (model->B3SOIPDshMod == 1) && (here->B3SOIPDrth0 != 0.0); if (here->B3SOIPDmode >= 0) { Gm = here->B3SOIPDgm; Gmbs = here->B3SOIPDgmbs; GmT = model->B3SOIPDtype * here->B3SOIPDgmT; FwdSum = Gm + Gmbs; RevSum = 0.0; cbgb = here->B3SOIPDcbgb; cbsb = here->B3SOIPDcbsb; cbdb = here->B3SOIPDcbdb; cbeb = here->B3SOIPDcbeb; cbT = model->B3SOIPDtype * here->B3SOIPDcbT; ceeb = here->B3SOIPDceeb; ceT = model->B3SOIPDtype * here->B3SOIPDceT; cggb = here->B3SOIPDcggb; cgsb = here->B3SOIPDcgsb; cgdb = here->B3SOIPDcgdb; cgT = model->B3SOIPDtype * here->B3SOIPDcgT; cdgb = here->B3SOIPDcdgb; cdsb = here->B3SOIPDcdsb; cddb = here->B3SOIPDcddb; cdeb = here->B3SOIPDcdeb; cdT = model->B3SOIPDtype * here->B3SOIPDcdT; cTt = here->pParam->B3SOIPDcth; gbbg = -here->B3SOIPDgbgs; gbbdp = -here->B3SOIPDgbds; gbbb = -here->B3SOIPDgbbs; gbbp = -here->B3SOIPDgbps; gbbT = -model->B3SOIPDtype * here->B3SOIPDgbT; gbbsp = - ( gbbg + gbbdp + gbbb + gbbp); gddpg = -here->B3SOIPDgjdg; gddpdp = -here->B3SOIPDgjdd; gddpb = -here->B3SOIPDgjdb; gddpT = -model->B3SOIPDtype * here->B3SOIPDgjdT; gddpsp = - ( gddpg + gddpdp + gddpb ); gsspg = -here->B3SOIPDgjsg; gsspdp = -here->B3SOIPDgjsd; gsspb = -here->B3SOIPDgjsb; gsspT = -model->B3SOIPDtype * here->B3SOIPDgjsT; gsspsp = - (gsspg + gsspdp + gsspb ); gppdp = 0; gppb = -here->B3SOIPDgbpbs; gppp = -here->B3SOIPDgbpps; gppT = -model->B3SOIPDtype * here->B3SOIPDgbpT; gTtg = here->B3SOIPDgtempg; gTtb = here->B3SOIPDgtempb; gTtdp = here->B3SOIPDgtempd; gTtt = here->B3SOIPDgtempT; gTtsp = - (gTtg + gTtb + gTtdp); sxpart = 0.6; dxpart = 0.4; } else { Gm = -here->B3SOIPDgm; Gmbs = -here->B3SOIPDgmbs; GmT = -model->B3SOIPDtype * here->B3SOIPDgmT; FwdSum = 0.0; RevSum = -Gm - Gmbs ; cdgb = - (here->B3SOIPDcdgb + here->B3SOIPDcggb + here->B3SOIPDcbgb); cdsb = - (here->B3SOIPDcddb + here->B3SOIPDcgdb + here->B3SOIPDcbdb); cddb = - (here->B3SOIPDcdsb + here->B3SOIPDcgsb + here->B3SOIPDcbsb); cdeb = - (here->B3SOIPDcdeb + here->B3SOIPDcbeb + here->B3SOIPDceeb); cdT = - model->B3SOIPDtype * (here->B3SOIPDcgT + here->B3SOIPDcbT + here->B3SOIPDcdT + here->B3SOIPDceT); ceeb = here->B3SOIPDceeb; ceT = model->B3SOIPDtype * here->B3SOIPDceT; cggb = here->B3SOIPDcggb; cgsb = here->B3SOIPDcgdb; cgdb = here->B3SOIPDcgsb; cgT = model->B3SOIPDtype * here->B3SOIPDcgT; cbgb = here->B3SOIPDcbgb; cbsb = here->B3SOIPDcbdb; cbdb = here->B3SOIPDcbsb; cbeb = here->B3SOIPDcbeb; cbT = model->B3SOIPDtype * here->B3SOIPDcbT; cTt = here->pParam->B3SOIPDcth; gbbg = -here->B3SOIPDgbgs; gbbb = -here->B3SOIPDgbbs; gbbp = -here->B3SOIPDgbps; gbbsp = -here->B3SOIPDgbds; gbbT = -model->B3SOIPDtype * here->B3SOIPDgbT; gbbdp = - ( gbbg + gbbsp + gbbb + gbbp); gddpg = -here->B3SOIPDgjsg; gddpsp = -here->B3SOIPDgjsd; gddpb = -here->B3SOIPDgjsb; gddpT = -model->B3SOIPDtype * here->B3SOIPDgjsT; gddpdp = - (gddpg + gddpsp + gddpb ); gsspg = -here->B3SOIPDgjdg; gsspsp = -here->B3SOIPDgjdd; gsspb = -here->B3SOIPDgjdb; gsspT = -model->B3SOIPDtype * here->B3SOIPDgjdT; gsspdp = - ( gsspg + gsspsp + gsspb ); gppb = -here->B3SOIPDgbpbs; gppp = -here->B3SOIPDgbpps; gppT = -model->B3SOIPDtype * here->B3SOIPDgbpT; gppdp = - (gppb + gppp); gTtt = here->B3SOIPDgtempT; gTtg = here->B3SOIPDgtempg; gTtb = here->B3SOIPDgtempb; gTtdp = here->B3SOIPDgtempd; gTtsp = - (gTtt + gTtg + gTtb + gTtdp); gTtg = here->B3SOIPDgtempg; gTtb = here->B3SOIPDgtempb; gTtsp = here->B3SOIPDgtempd; gTtt = here->B3SOIPDgtempT; gTtdp = - (gTtg + gTtb + gTtsp); sxpart = 0.6; sxpart = 0.4; dxpart = 0.6; } gdpr=here->B3SOIPDdrainConductance; gspr=here->B3SOIPDsourceConductance; gds= here->B3SOIPDgds; GSoverlapCap = here->B3SOIPDcgso; GDoverlapCap = here->B3SOIPDcgdo; GEoverlapCap = here->pParam->B3SOIPDcgeo; EDextrinsicCap = here->B3SOIPDgcde; ESextrinsicCap = here->B3SOIPDgcse; xcedb = -EDextrinsicCap * omega; xcdeb = (cdeb - EDextrinsicCap) * omega; xcddb = (cddb + GDoverlapCap + EDextrinsicCap) * omega; xceeb = (ceeb + GEoverlapCap + EDextrinsicCap + ESextrinsicCap) * omega; xcesb = -ESextrinsicCap * omega; xcssb = (GSoverlapCap + ESextrinsicCap - (cgsb + cbsb + cdsb)) * omega; xcseb = -(cbeb + cdeb + ceeb + ESextrinsicCap) * omega; xcegb = (- GEoverlapCap) * omega; xceT = ceT * omega; xcggb = (cggb + GDoverlapCap + GSoverlapCap + GEoverlapCap) * omega; xcgdb = (cgdb - GDoverlapCap ) * omega; xcgsb = (cgsb - GSoverlapCap) * omega; xcgeb = (- GEoverlapCap) * omega; xcgT = cgT * omega; xcdgb = (cdgb - GDoverlapCap) * omega; xcdsb = cdsb * omega; xcdT = cdT * omega; xcsgb = -(cggb + cbgb + cdgb + GSoverlapCap) * omega; xcsdb = -(cgdb + cbdb + cddb) * omega; xcsT = -(cgT + cbT + cdT + ceT) * omega; xcbgb = cbgb * omega; xcbdb = cbdb * omega; xcbsb = cbsb * omega; xcbeb = cbeb * omega; xcbT = cbT * omega; xcTt = cTt * omega; m = here->B3SOIPDm; *(here->B3SOIPDEdpPtr +1) += m * xcedb; *(here->B3SOIPDEspPtr +1) += m * xcesb; *(here->B3SOIPDDPePtr +1) += m * xcdeb; *(here->B3SOIPDSPePtr +1) += m * xcseb; *(here->B3SOIPDEgPtr +1) += m * xcegb; *(here->B3SOIPDGePtr +1) += m * xcgeb; *(here->B3SOIPDEePtr +1) += m * xceeb; *(here->B3SOIPDGgPtr +1) += m * xcggb; *(here->B3SOIPDGdpPtr +1) += m * xcgdb; *(here->B3SOIPDGspPtr +1) += m * xcgsb; *(here->B3SOIPDDPgPtr +1) += m * xcdgb; *(here->B3SOIPDDPdpPtr +1) += m * xcddb; *(here->B3SOIPDDPspPtr +1) += m * xcdsb; *(here->B3SOIPDSPgPtr +1) += m * xcsgb; *(here->B3SOIPDSPdpPtr +1) += m * xcsdb; *(here->B3SOIPDSPspPtr +1) += m * xcssb; *(here->B3SOIPDBePtr +1) += m * xcbeb; *(here->B3SOIPDBgPtr +1) += m * xcbgb; *(here->B3SOIPDBdpPtr +1) += m * xcbdb; *(here->B3SOIPDBspPtr +1) += m * xcbsb; *(here->B3SOIPDEbPtr +1) -= m * (xcegb + xceeb + xcedb + xcesb); *(here->B3SOIPDGbPtr +1) -= m * (xcggb + xcgdb + xcgsb + xcgeb); *(here->B3SOIPDDPbPtr +1) -= m * (xcdgb + xcddb + xcdsb + xcdeb); *(here->B3SOIPDSPbPtr +1) -= m * (xcsgb + xcsdb + xcssb + xcseb); *(here->B3SOIPDBbPtr +1) -= m * (xcbgb + xcbdb + xcbsb + xcbeb); if (selfheat) { *(here->B3SOIPDTemptempPtr + 1) += m * xcTt; *(here->B3SOIPDDPtempPtr + 1) += m * xcdT; *(here->B3SOIPDSPtempPtr + 1) += m * xcsT; *(here->B3SOIPDBtempPtr + 1) += m * xcbT; *(here->B3SOIPDEtempPtr + 1) += m * xceT; *(here->B3SOIPDGtempPtr + 1) += m * xcgT; } *(here->B3SOIPDEePtr) += 0.0; *(here->B3SOIPDDPgPtr) += m * (Gm + gddpg); *(here->B3SOIPDDPdpPtr) += m * (gdpr + gds + gddpdp + RevSum); *(here->B3SOIPDDPspPtr) -= m * (gds + FwdSum - gddpsp); *(here->B3SOIPDDPdPtr) -= m * gdpr; *(here->B3SOIPDSPgPtr) -= m * (Gm - gsspg); *(here->B3SOIPDSPdpPtr) -= m * (gds + RevSum - gsspdp); *(here->B3SOIPDSPspPtr) += m * (gspr + gds + FwdSum + gsspsp); *(here->B3SOIPDSPsPtr) -= m * gspr; *(here->B3SOIPDBePtr) += 0; *(here->B3SOIPDBgPtr) += m * gbbg; *(here->B3SOIPDBdpPtr) += m * gbbdp; *(here->B3SOIPDBspPtr) += m * gbbsp; *(here->B3SOIPDBbPtr) += m * gbbb; *(here->B3SOIPDEbPtr) += 0.0; *(here->B3SOIPDSPbPtr) -= m * Gmbs - gsspb; *(here->B3SOIPDDPbPtr) -= m * (-gddpb - Gmbs); if (selfheat) { *(here->B3SOIPDDPtempPtr) += m * (GmT + gddpT); *(here->B3SOIPDSPtempPtr) += m * (-GmT + gsspT); *(here->B3SOIPDBtempPtr) += m * gbbT; *(here->B3SOIPDTemptempPtr) += m * (gTtt + 1/here->pParam->B3SOIPDrth); *(here->B3SOIPDTempgPtr) += m * gTtg; *(here->B3SOIPDTempbPtr) += m * gTtb; *(here->B3SOIPDTempdpPtr) += m * gTtdp; *(here->B3SOIPDTempspPtr) += m * gTtsp; } *(here->B3SOIPDDdPtr) += m * gdpr; *(here->B3SOIPDDdpPtr) -= m * gdpr; *(here->B3SOIPDSsPtr) += m * gspr; *(here->B3SOIPDSspPtr) -= m * gspr; if (here->B3SOIPDbodyMod == 1) { (*(here->B3SOIPDBpPtr) -= m * gppp); (*(here->B3SOIPDPbPtr) += m * gppb); (*(here->B3SOIPDPpPtr) += m * gppp); } if (here->B3SOIPDdebugMod != 0) { *(here->B3SOIPDVbsPtr) += m * 1; *(here->B3SOIPDIdsPtr) += m * 1; *(here->B3SOIPDIcPtr) += m * 1; *(here->B3SOIPDIbsPtr) += m * 1; *(here->B3SOIPDIbdPtr) += m * 1; *(here->B3SOIPDIiiPtr) += m * 1; *(here->B3SOIPDIgidlPtr) += m * 1; *(here->B3SOIPDItunPtr) += m * 1; *(here->B3SOIPDIbpPtr) += m * 1; *(here->B3SOIPDCbgPtr) += m * 1; *(here->B3SOIPDCbbPtr) += m * 1; *(here->B3SOIPDCbdPtr) += m * 1; *(here->B3SOIPDQbfPtr) += m * 1; *(here->B3SOIPDQjsPtr) += m * 1; *(here->B3SOIPDQjdPtr) += m * 1; } } } return(OK); }