epilectrik
/
pyNgSpice
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity
Browse Source

bsim3v3.24: add OMP support 

model is still used by commercial (e.g. X-fab) libraries
pre-master-46
h_vogt 10 years ago
committed by rlar
parent
151fed6084
commit
4047a794e7
4 changed files with 248 additions and 0 deletions
Whitespace
Show all changes Ignore whitespace when comparing lines Ignore changes in amount of whitespace Ignore changes in whitespace at EOL
Split View
Diff Options
Show Stats Download Patch File Download Diff File
  1. 4
      src/spicelib/devices/bsim3v32/b3v32dest.c
  2. 159
      src/spicelib/devices/bsim3v32/b3v32ld.c
  3. 38
      src/spicelib/devices/bsim3v32/b3v32set.c
  4. 47
      src/spicelib/devices/bsim3v32/bsim3v32def.h

4
src/spicelib/devices/bsim3v32/b3v32dest.c
View File

@ -48,6 +48,10 @@ BSIM3v32destroy (GENmodel **inModel)
if(prev) FREE(prev);
}
if(oldmod) {
#ifdef USE_OMP
/* free just once for all models */
FREE(oldmod->BSIM3v32InstanceArray);
#endif
FREE(oldmod->BSIM3v32version);
FREE(oldmod);
}

159
src/spicelib/devices/bsim3v32/b3v32ld.c
View File

