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

started dual numbers in hicum2temp and a transfer struct in hicum2defs

This commit is contained in:
mariok 2020-05-12 17:14:17 +02:00
parent 479c2a19ea
commit 96c66dd645
7 changed files with 487 additions and 63 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
src/spicelib/devices/hicum2/Makefile.am
View File

@ -22,12 +22,13 @@ libhicum2_la_SOURCES = \
hicum2pzld.c \
hicum2setup.c \
hicum2soachk.c \
hicum2temp.c \
hicum2temp.hpp \
hicum2temp.cpp \
hicum2trunc.c
AM_CPPFLAGS = @AM_CPPFLAGS@ -I$(top_srcdir)/src/include -lstdc++ -std=c++11 -Og -Wall
AM_CFLAGS = -lstdc++ -I$(top_srcdir)/src/include -Og -Wall
AM_CXXFLAGS = -I$(top_srcdir)/src/include -lstdc++ -std=c++11 -Og -Wall
AM_CPPFLAGS = @AM_CPPFLAGS@ -I$(top_srcdir)/src/include -lstdc++ -std=c++11 -Og -Wall
AM_CFLAGS = -lstdc++ -I$(top_srcdir)/src/include -Og -Wall
AM_CXXFLAGS = -I$(top_srcdir)/src/include -lstdc++ -std=c++11 -Og -Wall
MAINTAINERCLEANFILES = Makefile.in

116
src/spicelib/devices/hicum2/hicum2defs.h
View File

