You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
 
 
 
 
 
 

135 lines
3.9 KiB

/**********
Copyright 1990 Regents of the University of California. All rights reserved.
Author: 1985 Thomas L. Quarles
Modifed: 2000 AlansFixes
**********/
#include "ngspice/config.h"
#include "ngspice/memory.h"
#include "ngspice/cktdefs.h"
#include "ngspice/devdefs.h"
#include "ngspice/sperror.h"
#include "ngspice/fteext.h"
#include "ngspice/ifsim.h"
#include "dev.h"
#ifdef XSPICE
/* gtri - add - wbk - 11/26/90 - add include for MIF global data */
#include "ngspice/mif.h"
/* gtri - end - wbk - 11/26/90 */
#endif
int
CKTinit(CKTcircuit **ckt) /* new circuit to create */
{
int i;
CKTcircuit *sckt = TMALLOC(CKTcircuit, 1);
*ckt = sckt;
if (sckt == NULL)
return(E_NOMEM);
/* gtri - begin - dynamically allocate the array of model lists */
/* CKThead used to be statically sized in CKTdefs.h, but has been changed */
/* to a ** pointer */
sckt->CKThead = TMALLOC(GENmodel *, DEVmaxnum);
if(sckt->CKThead == NULL) return(E_NOMEM);
/* gtri - end - dynamically allocate the array of model lists */
for (i = 0; i < DEVmaxnum; i++)
sckt->CKThead[i] = NULL;
sckt->CKTmaxEqNum = 1;
sckt->CKTnodes = NULL;
sckt->CKTlastNode = NULL;
sckt->CKTmatrix = NULL;
sckt->CKTgmin = 1e-12;
sckt->CKTgshunt=0;
sckt->CKTabstol = 1e-12;
sckt->CKTreltol = 1e-3;
sckt->CKTchgtol = 1e-14;
sckt->CKTvoltTol = 1e-6;
sckt->CKTtrtol = 7;
sckt->CKTbypass = 0;
sckt->CKTisSetup = 0;
#ifdef XSPICE
sckt->CKTadevFlag = 0; /* flag indicates A devices in the circuit */
#endif
sckt->CKTtranMaxIter = 10;
sckt->CKTdcMaxIter = 100;
sckt->CKTdcTrcvMaxIter = 50;
sckt->CKTintegrateMethod = TRAPEZOIDAL;
sckt->CKTorder = 1;
sckt->CKTmaxOrder = 2;
sckt->CKTpivotAbsTol = 1e-13;
sckt->CKTpivotRelTol = 1e-3;
sckt->CKTtemp = 300.15;
sckt->CKTnomTemp = 300.15;
sckt->CKTdefaultMosM = 1;
sckt->CKTdefaultMosL = 1e-4;
sckt->CKTdefaultMosW = 1e-4;
sckt->CKTdefaultMosAD = 0;
sckt->CKTdefaultMosAS = 0;
sckt->CKTsrcFact=1;
sckt->CKTdiagGmin=0;
/* PN: additions for circuit inventory */
sckt->CKTstat = TMALLOC(STATistics, 1);
if(sckt->CKTstat == NULL)
return(E_NOMEM);
sckt->CKTstat->STATdevNum = TMALLOC(STATdevList, DEVmaxnum);
if(sckt->CKTstat->STATdevNum == NULL)
return(E_NOMEM);
sckt->CKTtroubleNode = 0;
sckt->CKTtroubleElt = NULL;
sckt->CKTtimePoints = NULL;
if (sckt->CKTstat == NULL)
return E_NOMEM;
sckt->CKTnodeDamping = 0;
sckt->CKTabsDv = 0.5;
sckt->CKTrelDv = 2.0;
sckt->CKTvarHertz = 0;
sckt->DEVnameHash = nghash_init_pointer(100);
sckt->MODnameHash = nghash_init_pointer(100);
#ifdef XSPICE
/* gtri - begin - wbk - allocate/initialize substructs */
/* Allocate evt data structure */
sckt->evt = TMALLOC(Evt_Ckt_Data_t, 1);
if(! sckt->evt)
return(E_NOMEM);
/* Initialize options data */
sckt->evt->options.op_alternate = MIF_TRUE;
/* Allocate enh data structure */
sckt->enh = TMALLOC(Enh_Ckt_Data_t, 1);
if(! sckt->enh)
return(E_NOMEM);
/* Initialize non-zero, non-NULL data */
sckt->enh->breakpoint.current = 1.0e30;
sckt->enh->breakpoint.last = 1.0e30;
sckt->enh->ramp.ramptime = 0.0;
sckt->enh->conv_limit.enabled = MIF_TRUE;
sckt->enh->conv_limit.step = 0.25;
sckt->enh->conv_limit.abs_step = 0.1;
sckt->enh->rshunt_data.enabled = MIF_FALSE;
/* gtri - end - wbk - allocate/initialize substructs */
/* gtri - add - wbk - 01/12/91 - initialize g_mif_info */
g_mif_info.circuit.init = MIF_TRUE;
g_mif_info.circuit.anal_init = MIF_TRUE;
g_mif_info.circuit.anal_type = MIF_DC;
g_mif_info.instance = NULL;
g_mif_info.ckt = sckt;
g_mif_info.errmsg = NULL;
g_mif_info.auto_partial.global = MIF_FALSE;
g_mif_info.auto_partial.local = MIF_FALSE;
/* gtri - end - wbk - 01/12/91 */
#endif
return OK;
}