@ -31,12 +31,41 @@
#define DELTA_3 0.02
#define DELTA_4 0.02
#ifdef USE_OMP
int BSIM3v32LoadOMP(BSIM3v32instance *here, CKTcircuit *ckt);
void BSIM3v32LoadRhsMat(GENmodel *inModel, CKTcircuit *ckt);
#endif
int
BSIM3v32load (GENmodel *inModel, CKTcircuit *ckt)
{
#ifdef USE_OMP
int idx;
BSIM3v32model *model = (BSIM3v32model*)inModel;
int error = 0;
BSIM3v32instance **InstArray;
InstArray = model->BSIM3v32InstanceArray;
#pragma omp parallel for
for (idx = 0; idx < model->BSIM3v32InstCount; idx++) {
BSIM3v32instance *here = InstArray[idx];
int local_error = BSIM3v32LoadOMP(here, ckt);
if (local_error)
error = local_error;
}
BSIM3v32LoadRhsMat(inModel, ckt);
return error;
}
int BSIM3v32LoadOMP(BSIM3v32instance *here, CKTcircuit *ckt) {
BSIM3v32model *model = here->BSIM3v32modPtr;
#else
BSIM3v32model *model = (BSIM3v32model*)inModel;
BSIM3v32instance *here;
#endif
double SourceSatCurrent, DrainSatCurrent;
double ag0, qgd, qgs, qgb, von, cbhat, VgstNVt, ExpVgst;
double cdrain, cdhat, cdreq, ceqbd, ceqbs, ceqqb, ceqqd, ceqqg, ceq, geq;
@ -141,10 +170,12 @@ ChargeComputationNeeded =
((ckt->CKTmode & (MODEDCTRANCURVE | MODEAC | MODETRAN | MODEINITSMSIG)) ||
((ckt->CKTmode & MODETRANOP) && (ckt->CKTmode & MODEUIC)))
? 1 : 0;
#ifndef USE_OMP
for (; model != NULL; model = model->BSIM3v32nextModel)
{ for (here = model->BSIM3v32instances; here != NULL;
here = here->BSIM3v32nextInstance)
{
#endif
Check = 1;
ByPass = 0;
pParam = here->pParam;
@ -3255,6 +3286,15 @@ line900:
m = here->BSIM3v32m;
#ifdef USE_OMP
here->BSIM3v32rhsG = m * ceqqg;
here->BSIM3v32rhsB = m * (ceqbs + ceqbd + ceqqb);
here->BSIM3v32rhsD = m * (ceqbd - cdreq - ceqqd);
here->BSIM3v32rhsS = m * (cdreq + ceqbs + ceqqg
+ ceqqb + ceqqd);
if (here->BSIM3v32nqsMod)
here->BSIM3v32rhsQ = m * (cqcheq - cqdef);
#else
(*(ckt->CKTrhs + here->BSIM3v32gNode) -= m * ceqqg);
(*(ckt->CKTrhs + here->BSIM3v32bNode) -= m * (ceqbs + ceqbd + ceqqb));
(*(ckt->CKTrhs + here->BSIM3v32dNodePrime) += m * (ceqbd - cdreq - ceqqd));
@ -3262,12 +3302,66 @@ line900:
+ ceqqb + ceqqd));
if (here->BSIM3v32nqsMod)
*(ckt->CKTrhs + here->BSIM3v32qNode) += m * (cqcheq - cqdef);
#endif
/*
* load y matrix
*/
T1 = qdef * here->BSIM3v32gtau;
#ifdef USE_OMP
here->BSIM3v32DdPt = m * here->BSIM3v32drainConductance;
here->BSIM3v32GgPt = m * (gcggb - ggtg);
here->BSIM3v32SsPt = m * here->BSIM3v32sourceConductance;
here->BSIM3v32BbPt = m * (here->BSIM3v32gbd + here->BSIM3v32gbs
- gcbgb - gcbdb - gcbsb - here->BSIM3v32gbbs);
here->BSIM3v32DPdpPt = m * (here->BSIM3v32drainConductance
+ here->BSIM3v32gds + here->BSIM3v32gbd
+ RevSum + gcddb + dxpart * ggtd
+ T1 * ddxpart_dVd + gbdpdp);
here->BSIM3v32SPspPt = m * (here->BSIM3v32sourceConductance
+ here->BSIM3v32gds + here->BSIM3v32gbs
+ FwdSum + gcssb + sxpart * ggts
+ T1 * dsxpart_dVs + gbspsp);
here->BSIM3v32DdpPt = m * here->BSIM3v32drainConductance;
here->BSIM3v32GbPt = m * (gcggb + gcgdb + gcgsb + ggtb);
here->BSIM3v32GdpPt = m * (gcgdb - ggtd);
here->BSIM3v32GspPt = m * (gcgsb - ggts);
here->BSIM3v32SspPt = m * here->BSIM3v32sourceConductance;
here->BSIM3v32BgPt = m * (gcbgb - here->BSIM3v32gbgs);
here->BSIM3v32BdpPt = m * (gcbdb - here->BSIM3v32gbd + gbbdp);
here->BSIM3v32BspPt = m * (gcbsb - here->BSIM3v32gbs + gbbsp);
here->BSIM3v32DPdPt = m * here->BSIM3v32drainConductance;
here->BSIM3v32DPgPt = m * (Gm + gcdgb + dxpart * ggtg
+ T1 * ddxpart_dVg + gbdpg);
here->BSIM3v32DPbPt = m * (here->BSIM3v32gbd - Gmbs + gcdgb + gcddb
+ gcdsb - dxpart * ggtb
- T1 * ddxpart_dVb - gbdpb);
here->BSIM3v32DPspPt = m * (here->BSIM3v32gds + FwdSum - gcdsb
- dxpart * ggts - T1 * ddxpart_dVs - gbdpsp);
here->BSIM3v32SPgPt = m * (gcsgb - Gm + sxpart * ggtg
+ T1 * dsxpart_dVg + gbspg);
here->BSIM3v32SPsPt = m * here->BSIM3v32sourceConductance;
here->BSIM3v32SPbPt = m * (here->BSIM3v32gbs + Gmbs + gcsgb + gcsdb
+ gcssb - sxpart * ggtb
- T1 * dsxpart_dVb - gbspb);
here->BSIM3v32SPdpPt = m * (here->BSIM3v32gds + RevSum - gcsdb
- sxpart * ggtd - T1 * dsxpart_dVd - gbspdp);
if (here->BSIM3v32nqsMod)
{
here->BSIM3v32QqPt = m * (gqdef + here->BSIM3v32gtau);
here->BSIM3v32DPqPt = m * (dxpart * here->BSIM3v32gtau);
here->BSIM3v32SPqPt = m * (sxpart * here->BSIM3v32gtau);
here->BSIM3v32GqPt = m * here->BSIM3v32gtau;
here->BSIM3v32QgPt = m * (ggtg - gcqgb);
here->BSIM3v32QdpPt = m * (ggtd - gcqdb);
here->BSIM3v32QspPt = m * (ggts - gcqsb);
here->BSIM3v32QbPt = m * (ggtb - gcqbb);
}
#else
(*(here->BSIM3v32DdPtr) += m * here->BSIM3v32drainConductance);
(*(here->BSIM3v32GgPtr) += m * (gcggb - ggtg));
(*(here->BSIM3v32SsPtr) += m * here->BSIM3v32sourceConductance);
@ -3319,12 +3413,77 @@ line900:
*(here->BSIM3v32QspPtr) += m * (ggts - gcqsb);
*(here->BSIM3v32QbPtr) += m * (ggtb - gcqbb);
}
#endif
line1000: ;
#ifndef USE_OMP
} /* End of Mosfet Instance */
} /* End of Model Instance */
#endif
return(OK);
}
#ifdef USE_OMP
void BSIM3v32LoadRhsMat(GENmodel *inModel, CKTcircuit *ckt)
{
int InstCount, idx;
BSIM3v32instance **InstArray;
BSIM3v32instance *here;
BSIM3v32model *model = (BSIM3v32model*)inModel;
InstArray = model->BSIM3v32InstanceArray;
InstCount = model->BSIM3v32InstCount;
for (idx = 0; idx < InstCount; idx++) {
here = InstArray[idx];
/* Update b for Ax = b */
(*(ckt->CKTrhs + here->BSIM3v32gNode) -= here->BSIM3v32rhsG);
(*(ckt->CKTrhs + here->BSIM3v32bNode) -= here->BSIM3v32rhsB);
(*(ckt->CKTrhs + here->BSIM3v32dNodePrime) += here->BSIM3v32rhsD);
(*(ckt->CKTrhs + here->BSIM3v32sNodePrime) += here->BSIM3v32rhsS);
if (here->BSIM3v32nqsMod)
(*(ckt->CKTrhs + here->BSIM3v32qNode) += here->BSIM3v32rhsQ);
/* Update A for Ax = b */
(*(here->BSIM3v32DdPtr) += here->BSIM3v32DdPt);
(*(here->BSIM3v32GgPtr) += here->BSIM3v32GgPt);
(*(here->BSIM3v32SsPtr) += here->BSIM3v32SsPt);
(*(here->BSIM3v32BbPtr) += here->BSIM3v32BbPt);
(*(here->BSIM3v32DPdpPtr) += here->BSIM3v32DPdpPt);
(*(here->BSIM3v32SPspPtr) += here->BSIM3v32SPspPt);
(*(here->BSIM3v32DdpPtr) -= here->BSIM3v32DdpPt);
(*(here->BSIM3v32GbPtr) -= here->BSIM3v32GbPt);
(*(here->BSIM3v32GdpPtr) += here->BSIM3v32GdpPt);
(*(here->BSIM3v32GspPtr) += here->BSIM3v32GspPt);
(*(here->BSIM3v32SspPtr) -= here->BSIM3v32SspPt);
(*(here->BSIM3v32BgPtr) += here->BSIM3v32BgPt);
(*(here->BSIM3v32BdpPtr) += here->BSIM3v32BdpPt);
(*(here->BSIM3v32BspPtr) += here->BSIM3v32BspPt);
(*(here->BSIM3v32DPdPtr) -= here->BSIM3v32DPdPt);
(*(here->BSIM3v32DPgPtr) += here->BSIM3v32DPgPt);
(*(here->BSIM3v32DPbPtr) -= here->BSIM3v32DPbPt);
(*(here->BSIM3v32DPspPtr) -= here->BSIM3v32DPspPt);
(*(here->BSIM3v32SPgPtr) += here->BSIM3v32SPgPt);
(*(here->BSIM3v32SPsPtr) -= here->BSIM3v32SPsPt);
(*(here->BSIM3v32SPbPtr) -= here->BSIM3v32SPbPt);
(*(here->BSIM3v32SPdpPtr) -= here->BSIM3v32SPdpPt);
if (here->BSIM3v32nqsMod)
{
*(here->BSIM3v32QqPtr) += here->BSIM3v32QqPt;
*(here->BSIM3v32DPqPtr) += here->BSIM3v32DPqPt;
*(here->BSIM3v32SPqPtr) += here->BSIM3v32SPqPt;
*(here->BSIM3v32GqPtr) -= here->BSIM3v32GqPt;
*(here->BSIM3v32QgPtr) += here->BSIM3v32QgPt;
*(here->BSIM3v32QdpPtr) += here->BSIM3v32QdpPt;
*(here->BSIM3v32QspPtr) += here->BSIM3v32QspPt;
*(here->BSIM3v32QbPtr) += here->BSIM3v32QbPt;
}
}
}
#endif

