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

Some fixes for the Reliability Analysis

This commit is contained in:
Francesco Lannutti 2016-06-09 18:27:33 +02:00
parent a0bb42580c
commit a3534e2954
4 changed files with 140 additions and 71 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

198
src/spicelib/devices/bsim4/b4reliability.c
View File

@ -164,89 +164,142 @@ BSIM4reliability_internal (BSIM4instance *here, CKTcircuit *ckt, unsigned int mo
}
}
} else if (mode == 1) {
// In this mode, it doesn't matter if NOW the device is in stress or in recovery, since it's the latest timestep
if (here->relStruct->IsON == 1)
if (NowIsON)
{
// Calculate stress
delta = ckt->CKTtime - here->relStruct->time ;
// Update time and flag - Maybe Optional
here->relStruct->time = ckt->CKTtime ;
here->relStruct->IsON = 1 ;
// Calculate Aging - Giorgio Liatis' Model
ret = RELMODELcalculateAging ((GENinstance *)here, here->BSIM4modPtr->BSIM4modType, delta, 1e-12, 1) ;
if (ret == 1)
if (here->relStruct->IsON == 1)
{
return (E_INTERN) ;
// Until now, the device was ON - Calculate stress
delta = ckt->CKTtime - here->relStruct->time ;
// Update time and flag - Recovery begins
here->relStruct->time = ckt->CKTtime ;
here->relStruct->IsON = 1 ;
// Calculate Aging - Giorgio Liatis' Model
ret = RELMODELcalculateAging ((GENinstance *)here, here->BSIM4modPtr->BSIM4modType, delta, 1e-12, 1) ;
if (ret == 1)
{
return (E_INTERN) ;
}
// Update the semiperiod counter
here->relStruct->semiPeriods++ ;
} else if (here->relStruct->IsON == 0) {
// Until now, the device was OFF - Calculate recovery
delta = ckt->CKTtime - here->relStruct->time ;
// Update time and flag - Stress begins
here->relStruct->time = ckt->CKTtime ;
here->relStruct->IsON = 1 ;
// Calculate Aging - Giogio Liatis' Model
ret = RELMODELcalculateAging ((GENinstance *)here, here->BSIM4modPtr->BSIM4modType, delta, 1e-12, 0) ;
if (ret == 1)
{
return (E_INTERN) ;
}
// Update the semiperiod counter
here->relStruct->semiPeriods++ ;
} else {
fprintf (stderr, "Reliability Analysis Error\n") ;
}
// Update the semiperiod counter
here->relStruct->semiPeriods++ ;
} else if (here->relStruct->IsON == 0) {
// Calculate recovery
delta = ckt->CKTtime - here->relStruct->time ;
// Update time and flag - Maybe Optional
here->relStruct->time = ckt->CKTtime ;
here->relStruct->IsON = 0 ;
// Calculate Aging - Giorgio Liatis' Model
ret = RELMODELcalculateAging ((GENinstance *)here, here->BSIM4modPtr->BSIM4modType, delta, 1e-12, 0) ;
if (ret == 1)
{
return (E_INTERN) ;
}
// Update the semiperiod counter
here->relStruct->semiPeriods++ ;
} else {
fprintf (stderr, "Reliability Analysis Error\n") ;
if (here->relStruct->IsON == 1)
{
// Until now, the device was ON - Calculate stress
delta = ckt->CKTtime - here->relStruct->time ;
// Update time and flag - Recovery begins
here->relStruct->time = ckt->CKTtime ;
here->relStruct->IsON = 0 ;
// Calculate Aging - Giorgio Liatis' Model
ret = RELMODELcalculateAging ((GENinstance *)here, here->BSIM4modPtr->BSIM4modType, delta, 1e-12, 1) ;
if (ret == 1)
{
return (E_INTERN) ;
}
// Update the semiperiod counter
here->relStruct->semiPeriods++ ;
} else if (here->relStruct->IsON == 0) {
// Until now, the device was OFF - Calculate recovery
delta = ckt->CKTtime - here->relStruct->time ;
// Update time and flag - Stress begins
here->relStruct->time = ckt->CKTtime ;
here->relStruct->IsON = 0 ;
// Calculate Aging - Giogio Liatis' Model
ret = RELMODELcalculateAging ((GENinstance *)here, here->BSIM4modPtr->BSIM4modType, delta, 1e-12, 0) ;
if (ret == 1)
{
return (E_INTERN) ;
}
// Update the semiperiod counter
here->relStruct->semiPeriods++ ;
} else {
fprintf (stderr, "Reliability Analysis Error\n") ;
}
}
if (here->relStruct->deltaVth > 0) {
model->number_of_aged_instances++ ;
model->total_deltaVth += here->relStruct->deltaVth * 1000 ;
printf ("DEVICE OK!!!\tTime: %-.9gs\t\t", ckt->CKTtime) ;
printf ("DeltaVth: %-.9gmV\t\t", here->relStruct->deltaVth * 1000) ;
printf ("Device Name: %s\t\t", here->BSIM4name) ;
printf ("Device Type: %s\n\n", model->BSIM4modName) ;
} else if (here->relStruct->deltaVth > 0) {
printf ("\n\n\n\nWarning: PROBLEMATIC DEVICE!!!\tTime: %-.9gs\t\t", ckt->CKTtime) ;
printf ("DeltaVth: %-.9gmV\t\t", here->relStruct->deltaVth * 1000) ;
printf ("Device Name: %s\t\t", here->BSIM4name) ;
printf ("Device Type: %s\n\n\n\n\n\n", model->BSIM4modName) ;
} else {
printf ("THIS DEVICE IS OFF!!!\tTime: %-.9gs\t\t", ckt->CKTtime) ;
printf ("DeltaVth: %-.9gmV\t\t", here->relStruct->deltaVth * 1000) ;
printf ("Device Name: %s\t\t", here->BSIM4name) ;
printf ("Device Type: %s\n\n", model->BSIM4modName) ;
}
printf ("\tTime: %-.9gs\t\t", ckt->CKTtime) ;
printf ("DeltaVth: %-.9gmV\t\t", here->relStruct->deltaVth * 1000) ;
printf ("Device Name: %s\t\t", here->BSIM4name) ;
printf ("Device Type: %s\n\n", model->BSIM4modName) ;
/* Calculate fitting */
/* Count how many deltaVth we have */
unsigned int i ;
RELMODELrelList *current ;
i = 0 ;
current = here->relStruct->deltaVthList ;
while (current != NULL)
{
i++ ;
current = current->next ;
}
/* Assign list members to vectors */
double *timeFit, *deltaVthFit ;
timeFit = TMALLOC (double, i) ;
deltaVthFit = TMALLOC (double, i) ;
i = 0 ;
current = here->relStruct->deltaVthList ;
while (current != NULL)
{
timeFit [i] = current->time ;
deltaVthFit [i] = current->deltaVth ;
i++ ;
current = current->next ;
}
if (here->relStruct->semiPeriods > 1)
{
/* The model behavior is periodic - Use Fourier basis fitting */
double *deltaVthFit, f, factor_for_2pi, *fitting_matrix, target, *timeFit ;
RELMODELrelList *current ;
unsigned int columns, i, j, number_of_modes, number_of_periods, rows, size ;
/* Count how many deltaVth we have */
i = 0 ;
current = here->relStruct->deltaVthList ;
while (current != NULL)
{
i++ ;
current = current->next ;
}
/* Assign list members to vectors */
timeFit = TMALLOC (double, i) ;
deltaVthFit = TMALLOC (double, i) ;
i = 0 ;
current = here->relStruct->deltaVthList ;
while (current != NULL)
{
timeFit [i] = current->time ;
deltaVthFit [i] = current->deltaVth ;
i++ ;
current = current->next ;
}
/* Generate the fitting matrix */
double f, factor_for_2pi, *fitting_matrix, target ;
unsigned int columns, j, number_of_modes, number_of_periods, rows, size ;
number_of_periods = here->relStruct->semiPeriods / 2 ;
number_of_modes = 10 * (number_of_periods + 1) ;
@ -358,6 +411,11 @@ BSIM4reliability (GENmodel *inModel, CKTcircuit *ckt, unsigned int mode)
BSIM4reliability_internal (here, ckt, mode) ;
}
}
if (mode == 1) {
printf ("Number of aged BSIM4 PMOS instances: %u\n\twith a mean DeltaVth of: %-.9gmV\n\n\n",
model->number_of_aged_instances, model->total_deltaVth / model->number_of_aged_instances) ;
}
}
}

6
src/spicelib/devices/bsim4/b4set.c
View File

@ -2671,6 +2671,12 @@ do { if((here->ptr = SMPmakeElt(matrix,here->first,here->second))==(double *)NUL
#endif
}
#ifdef RELAN
model->number_of_aged_instances = 0 ;
model->total_deltaVth = 0 ;
#endif
}
#ifdef USE_OMP

2
src/spicelib/devices/bsim4/b4temp.c
View File

@ -1721,7 +1721,7 @@ int Size_Not_Found, i;
{
if (here->relStruct->IsON != -1)
{
here->BSIM4delvto = here->relStruct->deltaVth ;
here->BSIM4delvto = model->BSIM4type * here->relStruct->deltaVth ;
}
}
#endif

5
src/spicelib/devices/bsim4/bsim4def.h
View File

@ -2723,6 +2723,11 @@ typedef struct sBSIM4model
unsigned BSIM4pku0weGiven :1;
#ifdef RELAN
unsigned int number_of_aged_instances ;
double total_deltaVth ;
#endif
} BSIM4model;