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

Fixed temperature sweep.

This commit is contained in:
pnenzi 2003-09-20 18:05:27 +00:00
parent 92ba461127
commit bec68b1559
2 changed files with 32 additions and 19 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
ChangeLog
View File

@ -1,3 +1,9 @@
2003-09-18 Paolo Nenzi <p.nenzi@ieee.org>
* src/spicelib/frontend/dctrcv.c: Corrected temperature sweep.
Initial temp was the circuit temperature and not the start
value in the input line. There was a missing call to CKTtemp().
2003-08-24 Paolo Nenzi <p.nenzi@ieee.org> 2003-08-24 Paolo Nenzi <p.nenzi@ieee.org>
* src/maths/sparse/, src/include/complex.h: Added a function * src/maths/sparse/, src/include/complex.h: Added a function

45
src/spicelib/analysis/dctrcurv.c
View File

@ -5,7 +5,6 @@ Modified: 1999 Paolo Nenzi
**********/ **********/
#include "ngspice.h" #include "ngspice.h"
#include <stdio.h>
#include "vsrc/vsrcdefs.h" #include "vsrc/vsrcdefs.h"
#include "isrc/isrcdefs.h" #include "isrc/isrcdefs.h"
@ -88,7 +87,6 @@ DCtrCurv(CKTcircuit *ckt, int restart)
cv->TRCVvType[i] = rcode; cv->TRCVvType[i] = rcode;
here->RESresist = cv->TRCVvStart[i]; here->RESresist = cv->TRCVvStart[i];
here->RESresGiven = 1; here->RESresGiven = 1;
printf("** Resistor sweep is highly alpha code\n**Results may not be accurate.\n");
goto found; goto found;
} }
} }
@ -142,7 +140,7 @@ DCtrCurv(CKTcircuit *ckt, int restart)
cv->TRCVvSave[i]=ckt->CKTtemp; /* Saves the old circuit temperature */ cv->TRCVvSave[i]=ckt->CKTtemp; /* Saves the old circuit temperature */
cv->TRCVvType[i]=TEMP_CODE; /* Set the sweep type code */ cv->TRCVvType[i]=TEMP_CODE; /* Set the sweep type code */
ckt->CKTtemp = cv->TRCVvStart[i] + CONSTCtoK; /* Set the new circuit temp */ ckt->CKTtemp = cv->TRCVvStart[i] + CONSTCtoK; /* Set the new circuit temp */
printf("Temperature sweep is alpha code\nresults may not be accurate\n"); CKTtemp(ckt);
goto found; goto found;
} }
@ -158,22 +156,32 @@ found:;
error = CKTnames(ckt,&numNames,&nameList); error = CKTnames(ckt,&numNames,&nameList);
if(error) return(error); if(error) return(error);
if (cv->TRCVvType[i]==vcode) if (cv->TRCVvType[i]==vcode)
(*(SPfrontEnd->IFnewUid))((void *)ckt,&varUid,(IFuid )NULL, (*(SPfrontEnd->IFnewUid))((void *)ckt,&varUid,(IFuid )NULL,
"v-sweep", UID_OTHER, (void **)NULL); "v-sweep", UID_OTHER, (void **)NULL);
else if (cv->TRCVvType[i]==icode)
(*(SPfrontEnd->IFnewUid))((void *)ckt,&varUid,(IFuid )NULL, else {
"i-sweep", UID_OTHER, (void **)NULL); if (cv->TRCVvType[i]==icode)
else if (cv->TRCVvType[i]==TEMP_CODE) (*(SPfrontEnd->IFnewUid))((void *)ckt,&varUid,(IFuid )NULL,
(*(SPfrontEnd->IFnewUid))((void *)ckt,&varUid,(IFuid )NULL, "i-sweep", UID_OTHER, (void **)NULL);
"temp-sweep", UID_OTHER, (void **)NULL);
/* PN Resistance Sweep */ else {
else if (cv->TRCVvType[i]==rcode) if (cv->TRCVvType[i]==TEMP_CODE)
(*(SPfrontEnd->IFnewUid))((void *)ckt,&varUid,(IFuid )NULL, (*(SPfrontEnd->IFnewUid))((void *)ckt,&varUid,(IFuid )NULL,
"res-sweep", UID_OTHER, (void **)NULL); "temp-sweep", UID_OTHER, (void **)NULL);
else
(*(SPfrontEnd->IFnewUid))((void *)ckt,&varUid,(IFuid )NULL, else {
"?-sweep", UID_OTHER, (void **)NULL); if (cv->TRCVvType[i]==rcode)
(*(SPfrontEnd->IFnewUid))((void *)ckt,&varUid,(IFuid )NULL,
"res-sweep", UID_OTHER, (void **)NULL);
else
(*(SPfrontEnd->IFnewUid))((void *)ckt,&varUid,(IFuid )NULL,
"?-sweep", UID_OTHER, (void **)NULL);
} /* icode */
} /* TEMP_CODE */
} /* rcode*/
error = (*(SPfrontEnd->OUTpBeginPlot))((void *)ckt, error = (*(SPfrontEnd->OUTpBeginPlot))((void *)ckt,
(void*)ckt->CKTcurJob, ckt->CKTcurJob->JOBname, (void*)ckt->CKTcurJob, ckt->CKTcurJob->JOBname,
@ -218,8 +226,7 @@ resume:
} else if(cv->TRCVvType[i]==rcode) { /* resistance */ } else if(cv->TRCVvType[i]==rcode) { /* resistance */
if((((RESinstance*)(cv->TRCVvElt[i]))->RESresist)* if((((RESinstance*)(cv->TRCVvElt[i]))->RESresist)*
SIGN(1.,cv->TRCVvStep[i]) - SIGN(1.,cv->TRCVvStep[i]) -
SIGN(1.,cv->TRCVvStep[i]) * cv->TRCVvStop[i] > SIGN(1.,cv->TRCVvStep[i]) * cv->TRCVvStop[i] > 0.5 * fabs(cv->TRCVvStep[i]))
0.5 * fabs(cv->TRCVvStep[i]))
{ {
i++ ; i++ ;
firstTime=1; firstTime=1;