38
src/spicelib/devices/bsim3v32/b3v32set.c
View File

@ -38,6 +38,11 @@ CKTnode *tmp;
CKTnode *tmpNode;
IFuid tmpName;
#ifdef USE_OMP
int idx, InstCount;
BSIM3v32instance **InstArray;
#endif
/* loop through all the BSIM3v32 device models */
for( ; model != NULL; model = model->BSIM3v32nextModel )
{
@ -1091,6 +1096,39 @@ do { if((here->ptr = SMPmakeElt(matrix, here->first, here->second)) == NULL){\
}
}
#ifdef USE_OMP
InstCount = 0;
model = (BSIM3v32model*)inModel;
/* loop through all the BSIM3 device models
to count the number of instances */
for (; model != NULL; model = model->BSIM3v32nextModel)
{
/* loop through all the instances of the model */
for (here = model->BSIM3v32instances; here != NULL;
here = here->BSIM3v32nextInstance)
{
InstCount++;
}
}
InstArray = TMALLOC(BSIM3v32instance*, InstCount);
model = (BSIM3v32model*)inModel;
idx = 0;
for (; model != NULL; model = model->BSIM3v32nextModel)
{
/* loop through all the instances of the model */
for (here = model->BSIM3v32instances; here != NULL;
here = here->BSIM3v32nextInstance)
{
InstArray[idx] = here;
idx++;
}
/* set the array pointer and instance count into each model */
model->BSIM3v32InstCount = InstCount;
model->BSIM3v32InstanceArray = InstArray;
}
#endif
return(OK);
}

