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

first try dctrcurv.c

This commit is contained in:
rlar 2018-04-28 19:07:07 +02:00
parent 958e209866
commit d8611aede8
1 changed files with 34 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -28,6 +28,10 @@ Modified: 1999 Paolo Nenzi
static double actval, actdiff;
#endif
#ifdef USE_CUSPICE
#include "ngspice/CUSPICE/CUSPICE.h"
#endif
int
DCtrCurv(CKTcircuit *ckt, int restart)
@ -35,7 +39,11 @@ DCtrCurv(CKTcircuit *ckt, int restart)
TRCV *job = (TRCV *) ckt->CKTcurJob;
int i;
#ifdef USE_CUSPICE
int status;
#else
double *temp;
#endif
int converged;
int rcode;
int vcode;
@ -285,11 +293,17 @@ DCtrCurv(CKTcircuit *ckt, int restart)
DEVices[rcode]->DEVload(job->TRCVvElt[i]->GENmodPtr, ckt);
}
#ifdef USE_CUSPICE
status = cuCKTstatesCircularBuffer(ckt);
if (status != 0)
return (E_NOMEM);
#else
/* Rotate state vectors. */
temp = ckt->CKTstates[ckt->CKTmaxOrder + 1];
for (j = ckt->CKTmaxOrder; j >= 0; j--)
ckt->CKTstates[j + 1] = ckt->CKTstates[j];
ckt->CKTstate0 = temp;
#endif
/* do operation */
#ifdef XSPICE
@ -444,8 +458,14 @@ DCtrCurv(CKTcircuit *ckt, int restart)
if (firstTime) {
firstTime = 0;
#ifdef USE_CUSPICE
status = cuCKTstate01copy(ckt);
if (status != 0)
return (E_NOMEM);
#else
memcpy(ckt->CKTstate1, ckt->CKTstate0,
(size_t) ckt->CKTnumStates * sizeof(double));
#endif
}
i = 0;
@ -459,9 +479,21 @@ DCtrCurv(CKTcircuit *ckt, int restart)
((ISRCinstance*)(job->TRCVvElt[i]))->ISRCdcValue +=
job->TRCVvStep[i];
} else if (job->TRCVvType[i] == rcode) { /* resistance */
((RESinstance*)(job->TRCVvElt[i]))->RESresist +=
RESinstance *here = (RESinstance *)(job->TRCVvElt[i]);
here->RESresist +=
job->TRCVvStep[i];
RESupdate_conduct((RESinstance *)(job->TRCVvElt[i]), FALSE);
RESupdate_conduct(here, FALSE);
#ifdef USE_CUSPICE
RESmodel *model = RESmodPtr(here) ;
if (model->RESinitCUDA) {
model->RESparamCPU.RESconductArray[here->REScudaIndex] = here->RESconduct;
status = cuREStemp ((GENmodel *)model);
if (status != 0)
return E_NOMEM;
}
#endif
DEVices[rcode]->DEVload(job->TRCVvElt[i]->GENmodPtr, ckt);
} else if (job->TRCVvType[i] == TEMP_CODE) { /* temperature */
ckt->CKTtemp += job->TRCVvStep[i];