45 lines
1.1 KiB
C
45 lines
1.1 KiB
C
/**********
|
|
Copyright 1990 Regents of the University of California. All rights reserved.
|
|
Author: 1985 Thomas L. Quarles
|
|
**********/
|
|
|
|
/* NIreinit(ckt)
|
|
* Perform reinitialization necessary for the numerical iteration
|
|
* package - the matrix has now been fully accessed once, so we know
|
|
* how big it is, so allocate RHS vector
|
|
*/
|
|
|
|
#include "ngspice/ngspice.h"
|
|
#include "ngspice/cktdefs.h"
|
|
#include "ngspice/smpdefs.h"
|
|
#include "ngspice/sperror.h"
|
|
|
|
|
|
#define CKALLOC(ptr,size,type) if(( ckt->ptr =\
|
|
TMALLOC(type, size)) == NULL) return(E_NOMEM);
|
|
|
|
int
|
|
NIreinit( CKTcircuit *ckt)
|
|
{
|
|
int size;
|
|
#ifdef PREDICTOR
|
|
int i;
|
|
#endif
|
|
|
|
size = SMPmatSize(ckt->CKTmatrix);
|
|
CKALLOC(CKTrhs,size+1,double);
|
|
CKALLOC(CKTrhsOld,size+1,double);
|
|
CKALLOC(CKTrhsSpare,size+1,double);
|
|
CKALLOC(CKTirhs,size+1,double);
|
|
CKALLOC(CKTirhsOld,size+1,double);
|
|
CKALLOC(CKTirhsSpare,size+1,double);
|
|
#ifdef PREDICTOR
|
|
CKALLOC(CKTpred,size+1,double);
|
|
for( i=0;i<8;i++) {
|
|
CKALLOC(CKTsols[i],size+1,double);
|
|
}
|
|
#endif /* PREDICTOR */
|
|
ckt->CKTniState = NISHOULDREORDER | NIACSHOULDREORDER | NIPZSHOULDREORDER;
|
|
return(0);
|
|
}
|