47
src/spicelib/devices/bsim3v32/bsim3v32def.h
View File

@ -171,6 +171,48 @@ typedef struct sBSIM3v32instance
double *BSIM3v32SPqPtr;
double *BSIM3v32BqPtr;
#ifdef USE_OMP
/* per instance storage of results, to update matrix at a later stge */
double BSIM3v32rhsG;
double BSIM3v32rhsB;
double BSIM3v32rhsD;
double BSIM3v32rhsS;
double BSIM3v32rhsQ;
double BSIM3v32DdPt;
double BSIM3v32GgPt;
double BSIM3v32SsPt;
double BSIM3v32BbPt;
double BSIM3v32DPdpPt;
double BSIM3v32SPspPt;
double BSIM3v32DdpPt;
double BSIM3v32GbPt;
double BSIM3v32GdpPt;
double BSIM3v32GspPt;
double BSIM3v32SspPt;
double BSIM3v32BdpPt;
double BSIM3v32BspPt;
double BSIM3v32DPspPt;
double BSIM3v32DPdPt;
double BSIM3v32BgPt;
double BSIM3v32DPgPt;
double BSIM3v32SPgPt;
double BSIM3v32SPsPt;
double BSIM3v32DPbPt;
double BSIM3v32SPbPt;
double BSIM3v32SPdpPt;
double BSIM3v32QqPt;
double BSIM3v32QdpPt;
double BSIM3v32QgPt;
double BSIM3v32QspPt;
double BSIM3v32QbPt;
double BSIM3v32DPqPt;
double BSIM3v32GqPt;
double BSIM3v32SPqPt;
double BSIM3v32BqPt;
#endif
#define BSIM3v32vbd BSIM3v32states+ 0
#define BSIM3v32vbs BSIM3v32states+ 1
#define BSIM3v32vgs BSIM3v32states+ 2
@ -836,6 +878,11 @@ typedef struct sBSIM3v32model
struct bsim3v32SizeDependParam *pSizeDependParamKnot;
#ifdef USE_OMP
int BSIM3v32InstCount;
struct sBSIM3v32instance **BSIM3v32InstanceArray;
#endif
/* Flags */
unsigned BSIM3v32mobModGiven :1;
unsigned BSIM3v32binUnitGiven :1;

Write Preview
Loading…
Cancel
Save