@ -36,6 +36,11 @@ enum {
HICUMNSRCS
};
typedef struct sDualDouble {
double rpart;
double dpart;
} dual_double;
/* data needed to describe a single instance */
typedef struct sHICUMinstance {
@ -76,65 +81,62 @@ typedef struct sHICUMinstance {
double HICUMdtemp_sh; /* instance delta temperature because of self-heating */
double HICUMm; /* multiply factor for the hicum */
double HICUMmg;
double HICUMvt0;
double HICUMvt;
double HICUMdT;
double HICUMqtt0;
double HICUMln_qtt0;
double HICUMibeis_t;
double HICUMireis_t;
double HICUMibeps_t;
double HICUMireps_t;
double HICUMibcis_t;
double HICUMibcxs_t;
double HICUMitss_t;
double HICUMiscs_t;
double HICUMqp0_t;
double HICUMvlim_t;
double HICUMtef0_t;
double HICUMc10_t;
double HICUMrci0_t;
double HICUMvces_t;
double HICUMt0_t;
double HICUMthcs_t;
double HICUMfavl_t;
double HICUMqavl_t;
double HICUMkavl_t;
double HICUMrbi0_t;
double HICUMibets_t;
double HICUMabet_t;
double HICUMcjcx01_t;
double HICUMcjcx02_t;
double HICUMrcx_t;
double HICUMrbx_t;
double HICUMre_t;
double HICUMrbi;
double HICUMtsf_t;
double HICUMcscp0_t;
double HICUMvdsp_t;
double HICUMvptsp_t;
double HICUMahjei_t;
double HICUMhjei0_t;
double HICUMhf0_t;
double HICUMhfc_t;
double HICUMhfe_t;
double HICUMrth_t;
dual_double HICUMvt;
dual_double HICUMdT;
dual_double HICUMibeis_t;
dual_double HICUMireis_t;
dual_double HICUMibeps_t;
dual_double HICUMireps_t;
dual_double HICUMibcis_t;
dual_double HICUMibcxs_t;
dual_double HICUMitss_t;
dual_double HICUMiscs_t;
dual_double HICUMqp0_t;
dual_double HICUMvlim_t;
dual_double HICUMtef0_t;
dual_double HICUMc10_t;
dual_double HICUMrci0_t;
dual_double HICUMvces_t;
dual_double HICUMt0_t;
dual_double HICUMthcs_t;
dual_double HICUMfavl_t;
dual_double HICUMqavl_t;
dual_double HICUMkavl_t;
dual_double HICUMrbi0_t;
dual_double HICUMibets_t;
dual_double HICUMabet_t;
dual_double HICUMcjcx01_t;
dual_double HICUMcjcx02_t;
dual_double HICUMrcx_t;
dual_double HICUMrbx_t;
dual_double HICUMre_t;
dual_double HICUMrbi;
dual_double HICUMtsf_t;
dual_double HICUMcscp0_t;
dual_double HICUMvdsp_t;
dual_double HICUMvptsp_t;
dual_double HICUMahjei_t;
dual_double HICUMhjei0_t;
dual_double HICUMhf0_t;
dual_double HICUMhfc_t;
dual_double HICUMhfe_t;
dual_double HICUMrth_t;
double HICUMvdei_t;
double HICUMajei_t;
double HICUMcjci0_t;
double HICUMvdci_t;
double HICUMvptci_t;
double HICUMvdep_t;
double HICUMajep_t;
double HICUMcjep0_t;
double HICUMcjei0_t;
double HICUMvdcx_t;
double HICUMvptcx_t;
double HICUMcjs0_t;
double HICUMvds_t;
double HICUMvpts_t;
dual_double HICUMvdei_t;
dual_double HICUMajei_t;
dual_double HICUMcjci0_t;
dual_double HICUMvdci_t;
dual_double HICUMvptci_t;
dual_double HICUMvdep_t;
dual_double HICUMajep_t;
dual_double HICUMcjep0_t;
dual_double HICUMcjei0_t;
dual_double HICUMvdcx_t;
dual_double HICUMvptcx_t;
dual_double HICUMcjs0_t;
dual_double HICUMvds_t;
dual_double HICUMvpts_t;
double HICUMiavl;
double HICUMpterm;

2
src/spicelib/devices/hicum2/hicum2ext.h
View File

@ -20,7 +20,7 @@ extern int HICUMparam(int,IFvalue*,GENinstance*,IFvalue*);
extern int HICUMpzLoad(GENmodel*, CKTcircuit*, SPcomplex*);
extern int HICUMsetup(SMPmatrix*,GENmodel*,CKTcircuit*,int*);
extern int HICUMunsetup(GENmodel*,CKTcircuit*);
extern int HICUMtemp(GENmodel*,CKTcircuit*);
// extern int HICUMtemp(GENmodel*,CKTcircuit*); // moved to hicum2temp.hpp
extern int HICUMtrunc(GENmodel*,CKTcircuit*,double*);
extern int HICUMnoise(int,int,GENmodel*,CKTcircuit*,Ndata*,double*);
extern int HICUMsoaCheck(CKTcircuit *, GENmodel *);

1
src/spicelib/devices/hicum2/hicum2init.c
View File

@ -11,6 +11,7 @@
#include "hicum2ext.h"
#include "hicum2init.h"
#include "hicumL2.hpp"
#include "hicum2temp.hpp"
SPICEdev HICUMinfo = {

10
src/spicelib/devices/hicum2/hicum2temp.c
View File

@ -98,6 +98,7 @@ int hicum_thermal_update(HICUMmodel *inModel, HICUMinstance *inInstance)
double k10,k20,avs,vgb_t0,vge_t0,vgbe_t0,vgbe0,vgbc0,vgsc0;
double zetabci,zetabcxt,zetasct;
double k1,k2,dvg0,vge_t,vgb_t,vgbe_t,cratio_t,a;
double k1_dT,k2_dT,dvg0_dT,vge_t_dT,vgb_t_dT,vgbe_t_dT,cratio_t_dT,a_dT;
double Tnom, dT, zetatef, cjcx01, cjcx02, C_1;
double cjci0_t, vdci_t, vptci_t, cjep0_t, vdep_t, ajep_t, vdcx_t, vptcx_t, cscp0_t, vdsp_t, vptsp_t, cjs0_t, vds_t, vpts_t;
@ -128,14 +129,22 @@ int hicum_thermal_update(HICUMmodel *inModel, HICUMinstance *inInstance)
}
here->HICUMvt0 = Tnom * CONSTKoverQ;
here->HICUMvt = here->HICUMtemp * CONSTKoverQ;
here->HICUMvt_dT = CONSTKoverQ;
dT = here->HICUMtemp-Tnom;
here->HICUMqtt0 = here->HICUMtemp/Tnom;
here->HICUMqtt0_dT = 1/Tnom;
here->HICUMln_qtt0 = log(here->HICUMqtt0);
here->HICUMln_qtt0_dT = 1/here->HICUMtemp; // d(log(x/a))/dx = 1/x
k1 = model->HICUMf1vg*here->HICUMtemp*log(here->HICUMtemp);
k1_dT = model->HICUMf1vg*(log(here->HICUMtemp) + 1);
k2 = model->HICUMf2vg*here->HICUMtemp;
k2_dT = model->HICUMf2vg;
vgb_t = model->HICUMvgb+k1+k2;
vgb_t_dT = k1_dT+k2_dT;
vge_t = model->HICUMvge+k1+k2;
vge_t_dT = k1_dT+k2_dT;
vgbe_t = (vgb_t+vge_t)/2;
vgbe_t_dT = (vgb_t_dT+vge_t_dT)/2;
here->HICUMtVcrit = here->HICUMvt *
log(here->HICUMvt / (CONSTroot2*model->HICUMibeis*here->HICUMarea*here->HICUMm));
@ -350,4 +359,3 @@ int hicum_thermal_update(HICUMmodel *inModel, HICUMinstance *inInstance)
return(0);
}

399
src/spicelib/devices/hicum2/hicum2temp.cpp Normal file
View File

@ -0,0 +1,399 @@
/**********
Copyright 1990 Regents of the University of California. All rights reserved.
Author: 1985 Thomas L. Quarles
Model Author: 1990 Michael Schröter TU Dresden
Spice3 Implementation: 2019 Dietmar Warning
**********/
#include "cmath"
#ifndef M_PI
#define M_PI 3.1415926535897932384626433832795
#endif
#include <duals/dual>
#include "hicum2temp.hpp"
#include <functional>
#include <fenv.h> //trap NAN
//ngspice header files written in C
#ifdef __cplusplus
extern "C"
{
#endif
#include "ngspice/ngspice.h"
#include "ngspice/cktdefs.h"
#include "ngspice/smpdefs.h"
#include "hicum2defs.h"
#include "ngspice/const.h"
#include "ngspice/sperror.h"
#include "ngspice/ifsim.h"
#include "ngspice/suffix.h"
#ifdef __cplusplus
}
#endif
using namespace duals::literals;
#define TMAX 326.85
#define TMIN -100.0
#define LN_EXP_LIMIT 11.0
void TMPHICJ(double , double , double , double , double ,
double , double , double , double , double , double ,
double *, double *, double *);
// TEMPERATURE UPDATE OF JUNCTION CAPACITANCE RELATED PARAMETERS
// INPUT:
// mostly model parameters
// x : zero bias junction capacitance
// y : junction built-in potential
// z : grading co-efficient
// w : ratio of maximum to zero-bias value of capacitance or punch-through voltage
// is_al : condition factor to check what "w" stands for
// vgeff : band-gap voltage
// IMPLICIT INPUT:
// vt : thermal voltage
// vt0,qtt0,ln_qtt0,mg : other model variables
// OUTPUT:
// c_j_t : temperature update of "c_j"
// vd_t : temperature update of "vd0"
// w_t : temperature update of "w"
void TMPHICJ(duals::duald vt, double vt0, duals::duald qtt0, duals::duald ln_qtt0, double mg,
double c_j, double vd0, double z, double w, double is_al, double vgeff,
duals::duald *c_j_t, duals::duald *vd_t, duals::duald *w_t)
{
duals::duald vdj0,vdjt,vdt;
if (c_j > 0.0) {
vdj0 = 2*vt0*log(exp(vd0*0.5/vt0)-exp(-0.5*vd0/vt0));
vdjt = vdj0*qtt0+vgeff*(1-qtt0)-mg*vt*ln_qtt0;
vdt = vdjt+2*vt*log(0.5*(1+sqrt(1+4*exp(-vdjt/vt))));
*vd_t = vdt;
*c_j_t = c_j*exp(z*log(vd0/(*vd_t)));
if (is_al == 1) {
*w_t = w*(*vd_t)/vd0;
} else {
*w_t = w;
}
} else {
*c_j_t = c_j;
*vd_t = vd0;
*w_t = w;
}
}
void hicum_TMPHICJ(duals::duald vt, double vt0, duals::duald qtt0, duals::duald ln_qtt0, double mg,
double c_j, double vd0, double z, double w, double is_al, double vgeff,
double *c_j_t, double *vd_t, double *w_t,
double *c_j_t_dT, double *vd_t_dT, double *w_t_dT)
{
duals::duald c_j_t_result = 0;
duals::duald vd_t_result = 0;
duals::duald w_t_result = 0;
TMPHICJ(vt, vt0, qtt0, ln_qtt0, mg, c_j, vd0, z, w, is_al, vgeff, &c_j_t_result, &vd_t_result, &w_t_result);
*c_j_t = c_j_t_result.rpart();
*c_j_t_dT = c_j_t_result.dpart();
*vd_t = vd_t_result.rpart();
*vd_t_dT = vd_t_result.dpart();
*w_t = w_t_result.rpart();
*w_t_dT = w_t_result.dpart();
}
int iret, hicum_thermal_update(HICUMmodel *, HICUMinstance *);
int
HICUMtemp(GENmodel *inModel, CKTcircuit *ckt)
/* Pre-compute many useful parameters
*/
{
HICUMmodel *model = (HICUMmodel *)inModel;
HICUMinstance *here;
/* loop through all the bipolar models */
for( ; model != NULL; model = HICUMnextModel(model)) {
/* loop through all the instances of the model */
for (here = HICUMinstances(model); here != NULL ;
here=HICUMnextInstance(here)) {
if(!here->HICUMtempGiven) here->HICUMtemp = ckt->CKTtemp;
if(here->HICUMdtempGiven) here->HICUMtemp = here->HICUMtemp + here->HICUMdtemp;
iret = hicum_thermal_update(model, here);
}
}
return(OK);
}
int hicum_thermal_update(HICUMmodel *inModel, HICUMinstance *inInstance)
{
HICUMmodel *model = (HICUMmodel *)inModel;
HICUMinstance *here = (HICUMinstance *)inInstance;
double mg,k10,k20,avs,vgb_t0,vge_t0,vgbe_t0,vgbe0,vgbc0,vgsc0;
double zetabci,zetabcxt,zetasct;
duals::duald temp, dT, vt, qtt0, ln_qtt0;
duals::duald k1,k2,dvg0,vge_t,vgb_t,vgbe_t,cratio_t,a;
double Tnom, zetatef, cjcx01, cjcx02, C_1;
double cjci0_t, vdci_t, vptci_t, cjep0_t, vdep_t, ajep_t, vdcx_t, vptcx_t, cscp0_t, vdsp_t, vptsp_t, cjs0_t, vds_t, vpts_t;
Tnom = model->HICUMtnom;
k10 = model->HICUMf1vg*Tnom*log(Tnom);
k20 = model->HICUMf2vg*Tnom;
avs = model->HICUMalvs*Tnom;
vgb_t0 = model->HICUMvgb+k10+k20;
vge_t0 = model->HICUMvge+k10+k20;
vgbe_t0 = (vgb_t0+vge_t0)/2;
vgbe0 = (model->HICUMvgb+model->HICUMvge)/2;
vgbc0 = (model->HICUMvgb+model->HICUMvgc)/2;
vgsc0 = (model->HICUMvgs+model->HICUMvgc)/2;
mg = 3-model->HICUMf1vg/CONSTKoverQ;
zetabci = mg+1-model->HICUMzetaci;
zetabcxt= mg+1-model->HICUMzetacx;
zetasct = mg-1.5;
// Limit temperature to avoid FPEs in equations
if(here->HICUMtemp < TMIN + CONSTCtoK) {
here->HICUMtemp = TMIN + CONSTCtoK;
} else {
if (here->HICUMtemp > TMAX + CONSTCtoK) {
here->HICUMtemp = TMAX + CONSTCtoK;
}
}
temp = here->HICUMtemp+1_e;
vt = temp*CONSTKoverQ; // dual temperature
here->HICUMvt0 = Tnom * CONSTKoverQ;
here->HICUMvt.rpart = vt.rpart();
here->HICUMvt.dpart = vt.dpart();
dT = temp-Tnom;
qtt0 = temp/Tnom;
ln_qtt0 = log(qtt0);
k1 = model->HICUMf1vg*temp*log(temp);
k2 = model->HICUMf2vg*temp;
vgb_t = model->HICUMvgb+k1+k2;
vge_t = model->HICUMvge+k1+k2;
vgbe_t = (vgb_t+vge_t)/2;
here->HICUMtVcrit = here->HICUMvt.rpart *
log(here->HICUMvt.rpart / (CONSTroot2*model->HICUMibeis*here->HICUMarea*here->HICUMm));
//Internal b-e junction capacitance
// TMPHICJ(here->HICUMvt0,here->HICUMvt,here->HICUMqtt0,here->HICUMln_qtt0,here->HICUMmg,model->HICUMcjei0,model->HICUMvdei,model->HICUMzei,model->HICUMajei,1,vgbe0,&here->HICUMcjei0_t,&here->HICUMvdei_t,&here->HICUMajei_t);
hicum_TMPHICJ(vt, here->HICUMvt0, qtt0, ln_qtt0, mg,
model->HICUMcjei0, model->HICUMvdei, model->HICUMzei, model->HICUMajei, 1, vgbe0,
&here->HICUMcjei0_t.rpart, &here->HICUMvdei_t.rpart, &here->HICUMajei_t.rpart,
&here->HICUMcjei0_t.dpart, &here->HICUMvdei_t.dpart, &here->HICUMajei_t.dpart);
// if (model->HICUMflcomp == 0.0 || model->HICUMflcomp == 2.1) {
// double V_gT, r_VgVT, k;
// V_gT = 3.0*here->HICUMvt*here->HICUMln_qtt0 + model->HICUMvgb*(here->HICUMqtt0-1.0);
// r_VgVT = V_gT/here->HICUMvt;
// //Internal b-e diode saturation currents
// a = model->HICUMmcf*r_VgVT/model->HICUMmbei - model->HICUMalb*dT;
// here->HICUMibeis_t = model->HICUMibeis*exp(a);
// a = model->HICUMmcf*r_VgVT/model->HICUMmrei - model->HICUMalb*dT;
// here->HICUMireis_t = model->HICUMireis*exp(a);
// a = model->HICUMmcf*r_VgVT/model->HICUMmbep - model->HICUMalb*dT;
// //Peripheral b-e diode saturation currents
// here->HICUMibeps_t = model->HICUMibeps*exp(a);
// a = model->HICUMmcf*r_VgVT/model->HICUMmrep - model->HICUMalb*dT;
// here->HICUMireps_t = model->HICUMireps*exp(a);
// //Internal b-c diode saturation current
// a = r_VgVT/model->HICUMmbci;
// here->HICUMibcis_t = model->HICUMibcis*exp(a);
// //External b-c diode saturation currents
// a = r_VgVT/model->HICUMmbcx;
// here->HICUMibcxs_t = model->HICUMibcxs*exp(a);
// //Saturation transfer current for substrate transistor
// a = r_VgVT/model->HICUMmsf;
// here->HICUMitss_t = model->HICUMitss*exp(a);
// //Saturation current for c-s diode
// a = r_VgVT/model->HICUMmsc;
// here->HICUMiscs_t = model->HICUMiscs*exp(a);
// //Zero bias hole charge
// a = here->HICUMvdei_t/model->HICUMvdei;
// here->HICUMqp0_t = model->HICUMqp0*(1.0+0.5*model->HICUMzei*(1.0-a));
// //Voltage separating ohmic and saturation velocity regime
// a = model->HICUMvlim*(1.0-model->HICUMalvs*dT)*exp(model->HICUMzetaci*here->HICUMln_qtt0);
// k = (a-here->HICUMvt)/here->HICUMvt;
// if (k < LN_EXP_LIMIT) {
// here->HICUMvlim_t = here->HICUMvt + here->HICUMvt*log(1.0+exp(k));
// } else {
// here->HICUMvlim_t = a;
// }
// //Neutral emitter storage time
// a = 1.0+model->HICUMalb*dT;
// k = 0.5*(a+sqrt(a*a+0.01));
// here->HICUMtef0_t = model->HICUMtef0*here->HICUMqtt0/k;
// } else {
// //Internal b-e diode saturation currents
// here->HICUMibeis_t = model->HICUMibeis*exp(model->HICUMzetabet*here->HICUMln_qtt0+model->HICUMvge/here->HICUMvt*(here->HICUMqtt0-1));
// if (model->HICUMflcomp>=2.3) {
// here->HICUMireis_t = model->HICUMireis*exp(here->HICUMmg/model->HICUMmrei*here->HICUMln_qtt0+vgbe0/(model->HICUMmrei*here->HICUMvt)*(here->HICUMqtt0-1));
// } else {
// here->HICUMireis_t = model->HICUMireis*exp(0.5*here->HICUMmg*here->HICUMln_qtt0+0.5*vgbe0/here->HICUMvt*(here->HICUMqtt0-1));
// }
// //Peripheral b-e diode saturation currents
// here->HICUMibeps_t = model->HICUMibeps*exp(model->HICUMzetabet*here->HICUMln_qtt0+model->HICUMvge/here->HICUMvt*(here->HICUMqtt0-1));
// if (model->HICUMflcomp>=2.3) {
// here->HICUMireps_t = model->HICUMireps*exp(here->HICUMmg/model->HICUMmrep*here->HICUMln_qtt0+vgbe0/(model->HICUMmrep*here->HICUMvt)*(here->HICUMqtt0-1));
// } else {
// here->HICUMireps_t = model->HICUMireps*exp(0.5*here->HICUMmg*here->HICUMln_qtt0+0.5*vgbe0/here->HICUMvt*(here->HICUMqtt0-1));
// }
// //Internal b-c diode saturation currents
// here->HICUMibcis_t = model->HICUMibcis*exp(zetabci*here->HICUMln_qtt0+model->HICUMvgc/here->HICUMvt*(here->HICUMqtt0-1));
// //External b-c diode saturation currents
// here->HICUMibcxs_t = model->HICUMibcxs*exp(zetabcxt*here->HICUMln_qtt0+model->HICUMvgc/here->HICUMvt*(here->HICUMqtt0-1));
// //Saturation transfer current for substrate transistor
// here->HICUMitss_t = model->HICUMitss*exp(zetasct*here->HICUMln_qtt0+model->HICUMvgc/here->HICUMvt*(here->HICUMqtt0-1));
// //Saturation current for c-s diode
// here->HICUMiscs_t = model->HICUMiscs*exp(zetasct*here->HICUMln_qtt0+model->HICUMvgs/here->HICUMvt*(here->HICUMqtt0-1));
// //Zero bias hole charge
// a = exp(model->HICUMzei*log(here->HICUMvdei_t/model->HICUMvdei));
// here->HICUMqp0_t = model->HICUMqp0*(2.0-a);
// //Voltage separating ohmic and saturation velocity regime
// here->HICUMvlim_t = model->HICUMvlim*exp((model->HICUMzetaci-avs)*here->HICUMln_qtt0);
// //Neutral emitter storage time
// if (model->HICUMflcomp >= 2.3) {
// here->HICUMtef0_t = model->HICUMtef0;
// } else {
// zetatef = model->HICUMzetabet-model->HICUMzetact-0.5;
// dvg0 = model->HICUMvgb-model->HICUMvge;
// here->HICUMtef0_t = model->HICUMtef0*exp(zetatef*here->HICUMln_qtt0-dvg0/here->HICUMvt*(here->HICUMqtt0-1));
// }
// }
// //GICCR prefactor
// here->HICUMc10_t = model->HICUMc10*exp(model->HICUMzetact*here->HICUMln_qtt0+model->HICUMvgb/here->HICUMvt*(here->HICUMqtt0-1));
// // Low-field internal collector resistance
// here->HICUMrci0_t = model->HICUMrci0*exp(model->HICUMzetaci*here->HICUMln_qtt0);
// //Voltage separating ohmic and saturation velocity regime
// //vlim_t = model->HICUMvlim*exp((model->HICUMzetaci-avs)*here->HICUMln_qtt0);
// //Internal c-e saturation voltage
// here->HICUMvces_t = model->HICUMvces*(1+model->HICUMalces*dT);
// //Internal b-c diode saturation current
// //ibcis_t = model->HICUMibcis*exp(zetabci*here->HICUMln_qtt0+model->HICUMvgc/here->HICUMvt*(here->HICUMqtt0-1));
// //Internal b-c junction capacitance
// TMPHICJ(here->HICUMvt0,here->HICUMvt,here->HICUMqtt0,here->HICUMln_qtt0,here->HICUMmg,model->HICUMcjci0,model->HICUMvdci,model->HICUMzci,model->HICUMvptci,0,vgbc0,&cjci0_t,&vdci_t,&vptci_t);
// here->HICUMcjci0_t = cjci0_t;
// here->HICUMvdci_t = vdci_t;
// here->HICUMvptci_t = vptci_t;
// //Low-current forward transit time
// here->HICUMt0_t = model->HICUMt0*(1+model->HICUMalt0*dT+model->HICUMkt0*dT*dT);
// //Saturation time constant at high current densities
// here->HICUMthcs_t = model->HICUMthcs*exp((model->HICUMzetaci-1)*here->HICUMln_qtt0);
// //Avalanche current factors
// here->HICUMfavl_t = model->HICUMfavl*exp(model->HICUMalfav*dT);
// here->HICUMqavl_t = model->HICUMqavl*exp(model->HICUMalqav*dT);
// here->HICUMkavl_t = model->HICUMkavl*exp(model->HICUMalkav*dT);
// //Zero bias internal base resistance
// here->HICUMrbi0_t = model->HICUMrbi0*exp(model->HICUMzetarbi*here->HICUMln_qtt0);
// //Peripheral b-e junction capacitance
// TMPHICJ(here->HICUMvt0,here->HICUMvt,here->HICUMqtt0,here->HICUMln_qtt0,here->HICUMmg,model->HICUMcjep0,model->HICUMvdep,model->HICUMzep,model->HICUMajep,1,vgbe0,&cjep0_t,&vdep_t,&ajep_t);
// here->HICUMcjep0_t = cjep0_t;
// here->HICUMvdep_t = vdep_t;
// here->HICUMajep_t = ajep_t;
// //Tunneling current factors
// if (model->HICUMibets > 0) { // HICTUN_T
// double a_eg,ab,aa;
// ab = 1.0;
// aa = 1.0;
// a_eg=vgbe_t0/vgbe_t;
// if(model->HICUMtunode==1 && model->HICUMcjep0 > 0.0 && model->HICUMvdep >0.0) {
// ab = (here->HICUMcjep0_t/model->HICUMcjep0)*sqrt(a_eg)*vdep_t*vdep_t/(model->HICUMvdep*model->HICUMvdep);
// aa = (model->HICUMvdep/vdep_t)*(model->HICUMcjep0/here->HICUMcjep0_t)*pow(a_eg,-1.5);
// } else if (model->HICUMtunode==0 && model->HICUMcjei0 > 0.0 && model->HICUMvdei >0.0) {
// ab = (here->HICUMcjei0_t/model->HICUMcjei0)*sqrt(a_eg)*here->HICUMvdei_t*here->HICUMvdei_t/(model->HICUMvdei*model->HICUMvdei);
// aa = (model->HICUMvdei/here->HICUMvdei_t)*(model->HICUMcjei0/here->HICUMcjei0_t)*pow(a_eg,-1.5);
// }
// here->HICUMibets_t = model->HICUMibets*ab;
// here->HICUMabet_t = model->HICUMabet*aa;
// } else {
// here->HICUMibets_t = 0;
// here->HICUMabet_t = 1;
// }
// //Depletion capacitance splitting at b-c junction
// //Capacitances at peripheral and external base node
// C_1 = (1.0-model->HICUMfbcpar)*(model->HICUMcjcx0+model->HICUMcbcpar);
// if (C_1 >= model->HICUMcbcpar) {
// cjcx01 = C_1-model->HICUMcbcpar;
// cjcx02 = model->HICUMcjcx0-cjcx01;
// } else {
// cjcx01 = 0.0;
// cjcx02 = model->HICUMcjcx0;
// }
// //Temperature mapping for tunneling current is done inside HICTUN
// TMPHICJ(here->HICUMvt0,here->HICUMvt,here->HICUMqtt0,here->HICUMln_qtt0,here->HICUMmg,1.0,model->HICUMvdcx,model->HICUMzcx,model->HICUMvptcx,0,vgbc0,&cratio_t,&vdcx_t,&vptcx_t);
// here->HICUMcjcx01_t=cratio_t*cjcx01;
// here->HICUMcjcx02_t=cratio_t*cjcx02;
// here->HICUMvdcx_t = vdcx_t;
// here->HICUMvptcx_t = vptcx_t;
// //External b-c diode saturation currents
// //ibcxs_t = model->HICUMibcxs*exp(zetabcxt*here->HICUMln_qtt0+model->HICUMvgc/here->HICUMvt*(qtt0-1));
// //Constant external series resistances
// here->HICUMrcx_t = model->HICUMrcx*exp(model->HICUMzetarcx*here->HICUMln_qtt0);
// here->HICUMrbx_t = model->HICUMrbx*exp(model->HICUMzetarbx*here->HICUMln_qtt0);
// here->HICUMre_t = model->HICUMre*exp(model->HICUMzetare*here->HICUMln_qtt0);
// //Forward transit time in substrate transistor
// here->HICUMtsf_t = model->HICUMtsf*exp((model->HICUMzetacx-1.0)*here->HICUMln_qtt0);
// //Capacitance for c-s junction
// TMPHICJ(here->HICUMvt0,here->HICUMvt,here->HICUMqtt0,here->HICUMln_qtt0,here->HICUMmg,model->HICUMcjs0,model->HICUMvds,model->HICUMzs,model->HICUMvpts,0,vgsc0,&cjs0_t,&vds_t,&vpts_t);
// here->HICUMcjs0_t = cjs0_t;
// here->HICUMvds_t = vds_t;
// here->HICUMvpts_t = vpts_t;
// /*Peripheral s-c capacitance
// * Note, thermal update only required for model->HICUMvds > 0
// * Save computional effort otherwise
// */
// if (model->HICUMvdsp > 0) {
// TMPHICJ(here->HICUMvt0,here->HICUMvt,here->HICUMqtt0,here->HICUMln_qtt0,here->HICUMmg,model->HICUMcscp0,model->HICUMvdsp,model->HICUMzsp,model->HICUMvptsp,0,vgsc0,&cscp0_t,&vdsp_t,&vptsp_t);
// here->HICUMcscp0_t = cscp0_t;
// here->HICUMvdsp_t = vdsp_t;
// here->HICUMvptsp_t = vptsp_t;
// } else {
// // Avoid uninitialized variables
// here->HICUMcscp0_t = model->HICUMcscp0;
// here->HICUMvdsp_t = model->HICUMvdsp;
// here->HICUMvptsp_t = model->HICUMvptsp;
// }
// here->HICUMahjei_t = model->HICUMahjei*exp(model->HICUMzetahjei*here->HICUMln_qtt0);
// here->HICUMhjei0_t = model->HICUMhjei*exp(model->HICUMdvgbe/here->HICUMvt*(exp(model->HICUMzetavgbe*log(here->HICUMqtt0))-1));
// here->HICUMhf0_t = model->HICUMhf0*exp(model->HICUMdvgbe/here->HICUMvt*(here->HICUMqtt0-1));
// if (model->HICUMflcomp >= 2.3) {
// here->HICUMhfe_t = model->HICUMhfe*exp((model->HICUMvgb-model->HICUMvge)/here->HICUMvt*(here->HICUMqtt0-1));
// here->HICUMhfc_t = model->HICUMhfc*exp((model->HICUMvgb-model->HICUMvgc)/here->HICUMvt*(here->HICUMqtt0-1));
// } else {
// here->HICUMhfe_t = model->HICUMhfe;
// here->HICUMhfc_t = model->HICUMhfc;
// }
// here->HICUMrth_t = model->HICUMrth*exp(model->HICUMzetarth*here->HICUMln_qtt0)*(1+model->HICUMalrth*dT);
return(0);
}

13
src/spicelib/devices/hicum2/hicum2temp.hpp Normal file
View File

@ -0,0 +1,13 @@
#ifndef hicumL2_temp
#define hicumL2_temp
#include "hicum2defs.h"
#ifdef __cplusplus
extern "C" {
#endif
int HICUMtemp(GENmodel *inModel, CKTcircuit *ckt);
#ifdef __cplusplus
}
#endif
#endif /* hicumL2_temp */