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

bsim6.va, macro for lexp to avoid derivatives 

fixme, adsmXml can't do
  a ? b : c ? d : e
needs parentheses
  a ? b : (c ? d : e)
bug in parser
This commit is contained in:
rlar 2017-05-26 20:24:32 +02:00 committed by Holger Vogt
parent 37c5813d68
commit 4f46cd20ee
1 changed files with 55 additions and 65 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

120
src/spicelib/devices/adms/bsim6/admsva/bsim6.va
View File

@ -624,19 +624,9 @@ electrical d, g, s, b;
`endif
// Clamped exponential function
analog function real lexp;
input x;
real x;
begin
if(x > `EXPL_THRESHOLD) begin
lexp = `MAX_EXPL * (1.0+(x)-`EXPL_THRESHOLD);
end else if(x < -`EXPL_THRESHOLD) begin
lexp = `MIN_EXPL;
end else begin
lexp = exp(x);
end
end
endfunction
`define lexp(x) \
(((x) > `EXPL_THRESHOLD) ? (`MAX_EXPL * (1.0+(x)-`EXPL_THRESHOLD)) : \
(((x) < -`EXPL_THRESHOLD) ? (`MIN_EXPL) : exp(x)))
// Clamped log function
analog function real lln;
@ -2747,7 +2737,7 @@ DevTemp = $temperature + DTEMP;
Eg0 = BG0SUB - TBGASUB * Tnom * Tnom / (Tnom + TBGBSUB);
T1y = (DevTemp / Tnom) * sqrt(DevTemp / Tnom);
ni = NI0SUB * T1y * lexp(Eg / (2.0 * Vtm0) - Eg / (2.0 * Vtm));
ni = NI0SUB * T1y * `lexp(Eg / (2.0 * Vtm0) - Eg / (2.0 * Vtm));
`ifdef __SHMOD__
if ((SHMOD != 0) && (RTH0 > 0)) begin
T0 = lln(NDEP_i/ni);
@ -2832,7 +2822,7 @@ DevTemp = $temperature + DTEMP;
BETA0_t = BETA0_i * pow(TRatio, IIT_i);
BGIDL_t = BGIDL_i * `hypsmooth(1.0 + TGIDL_i * delTemp - 1.0E-6, 1.0E-3);
BGISL_t = BGISL_i * `hypsmooth(1.0 + TGIDL_i * delTemp - 1.0E-6, 1.0E-3);
igtemp = lexp(IGT_i * lln(TRatio)); //pow(TRatio, IGT_i);
igtemp = `lexp(IGT_i * lln(TRatio)); //pow(TRatio, IGT_i);
//***Diode Model temperature Code Start***//
CJS_t = CJS * `hypsmooth(1.0 + TCJ * delTemp - 1.0E-6, 1.0E-3);
@ -2850,20 +2840,20 @@ DevTemp = $temperature + DTEMP;
T0y = Eg0 / Vtm0 - Eg / Vtm;
T1y = lln(TRatio);
T3y = lexp((T0y + XTIS * T1y) / NJS);
T3y = `lexp((T0y + XTIS * T1y) / NJS);
JSS_t = JSS * T3y;
JSWS_t = JSWS * T3y;
JSWGS_t = JSWGS * T3y;
T3y = lexp((T0y + XTID * T1y) / NJD);
T3y = `lexp((T0y + XTID * T1y) / NJD);
JSD_t = JSD * T3y;
JSWD_t = JSWD * T3y;
JSWGD_t = JSWGD * T3y;
JTSS_t = JTSS * lexp(Eg0 * XTSS * (TRatio - 1.0) / Vtm);
JTSSWS_t = JTSSWS * lexp(Eg0 * XTSSWS * (TRatio - 1.0) / Vtm);
JTSSWGS_t = JTSSWGS * (sqrt(JTWEFF / Weffcj) + 1.0) * lexp(Eg0 * XTSSWGS * (TRatio - 1) / Vtm);
JTSD_t = JTSD * lexp(Eg0 * XTSD * (TRatio - 1.0) / Vtm);
JTSSWD_t = JTSSWD * lexp(Eg0 * XTSSWD * (TRatio - 1.0) / Vtm);
JTSSWGD_t = JTSSWGD * (sqrt(JTWEFF / Weffcj) + 1.0) * lexp(Eg0 * XTSSWGD * (TRatio - 1) / Vtm);
JTSS_t = JTSS * `lexp(Eg0 * XTSS * (TRatio - 1.0) / Vtm);
JTSSWS_t = JTSSWS * `lexp(Eg0 * XTSSWS * (TRatio - 1.0) / Vtm);
JTSSWGS_t = JTSSWGS * (sqrt(JTWEFF / Weffcj) + 1.0) * `lexp(Eg0 * XTSSWGS * (TRatio - 1) / Vtm);
JTSD_t = JTSD * `lexp(Eg0 * XTSD * (TRatio - 1.0) / Vtm);
JTSSWD_t = JTSSWD * `lexp(Eg0 * XTSSWD * (TRatio - 1.0) / Vtm);
JTSSWGD_t = JTSSWGD * (sqrt(JTWEFF / Weffcj) + 1.0) * `lexp(Eg0 * XTSSWGD * (TRatio - 1) / Vtm);
//All NJT*'s smoothed to 0.01 to prevent divide by zero / negative values
NJTS_t = `hypsmooth(NJTS * (1.0 + TNJTS * (TRatio-1.0)) - 0.01, 1.0E-3) + 0.01;
@ -2925,16 +2915,16 @@ DevTemp = $temperature + DTEMP;
Isbs = ASeff * JSS_t + PSeff * JSWS_t + Weffcj * NF * JSWGS_t;
if(Isbs > 0.0) begin
Nvtms = Vtm * NJS;
XExpBVS = lexp(-BVS / Nvtms) * XJBVS;
XExpBVS = `lexp(-BVS / Nvtms) * XJBVS;
T2y = max(IJTHSFWD / Isbs, 10.0);
Tb = 1.0 + T2y - XExpBVS;
VjsmFwd = Nvtms * lln(0.5 * (Tb + sqrt(Tb * Tb + 4 * XExpBVS)));
T0y = lexp(VjsmFwd / Nvtms);
T0y = `lexp(VjsmFwd / Nvtms);
IVjsmFwd = Isbs * (T0y - XExpBVS / T0y + XExpBVS - 1.0);
SslpFwd = Isbs * (T0y + XExpBVS / T0y) / Nvtms;
T2y = `hypsmooth(IJTHSREV / Isbs - 10.0, 1.0E-3) + 10.0;
VjsmRev = -BVS - Nvtms * lln((T2y - 1.0) / XJBVS);
T1y = XJBVS * lexp(-(BVS + VjsmRev) / Nvtms);
T1y = XJBVS * `lexp(-(BVS + VjsmRev) / Nvtms);
IVjsmRev = Isbs * (1.0 + T1y);
SslpRev = -Isbs * T1y / Nvtms;
end else begin
@ -2952,16 +2942,16 @@ DevTemp = $temperature + DTEMP;
Isbd = ADeff * JSD_t + PDeff * JSWD_t + Weffcj * NF * JSWGD_t;
if(Isbd > 0.0) begin
Nvtmd = Vtm * NJD;
XExpBVD = lexp(-BVD / Nvtmd) * XJBVD;
XExpBVD = `lexp(-BVD / Nvtmd) * XJBVD;
T2y = max(IJTHDFWD / Isbd, 10.0);
Tb = 1.0 + T2y - XExpBVD;
VjdmFwd = Nvtmd * lln(0.5 * (Tb + sqrt(Tb * Tb + 4 * XExpBVD)));
T0y = lexp(VjdmFwd / Nvtmd);
T0y = `lexp(VjdmFwd / Nvtmd);
IVjdmFwd = Isbd * (T0y - XExpBVD / T0y + XExpBVD - 1.0);
DslpFwd = Isbd * (T0y + XExpBVD / T0y) / Nvtmd;
T2y = `hypsmooth(IJTHDREV / Isbd - 10.0, 1.0E-3) + 10.0;
VjdmRev = -BVD - Nvtmd * lln((T2y - 1.0) / XJBVD);
T1y = XJBVD * lexp(-(BVD + VjdmRev) / Nvtmd);
T1y = XJBVD * `lexp(-(BVD + VjdmRev) / Nvtmd);
IVjdmRev = Isbd * (1.0 + T1y);
DslpRev = -Isbd * T1y / Nvtmd;
end else begin
@ -3423,7 +3413,7 @@ DevTemp = $temperature + DTEMP;
Moc = Moc * MdL;
// Calculate Vadits -- Ref BSIM4
T1 = lexp(PDITSD_i * Vds);
T1 = `lexp(PDITSD_i * Vds);
if (PDITS_i > 0.0) begin
T2y = 1.0 + PDITSL * Leff;
VaDITS = (1.0 + T2y * T1) / PDITS_i;
@ -3566,7 +3556,7 @@ DevTemp = $temperature + DTEMP;
Iii = 0.0;
else if (diffVds > BETA0_t / `EXPL_THRESHOLD) begin
T1 = - BETA0_t / diffVds;
Iii = ALPHA0_i * diffVds * ids * lexp(T1)/Mscbe;
Iii = ALPHA0_i * diffVds * ids * `lexp(T1)/Mscbe;
end else begin
Iii = ALPHA0_i * diffVds * ids * `MIN_EXPL/Mscbe;
end
@ -3587,21 +3577,21 @@ DevTemp = $temperature + DTEMP;
// Igbinv
if (IGBMOD != 0) begin
T1 = Voxmacc / NIGBACC_i / Vt; //representative of Qb in acc
Vaux_Igbacc = NIGBACC_i * Vt * lln(1 + lexp(-T1));
Vaux_Igbacc = NIGBACC_i * Vt * lln(1 + `lexp(-T1));
T2 = AIGBACC_i - BIGBACC_i * Voxmacc;
T3 = 1.0 + CIGBACC_i * Voxmacc;
T4 = -7.45669e11 * TOXE * T2 * T3;
T5 = lexp (T4);
T5 = `lexp(T4);
T6 = 4.97232e-7;
igbacc = NF * Weff * Leff * T6 * ToxRatio * Vg * Vaux_Igbacc * T5;
igbacc = igbacc * igtemp;
T1 = (Voxminv - EIGBINV_i) / NIGBINV_i / Vt;
Vaux_Igbinv = NIGBINV_i * Vt * lln(1.0 + lexp(T1));
Vaux_Igbinv = NIGBINV_i * Vt * lln(1.0 + `lexp(T1));
T2 = AIGBINV_i - BIGBINV_i * Voxminv;
T3 = 1.0 + CIGBINV_i * Voxminv;
T4 = -9.82222e11 * TOXE * T2 * T3;
T5 = lexp (T4);
T5 = `lexp(T4);
T6 = 3.75956e-7;
igbinv = NF * Weff * Leff * T6 * ToxRatio * Vg * Vaux_Igbinv * T5;
igbinv = igbinv * igtemp;
@ -3613,13 +3603,13 @@ DevTemp = $temperature + DTEMP;
T1 = AIGC_i - BIGC_i * Voxminv;
T2 = 1.0 + CIGC_i * Voxminv;
T3 = Bechvb * T1 * T2;
T4 = nq * nVt * (qs + qdeff) * lexp(T3);
T4 = nq * nVt * (qs + qdeff) * `lexp(T3);
igc0 = NF * Aechvb * T4 * (Vg + 0.5 *Vdsx - 0.5*(Vs+Vd)) * igtemp;
// Gate-current partitioning
Vdseffx = sqrt(Vdseff*Vdseff + 0.01) - 0.1;
T1 = PIGCD_i * Vdseffx;
T1_exp = lexp(-T1);
T1_exp = `lexp(-T1);
T3 = T1 + T1_exp -1.0 + 1.0E-4;
T4 = 1.0 - (T1 + 1.0) * T1_exp + 1.0E-4;
T5 = T1 * T1 + 2.0E-4;
@ -3637,7 +3627,7 @@ DevTemp = $temperature + DTEMP;
T1 = AIGS_i - BIGS_i * Vgs_eff;
T2 = 1.0 + CIGS_i * Vgs_eff;
T3 = BechvbEdge * T1 * T2;
T4 = lexp(T3);
T4 = `lexp(T3);
igs_mult = igtemp * NF * AechvbEdge * DLCIG_i;
igs = igs_mult * Vgs_noswap * Vgs_eff * T4;
@ -3648,7 +3638,7 @@ DevTemp = $temperature + DTEMP;
T1 = AIGD_i - BIGD_i * Vgd_eff;
T2 = 1.0 + CIGD_i * Vgd_eff;
T3 = BechvbEdge * T1 * T2;
T4 = lexp(T3);
T4 = `lexp(T3);
igd_mult = igtemp * NF * AechvbEdge * DLCIGD_i;
igd = igd_mult * Vgd_noswap * Vgd_eff * T4;
end
@ -3678,7 +3668,7 @@ DevTemp = $temperature + DTEMP;
T5 = `hypsmooth(T3/T4, 1.0E-6) - 1.0E-6;
end else
T5 = 1.0;
T6 = AGIDL_i * Weff * T1 * lexp(-T2) * T5;
T6 = AGIDL_i * Weff * T1 * `lexp(-T2) * T5;
end
igidl = T6;
@ -3695,7 +3685,7 @@ DevTemp = $temperature + DTEMP;
T5 = `hypsmooth(T3/T4, 1.0E-6) - 1.0E-6;
end else
T5 = 1.0;
T6 = AGISL_i * Weff * T1 * lexp(-T2) * T5;
T6 = AGISL_i * Weff * T1 * `lexp(-T2) * T5;
end
igisl = T6;
end // GIDLMOD
@ -3708,14 +3698,14 @@ DevTemp = $temperature + DTEMP;
if(Isbs > 0.0) begin
if (Vbs_jct < VjsmRev) begin
T0 = Vbs_jct / Nvtms;
T1 = lexp(T0) - 1.0;
T1 = `lexp(T0) - 1.0;
T2 = IVjsmRev + SslpRev * (Vbs_jct - VjsmRev);
Ibs = T1 * T2;
end else if (Vbs_jct <= VjsmFwd) begin
T0 = Vbs_jct / Nvtms;
T1 = (BVS + Vbs_jct) / Nvtms;
T2 = lexp(-T1);
Ibs = Isbs * (lexp(T0) + XExpBVS - 1.0 - XJBVS * T2);
T2 = `lexp(-T1);
Ibs = Isbs * (`lexp(T0) + XExpBVS - 1.0 - XJBVS * T2);
end else
Ibs = IVjsmFwd + SslpFwd * (Vbs_jct - VjsmFwd);
end else
@ -3725,33 +3715,33 @@ DevTemp = $temperature + DTEMP;
if(JTSS_t > 0.0) begin
if((VTSS - Vbs_jct) < (VTSS * 1.0e-3)) begin
T0 = -Vbs_jct / Vtm0 / NJTS_t;
T1 = lexp(T0 * 1.0e+3) - 1.0;
T1 = `lexp(T0 * 1.0e+3) - 1.0;
Ibs = Ibs - ASeff * JTSS_t * T1;
end else begin
T0 = -Vbs_jct / Vtm0 / NJTS_t;
T1 = lexp(T0 * VTSS / (VTSS - Vbs_jct)) - 1.0;
T1 = `lexp(T0 * VTSS / (VTSS - Vbs_jct)) - 1.0;
Ibs = Ibs - ASeff * JTSS_t * T1;
end
end
if(JTSSWS_t > 0.0) begin
if((VTSSWS - Vbs_jct) < (VTSSWS * 1.0e-3)) begin
T0 = -Vbs_jct / Vtm0 / NJTSSW_t;
T1 = lexp(T0 * 1.0e+3) - 1.0;
T1 = `lexp(T0 * 1.0e+3) - 1.0;
Ibs = Ibs - PSeff * JTSSWS_t * T1;
end else begin
T0 = -Vbs_jct / Vtm0 / NJTSSW_t;
T1 = lexp(T0 * VTSSWS / (VTSSWS - Vbs_jct)) - 1.0;
T1 = `lexp(T0 * VTSSWS / (VTSSWS - Vbs_jct)) - 1.0;
Ibs = Ibs - PSeff * JTSSWS_t * T1;
end
end
if(JTSSWGS_t > 0.0) begin
if((VTSSWGS - Vbs_jct) < (VTSSWGS * 1.0e-3)) begin
T0 = -Vbs_jct / Vtm0 / NJTSSWG_t;
T1 = lexp(T0 * 1.0e+3) - 1.0;
T1 = `lexp(T0 * 1.0e+3) - 1.0;
Ibs = Ibs - Weffcj * NF * JTSSWGS_t * T1;
end else begin
T0 = -Vbs_jct / Vtm0 / NJTSSWG_t;
T1 = lexp(T0 * VTSSWGS / (VTSSWGS - Vbs_jct)) - 1.0;
T1 = `lexp(T0 * VTSSWGS / (VTSSWGS - Vbs_jct)) - 1.0;
Ibs = Ibs - Weffcj * NF * JTSSWGS_t * T1;
end
end
@ -3761,14 +3751,14 @@ DevTemp = $temperature + DTEMP;
if(Isbd > 0.0) begin
if (Vbd_jct < VjdmRev) begin
T0 = Vbd_jct / Nvtmd;
T1 = lexp(T0) - 1.0;
T1 = `lexp(T0) - 1.0;
T2 = IVjdmRev + DslpRev * (Vbd_jct - VjdmRev);
Ibd = T1 * T2;
end else if (Vbd_jct <= VjdmFwd) begin
T0 = Vbd_jct / Nvtmd;
T1 = (BVD + Vbd_jct) / Nvtmd;
T2 = lexp(-T1);
Ibd = Isbd * (lexp(T0) + XExpBVD - 1.0 - XJBVD * T2);
T2 = `lexp(-T1);
Ibd = Isbd * (`lexp(T0) + XExpBVD - 1.0 - XJBVD * T2);
end else
Ibd = IVjdmFwd + DslpFwd * (Vbd_jct - VjdmFwd);
end else
@ -3778,33 +3768,33 @@ DevTemp = $temperature + DTEMP;
if(JTSD_t > 0.0) begin
if((VTSD - Vbd_jct) < (VTSD * 1.0e-3)) begin
T0 = -Vbd_jct / Vtm0 / NJTSD_t;
T1 = lexp(T0 * 1.0e+3) - 1.0;
T1 = `lexp(T0 * 1.0e+3) - 1.0;
Ibd = Ibd - ADeff * JTSD_t * T1;
end else begin
T0 = -Vbd_jct / Vtm0 / NJTSD_t;
T1 = lexp(T0 * VTSD/ (VTSD - Vbd_jct)) - 1.0;
T1 = `lexp(T0 * VTSD/ (VTSD - Vbd_jct)) - 1.0;
Ibd = Ibd - ADeff * JTSD_t * T1;
end
end
if(JTSSWD_t > 0.0) begin
if((VTSSWD - Vbd_jct) < (VTSSWD * 1.0e-3)) begin
T0 = -Vbd_jct / Vtm0 / NJTSSWD_t;
T1 = lexp(T0 * 1.0e+3) - 1.0;
T1 = `lexp(T0 * 1.0e+3) - 1.0;
Ibd = Ibd - PDeff * JTSSWD_t * T1;
end else begin
T0 = -Vbd_jct / Vtm0 / NJTSSWD_t;
T1 = lexp(T0 * VTSSWD / (VTSSWD - Vbd_jct)) - 1.0;
T1 = `lexp(T0 * VTSSWD / (VTSSWD - Vbd_jct)) - 1.0;
Ibd = Ibd - PDeff * JTSSWD_t * T1;
end
end
if(JTSSWGD_t > 0.0) begin
if((VTSSWGD - Vbd_jct) < (VTSSWGD * 1.0e-3)) begin
T0 = -Vbd_jct / Vtm0 / NJTSSWGD_t;
T1 = lexp(T0 * 1.0e+3) - 1.0;
T1 = `lexp(T0 * 1.0e+3) - 1.0;
Ibd = Ibd - Weffcj * NF * JTSSWGD_t * T1;
end else begin
T0 = -Vbd_jct / Vtm0 / NJTSSWGD_t;
T1 = lexp(T0 * VTSSWGD / (VTSSWGD - Vbd_jct)) - 1.0;
T1 = `lexp(T0 * VTSSWGD / (VTSSWGD - Vbd_jct)) - 1.0;
Ibd = Ibd - Weffcj * NF * JTSSWGD_t * T1;
end
end
@ -3816,7 +3806,7 @@ DevTemp = $temperature + DTEMP;
if (T1<0.9) begin
arg = 1.0 - T1;
if (MJS == 0.5) sarg = 1.0 / sqrt(arg);
else sarg = lexp(-MJS * lln(arg));
else sarg = `lexp(-MJS * lln(arg));
Qbs = PBS_t * Czbs * (1.0 - arg * sarg) / (1.0 - MJS);
end else begin
T2 = pow(0.1,-MJS);
@ -3834,7 +3824,7 @@ DevTemp = $temperature + DTEMP;
if (T1<0.9) begin
arg = 1.0 - T1;
if (MJSWS == 0.5) sarg = 1.0 / sqrt(arg);
else sarg = lexp(-MJSWS * lln(arg));
else sarg = `lexp(-MJSWS * lln(arg));
Qbs = Qbs + PBSWS_t * Czbssw * (1.0 - arg * sarg) / (1.0 - MJSWS);
end else begin
T2 = pow(0.1,-MJSWS);
@ -3850,7 +3840,7 @@ DevTemp = $temperature + DTEMP;
if (T1<0.9) begin
arg = 1.0 - T1;
if (MJSWGS == 0.5) sarg = 1.0 / sqrt(arg);
else sarg = lexp(-MJSWGS * lln(arg));
else sarg = `lexp(-MJSWGS * lln(arg));
Qbs = Qbs + PBSWGS_t * Czbsswg * (1.0 - arg * sarg) / (1.0 - MJSWGS);
end else begin
T2 = pow(0.1,-MJSWGS);
@ -3867,7 +3857,7 @@ DevTemp = $temperature + DTEMP;
if (T1<0.9) begin
arg = 1.0 - T1;
if (MJD == 0.5) sarg = 1.0 / sqrt(arg);
else sarg = lexp(-MJD * lln(arg));
else sarg = `lexp(-MJD * lln(arg));
Qbd = PBD_t * Czbd * (1.0 - arg * sarg) / (1.0 - MJD);
end else begin
T2 = pow(0.1,-MJD);
@ -3885,7 +3875,7 @@ DevTemp = $temperature + DTEMP;
if (T1<0.9) begin
arg = 1.0 - T1;
if (MJSWD == 0.5) sarg = 1.0 / sqrt(arg);
else sarg = lexp(-MJSWD * lln(arg));
else sarg = `lexp(-MJSWD * lln(arg));
Qbd = Qbd + PBSWD_t * Czbdsw * (1.0 - arg * sarg) / (1.0 - MJSWD);
end else begin
T2 = pow(0.1,-MJSWD);
@ -3901,7 +3891,7 @@ DevTemp = $temperature + DTEMP;
if (T1<0.9) begin
arg = 1.0 - T1;
if (MJSWGD == 0.5) sarg = 1.0 / sqrt(arg);
else sarg = lexp(-MJSWGD * lln(arg));
else sarg = `lexp(-MJSWGD * lln(arg));
Qbd = Qbd + PBSWGD_t * Czbdswg * (1.0 - arg * sarg) / (1.0 - MJSWGD);
end else begin
T2 = pow(0.1,-MJSWGD);