diff --git a/src/spicelib/devices/dio/dio.c b/src/spicelib/devices/dio/dio.c index 2ebeb2e41..4620efc71 100644 --- a/src/spicelib/devices/dio/dio.c +++ b/src/spicelib/devices/dio/dio.c @@ -122,6 +122,7 @@ IFparm DIOmPTable[] = { /* model parameters */ OPU( "cond", DIO_MOD_COND,IF_REAL, "Ohmic conductance"), IOP( "isr", DIO_MOD_ISR, IF_REAL, "Recombination saturation current"), IOP( "nr", DIO_MOD_NR, IF_REAL, "Recombination current emission coefficient"), + IOP( "vp", DIO_MOD_VP, IF_REAL, "Soft reverse recovery parameter"), /* SOA parameters */ IOPX( "fv_max", DIO_MOD_FV_MAX, IF_REAL, "maximum voltage in forward direction"), diff --git a/src/spicelib/devices/dio/dioacld.c b/src/spicelib/devices/dio/dioacld.c index 62ad65efd..56469d2b6 100644 --- a/src/spicelib/devices/dio/dioacld.c +++ b/src/spicelib/devices/dio/dioacld.c @@ -91,6 +91,20 @@ DIOacLoad(GENmodel *inModel, CKTcircuit *ckt) (*(here->DIOnegTempPtr) += -dIdioSw_dT); } } + if ((here->DIOqpNode > 0) && (model->DIOsoftRevRecParam!=0) && (here->DIOtTransitTime!=0)) { + double gdres = *(ckt->CKTstate0 + here->DIOresConduct); + /* QP subcircuit */ + double fac = here->DIOtTransitTime / model->DIOsoftRevRecParam; + double dcrrdvd = fac*gdres; + *(here->DIOqpQpPtr) += 1/model->DIOsoftRevRecParam; + *(here->DIOqpQpPtr + 1) += here->DIOtTransitTime * ckt->CKTomega; + *(here->DIOqpPosPrimePtr) += -dcrrdvd; + *(here->DIOqpNegPtr) += dcrrdvd; + /* Gain of VCVS (1-vp)/tau * j*omega*tau = (1-vp) * j*omega */ + double xgain = (1 - model->DIOsoftRevRecParam) * ckt->CKTomega; + *(here->DIOposPrimeQpPtr + 1) += xgain; + *(here->DIOnegQpPtr + 1) += -xgain; + } } } return(OK); diff --git a/src/spicelib/devices/dio/dioask.c b/src/spicelib/devices/dio/dioask.c index 67009b889..dc96ddc14 100644 --- a/src/spicelib/devices/dio/dioask.c +++ b/src/spicelib/devices/dio/dioask.c @@ -78,6 +78,8 @@ DIOask (CKTcircuit *ckt, GENinstance *inst, int which, IFvalue *value, return(OK); case DIO_CAP: value->rValue = here->DIOcap; + if ((here->DIOqpNode > 0) && (here->DIOtTransitTime!=0)) + value->rValue += here->DIOtTransitTime * *(ckt->CKTstate0+here->DIOconduct); return(OK); case DIO_CHARGE: value->rValue = *(ckt->CKTstate0+here->DIOcapCharge); diff --git a/src/spicelib/devices/dio/diobindCSC.c b/src/spicelib/devices/dio/diobindCSC.c index 29cee9cff..89994448e 100644 --- a/src/spicelib/devices/dio/diobindCSC.c +++ b/src/spicelib/devices/dio/diobindCSC.c @@ -54,6 +54,14 @@ DIObindCSC (GENmodel *inModel, CKTcircuit *ckt) CREATE_KLU_BINDING_TABLE(DIOposSwPrimeTempPtr, DIOposSwPrimeTempBinding, DIOposSwPrimeNode, DIOtempNode); } } + /* rev-rec */ + if (model->DIOsoftRevRecParamGiven && model->DIOsoftRevRecParam!=0 && model->DIOtransitTime!=0) { + CREATE_KLU_BINDING_TABLE(DIOqpQpPtr , DIOqpQpBinding , DIOqpNode, DIOqpNode); + CREATE_KLU_BINDING_TABLE(DIOqpPosPrimePtr, DIOqpPosPrimeBinding, DIOqpNode, DIOposPrimeNode); + CREATE_KLU_BINDING_TABLE(DIOqpNegPtr , DIOqpNegBinding , DIOqpNode, DIOnegNode); + CREATE_KLU_BINDING_TABLE(DIOposPrimeQpPtr, DIOposPrimeQpBinding, DIOposPrimeNode, DIOqpNode); + CREATE_KLU_BINDING_TABLE(DIOnegQpPtr , DIOnegQpBinding , DIOnegNode, DIOqpNode); + } } } @@ -103,6 +111,14 @@ DIObindCSCComplex (GENmodel *inModel, CKTcircuit *ckt) CONVERT_KLU_BINDING_TABLE_TO_COMPLEX(DIOposSwPrimeTempPtr, DIOposSwPrimeTempBinding, DIOposSwPrimeNode, DIOtempNode); } } + /* rev-rec */ + if (model->DIOsoftRevRecParamGiven && model->DIOsoftRevRecParam!=0 && model->DIOtransitTime!=0) { + CONVERT_KLU_BINDING_TABLE_TO_COMPLEX(DIOqpQpPtr , DIOqpQpBinding , DIOqpNode, DIOqpNode); + CONVERT_KLU_BINDING_TABLE_TO_COMPLEX(DIOqpPosPrimePtr, DIOqpPosPrimeBinding, DIOqpNode, DIOposPrimeNode); + CONVERT_KLU_BINDING_TABLE_TO_COMPLEX(DIOqpNegPtr , DIOqpNegBinding , DIOqpNode, DIOnegNode); + CONVERT_KLU_BINDING_TABLE_TO_COMPLEX(DIOposPrimeQpPtr, DIOposPrimeQpBinding, DIOposPrimeNode, DIOqpNode); + CONVERT_KLU_BINDING_TABLE_TO_COMPLEX(DIOnegQpPtr , DIOnegQpBinding , DIOnegNode, DIOqpNode); + } } } @@ -152,6 +168,14 @@ DIObindCSCComplexToReal (GENmodel *inModel, CKTcircuit *ckt) CONVERT_KLU_BINDING_TABLE_TO_REAL(DIOposSwPrimeTempPtr, DIOposSwPrimeTempBinding, DIOposSwPrimeNode, DIOtempNode); } } + /* rev-rec */ + if (model->DIOsoftRevRecParamGiven && model->DIOsoftRevRecParam!=0 && model->DIOtransitTime!=0) { + CONVERT_KLU_BINDING_TABLE_TO_REAL(DIOqpQpPtr , DIOqpQpBinding , DIOqpNode, DIOqpNode); + CONVERT_KLU_BINDING_TABLE_TO_REAL(DIOqpPosPrimePtr, DIOqpPosPrimeBinding, DIOqpNode, DIOposPrimeNode); + CONVERT_KLU_BINDING_TABLE_TO_REAL(DIOqpNegPtr , DIOqpNegBinding , DIOqpNode, DIOnegNode); + CONVERT_KLU_BINDING_TABLE_TO_REAL(DIOposPrimeQpPtr, DIOposPrimeQpBinding, DIOposPrimeNode, DIOqpNode); + CONVERT_KLU_BINDING_TABLE_TO_REAL(DIOnegQpPtr , DIOnegQpBinding , DIOnegNode, DIOqpNode); + } } } diff --git a/src/spicelib/devices/dio/diodefs.h b/src/spicelib/devices/dio/diodefs.h index 298f330b9..c2637310b 100644 --- a/src/spicelib/devices/dio/diodefs.h +++ b/src/spicelib/devices/dio/diodefs.h @@ -44,6 +44,7 @@ typedef struct sDIOinstance { const int DIOtempNode; /* number of the temperature node of the diode */ int DIOposPrimeNode; /* number of positive prime node of diode */ int DIOposSwPrimeNode; /* number of positive prime node of diode sidewall */ + int DIOqpNode; /* number of soft recovery charge node */ double *DIOposPosPrimePtr; /* pointer to sparse matrix at * (positive,positive prime) */ @@ -83,6 +84,13 @@ typedef struct sDIOinstance { double *DIOtempPosSwPrimePtr; double *DIOposSwPrimeTempPtr; + /* rev-rec */ + double *DIOqpQpPtr; + double *DIOqpPosPrimePtr; + double *DIOqpNegPtr; + double *DIOposPrimeQpPtr; + double *DIOnegQpPtr; + double DIOcap; /* stores the diode capacitance */ double DIOcapSW; /* stores the diode Sw capacitance */ @@ -229,13 +237,19 @@ typedef struct sDIOinstance { /* self heating */ BindElement *DIOtempPosBinding; BindElement *DIOtempPosPrimeBinding; - BindElement *DIOtempNegBinding; + BindElement *DIOtempNegBinding; BindElement *DIOtempTempBinding; BindElement *DIOposTempBinding; BindElement *DIOposPrimeTempBinding; BindElement *DIOnegTempBinding; BindElement *DIOtempPosSwPrimeBinding; BindElement *DIOposSwPrimeTempBinding; + /* rev-rec */ + BindElement *DIOqpQpBinding; + BindElement *DIOqpPosPrimeBinding; + BindElement *DIOqpNegBinding; + BindElement *DIOposPrimeQpBinding; + BindElement *DIOnegQpBinding; #endif } DIOinstance ; @@ -244,7 +258,6 @@ typedef struct sDIOinstance { #define DIOsenCeq DIOsens + 3 /* stores the perturbed values of ceq */ #define DIOdphidp DIOsens + 6 - #define DIOvoltage DIOstate #define DIOcurrent DIOstate+1 #define DIOconduct DIOstate+2 @@ -256,18 +269,25 @@ typedef struct sDIOinstance { #define DIOcapChargeSW DIOstate+8 #define DIOcapCurrentSW DIOstate+9 -#define DIOqth DIOstate+10 /* thermal capacitor charge */ -#define DIOcqth DIOstate+11 /* thermal capacitor current */ - -#define DIOdeltemp DIOstate+12 /* thermal voltage over rth0 */ +#define DIOqth DIOstate+10 /* thermal capacitor charge */ +#define DIOcqth DIOstate+11 /* thermal capacitor current */ +#define DIOdeltemp DIOstate+12 /* thermal voltage over rth0 */ #define DIOdIdio_dT DIOstate+13 #define DIOdIdioSW_dT DIOstate+14 +/* rev-rec */ +#define DIOsrcapCharge DIOstate+15 +#define DIOsrcapCurrent DIOstate+16 +#define DIOqp DIOstate+17 +#define DIOresCurrent DIOstate+18 +#define DIOresConduct DIOstate+19 +#define DIOcqcsr DIOstate+20 +#define DIOgqcsr DIOstate+21 -#define DIOnumStates 15 +#define DIOnumStates 22 -#define DIOsensxp DIOstate+15 /* charge sensitivities and their derivatives. - * +16 for the derivatives - pointer to the - * beginning of the array */ +#define DIOsensxp DIOstate+21 /* charge sensitivities and their derivatives. + * +22 for the derivatives - pointer to the + * beginning of the array */ #define DIOnumSenStates 2 @@ -342,6 +362,7 @@ typedef struct sDIOmodel { /* model structure for a diode */ unsigned DIOte_maxGiven : 1; unsigned DIOrecSatCurGiven : 1; unsigned DIOrecEmissionCoeffGiven : 1; + unsigned DIOsoftRevRecParamGiven : 1; unsigned DIOrth0Given :1; unsigned DIOcth0Given :1; @@ -418,6 +439,7 @@ typedef struct sDIOmodel { /* model structure for a diode */ double DIOte_max; /* maximum temperature */ double DIOrecSatCur; /* Recombination saturation current */ double DIOrecEmissionCoeff; /* Recombination emission coefficient */ + double DIOsoftRevRecParam; /* Soft reverse recovery parameter */ double DIOrth0; double DIOcth0; @@ -526,6 +548,7 @@ enum { DIO_MOD_PD_MAX, DIO_MOD_ISR, DIO_MOD_NR, + DIO_MOD_VP, DIO_MOD_RTH0, DIO_MOD_CTH0, diff --git a/src/spicelib/devices/dio/dioload.c b/src/spicelib/devices/dio/dioload.c index 6f53966eb..513aa42dd 100644 --- a/src/spicelib/devices/dio/dioload.c +++ b/src/spicelib/devices/dio/dioload.c @@ -26,6 +26,7 @@ DIOload(GENmodel *inModel, CKTcircuit *ckt) double argsw; double capd, capdsw=0.0; double cd, cdb, cdsw, cdb_dT, cdsw_dT; + double cdres, gdres; double cdeq; double cdhat, cdhatsw=0.0; double ceq; @@ -43,6 +44,7 @@ DIOload(GENmodel *inModel, CKTcircuit *ckt) double evd; double evrev; double gd, gdb, gdsw, gen_fac, gen_fac_vd; + double capsr, gqcsr, cqcsr; double t1, evd_rec, cdb_rec, gdb_rec, cdb_rec_dT; double geq; double gspr; /* area-scaled conductance */ @@ -52,6 +54,7 @@ DIOload(GENmodel *inModel, CKTcircuit *ckt) double tol; /* temporary for tolerence calculations */ #endif double vd, vdsw=0.0; /* current diode voltage */ + double vqp; double vdtemp; double vt; /* K t / Q */ double vte, vtesw, vtetun, vtebrk, vterec; @@ -74,6 +77,7 @@ DIOload(GENmodel *inModel, CKTcircuit *ckt) here=DIOnextInstance(here)) { int selfheat = ((here->DIOtempNode > 0) && (here->DIOthermal) && (model->DIOrth0Given)); + int revrec = ((here->DIOqpNode > 0) && (model->DIOsoftRevRecParam!=0) && (here->DIOtTransitTime!=0)); /* * this routine loads diodes for dc and transient analyses. @@ -120,10 +124,12 @@ DIOload(GENmodel *inModel, CKTcircuit *ckt) vd = *(ckt->CKTstate1 + here->DIOvoltage); if (model->DIOresistSWGiven) vdsw = *(ckt->CKTstate1 + here->DIOvoltageSW); delTemp = *(ckt->CKTstate1 + here->DIOdeltemp); + vqp = *(ckt->CKTstate1 + here->DIOqp); } else{ vd = *(ckt->CKTstate0 + here->DIOvoltage); if (model->DIOresistSWGiven) vdsw = *(ckt->CKTstate0 + here->DIOvoltageSW); delTemp = *(ckt->CKTstate0 + here->DIOdeltemp); + vqp = *(ckt->CKTstate0 + here->DIOqp); } #ifdef SENSDEBUG @@ -137,24 +143,30 @@ DIOload(GENmodel *inModel, CKTcircuit *ckt) vd= *(ckt->CKTstate0 + here->DIOvoltage); if (model->DIOresistSWGiven) vdsw = *(ckt->CKTstate0 + here->DIOvoltageSW); delTemp = *(ckt->CKTstate0 + here->DIOdeltemp); + vqp= *(ckt->CKTstate0 + here->DIOqp); } else if (ckt->CKTmode & MODEINITTRAN) { vd= *(ckt->CKTstate1 + here->DIOvoltage); if (model->DIOresistSWGiven) vdsw = *(ckt->CKTstate1 + here->DIOvoltageSW); delTemp = *(ckt->CKTstate1 + here->DIOdeltemp); + vqp= *(ckt->CKTstate1 + here->DIOqp); } else if ( (ckt->CKTmode & MODEINITJCT) && (ckt->CKTmode & MODETRANOP) && (ckt->CKTmode & MODEUIC) ) { vd=here->DIOinitCond; if (model->DIOresistSWGiven) vdsw = here->DIOinitCond; + vqp=0; } else if ( (ckt->CKTmode & MODEINITJCT) && here->DIOoff) { vd=vdsw=0; delTemp = 0.0; + vqp=0; } else if ( ckt->CKTmode & MODEINITJCT) { vd=here->DIOtVcrit; vdsw=here->DIOtVcritSW; delTemp = 0.0; + vqp=0; } else if ( ckt->CKTmode & MODEINITFIX && here->DIOoff) { vd=vdsw=0; delTemp = 0.0; + vqp=0; } else { #ifndef PREDICTOR if (ckt->CKTmode & MODEINITPRED) { @@ -177,12 +189,22 @@ DIOload(GENmodel *inModel, CKTcircuit *ckt) *(ckt->CKTstate0 + here->DIOdIdioSW_dT) = *(ckt->CKTstate1 + here->DIOdIdioSW_dT); } + vqp = DEVpred(ckt,here->DIOqp); + *(ckt->CKTstate0 + here->DIOresCurrent) = + *(ckt->CKTstate1 + here->DIOresCurrent); + *(ckt->CKTstate0 + here->DIOresConduct) = + *(ckt->CKTstate1 + here->DIOresConduct); + *(ckt->CKTstate0 + here->DIOcqcsr) = + *(ckt->CKTstate1 + here->DIOcqcsr); + *(ckt->CKTstate0 + here->DIOgqcsr) = + *(ckt->CKTstate1 + here->DIOgqcsr); } else { #endif /* PREDICTOR */ vd = *(ckt->CKTrhsOld+here->DIOposPrimeNode)- *(ckt->CKTrhsOld + here->DIOnegNode); if (model->DIOresistSWGiven) vdsw = *(ckt->CKTrhsOld+here->DIOposSwPrimeNode)- *(ckt->CKTrhsOld + here->DIOnegNode); + if (selfheat) delTemp = *(ckt->CKTrhsOld + here->DIOtempNode); else @@ -192,6 +214,7 @@ DIOload(GENmodel *inModel, CKTcircuit *ckt) *(ckt->CKTstate1+here->DIOqth) = *(ckt->CKTstate0+here->DIOqth); } + vqp = *(ckt->CKTrhsOld+here->DIOqpNode); #ifndef PREDICTOR } #endif /* PREDICTOR */ @@ -238,6 +261,11 @@ DIOload(GENmodel *inModel, CKTcircuit *ckt) gdsw= *(ckt->CKTstate0 + here->DIOconductSW); dIdioSw_dT= *(ckt->CKTstate0 + here->DIOdIdioSW_dT); } + vqp= *(ckt->CKTstate0 + here->DIOqp); + cdres= *(ckt->CKTstate0 + here->DIOresCurrent); + gdres= *(ckt->CKTstate0 + here->DIOresConduct); + cqcsr= *(ckt->CKTstate0 + here->DIOcqcsr); + gqcsr= *(ckt->CKTstate0 + here->DIOgqcsr); goto load; } } @@ -304,7 +332,7 @@ next1: if (model->DIOsatSWCurGiven) { /* sidewall current */ double vds; - if (model->DIOresistSWGiven) + if (model->DIOresistSWGiven) vds = vdsw; /* sidewall voltage used */ else vds = vd; /* common voltage used */ @@ -470,7 +498,6 @@ next1: gdb = ((1+sqrt_ikx)*gdb + cdb*gdb/(2*sqrt_ikx*ikr_area_m))/(1+2*sqrt_ikx - cdb/ikr_area_m); cdb = cdb/(1+sqrt_ikx); } - } if ( (model->DIOforwardSWKneeCurrentGiven) && (cdsw > 1.0e-18) ) { @@ -493,6 +520,11 @@ next1: dIdioSw_dT = cdsw_dT; } + cdres = cd; + gdres = gd; + cqcsr = 0; + gqcsr = 0; + if ((ckt->CKTmode & (MODEDCTRANCURVE | MODETRAN | MODEAC | MODEINITSMSIG)) || ((ckt->CKTmode & MODETRANOP) && (ckt->CKTmode & MODEUIC))) { /* @@ -528,22 +560,47 @@ next1: deplcapSW = czof2SW*(here->DIOtF3SW+model->DIOgradingSWCoeff*vdx/here->DIOtJctSWPot); } - diffcharge = here->DIOtTransitTime*cd; - diffcap = here->DIOtTransitTime*gd; - if (!model->DIOresistSWGiven) { + if (revrec) { + /* + soft recovery with TT!=0 + add only depletion capacitance. + */ *(ckt->CKTstate0 + here->DIOcapCharge) = - diffcharge + deplcharge + deplchargeSW + (here->DIOcmetal + here->DIOcpoly)*vd; - capd = diffcap + deplcap + deplcapSW + here->DIOcmetal + here->DIOcpoly; + deplcharge + deplchargeSW + (here->DIOcmetal + here->DIOcpoly)*vd; + + capd = deplcap + deplcapSW + here->DIOcmetal + here->DIOcpoly; here->DIOcap = capd; + /* + DIOcap is now equal only to depletion capacitance + overlap capacitance. + Diffusion capacitance is modelled via Qp so there is no clear way to define it. + Situation is similar to the one when we have an NQS model for the charge. + */ + + /* Now prepare the charge for the capacitor connected to the QP node */ + *(ckt->CKTstate0 + here->DIOsrcapCharge) = here->DIOtTransitTime * vqp; + capsr = here->DIOtTransitTime; } else { - *(ckt->CKTstate0 + here->DIOcapCharge) = - diffcharge + deplcharge + (here->DIOcmetal + here->DIOcpoly)*vd; - capd = diffcap + deplcap + here->DIOcmetal + here->DIOcpoly; - here->DIOcap = capd; - *(ckt->CKTstate0 + here->DIOcapChargeSW) = - deplcapSW; - capdsw = deplcapSW; - here->DIOcapSW = capdsw; + /* no soft recovery of soft recovery with TT=0 (i.e. no soft recovery due to TT=0) */ + diffcharge = here->DIOtTransitTime*cd; + diffcap = here->DIOtTransitTime*gd; + if (!model->DIOresistSWGiven) { + *(ckt->CKTstate0 + here->DIOcapCharge) = + diffcharge + deplcharge + deplchargeSW + (here->DIOcmetal + here->DIOcpoly)*vd; + capd = diffcap + deplcap + deplcapSW + here->DIOcmetal + here->DIOcpoly; + here->DIOcap = capd; + } else { + *(ckt->CKTstate0 + here->DIOcapCharge) = + diffcharge + deplcharge + (here->DIOcmetal + here->DIOcpoly)*vd; + capd = diffcap + deplcap + here->DIOcmetal + here->DIOcpoly; + here->DIOcap = capd; + *(ckt->CKTstate0 + here->DIOcapChargeSW) = + deplcapSW; + capdsw = deplcapSW; + here->DIOcapSW = capdsw; + } + + *(ckt->CKTstate0 + here->DIOsrcapCharge) = 0; + capsr = 0; } /* * store small-signal parameters @@ -564,6 +621,10 @@ next1: *(ckt->CKTstate0 + here->DIOconductSW) = gdsw; *(ckt->CKTstate0 + here->DIOdIdioSW_dT) = dIdioSw_dT; } + *(ckt->CKTstate0 + here->DIOresCurrent) = cdres; + *(ckt->CKTstate0 + here->DIOresConduct) = gdres; + *(ckt->CKTstate0 + here->DIOcqcsr) = cqcsr; + *(ckt->CKTstate0 + here->DIOgqcsr) = gqcsr; #ifdef SENSDEBUG printf("storing small signal parameters\n"); printf("cd = %.7e,vd = %.7e\n",cd,vd); @@ -580,6 +641,7 @@ next1: *(ckt->CKTstate0 + here->DIOcurrent) = cd; if (model->DIOresistSWGiven) *(ckt->CKTstate0 + here->DIOcurrentSW) = cdsw; + *(ckt->CKTstate0 + here->DIOresCurrent) = cdres; #ifdef SENSDEBUG printf("storing parameters for transient sensitivity\n" ); @@ -613,6 +675,21 @@ next1: *(ckt->CKTstate1 + here->DIOcapCurrentSW) = *(ckt->CKTstate0 + here->DIOcapCurrentSW); } + if (revrec) { + /* soft recovery subcircuit */ + if (ckt->CKTmode & MODEINITTRAN) { + *(ckt->CKTstate1 + here->DIOsrcapCharge) = + *(ckt->CKTstate0 + here->DIOsrcapCharge); + } + error = NIintegrate(ckt,&geq,&ceq,capsr,here->DIOsrcapCharge); + if(error) return(error); + gqcsr = geq; + cqcsr = *(ckt->CKTstate0 + here->DIOsrcapCurrent); + if (ckt->CKTmode & MODEINITTRAN) { + *(ckt->CKTstate1 + here->DIOsrcapCurrent) = + *(ckt->CKTstate0 + here->DIOsrcapCurrent); + } + } if (selfheat) { error = NIintegrate(ckt, &gcTt, &ceqqth, model->DIOcth0, here->DIOqth); @@ -654,6 +731,11 @@ next2: *(ckt->CKTstate0 + here->DIOvoltage) = vd; *(ckt->CKTstate0 + here->DIOconductSW) = gdsw; *(ckt->CKTstate0 + here->DIOdIdioSW_dT) = dIdioSw_dT; } + *(ckt->CKTstate0 + here->DIOqp) = vqp; + *(ckt->CKTstate0 + here->DIOresCurrent) = cdres; + *(ckt->CKTstate0 + here->DIOresConduct) = gdres; + *(ckt->CKTstate0 + here->DIOcqcsr) = cqcsr; + *(ckt->CKTstate0 + here->DIOgqcsr) = gqcsr; if(SenCond) continue; #ifndef NOBYPASS @@ -754,6 +836,29 @@ next2: *(ckt->CKTstate0 + here->DIOvoltage) = vd; (*(here->DIOnegTempPtr) += -dIdioSw_dT); } } + + if (revrec) { + double fac, ceqrr, dcrrdvd, grr; + double gain, ceqrrd, geqrrd; + /* QP subcircuit */ + fac = here->DIOtTransitTime / model->DIOsoftRevRecParam; + dcrrdvd = fac*gdres; + ceqrr = -fac*cdres + cqcsr + dcrrdvd*vd - gqcsr*vqp; + grr = 1/model->DIOsoftRevRecParam; + *(ckt->CKTrhs + here->DIOqpNode) -= ceqrr; + *(here->DIOqpQpPtr) += grr + gqcsr; + *(here->DIOqpPosPrimePtr) += -dcrrdvd; + *(here->DIOqpNegPtr) += dcrrdvd; + /* Contribution to diode current */ + gain = (1 - model->DIOsoftRevRecParam) / here->DIOtTransitTime; + /* Linear contribution -(1-vp)/tau*ddt(Qp) */ + geqrrd = gain*gqcsr; + ceqrrd = gain*cqcsr - geqrrd*vqp; + *(ckt->CKTrhs + here->DIOposPrimeNode) -= ceqrrd; + *(ckt->CKTrhs + here->DIOnegNode) += ceqrrd; + *(here->DIOposPrimeQpPtr) += geqrrd; + *(here->DIOnegQpPtr) += -geqrrd; + } } } return(OK); diff --git a/src/spicelib/devices/dio/diomask.c b/src/spicelib/devices/dio/diomask.c index 1f4e11ce5..b4f35d4af 100644 --- a/src/spicelib/devices/dio/diomask.c +++ b/src/spicelib/devices/dio/diomask.c @@ -199,6 +199,9 @@ DIOmAsk (CKTcircuit *ckt, GENmodel *inModel, int which, IFvalue *value) case DIO_MOD_NR: value->rValue = model->DIOrecEmissionCoeff; return(OK); + case DIO_MOD_VP: + value->rValue = model->DIOsoftRevRecParam; + return(OK); case DIO_MOD_RTH0: value->rValue = model->DIOrth0; return(OK); diff --git a/src/spicelib/devices/dio/diompar.c b/src/spicelib/devices/dio/diompar.c index 498421160..164d0186b 100644 --- a/src/spicelib/devices/dio/diompar.c +++ b/src/spicelib/devices/dio/diompar.c @@ -245,6 +245,10 @@ DIOmParam(int param, IFvalue *value, GENmodel *inModel) model->DIOrecEmissionCoeff = value->rValue; model->DIOrecEmissionCoeffGiven = TRUE; break; + case DIO_MOD_VP: + model->DIOsoftRevRecParam = value->rValue; + model->DIOsoftRevRecParamGiven = TRUE; + break; case DIO_MOD_RTH0: model->DIOrth0 = value->rValue; model->DIOrth0Given = TRUE; diff --git a/src/spicelib/devices/dio/diosetup.c b/src/spicelib/devices/dio/diosetup.c index df9555be8..2cccd6144 100644 --- a/src/spicelib/devices/dio/diosetup.c +++ b/src/spicelib/devices/dio/diosetup.c @@ -222,6 +222,9 @@ DIOsetup(SMPmatrix *matrix, GENmodel *inModel, CKTcircuit *ckt, int *states) if(!model->DIOrecSatCurGiven) { model->DIOrecSatCur = 1e-14; } + if (!model->DIOsoftRevRecParamGiven) { + model->DIOsoftRevRecParam = 0.0; + } /* set lower limit of saturation current */ if (model->DIOsatCur < ckt->CKTepsmin) @@ -412,6 +415,18 @@ DIOsetup(SMPmatrix *matrix, GENmodel *inModel, CKTcircuit *ckt, int *states) } } + /* rev-rec */ + if (model->DIOsoftRevRecParamGiven && model->DIOsoftRevRecParam!=0 && model->DIOtransitTime!=0) { + if(here->DIOqpNode == 0) { + error = CKTmkVolt(ckt, &tmp, here->DIOname, "qp"); + if(error) return(error); + here->DIOqpNode = tmp->number; + } + } else { + here->DIOqpNode = 0; + } + + int selfheat = ((here->DIOtempNode > 0) && (here->DIOthermal) && (model->DIOrth0Given)); /* macro to make elements with built in test for out of memory */ @@ -450,7 +465,15 @@ do { if((here->ptr = SMPmakeElt(matrix, here->first, here->second)) == NULL){\ TSTALLOC(DIOposSwPrimeTempPtr, DIOposSwPrimeNode, DIOtempNode); } } - + + /* rev-rec */ + if (model->DIOsoftRevRecParamGiven && model->DIOsoftRevRecParam!=0 && model->DIOtransitTime!=0) { + TSTALLOC(DIOqpQpPtr , DIOqpNode, DIOqpNode); + TSTALLOC(DIOqpPosPrimePtr, DIOqpNode, DIOposPrimeNode); + TSTALLOC(DIOqpNegPtr , DIOqpNode, DIOnegNode); + TSTALLOC(DIOposPrimeQpPtr, DIOposPrimeNode, DIOqpNode); + TSTALLOC(DIOnegQpPtr, DIOnegNode, DIOqpNode); + } } } return(OK); @@ -474,7 +497,7 @@ DIOunsetup( if (here->DIOposPrimeNode > 0 && here->DIOposPrimeNode != here->DIOposNode) CKTdltNNum(ckt, here->DIOposPrimeNode); - here->DIOposPrimeNode = 0; + here->DIOposPrimeNode = 0; if(model->DIOresistSWGiven) { /* separate sidewall */ @@ -484,6 +507,11 @@ DIOunsetup( here->DIOposSwPrimeNode = 0; } + /* rev-rec */ + if (here->DIOqpNode > 0) + CKTdltNNum(ckt, here->DIOqpNode); + here->DIOqpNode = 0; + } } return OK; diff --git a/src/spicelib/devices/dio/diotrunc.c b/src/spicelib/devices/dio/diotrunc.c index b940528be..04070a1d8 100644 --- a/src/spicelib/devices/dio/diotrunc.c +++ b/src/spicelib/devices/dio/diotrunc.c @@ -22,6 +22,8 @@ DIOtrunc(GENmodel *inModel, CKTcircuit *ckt, double *timeStep) for(here=DIOinstances(model);here!=NULL;here = DIOnextInstance(here)){ CKTterr(here->DIOcapCharge,ckt,timeStep); if (model->DIOresistSWGiven) CKTterr(here->DIOcapChargeSW,ckt,timeStep); + if (model->DIOsoftRevRecParam!=0 && here->DIOtTransitTime!=0) + CKTterr(here->DIOsrcapCharge,ckt,timeStep); } } return(OK);