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

Added MOS1 to the "new_kirchhoff" KCL verification

This commit is contained in:
Francesco Lannutti 2013-01-17 00:02:54 +01:00
parent d98e03c194
commit a1eeb1f278
1 changed files with 34 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

40
src/spicelib/devices/mos1/mos1load.c
View File

@ -33,13 +33,13 @@ MOS1load(GENmodel *inModel, CKTcircuit *ckt)
double cbhat; double cbhat;
double cdhat; double cdhat;
double cdrain; double cdrain;
double cdreq; double cdreq, cdreq_fvk ;
double ceq; double ceq;
double ceqbd; double ceqbd, ceqbd_fvk ;
double ceqbs; double ceqbs, ceqbs_fvk ;
double ceqgb; double ceqgb, ceqgb_fvk ;
double ceqgd; double ceqgd, ceqgd_fvk ;
double ceqgs; double ceqgs, ceqgs_fvk ;
double delvbd; double delvbd;
double delvbs; double delvbs;
double delvds; double delvds;
@ -847,10 +847,13 @@ next1: if(vbs <= -3*vt) {
*/ */
gcgs=0; gcgs=0;
ceqgs=0; ceqgs=0;
ceqgs_fvk = 0 ;
gcgd=0; gcgd=0;
ceqgd=0; ceqgd=0;
ceqgd_fvk = 0 ;
gcgb=0; gcgb=0;
ceqgb=0; ceqgb=0;
ceqgb_fvk = 0 ;
} else { } else {
if(capgs == 0) *(ckt->CKTstate0 + here->MOS1cqgs) =0; if(capgs == 0) *(ckt->CKTstate0 + here->MOS1cqgs) =0;
if(capgd == 0) *(ckt->CKTstate0 + here->MOS1cqgd) =0; if(capgd == 0) *(ckt->CKTstate0 + here->MOS1cqgd) =0;
@ -867,10 +870,13 @@ next1: if(vbs <= -3*vt) {
if(error) return(error); if(error) return(error);
ceqgs=ceqgs-gcgs*vgs+ckt->CKTag[0]* ceqgs=ceqgs-gcgs*vgs+ckt->CKTag[0]*
*(ckt->CKTstate0 + here->MOS1qgs); *(ckt->CKTstate0 + here->MOS1qgs);
ceqgs_fvk = *(ckt->CKTstate0 + here->MOS1cqgs) ;
ceqgd=ceqgd-gcgd*vgd+ckt->CKTag[0]* ceqgd=ceqgd-gcgd*vgd+ckt->CKTag[0]*
*(ckt->CKTstate0 + here->MOS1qgd); *(ckt->CKTstate0 + here->MOS1qgd);
ceqgd_fvk = *(ckt->CKTstate0 + here->MOS1cqgd) ;
ceqgb=ceqgb-gcgb*vgb+ckt->CKTag[0]* ceqgb=ceqgb-gcgb*vgb+ckt->CKTag[0]*
*(ckt->CKTstate0 + here->MOS1qgb); *(ckt->CKTstate0 + here->MOS1qgb);
ceqgb_fvk = *(ckt->CKTstate0 + here->MOS1cqgb) ;
} }
/* /*
* store charge storage info for meyer's cap in lx table * store charge storage info for meyer's cap in lx table
@ -881,18 +887,22 @@ next1: if(vbs <= -3*vt) {
*/ */
ceqbs = model->MOS1type * ceqbs = model->MOS1type *
(here->MOS1cbs-(here->MOS1gbs)*vbs); (here->MOS1cbs-(here->MOS1gbs)*vbs);
ceqbs_fvk = model->MOS1type * here->MOS1cbs ;
ceqbd = model->MOS1type * ceqbd = model->MOS1type *
(here->MOS1cbd-(here->MOS1gbd)*vbd); (here->MOS1cbd-(here->MOS1gbd)*vbd);
ceqbd_fvk = model->MOS1type * here->MOS1cbd ;
if (here->MOS1mode >= 0) { if (here->MOS1mode >= 0) {
xnrm=1; xnrm=1;
xrev=0; xrev=0;
cdreq=model->MOS1type*(cdrain-here->MOS1gds*vds- cdreq=model->MOS1type*(cdrain-here->MOS1gds*vds-
here->MOS1gm*vgs-here->MOS1gmbs*vbs); here->MOS1gm*vgs-here->MOS1gmbs*vbs);
cdreq_fvk = model->MOS1type * cdrain ;
} else { } else {
xnrm=0; xnrm=0;
xrev=1; xrev=1;
cdreq = -(model->MOS1type)*(cdrain-here->MOS1gds*(-vds)- cdreq = -(model->MOS1type)*(cdrain-here->MOS1gds*(-vds)-
here->MOS1gm*vgd-here->MOS1gmbs*vbd); here->MOS1gm*vgd-here->MOS1gmbs*vbd);
cdreq_fvk = - model->MOS1type * cdrain ;
} }
*(ckt->CKTrhs + here->MOS1gNode) -= *(ckt->CKTrhs + here->MOS1gNode) -=
(model->MOS1type * (ceqgs + ceqgb + ceqgd)); (model->MOS1type * (ceqgs + ceqgb + ceqgd));
@ -934,6 +944,24 @@ next1: if(vbs <= -3*vt) {
*(here->MOS1SPbPtr) += (-here->MOS1gbs-(xnrm-xrev)*here->MOS1gmbs); *(here->MOS1SPbPtr) += (-here->MOS1gbs-(xnrm-xrev)*here->MOS1gmbs);
*(here->MOS1SPdpPtr) += (-here->MOS1gds-xrev* *(here->MOS1SPdpPtr) += (-here->MOS1gds-xrev*
(here->MOS1gm+here->MOS1gmbs)); (here->MOS1gm+here->MOS1gmbs));
/* KCL verification - Linear-Dynamic Part */
*(ckt->CKTfvk+here->MOS1gNode) += model->MOS1type * (ceqgs_fvk + ceqgb_fvk + ceqgd_fvk) ;
*(ckt->CKTfvk+here->MOS1bNode) += ceqbs_fvk + ceqbd_fvk - model->MOS1type * ceqgb_fvk ;
*(ckt->CKTfvk+here->MOS1dNodePrime) -= (ceqbd_fvk - cdreq_fvk + model->MOS1type * ceqgd_fvk) ;
*(ckt->CKTfvk+here->MOS1sNodePrime) -= (cdreq_fvk + ceqbs_fvk + model->MOS1type * ceqgs_fvk) ;
/* KCL verification - Linear-Static Part */
*(ckt->CKTfvk+here->MOS1dNode) += here->MOS1drainConductance * *(ckt->CKTrhsOld+here->MOS1dNode) ;
*(ckt->CKTfvk+here->MOS1sNode) += here->MOS1sourceConductance * *(ckt->CKTrhsOld+here->MOS1sNode) ;
*(ckt->CKTfvk+here->MOS1dNodePrime) += here->MOS1drainConductance * *(ckt->CKTrhsOld+here->MOS1dNodePrime) ;
*(ckt->CKTfvk+here->MOS1sNodePrime) += here->MOS1sourceConductance * *(ckt->CKTrhsOld+here->MOS1sNodePrime) ;
*(ckt->CKTfvk+here->MOS1dNode) -= here->MOS1drainConductance * *(ckt->CKTrhsOld+here->MOS1dNodePrime) ;
*(ckt->CKTfvk+here->MOS1sNode) -= here->MOS1sourceConductance * *(ckt->CKTrhsOld+here->MOS1sNodePrime) ;
*(ckt->CKTfvk+here->MOS1dNodePrime) -= here->MOS1drainConductance * *(ckt->CKTrhsOld+here->MOS1dNode) ;
*(ckt->CKTfvk+here->MOS1sNodePrime) -= here->MOS1sourceConductance * *(ckt->CKTrhsOld+here->MOS1sNode) ;
} }
} }
return(OK); return(OK);