int
NUMDsetup(SMPmatrix *matrix, GENmodel *inModel, CKTcircuit *ckt, int *states)
/*
* load the structure with those pointers needed later for fast matrix
* loading
*/
{
register NUMDmodel *model = (NUMDmodel *) inModel;
register NUMDinstance *inst;
METHcard *methods;
MODLcard *models;
OPTNcard *options;
OUTPcard *outputs;
char *icFileName = NULL;
size_t nameLen;
int error;
int xMeshSize;
ONEdevice *pDevice;
ONEcoord *xCoordList = NULL;
ONEdomain *domainList = NULL;
ONEmaterial *pM, *pMaterial = NULL, *materialList = NULL;
DOPprofile *profileList = NULL;
DOPtable *dopTableList = NULL;
double startTime;
/* loop through all the models */
for (; model != NULL; model = model->NUMDnextModel) {
if (!model->NUMDpInfo) {
TSCALLOC(model->NUMDpInfo, 1, ONEtranInfo);
}
methods = model->NUMDmethods;
if (!methods) {
TSCALLOC(methods, 1, METHcard);
model->NUMDmethods = methods;
}
models = model->NUMDmodels;
if (!models) {
TSCALLOC(models, 1, MODLcard);
model->NUMDmodels = models;
}
options = model->NUMDoptions;
if (!options) {
TSCALLOC(options, 1, OPTNcard);
model->NUMDoptions = options;
}
outputs = model->NUMDoutputs;
if (!outputs) {
TSCALLOC(outputs, 1, OUTPcard);
model->NUMDoutputs = outputs;
}
if (!methods->METHvoltPredGiven) {
methods->METHvoltPred = FALSE;
}
if (!methods->METHmobDerivGiven) {
methods->METHmobDeriv = TRUE;
}
if (!methods->METHoneCarrierGiven) {
methods->METHoneCarrier = FALSE;
}
if (!methods->METHacAnalysisMethodGiven) {
methods->METHacAnalysisMethod = SOR;
}
if (!methods->METHdabstolGiven) {
methods->METHdabstol = DABSTOL1D;
}
if (!methods->METHdreltolGiven) {
methods->METHdreltol = ckt->CKTreltol;
}
if (!methods->METHitLimGiven) {
methods->METHitLim = 20;
}
if (!methods->METHomegaGiven || methods->METHomega <= 0.0) {
methods->METHomega = 2.0 * M_PI /* radians/sec */ ;
}
if (!options->OPTNdefaGiven || options->OPTNdefa <= 0.0) {
options->OPTNdefa = 1.0e4 /* cm^2 */ ;
}
if (!options->OPTNdeviceTypeGiven) {
options->OPTNdeviceType = OPTN_DIODE;
}
if (!options->OPTNicFileGiven) {
options->OPTNicFile = NULL;
options->OPTNunique = FALSE; /* Can't form a unique name. */
}
if (!options->OPTNuniqueGiven) {
options->OPTNunique = FALSE;
}
/* Set up the rest of the card lists */
if ((error = MODLsetup(model->NUMDmodels)) != 0)
return (error);
BandGapNarrowing = models->MODLbandGapNarrowing;
ConcDepLifetime = models->MODLconcDepLifetime;
TempDepMobility = models->MODLtempDepMobility;
ConcDepMobility = models->MODLconcDepMobility;
if ((error = OUTPsetup(model->NUMDoutputs)) != 0)
return (error);
if ((error = MATLsetup(model->NUMDmaterials, &materialList)) != 0)
return (error);
if ((error = MOBsetup(model->NUMDmobility, materialList)) != 0)
return (error);
if ((error = MESHsetup('x', model->NUMDxMeshes, &xCoordList, &xMeshSize)) != 0)
return (error);
if ((error = DOMNsetup(model->NUMDdomains, &domainList,
xCoordList, NULL, materialList)) != 0)
return (error);
if ((error = BDRYsetup(model->NUMDboundaries,
xCoordList, NULL, domainList)) != 0)
return (error);
if ((error = CONTsetup(model->NUMDcontacts, NULL)) != 0)
return (error);
if ((error = DOPsetup(model->NUMDdopings, &profileList,
&dopTableList, xCoordList, NULL)) != 0)
return (error);
model->NUMDmatlInfo = materialList;
model->NUMDprofiles = profileList;
model->NUMDdopTables = dopTableList;
/* loop through all the instances of the model */
for (inst = model->NUMDinstances; inst != NULL;
inst = inst->NUMDnextInstance) {
startTime = SPfrontEnd->IFseconds();
if ((!inst->NUMDprintGiven)) {
inst->NUMDprint = 0;
} else if (inst->NUMDprint <= 0) {
inst->NUMDprint = 1;
}
if ((!inst->NUMDicFileGiven)) {
if (options->OPTNunique) {
nameLen = strlen(options->OPTNicFile) + strlen(inst->NUMDname) + 1;
TSCALLOC(icFileName, nameLen+1, char);
sprintf(icFileName, "%s.%s", options->OPTNicFile, inst->NUMDname);
icFileName[nameLen] = '\0';
inst->NUMDicFile = icFileName;
} else if (options->OPTNicFile != NULL) {
nameLen = strlen(options->OPTNicFile);
TSCALLOC(icFileName, nameLen+1, char);
icFileName = strcpy(icFileName, options->OPTNicFile);
inst->NUMDicFile = icFileName;
} else {
int
NUMDsetup(SMPmatrix *matrix, GENmodel *inModel, CKTcircuit *ckt, int *states)
/*
* load the structure with those pointers needed later for fast matrix
* loading
*/
{
register NUMDmodel *model = (NUMDmodel *) inModel;
register NUMDinstance *inst;
METHcard *methods;
MODLcard *models;
OPTNcard *options;
OUTPcard *outputs;
int error;
int xMeshSize;
ONEdevice *pDevice;
ONEcoord *xCoordList = NULL;
ONEdomain *domainList = NULL;
ONEmaterial *pM, *pMaterial = NULL, *materialList = NULL;
DOPprofile *profileList = NULL;
DOPtable *dopTableList = NULL;
double startTime;
/* loop through all the models */
for (; model != NULL; model = model->NUMDnextModel) {
if (!model->NUMDpInfo) {
TSCALLOC(model->NUMDpInfo, 1, ONEtranInfo);
}
methods = model->NUMDmethods;
if (!methods) {
TSCALLOC(methods, 1, METHcard);
model->NUMDmethods = methods;
}
models = model->NUMDmodels;
if (!models) {
TSCALLOC(models, 1, MODLcard);
model->NUMDmodels = models;
}
options = model->NUMDoptions;
if (!options) {
TSCALLOC(options, 1, OPTNcard);
model->NUMDoptions = options;
}
outputs = model->NUMDoutputs;
if (!outputs) {
TSCALLOC(outputs, 1, OUTPcard);
model->NUMDoutputs = outputs;
}
if (!methods->METHvoltPredGiven) {
methods->METHvoltPred = FALSE;
}
if (!methods->METHmobDerivGiven) {
methods->METHmobDeriv = TRUE;
}
if (!methods->METHoneCarrierGiven) {
methods->METHoneCarrier = FALSE;
}
if (!methods->METHacAnalysisMethodGiven) {
methods->METHacAnalysisMethod = SOR;
}
if (!methods->METHdabstolGiven) {
methods->METHdabstol = DABSTOL1D;
}
if (!methods->METHdreltolGiven) {
methods->METHdreltol = ckt->CKTreltol;
}
if (!methods->METHitLimGiven) {
methods->METHitLim = 20;
}
if (!methods->METHomegaGiven || methods->METHomega <= 0.0) {
methods->METHomega = 2.0 * M_PI /* radians/sec */ ;
}
if (!options->OPTNdefaGiven || options->OPTNdefa <= 0.0) {
options->OPTNdefa = 1.0e4 /* cm^2 */ ;
}
if (!options->OPTNdeviceTypeGiven) {
options->OPTNdeviceType = OPTN_DIODE;
}
if (!options->OPTNicFileGiven) {
options->OPTNicFile = NULL;
options->OPTNunique = FALSE; /* Can't form a unique name. */
}
if (!options->OPTNuniqueGiven) {
options->OPTNunique = FALSE;
}
/* Set up the rest of the card lists */
if ((error = MODLsetup(model->NUMDmodels)) != 0)
return (error);
BandGapNarrowing = models->MODLbandGapNarrowing;
ConcDepLifetime = models->MODLconcDepLifetime;
TempDepMobility = models->MODLtempDepMobility;
ConcDepMobility = models->MODLconcDepMobility;
if ((error = OUTPsetup(model->NUMDoutputs)) != 0)
return (error);
if ((error = MATLsetup(model->NUMDmaterials, &materialList)) != 0)
return (error);
if ((error = MOBsetup(model->NUMDmobility, materialList)) != 0)
return (error);
if ((error = MESHsetup('x', model->NUMDxMeshes, &xCoordList, &xMeshSize)) != 0)
return (error);
if ((error = DOMNsetup(model->NUMDdomains, &domainList,
xCoordList, NULL, materialList)) != 0)
return (error);
if ((error = BDRYsetup(model->NUMDboundaries,
xCoordList, NULL, domainList)) != 0)
return (error);
if ((error = CONTsetup(model->NUMDcontacts, NULL)) != 0)
return (error);
if ((error = DOPsetup(model->NUMDdopings, &profileList,
&dopTableList, xCoordList, NULL)) != 0)
return (error);
model->NUMDmatlInfo = materialList;
model->NUMDprofiles = profileList;
model->NUMDdopTables = dopTableList;
/* loop through all the instances of the model */
for (inst = model->NUMDinstances; inst != NULL;
inst = inst->NUMDnextInstance) {
startTime = SPfrontEnd->IFseconds();
if ((!inst->NUMDprintGiven)) {
inst->NUMDprint = 0;
} else if (inst->NUMDprint <= 0) {
inst->NUMDprint = 1;
}
if ((!inst->NUMDicFileGiven)) {
if (options->OPTNunique) {
inst->NUMDicFile = tprintf("%s.%s", options->OPTNicFile, inst->NUMDname);
} else if (options->OPTNicFile != NULL) {
inst->NUMDicFile = tprintf("%s", options->OPTNicFile);
} else {
inst->NUMDicFile = NULL;
}
}
inst->NUMDstate = *states;
*states += NUMDnumStates;
if (!inst->NUMDpDevice) {
/* Assign the mesh info to each instance. */
TSCALLOC(pDevice, 1, ONEdevice);
TSCALLOC(pDevice->pStats, 1, ONEstats);
pDevice->name = inst->NUMDname;
pDevice->solverType = SLV_NONE;
pDevice->numNodes = xMeshSize;
pDevice->abstol = methods->METHdabstol;
pDevice->reltol = methods->METHdreltol;
pDevice->rhsImag = NULL;
TSCALLOC(pDevice->elemArray, pDevice->numNodes, ONEelem *);
/* Create a copy of material data that can change with temperature. */
pDevice->pMaterials = NULL;
for (pM = materialList; pM != NULL; pM = pM->next) {
if (pDevice->pMaterials == NULL) {
TSCALLOC(pMaterial, 1, ONEmaterial);
pDevice->pMaterials = pMaterial;
} else {
TSCALLOC(pMaterial->next, 1, ONEmaterial);
pMaterial = pMaterial->next;
}
/* Copy everything, then fix the incorrect pointer. */
memcpy(pMaterial, pM, sizeof(ONEmaterial));
pMaterial->next = NULL;
}
/* generate the mesh structure for the device */
ONEbuildMesh(pDevice, xCoordList, domainList, pDevice->pMaterials);
/* store the device info in the instance */
inst->NUMDpDevice = pDevice;
}
/* Now update the state pointers. */
ONEgetStatePointers(inst->NUMDpDevice, states);
/* Wipe out statistics from previous runs (if any). */
memset(inst->NUMDpDevice->pStats, 0, sizeof(ONEstats));
inst->NUMDpDevice->pStats->totalTime[STAT_SETUP] +=
SPfrontEnd->IFseconds() - startTime;
/* macro to make elements with built in test for out of memory */
#define TSTALLOC(ptr,first,second) \
do { if ((inst->ptr = SMPmakeElt(matrix, inst->first, inst->second)) == NULL){\
return(E_NOMEM);\
} } while(0)
TSTALLOC(NUMDposPosPtr, NUMDposNode, NUMDposNode);
TSTALLOC(NUMDnegNegPtr, NUMDnegNode, NUMDnegNode);
TSTALLOC(NUMDnegPosPtr, NUMDnegNode, NUMDposNode);
TSTALLOC(NUMDposNegPtr, NUMDposNode, NUMDnegNode);
}
/* Clean up lists */
killCoordInfo(xCoordList);
killDomainInfo(domainList);
}
int
NBJT2setup(SMPmatrix *matrix, GENmodel *inModel, CKTcircuit *ckt, int *states)
/*
* load the structure with those pointers needed later for fast matrix
* loading
*/
{
register NBJT2model *model = (NBJT2model *) inModel;
register NBJT2instance *inst;
METHcard *methods;
MODLcard *models;
OPTNcard *options;
OUTPcard *outputs;
int error, xIndex;
int xMeshSize, yMeshSize;
TWOdevice *pDevice;
TWOcoord *xCoordList = NULL;
TWOcoord *yCoordList = NULL;
TWOdomain *domainList = NULL;
TWOelectrode *electrodeList = NULL;
TWOmaterial *pM, *pMaterial = NULL, *materialList = NULL;
DOPprofile *profileList = NULL;
DOPtable *dopTableList = NULL;
double startTime;
/* loop through all the models */
for (; model != NULL; model = model->NBJT2nextModel) {
if (!model->NBJT2pInfo) {
TSCALLOC(model->NBJT2pInfo, 1, TWOtranInfo);
}
methods = model->NBJT2methods;
if (!methods) {
TSCALLOC(methods, 1, METHcard);
model->NBJT2methods = methods;
}
models = model->NBJT2models;
if (!models) {
TSCALLOC(models, 1, MODLcard);
model->NBJT2models = models;
}
options = model->NBJT2options;
if (!options) {
TSCALLOC(options, 1, OPTNcard);
model->NBJT2options = options;
}
outputs = model->NBJT2outputs;
if (!outputs) {
TSCALLOC(outputs, 1, OUTPcard);
model->NBJT2outputs = outputs;
}
if (!methods->METHvoltPredGiven) {
methods->METHvoltPred = FALSE;
}
if (!methods->METHmobDerivGiven) {
methods->METHmobDeriv = TRUE;
}
if (!methods->METHoneCarrierGiven) {
methods->METHoneCarrier = FALSE;
}
if (!methods->METHacAnalysisMethodGiven) {
methods->METHacAnalysisMethod = SOR;
}
if (!methods->METHdabstolGiven) {
methods->METHdabstol = DABSTOL2D;
}
if (!methods->METHdreltolGiven) {
methods->METHdreltol = ckt->CKTreltol;
}
if (!methods->METHitLimGiven) {
methods->METHitLim = 50;
}
if (!methods->METHomegaGiven || methods->METHomega <= 0.0) {
methods->METHomega = 2.0 * M_PI /* radians/sec */ ;
}
if (!options->OPTNdefaGiven || options->OPTNdefa <= 0.0) {
options->OPTNdefa = 1.0e4 /* cm^2 */ ;
}
if (!options->OPTNdeflGiven || options->OPTNdefl <= 0.0) {
options->OPTNdefl = 1.0e2 /* cm */ ;
}
if (!options->OPTNdefwGiven && options->OPTNdefaGiven) {
options->OPTNdefw = options->OPTNdefa / options->OPTNdefl;
} else if (!options->OPTNdefwGiven || options->OPTNdefw <= 0.0) {
options->OPTNdefw = 1.0e2 /* cm */ ;
}
if (!options->OPTNdeviceTypeGiven) {
options->OPTNdeviceType = OPTN_BIPOLAR;
}
if (!options->OPTNicFileGiven) {
options->OPTNicFile = NULL;
options->OPTNunique = FALSE; /* Can't form a unique name. */
}
if (!options->OPTNuniqueGiven) {
options->OPTNunique = FALSE;
}
OneCarrier = methods->METHoneCarrier;
/* Set up the rest of the card lists */
if ((error = MODLsetup(model->NBJT2models)) != 0)
return (error);
BandGapNarrowing = models->MODLbandGapNarrowing;
ConcDepLifetime = models->MODLconcDepLifetime;
TempDepMobility = models->MODLtempDepMobility;
ConcDepMobility = models->MODLconcDepMobility;
SurfaceMobility = models->MODLsurfaceMobility;
if ((error = OUTPsetup(model->NBJT2outputs)) != 0)
return (error);
if ((error = MATLsetup(model->NBJT2materials, &materialList)) != 0)
return (error);
if ((error = MOBsetup(model->NBJT2mobility, materialList)) != 0)
return (error);
if ((error = MESHsetup('x', model->NBJT2xMeshes, &xCoordList, &xMeshSize)) != 0)
return (error);
if ((error = MESHsetup('y', model->NBJT2yMeshes, &yCoordList, &yMeshSize)) != 0)
return (error);
if ((error = DOMNsetup(model->NBJT2domains, &domainList,
xCoordList, yCoordList, materialList)) != 0)
return (error);
if ((error = BDRYsetup(model->NBJT2boundaries,
xCoordList, yCoordList, domainList)) != 0)
return (error);
if ((error = ELCTsetup(model->NBJT2electrodes, &electrodeList,
xCoordList, yCoordList)) != 0)
return (error);
/* Make sure electrodes are OK. */
checkElectrodes(electrodeList, 3); /* NBJT2 has 3 electrodes */
if ((error = CONTsetup(model->NBJT2contacts, electrodeList)) != 0)
return (error);
if ((error = DOPsetup(model->NBJT2dopings, &profileList,
&dopTableList, xCoordList, yCoordList)) != 0)
return (error);
model->NBJT2matlInfo = materialList;
model->NBJT2profiles = profileList;
model->NBJT2dopTables = dopTableList;
/* loop through all the instances of the model */
for (inst = model->NBJT2instances; inst != NULL;
inst = inst->NBJT2nextInstance) {
startTime = SPfrontEnd->IFseconds();
if (!inst->NBJT2printGiven) {
inst->NBJT2print = 0;
} else if (inst->NBJT2print <= 0) {
inst->NBJT2print = 1;
}
if (!inst->NBJT2icFileGiven) {
if (options->OPTNunique) {
inst->NBJT2icFile = tprintf("%s.%s", options->OPTNicFile, inst->NBJT2name);
} else if (options->OPTNicFile != NULL) {
inst->NBJT2icFile = tprintf("%s", options->OPTNicFile);
} else {
inst->NBJT2icFile = NULL;
}
}
inst->NBJT2state = *states;
*states += NBJT2numStates;
if (!inst->NBJT2pDevice) {
/* Assign the mesh and profile info to each instance. */
TSCALLOC(pDevice, 1, TWOdevice);
TSCALLOC(pDevice->pStats, 1, TWOstats);
pDevice->name = inst->NBJT2name;
pDevice->solverType = SLV_NONE;
pDevice->numXNodes = xMeshSize;
pDevice->numYNodes = yMeshSize;
pDevice->xScale = MESHmkArray(xCoordList, xMeshSize);
pDevice->yScale = MESHmkArray(yCoordList, yMeshSize);
pDevice->abstol = methods->METHdabstol;
pDevice->reltol = methods->METHdreltol;
pDevice->rhsImag = NULL;
TSCALLOC(pDevice->elemArray, pDevice->numXNodes, TWOelem **);
for (xIndex = 1; xIndex < pDevice->numXNodes; xIndex++) {
TSCALLOC(pDevice->elemArray[xIndex], pDevice->numYNodes, TWOelem *);
}
/* Create a copy of material data that can change with temperature. */
pDevice->pMaterials = NULL;
for (pM = materialList; pM != NULL; pM = pM->next) {
if (pDevice->pMaterials == NULL) {
TSCALLOC(pMaterial, 1, TWOmaterial);
pDevice->pMaterials = pMaterial;
} else {
TSCALLOC(pMaterial->next, 1, TWOmaterial);
pMaterial = pMaterial->next;
}
/* Copy everything, then fix the incorrect pointer. */
bcopy(pM, pMaterial, sizeof(TWOmaterial));
pMaterial->next = NULL;
}
/* Generate the mesh structure for the device. */
TWObuildMesh(pDevice, domainList, electrodeList, pDevice->pMaterials);
/* Store the device info in the instance. */
inst->NBJT2pDevice = pDevice;
}
/* Now update the state pointers. */
TWOgetStatePointers(inst->NBJT2pDevice, states);
/* Wipe out statistics from previous runs (if any). */
bzero(inst->NBJT2pDevice->pStats, sizeof(TWOstats));
inst->NBJT2pDevice->pStats->totalTime[STAT_SETUP] +=
SPfrontEnd->IFseconds() - startTime;
/* macro to make elements with built in test for out of memory */
#define TSTALLOC(ptr,first,second) \
do { if ((inst->ptr = SMPmakeElt(matrix, inst->first, inst->second)) == NULL){\
return(E_NOMEM);\
} } while(0)
TSTALLOC(NBJT2colColPtr, NBJT2colNode, NBJT2colNode);
TSTALLOC(NBJT2colBasePtr, NBJT2colNode, NBJT2baseNode);
TSTALLOC(NBJT2colEmitPtr, NBJT2colNode, NBJT2emitNode);
TSTALLOC(NBJT2baseColPtr, NBJT2baseNode, NBJT2colNode);
TSTALLOC(NBJT2baseBasePtr, NBJT2baseNode, NBJT2baseNode);
TSTALLOC(NBJT2baseEmitPtr, NBJT2baseNode, NBJT2emitNode);
TSTALLOC(NBJT2emitColPtr, NBJT2emitNode, NBJT2colNode);
TSTALLOC(NBJT2emitBasePtr, NBJT2emitNode, NBJT2baseNode);
TSTALLOC(NBJT2emitEmitPtr, NBJT2emitNode, NBJT2emitNode);
}
/* Clean up lists */
killCoordInfo(xCoordList);
killCoordInfo(yCoordList);
killDomainInfo(domainList);
killElectrodeInfo(electrodeList);
}
