int
NUMD2acLoad(GENmodel *inModel, CKTcircuit *ckt)
{
register NUMD2model *model = (NUMD2model *) inModel;
register NUMD2instance *inst;
SPcomplex y;
double startTime;
/* loop through all the diode models */
for (; model != NULL; model = NUMD2nextModel(model)) {
FieldDepMobility = model->NUMD2models->MODLfieldDepMobility;
TransDepMobility = model->NUMD2models->MODLtransDepMobility;
SurfaceMobility = model->NUMD2models->MODLsurfaceMobility;
Srh = model->NUMD2models->MODLsrh;
Auger = model->NUMD2models->MODLauger;
AvalancheGen = model->NUMD2models->MODLavalancheGen;
OneCarrier = model->NUMD2methods->METHoneCarrier;
AcAnalysisMethod = model->NUMD2methods->METHacAnalysisMethod;
MobDeriv = model->NUMD2methods->METHmobDeriv;
TWOacDebug = model->NUMD2outputs->OUTPacDebug;
for (inst = NUMD2instances(model); inst != NULL;
inst = NUMD2nextInstance(inst)) {
startTime = SPfrontEnd->IFseconds();
/* Get Temp.-Dep. Global Parameters */
GLOBgetGlobals(&(inst->NUMD2globals));
model->NUMD2methods->METHacAnalysisMethod =
NUMD2admittance(inst->NUMD2pDevice, ckt->CKTomega, &y);
*(inst->NUMD2posPosPtr) += y.real;
*(inst->NUMD2posPosPtr + 1) += y.imag;
*(inst->NUMD2negNegPtr) += y.real;
*(inst->NUMD2negNegPtr + 1) += y.imag;
*(inst->NUMD2negPosPtr) -= y.real;
*(inst->NUMD2negPosPtr + 1) -= y.imag;
*(inst->NUMD2posNegPtr) -= y.real;
*(inst->NUMD2posNegPtr + 1) -= y.imag;
if (ckt->CKTomega != 0.0) {
inst->NUMD2c11 = y.imag / ckt->CKTomega;
} else {
inst->NUMD2c11 = 0.0; /* XXX What else can be done?! */
}
inst->NUMD2y11r = y.real;
inst->NUMD2y11i = y.imag;
inst->NUMD2smSigAvail = TRUE;
inst->NUMD2pDevice->pStats->totalTime[STAT_AC] +=
SPfrontEnd->IFseconds() - startTime;
}
}
return (OK);
}
int
NUMDpzLoad(GENmodel *inModel, CKTcircuit *ckt, SPcomplex *s)
{
register NUMDmodel *model = (NUMDmodel *) inModel;
register NUMDinstance *inst;
SPcomplex y;
double startTime;
NG_IGNORE(ckt);
/* loop through all the diode models */
for (; model != NULL; model = model->NUMDnextModel) {
FieldDepMobility = model->NUMDmodels->MODLfieldDepMobility;
Srh = model->NUMDmodels->MODLsrh;
Auger = model->NUMDmodels->MODLauger;
AvalancheGen = model->NUMDmodels->MODLavalancheGen;
AcAnalysisMethod = model->NUMDmethods->METHacAnalysisMethod;
MobDeriv = model->NUMDmethods->METHmobDeriv;
ONEacDebug = model->NUMDoutputs->OUTPacDebug;
for (inst = model->NUMDinstances; inst != NULL;
inst = inst->NUMDnextInstance) {
startTime = SPfrontEnd->IFseconds();
/* Get Temp.-Dep. Global Parameters */
GLOBgetGlobals(&(inst->NUMDglobals));
NUMDys(inst->NUMDpDevice, s, &y);
*(inst->NUMDposPosPtr) += y.real;
*(inst->NUMDposPosPtr + 1) += y.imag;
*(inst->NUMDnegNegPtr) += y.real;
*(inst->NUMDnegNegPtr + 1) += y.imag;
*(inst->NUMDnegPosPtr) -= y.real;
*(inst->NUMDnegPosPtr + 1) -= y.imag;
*(inst->NUMDposNegPtr) -= y.real;
*(inst->NUMDposNegPtr + 1) -= y.imag;
inst->NUMDpDevice->pStats->totalTime[STAT_AC] +=
SPfrontEnd->IFseconds() - startTime;
}
}
return (OK);
}
int
NBJTpzLoad(GENmodel *inModel, CKTcircuit *ckt, SPcomplex *s)
{
register NBJTmodel *model = (NBJTmodel *) inModel;
register NBJTinstance *inst;
SPcomplex yIeVce, yIeVbe;
SPcomplex yIcVce, yIcVbe;
double startTime;
NG_IGNORE(ckt);
for (; model != NULL; model = NBJTnextModel(model)) {
FieldDepMobility = model->NBJTmodels->MODLfieldDepMobility;
Srh = model->NBJTmodels->MODLsrh;
Auger = model->NBJTmodels->MODLauger;
AvalancheGen = model->NBJTmodels->MODLavalancheGen;
AcAnalysisMethod = model->NBJTmethods->METHacAnalysisMethod;
MobDeriv = model->NBJTmethods->METHmobDeriv;
ONEacDebug = model->NBJToutputs->OUTPacDebug;
for (inst = NBJTinstances(model); inst != NULL;
inst = NBJTnextInstance(inst)) {
startTime = SPfrontEnd->IFseconds();
/* Get Temp.-Dep. Global Parameters */
GLOBgetGlobals(&(inst->NBJTglobals));
NBJTys(inst->NBJTpDevice, s,
&yIeVce, &yIcVce, &yIeVbe, &yIcVbe);
if (ONEacDebug) {
fprintf(stdout, "BJT admittances: %s:%s at s = % .5g, % .5g\n",
model->NBJTmodName, inst->NBJTname, s->real, s->imag);
fprintf(stdout, "Ycc: % .5g,% .5g\n",
yIcVce.real, yIcVce.imag);
fprintf(stdout, "Ycb: % .5g,% .5g\n",
yIcVbe.real, yIcVbe.imag);
fprintf(stdout, "Ybc: % .5g,% .5g\n",
yIeVce.real - yIcVce.real, yIeVce.imag - yIcVce.imag);
fprintf(stdout, "Ybb: % .5g,% .5g\n",
yIeVbe.real - yIcVbe.real, yIeVbe.imag - yIcVbe.imag);
}
*(inst->NBJTcolColPtr) += yIcVce.real;
*(inst->NBJTcolColPtr + 1) += yIcVce.imag;
*(inst->NBJTcolBasePtr) += yIcVbe.real;
*(inst->NBJTcolBasePtr + 1) += yIcVbe.imag;
*(inst->NBJTcolEmitPtr) -= yIcVbe.real + yIcVce.real;
*(inst->NBJTcolEmitPtr + 1) -= yIcVbe.imag + yIcVce.imag;
*(inst->NBJTbaseColPtr) -= yIcVce.real - yIeVce.real;
*(inst->NBJTbaseColPtr + 1) -= yIcVce.imag - yIeVce.imag;
*(inst->NBJTbaseBasePtr) -= yIcVbe.real - yIeVbe.real;
*(inst->NBJTbaseBasePtr + 1) -= yIcVbe.imag - yIeVbe.imag;
*(inst->NBJTbaseEmitPtr) += yIcVbe.real + yIcVce.real - yIeVbe.real - yIeVce.real;
*(inst->NBJTbaseEmitPtr + 1) += yIcVbe.imag + yIcVce.imag - yIeVbe.imag - yIeVce.imag;
*(inst->NBJTemitColPtr) -= yIeVce.real;
*(inst->NBJTemitColPtr + 1) -= yIeVce.imag;
*(inst->NBJTemitBasePtr) -= yIeVbe.real;
*(inst->NBJTemitBasePtr + 1) -= yIeVbe.imag;
*(inst->NBJTemitEmitPtr) += yIeVbe.real + yIeVce.real;
*(inst->NBJTemitEmitPtr + 1) += yIeVbe.imag + yIeVce.imag;
inst->NBJTpDevice->pStats->totalTime[STAT_AC] +=
SPfrontEnd->IFseconds() - startTime;
}
}
return (OK);
}
int
NBJTload(GENmodel *inModel, CKTcircuit *ckt)
{
register NBJTmodel *model = (NBJTmodel *) inModel;
register NBJTinstance *inst;
register ONEdevice *pDevice;
double startTime, startTime2, totalTime, totalTime2;
double tol;
double ic, ie;
double iceq, ieeq;
double ichat = 0.0, iehat = 0.0;
double delVce, delVbe;
double vce, vbe /*, vbc*/;
double dIeDVce, dIeDVbe;
double dIcDVce, dIcDVbe;
double xfact;
int icheck;
int icheck1;
int i;
double deltaNorm[7];
int devConverged = 0;
int numDevNonCon;
int deviceType;
int doInitSolve;
int doVoltPred;
char *initStateName;
/* loop through all the models */
for (; model != NULL; model = NBJTnextModel(model)) {
FieldDepMobility = model->NBJTmodels->MODLfieldDepMobility;
Srh = model->NBJTmodels->MODLsrh;
Auger = model->NBJTmodels->MODLauger;
AvalancheGen = model->NBJTmodels->MODLavalancheGen;
MobDeriv = model->NBJTmethods->METHmobDeriv;
MaxIterations = model->NBJTmethods->METHitLim;
ONEdcDebug = model->NBJToutputs->OUTPdcDebug;
ONEtranDebug = model->NBJToutputs->OUTPtranDebug;
ONEacDebug = model->NBJToutputs->OUTPacDebug;
deviceType = model->NBJToptions->OPTNdeviceType;
doVoltPred = model->NBJTmethods->METHvoltPred;
if (ckt->CKTmode & MODEINITPRED) {
/* compute normalized deltas and predictor coeff */
if (!(ckt->CKTmode & MODEDCTRANCURVE)) {
model->NBJTpInfo->order = ckt->CKTorder;
model->NBJTpInfo->method = ckt->CKTintegrateMethod;
for (i = 0; i <= ckt->CKTmaxOrder; i++) {
deltaNorm[i] = ckt->CKTdeltaOld[i] / TNorm;
}
computeIntegCoeff(ckt->CKTintegrateMethod, ckt->CKTorder,
model->NBJTpInfo->intCoeff, deltaNorm);
computePredCoeff(ckt->CKTintegrateMethod, ckt->CKTorder,
model->NBJTpInfo->predCoeff, deltaNorm);
}
} else if (ckt->CKTmode & MODEINITTRAN) {
model->NBJTpInfo->order = ckt->CKTorder;
model->NBJTpInfo->method = ckt->CKTintegrateMethod;
for (i = 0; i <= ckt->CKTmaxOrder; i++) {
deltaNorm[i] = ckt->CKTdeltaOld[i] / TNorm;
}
computeIntegCoeff(ckt->CKTintegrateMethod, ckt->CKTorder,
model->NBJTpInfo->intCoeff, deltaNorm);
}
/* loop through all the instances of the model */
for (inst = NBJTinstances(model); inst != NULL;
inst = NBJTnextInstance(inst)) {
pDevice = inst->NBJTpDevice;
totalTime = 0.0;
startTime = SPfrontEnd->IFseconds();
/* Get Temp.-Dep. Global Parameters */
GLOBgetGlobals(&(inst->NBJTglobals));
/*
* initialization
*/
pDevice->devStates = ckt->CKTstates;
icheck = 1;
doInitSolve = FALSE;
initStateName = NULL;
if (ckt->CKTmode & MODEINITSMSIG) {
vbe = *(ckt->CKTstate0 + inst->NBJTvbe);
vce = *(ckt->CKTstate0 + inst->NBJTvce);
delVbe = 0.0;
delVce = 0.0;
NBJTsetBCs(pDevice, vce, vbe);
} else if (ckt->CKTmode & MODEINITTRAN) {
*(ckt->CKTstate0 + inst->NBJTvbe) =
*(ckt->CKTstate1 + inst->NBJTvbe);
*(ckt->CKTstate0 + inst->NBJTvce) =
*(ckt->CKTstate1 + inst->NBJTvce);
vbe = *(ckt->CKTstate1 + inst->NBJTvbe);
vce = *(ckt->CKTstate1 + inst->NBJTvce);
ONEsaveState(pDevice);
delVbe = 0.0;
delVce = 0.0;
} else if ((ckt->CKTmode & MODEINITJCT) &&
(ckt->CKTmode & MODETRANOP) && (ckt->CKTmode & MODEUIC)) {
doInitSolve = TRUE;
initStateName = inst->NBJTicFile;
vbe = 0.0;
vce = 0.0;
delVbe = vbe;
delVce = vce;
} else if ((ckt->CKTmode & MODEINITJCT) && (inst->NBJToff == 0)) {
doInitSolve = TRUE;
initStateName = inst->NBJTicFile;
vbe = inst->NBJTtype * 0.6;
vce = inst->NBJTtype * 1.0;
delVbe = vbe;
delVce = vce;
} else if (ckt->CKTmode & MODEINITJCT) {
doInitSolve = TRUE;
vbe = 0.0;
vce = 0.0;
delVbe = vbe;
delVce = vce;
} else if ((ckt->CKTmode & MODEINITFIX) && inst->NBJToff) {
vbe = 0.0;
vce = 0.0;
delVbe = vbe;
delVce = vce;
} else {
if (ckt->CKTmode & MODEINITPRED) {
*(ckt->CKTstate0 + inst->NBJTvbe) =
*(ckt->CKTstate1 + inst->NBJTvbe);
*(ckt->CKTstate0 + inst->NBJTvce) =
*(ckt->CKTstate1 + inst->NBJTvce);
*(ckt->CKTstate0 + inst->NBJTic) =
*(ckt->CKTstate1 + inst->NBJTic);
*(ckt->CKTstate0 + inst->NBJTie) =
*(ckt->CKTstate1 + inst->NBJTie);
*(ckt->CKTstate0 + inst->NBJTdIeDVce) =
*(ckt->CKTstate1 + inst->NBJTdIeDVce);
*(ckt->CKTstate0 + inst->NBJTdIeDVbe) =
*(ckt->CKTstate1 + inst->NBJTdIeDVbe);
*(ckt->CKTstate0 + inst->NBJTdIcDVce) =
*(ckt->CKTstate1 + inst->NBJTdIcDVce);
*(ckt->CKTstate0 + inst->NBJTdIcDVbe) =
*(ckt->CKTstate1 + inst->NBJTdIcDVbe);
if (!(ckt->CKTmode & MODEDCTRANCURVE)) {
/* no linear prediction on device voltages */
vbe = *(ckt->CKTstate1 + inst->NBJTvbe);
vce = *(ckt->CKTstate1 + inst->NBJTvce);
ONEpredict(pDevice, model->NBJTpInfo);
} else {
if (doVoltPred) {
/* linear prediction */
xfact=ckt->CKTdelta/ckt->CKTdeltaOld[1];
vbe = (1+xfact) * (*(ckt->CKTstate1 + inst->NBJTvbe))
- (xfact) * (*(ckt->CKTstate2 + inst->NBJTvbe));
vce = (1+xfact) * (*(ckt->CKTstate1 + inst->NBJTvce))
- (xfact) * (*(ckt->CKTstate2 + inst->NBJTvce));
} else {
vbe = *(ckt->CKTstate1 + inst->NBJTvbe);
vce = *(ckt->CKTstate1 + inst->NBJTvce);
}
}
} else {
/*
* compute new nonlinear branch voltages
*/
vbe = *(ckt->CKTrhsOld + inst->NBJTbaseNode) -
*(ckt->CKTrhsOld + inst->NBJTemitNode);
vce = *(ckt->CKTrhsOld + inst->NBJTcolNode) -
*(ckt->CKTrhsOld + inst->NBJTemitNode);
}
delVbe = vbe - *(ckt->CKTstate0 + inst->NBJTvbe);
delVce = vce - *(ckt->CKTstate0 + inst->NBJTvce);
ichat = *(ckt->CKTstate0 + inst->NBJTic) -
*(ckt->CKTstate0 + inst->NBJTdIcDVbe) * delVbe -
*(ckt->CKTstate0 + inst->NBJTdIcDVce) * delVce;
iehat = *(ckt->CKTstate0 + inst->NBJTie) -
*(ckt->CKTstate0 + inst->NBJTdIeDVbe) * delVbe -
*(ckt->CKTstate0 + inst->NBJTdIeDVce) * delVce;
#ifndef NOBYPASS
/*
* bypass if solution has not changed
*/
/*
* the following collections of if's would be just one if the average
* compiler could handle it, but many find the expression too
* complicated, thus the split.
*/
if ((ckt->CKTbypass) && pDevice->converged &&
(!(ckt->CKTmode & MODEINITPRED)) &&
(fabs(delVbe) < (ckt->CKTreltol * MAX(fabs(vbe),
fabs(*(ckt->CKTstate0 + inst->NBJTvbe))) +
ckt->CKTvoltTol)))
if ((fabs(delVce) < ckt->CKTreltol * MAX(fabs(vce),
fabs(*(ckt->CKTstate0 + inst->NBJTvce))) +
ckt->CKTvoltTol))
if ((fabs(ichat - *(ckt->CKTstate0 + inst->NBJTic)) <
ckt->CKTreltol * MAX(fabs(ichat),
fabs(*(ckt->CKTstate0 + inst->NBJTic))) +
ckt->CKTabstol))
if ((fabs(iehat - *(ckt->CKTstate0 + inst->NBJTie)) <
ckt->CKTreltol * MAX(fabs(iehat),
fabs(*(ckt->CKTstate0 + inst->NBJTie))) +
ckt->CKTabstol)) {
/*
* bypassing....
*/
vbe = *(ckt->CKTstate0 + inst->NBJTvbe);
vce = *(ckt->CKTstate0 + inst->NBJTvce);
ic = *(ckt->CKTstate0 + inst->NBJTic);
ie = *(ckt->CKTstate0 + inst->NBJTie);
dIeDVce = *(ckt->CKTstate0 + inst->NBJTdIeDVce);
dIeDVbe = *(ckt->CKTstate0 + inst->NBJTdIeDVbe);
dIcDVce = *(ckt->CKTstate0 + inst->NBJTdIcDVce);
dIcDVbe = *(ckt->CKTstate0 + inst->NBJTdIcDVbe);
goto load;
}
#endif /* NOBYPASS */
/*
* limit nonlinear branch voltages
*/
icheck1 = 1;
vbe = inst->NBJTtype * limitVbe(inst->NBJTtype * vbe,
inst->NBJTtype * *(ckt->CKTstate0 + inst->NBJTvbe), &icheck);
/*
vbc = vbe - vce;
vbc = inst->NBJTtype * limitVbe(inst->NBJTtype * vbc,
inst->NBJTtype * (*(ckt->CKTstate0 + inst->NBJTvbe)
- *(ckt->CKTstate0 + inst->NBJTvce)), &icheck1);
vce = vbe - vbc;
*/
vce = inst->NBJTtype * limitVce(inst->NBJTtype * vce,
inst->NBJTtype * *(ckt->CKTstate0 + inst->NBJTvce), &icheck1);
if (icheck1 == 1)
icheck = 1;
delVbe = vbe - *(ckt->CKTstate0 + inst->NBJTvbe);
delVce = vce - *(ckt->CKTstate0 + inst->NBJTvce);
}
if (doInitSolve) {
if (ONEdcDebug) {
printVoltages(stdout,
model->NBJTmodName, inst->NBJTname,
deviceType, 2, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0);
}
startTime2 = SPfrontEnd->IFseconds();
ONEequilSolve(pDevice);
totalTime2 = SPfrontEnd->IFseconds() - startTime2;
pDevice->pStats->totalTime[STAT_SETUP] += totalTime2;
pDevice->pStats->totalTime[STAT_DC] -= totalTime2;
ONEbiasSolve(pDevice, MaxIterations, FALSE, NULL);
*(ckt->CKTstate0 + inst->NBJTvbe) = 0.0;
*(ckt->CKTstate0 + inst->NBJTvce) = 0.0;
if (initStateName != NULL) {
if (ONEreadState(pDevice, initStateName, 2, &vce, &vbe ) < 0) {
fprintf(stderr,
"NBJTload: trouble reading state-file %s\n", initStateName);
} else {
NBJTsetBCs(pDevice, vce, vbe);
delVce = delVbe = 0.0;
}
}
}
/*
* determine dc current and derivatives using the numerical routines
*/
if (ckt->CKTmode & (MODEDCOP | MODETRANOP | MODEDCTRANCURVE | MODEINITSMSIG)) {
numDevNonCon = 0;
inst->NBJTc11 = inst->NBJTy11r = inst->NBJTy11i = 0.0;
inst->NBJTc12 = inst->NBJTy12r = inst->NBJTy12i = 0.0;
inst->NBJTc21 = inst->NBJTy21r = inst->NBJTy21i = 0.0;
inst->NBJTc22 = inst->NBJTy22r = inst->NBJTy22i = 0.0;
inst->NBJTsmSigAvail = FALSE;
devNonCon:
NBJTproject(pDevice, delVce, delVbe, vbe);
if (ONEdcDebug) {
printVoltages(stdout,
model->NBJTmodName, inst->NBJTname,
deviceType, 2, vce, delVce, vbe, delVbe, 0.0, 0.0);
}
ONEbiasSolve(pDevice, MaxIterations, FALSE, model->NBJTpInfo);
devConverged = pDevice->converged;
if (devConverged && finite(pDevice->rhsNorm)) {
/* compute the currents */
NBJTcurrent(pDevice, FALSE, NULL,
&ie, &ic);
NBJTconductance(pDevice, FALSE, NULL,
&dIeDVce, &dIcDVce, &dIeDVbe, &dIcDVbe);
/*
* Add Gmin terms to everything in case we are operating at
* abnormally low current levels
*/
ie += ckt->CKTgmin * (vce + vbe);
dIeDVce += ckt->CKTgmin;
dIeDVbe += ckt->CKTgmin;
ic += ckt->CKTgmin * (2.0 * vce - vbe);
dIcDVce += 2.0 * ckt->CKTgmin;
dIcDVbe -= ckt->CKTgmin;
} else {
/* reduce the voltage step until converged */
/* restore the boundary potential to previous value */
NBJTsetBCs(pDevice, vce - delVce, vbe - delVbe);
ONEstoreInitialGuess(pDevice);
ONEresetJacobian(pDevice);
delVbe *= 0.5;
delVce *= 0.5;
vbe = delVbe + *(ckt->CKTstate0 + inst->NBJTvbe);
vce = delVce + *(ckt->CKTstate0 + inst->NBJTvce);
numDevNonCon++;
icheck = 1;
if (numDevNonCon > 10) {
printVoltages(stderr,
model->NBJTmodName, inst->NBJTname,
deviceType, 2, vce, delVce, vbe, delVbe, 0.0, 0.0);
fprintf(stderr,
"*** Non-convergence during load ***\n");
totalTime += SPfrontEnd->IFseconds() - startTime;
pDevice->pStats->totalTime[STAT_DC] += totalTime;
ckt->CKTtroubleElt = (GENinstance *) inst;
return (E_BADMATRIX);
} else {
goto devNonCon;
}
}
}
if ((ckt->CKTmode & (MODETRAN | MODEAC)) ||
((ckt->CKTmode & MODETRANOP) && (ckt->CKTmode & MODEUIC)) ||
(ckt->CKTmode & MODEINITSMSIG)) {
/*
* store small-signal parameters
*/
if ((!(ckt->CKTmode & MODETRANOP)) ||
(!(ckt->CKTmode & MODEUIC))) {
if (ckt->CKTmode & MODEINITSMSIG) {
totalTime += SPfrontEnd->IFseconds() - startTime;
pDevice->pStats->totalTime[STAT_DC] += totalTime;
startTime2 = SPfrontEnd->IFseconds();
NBJTinitSmSig(inst);
pDevice->pStats->totalTime[STAT_AC] +=
SPfrontEnd->IFseconds() - startTime2;
continue;
} else {
inst->NBJTsmSigAvail = FALSE;
}
/*
* transient analysis
*/
if (ckt->CKTmode & MODEINITPRED) {
NBJTsetBCs(pDevice, vce, vbe);
ONEstoreInitialGuess(pDevice);
} else {
NBJTupdate(pDevice, delVce, delVbe, vbe, TRUE);
}
if (ONEtranDebug) {
printVoltages(stdout,
model->NBJTmodName, inst->NBJTname,
deviceType, 2, vce, delVce, vbe, delVbe, 0.0, 0.0);
}
ONEbiasSolve(pDevice, 0, TRUE, model->NBJTpInfo);
if (!finite(pDevice->rhsNorm)) {
/* Timestep took us to never-never land. */
totalTime += SPfrontEnd->IFseconds() - startTime;
pDevice->pStats->totalTime[STAT_TRAN] += totalTime;
ckt->CKTtroubleElt = (GENinstance *) inst;
return (E_BADMATRIX);
}
devConverged = ONEdeviceConverged(pDevice);
pDevice->converged = devConverged;
/* compute the currents */
NBJTcurrent(pDevice, TRUE,
model->NBJTpInfo->intCoeff, &ie, &ic);
NBJTconductance(pDevice, TRUE,
model->NBJTpInfo->intCoeff,
&dIeDVce, &dIcDVce, &dIeDVbe, &dIcDVbe);
/*
* Add Gmin terms to everything in case we are operating at
* abnormally low current levels
*/
ie += ckt->CKTgmin * (vce + vbe);
dIeDVce += ckt->CKTgmin;
dIeDVbe += ckt->CKTgmin;
ic += ckt->CKTgmin * (2.0 * vce - vbe);
dIcDVce += 2.0 * ckt->CKTgmin;
dIcDVbe -= ckt->CKTgmin;
}
}
/*
* check convergence
*/
if ((!(ckt->CKTmode & MODEINITFIX)) || (!(inst->NBJToff))) {
if (icheck == 1 || !devConverged) {
ckt->CKTnoncon++;
ckt->CKTtroubleElt = (GENinstance *) inst;
} else {
tol = ckt->CKTreltol * MAX(fabs(ichat), fabs(ic)) + ckt->CKTabstol;
if (fabs(ichat - ic) > tol) {
ckt->CKTnoncon++;
ckt->CKTtroubleElt = (GENinstance *) inst;
} else {
tol = ckt->CKTreltol * MAX(fabs(iehat), fabs(ie)) +
ckt->CKTabstol;
if (fabs(iehat - ie) > tol) {
ckt->CKTnoncon++;
ckt->CKTtroubleElt = (GENinstance *) inst;
}
}
}
}
*(ckt->CKTstate0 + inst->NBJTvbe) = vbe;
*(ckt->CKTstate0 + inst->NBJTvce) = vce;
*(ckt->CKTstate0 + inst->NBJTic) = ic;
*(ckt->CKTstate0 + inst->NBJTie) = ie;
*(ckt->CKTstate0 + inst->NBJTdIeDVce) = dIeDVce;
*(ckt->CKTstate0 + inst->NBJTdIeDVbe) = dIeDVbe;
*(ckt->CKTstate0 + inst->NBJTdIcDVce) = dIcDVce;
*(ckt->CKTstate0 + inst->NBJTdIcDVbe) = dIcDVbe;
load:
/*
* load current excitation vector
*/
iceq = ic - dIcDVce * vce - dIcDVbe * vbe;
ieeq = ie - dIeDVce * vce - dIeDVbe * vbe;
*(ckt->CKTrhs + inst->NBJTcolNode) -= iceq;
*(ckt->CKTrhs + inst->NBJTbaseNode) -= ieeq - iceq;
*(ckt->CKTrhs + inst->NBJTemitNode) += ieeq;
/*
* load y matrix
*/
*(inst->NBJTcolColPtr) += dIcDVce;
*(inst->NBJTcolBasePtr) += dIcDVbe;
*(inst->NBJTcolEmitPtr) -= dIcDVbe + dIcDVce;
*(inst->NBJTbaseColPtr) -= dIcDVce - dIeDVce;
*(inst->NBJTbaseBasePtr) -= dIcDVbe - dIeDVbe;
*(inst->NBJTbaseEmitPtr) += dIcDVbe + dIcDVce - dIeDVbe - dIeDVce;
*(inst->NBJTemitColPtr) -= dIeDVce;
*(inst->NBJTemitBasePtr) -= dIeDVbe;
*(inst->NBJTemitEmitPtr) += dIeDVbe + dIeDVce;
totalTime += SPfrontEnd->IFseconds() - startTime;
if (ckt->CKTmode & MODETRAN) {
pDevice->pStats->totalTime[STAT_TRAN] += totalTime;
} else {
pDevice->pStats->totalTime[STAT_DC] += totalTime;
}
}
}
return (OK);
}
int
NBJTacLoad(GENmodel *inModel, CKTcircuit *ckt)
{
register NBJTmodel *model = (NBJTmodel *) inModel;
register NBJTinstance *inst;
SPcomplex yIeVce, yIeVbe;
SPcomplex yIcVce, yIcVbe;
double startTime;
for (; model != NULL; model = model->NBJTnextModel) {
FieldDepMobility = model->NBJTmodels->MODLfieldDepMobility;
Srh = model->NBJTmodels->MODLsrh;
Auger = model->NBJTmodels->MODLauger;
AvalancheGen = model->NBJTmodels->MODLavalancheGen;
AcAnalysisMethod = model->NBJTmethods->METHacAnalysisMethod;
MobDeriv = model->NBJTmethods->METHmobDeriv;
ONEacDebug = model->NBJToutputs->OUTPacDebug;
for (inst = model->NBJTinstances; inst != NULL;
inst = inst->NBJTnextInstance) {
startTime = SPfrontEnd->IFseconds();
/* Get Temp.-Dep. Global Parameters */
GLOBgetGlobals(&(inst->NBJTglobals));
model->NBJTmethods->METHacAnalysisMethod =
NBJTadmittance(inst->NBJTpDevice, ckt->CKTomega,
&yIeVce, &yIcVce, &yIeVbe, &yIcVbe);
*(inst->NBJTcolColPtr) += yIcVce.real;
*(inst->NBJTcolColPtr + 1) += yIcVce.imag;
*(inst->NBJTcolBasePtr) += yIcVbe.real;
*(inst->NBJTcolBasePtr + 1) += yIcVbe.imag;
*(inst->NBJTcolEmitPtr) -= yIcVbe.real + yIcVce.real;
*(inst->NBJTcolEmitPtr + 1) -= yIcVbe.imag + yIcVce.imag;
*(inst->NBJTbaseColPtr) -= yIcVce.real - yIeVce.real;
*(inst->NBJTbaseColPtr + 1) -= yIcVce.imag - yIeVce.imag;
*(inst->NBJTbaseBasePtr) -= yIcVbe.real - yIeVbe.real;
*(inst->NBJTbaseBasePtr + 1) -= yIcVbe.imag - yIeVbe.imag;
*(inst->NBJTbaseEmitPtr) += yIcVbe.real + yIcVce.real - yIeVbe.real - yIeVce.real;
*(inst->NBJTbaseEmitPtr + 1) += yIcVbe.imag + yIcVce.imag - yIeVbe.imag - yIeVce.imag;
*(inst->NBJTemitColPtr) -= yIeVce.real;
*(inst->NBJTemitColPtr + 1) -= yIeVce.imag;
*(inst->NBJTemitBasePtr) -= yIeVbe.real;
*(inst->NBJTemitBasePtr + 1) -= yIeVbe.imag;
*(inst->NBJTemitEmitPtr) += yIeVbe.real + yIeVce.real;
*(inst->NBJTemitEmitPtr + 1) += yIeVbe.imag + yIeVce.imag;
if (ckt->CKTomega != 0.0) {
inst->NBJTc11 = yIcVce.imag / ckt->CKTomega;
inst->NBJTc12 = yIcVbe.imag / ckt->CKTomega;
inst->NBJTc21 = (yIeVce.imag - yIcVce.imag) / ckt->CKTomega;
inst->NBJTc22 = (yIeVbe.imag - yIcVbe.imag) / ckt->CKTomega;
} else {
inst->NBJTc11 = 0.0; /* XXX What else can be done?! */
inst->NBJTc12 = 0.0; /* XXX What else can be done?! */
inst->NBJTc21 = 0.0; /* XXX What else can be done?! */
inst->NBJTc22 = 0.0; /* XXX What else can be done?! */
}
inst->NBJTy11r = yIcVce.real;
inst->NBJTy11i = yIcVce.imag;
inst->NBJTy12r = yIcVbe.real;
inst->NBJTy12i = yIcVbe.imag;
inst->NBJTy21r = yIeVce.real - yIcVce.real;
inst->NBJTy21i = yIeVce.imag - yIcVce.imag;
inst->NBJTy22r = yIeVbe.real - yIcVbe.real;
inst->NBJTy22i = yIeVbe.imag - yIcVbe.imag;
inst->NBJTsmSigAvail = TRUE;
inst->NBJTpDevice->pStats->totalTime[STAT_AC] +=
SPfrontEnd->IFseconds() - startTime;
}
}
return (OK);
}