Ejemplo n.º 1
0
int
LTRAsetup(SMPmatrix *matrix, GENmodel *inModel, CKTcircuit *ckt, int *state)
/*
 * load the transmission line structure with those pointers needed later for
 * fast matrix loading
 */
{
  LTRAmodel *model = (LTRAmodel *) inModel;
  LTRAinstance *here;
  int error;
  CKTnode *tmp;

  NG_IGNORE(state);

  /* loop through all the transmission line models */
  for (; model != NULL; model = model->LTRAnextModel) {

    if (!model->LTRAnlGiven) {
      model->LTRAnl = .25;
    }
    if (!model->LTRAfGiven) {
      model->LTRAf = 1e9;
    }
    if (!model->LTRAreltolGiven) {
      model->LTRAreltol = 1;
    }
    if (!model->LTRAabstolGiven) {
      model->LTRAabstol = 1;
    }
    if (!model->LTRAresistGiven) {
      SPfrontEnd->IFerror (ERR_WARNING,
	  "%s: lossy line series resistance not given, assumed zero",
	  &(model->LTRAmodName));
      model->LTRAresist = 0.0;
      /* return(E_BADPARM); */
    }
    if (model->LTRAstLineReltol == 0.0)
      model->LTRAstLineReltol = ckt->CKTreltol;
    if (model->LTRAstLineAbstol == 0.0)
      model->LTRAstLineAbstol = ckt->CKTabstol;
    /* LTRAchopReltol and LTRAchopAbstol default zero */

    if ((model->LTRAhowToInterp != LTRA_MOD_LININTERP) &&
	(model->LTRAhowToInterp != LTRA_MOD_QUADINTERP) &&
	(model->LTRAhowToInterp != LTRA_MOD_MIXEDINTERP)) {

      /*
       * SPfrontEnd->IFerror (ERR_FATAL, "%s: have to specify one of
       * lininterp, quadinterp or mixedinterp", &(model->LTRAmodName));
       * return(E_BADPARM);
       */
      if (ckt->CKTtryToCompact) {
	model->LTRAhowToInterp = LTRA_MOD_LININTERP;
	SPfrontEnd->IFerror (ERR_WARNING,
	    "%s: using linear interpolation because trytocompact option specified",
	    &(model->LTRAmodName));
      } else {
	model->LTRAhowToInterp = LTRA_MOD_QUADINTERP;
      }
    }
    if ((model->LTRAstepLimit != LTRA_MOD_NOSTEPLIMIT))
      model->LTRAstepLimit = LTRA_MOD_STEPLIMIT;
    if ((model->LTRAlteConType != LTRA_MOD_FULLCONTROL) &&
	(model->LTRAlteConType != LTRA_MOD_HALFCONTROL))
      model->LTRAlteConType = LTRA_MOD_NOCONTROL;

    if (!model->LTRAconductGiven) {
      /*
       * SPfrontEnd->IFerror (ERR_WARNING, "%s: lossy line parallel
       * conductance not given, assumed zero", &(model->LTRAmodName));
       */
      model->LTRAconduct = 0.0;
      /* return(E_BADPARM); */
    }
    if (!model->LTRAinductGiven) {
      SPfrontEnd->IFerror (ERR_WARNING,
	  "%s: lossy line series inductance not given, assumed zero",
	  &(model->LTRAmodName));
      model->LTRAinduct = 0.0;
      /* return(E_BADPARM); */
    }
    if (!model->LTRAcapacGiven) {
      SPfrontEnd->IFerror (ERR_FATAL,
	  "%s: lossy line parallel capacitance not given, assumed zero",
	  &(model->LTRAmodName));
      model->LTRAcapac = 0.0;
      /* return(E_BADPARM); */
    }
    if (!model->LTRAlengthGiven) {
      SPfrontEnd->IFerror (ERR_FATAL,
	  "%s: lossy line length must be given",
	  &(model->LTRAmodName));
      return (E_BADPARM);
    }
    if ((model->LTRAresist == 0) && (model->LTRAconduct == 0) &&
	(model->LTRAcapac != 0) && (model->LTRAinduct != 0)) {
      model->LTRAspecialCase = LTRA_MOD_LC;
#ifdef LTRADEBUG
      SPfrontEnd->IFerror (ERR_INFO,
	  "%s: lossless line",
	  &(model->LTRAmodName));
#endif
    }
    if ((model->LTRAresist != 0) && (model->LTRAconduct == 0) &&
	(model->LTRAcapac != 0) && (model->LTRAinduct != 0)) {
      model->LTRAspecialCase = LTRA_MOD_RLC;
#ifdef LTRADEBUG
      SPfrontEnd->IFerror (ERR_INFO,
	  "%s: RLC line",
	  &(model->LTRAmodName));
#endif
    }
    if ((model->LTRAresist != 0) && (model->LTRAconduct == 0) &&
	(model->LTRAcapac != 0) && (model->LTRAinduct == 0)) {
      model->LTRAspecialCase = LTRA_MOD_RC;
#ifdef LTRADEBUG
      SPfrontEnd->IFerror (ERR_INFO,
	  "%s: RC line",
	  &(model->LTRAmodName));
#endif
    }
    if ((model->LTRAresist != 0) && (model->LTRAconduct == 0) &&
	(model->LTRAcapac == 0) && (model->LTRAinduct != 0)) {
      model->LTRAspecialCase = LTRA_MOD_RL;
      SPfrontEnd->IFerror (ERR_FATAL,
	  "%s: RL line not supported yet",
	  &(model->LTRAmodName));
      return (E_BADPARM);
#ifdef LTRADEBUG
#endif
    }
    if ((model->LTRAresist != 0) && (model->LTRAconduct != 0) &&
	(model->LTRAcapac == 0) && (model->LTRAinduct == 0)) {
      model->LTRAspecialCase = LTRA_MOD_RG;
#ifdef LTRADEBUG
      SPfrontEnd->IFerror (ERR_INFO,
	  "%s: RG line",
	  &(model->LTRAmodName));
#endif
    }
    if ((model->LTRAconduct != 0) && ((model->LTRAcapac != 0) ||
	    (model->LTRAinduct != 0))) {
      model->LTRAspecialCase = LTRA_MOD_LTRA;
      SPfrontEnd->IFerror (ERR_FATAL,
	  "%s: Nonzero G (except RG) line not supported yet",
	  &(model->LTRAmodName));
      return (E_BADPARM);
#ifdef LTRADEBUG
#endif
    }
    if ((model->LTRAresist == 0.0 ? 0 : 1) + (model->LTRAconduct
	    == 0.0 ? 0 : 1) + (model->LTRAinduct == 0.0 ? 0 : 1) +
	(model->LTRAcapac == 0.0 ? 0 : 1) <= 1) {
      SPfrontEnd->IFerror (ERR_FATAL,
	  "%s: At least two of R,L,G,C must be specified and nonzero",
	  &(model->LTRAmodName));
      return (E_BADPARM);
    }
    /* loop through all the instances of the model */
    for (here = model->LTRAinstances; here != NULL;
	here = here->LTRAnextInstance) {

      if (here->LTRAbrEq1 == 0) {
	error = CKTmkVolt(ckt, &tmp, here->LTRAname, "i1");
	if (error)
	  return (error);
	here->LTRAbrEq1 = tmp->number;
      }
      if (here->LTRAbrEq2 == 0) {
	error = CKTmkVolt(ckt, &tmp, here->LTRAname, "i2");
	if (error)
	  return (error);
	here->LTRAbrEq2 = tmp->number;
      }
      /* macro to make elements with built in test for out of memory */
#define TSTALLOC(ptr,first,second) \
do { if((here->ptr = SMPmakeElt(matrix, here->first, here->second)) == NULL){\
    return(E_NOMEM);\
} } while(0)

      TSTALLOC(LTRAibr1Pos1Ptr, LTRAbrEq1, LTRAposNode1);
      TSTALLOC(LTRAibr1Neg1Ptr, LTRAbrEq1, LTRAnegNode1);
      TSTALLOC(LTRAibr1Pos2Ptr, LTRAbrEq1, LTRAposNode2);
      TSTALLOC(LTRAibr1Neg2Ptr, LTRAbrEq1, LTRAnegNode2);
      TSTALLOC(LTRAibr1Ibr1Ptr, LTRAbrEq1, LTRAbrEq1);
      TSTALLOC(LTRAibr1Ibr2Ptr, LTRAbrEq1, LTRAbrEq2);
      TSTALLOC(LTRAibr2Pos1Ptr, LTRAbrEq2, LTRAposNode1);
      TSTALLOC(LTRAibr2Neg1Ptr, LTRAbrEq2, LTRAnegNode1);
      TSTALLOC(LTRAibr2Pos2Ptr, LTRAbrEq2, LTRAposNode2);
      TSTALLOC(LTRAibr2Neg2Ptr, LTRAbrEq2, LTRAnegNode2);
      TSTALLOC(LTRAibr2Ibr1Ptr, LTRAbrEq2, LTRAbrEq1);
      TSTALLOC(LTRAibr2Ibr2Ptr, LTRAbrEq2, LTRAbrEq2);
      TSTALLOC(LTRApos1Ibr1Ptr, LTRAposNode1, LTRAbrEq1);
      TSTALLOC(LTRAneg1Ibr1Ptr, LTRAnegNode1, LTRAbrEq1);
      TSTALLOC(LTRApos2Ibr2Ptr, LTRAposNode2, LTRAbrEq2);
      TSTALLOC(LTRAneg2Ibr2Ptr, LTRAnegNode2, LTRAbrEq2);
      /*
       * the following are done so that SMPpreOrder does not screw up on
       * occasion - for example, when one end of the lossy line is hanging
       */
      TSTALLOC(LTRApos1Pos1Ptr, LTRAposNode1, LTRAposNode1);
      TSTALLOC(LTRAneg1Neg1Ptr, LTRAnegNode1, LTRAnegNode1);
      TSTALLOC(LTRApos2Pos2Ptr, LTRAposNode2, LTRAposNode2);
      TSTALLOC(LTRAneg2Neg2Ptr, LTRAnegNode2, LTRAnegNode2);
    }
  }
  return (OK);
}
Ejemplo n.º 2
0
int hicum0setup (SMPmatrix *matrix, GENmodel *inModel, CKTcircuit *ckt, int *states)
     /* load the hicum0 device structure with those pointers needed later 
      * for fast matrix loading 
      */
{
  register hicum0model *model = (hicum0model*)inModel;
  register hicum0instance *here;
  
  /*  loop through all the hicum0 device models */
  for ( ;model != NULL ;model = model->hicum0nextModel )
  {
    if(model->is_Given == FALSE) model->is=1.0e-16;
    if(model->mcf_Given == FALSE) model->mcf=1.00;
    if(model->mcr_Given == FALSE) model->mcr=1.00;
    if(model->vef_Given == FALSE) model->vef=1.0e6;
    if(model->iqf_Given == FALSE) model->iqf=1.0e6;
    if(model->iqr_Given == FALSE) model->iqr=1.0e6;
    if(model->iqfh_Given == FALSE) model->iqfh=1.0e6;
    if(model->tfh_Given == FALSE) model->tfh=0.0;
    if(model->ibes_Given == FALSE) model->ibes=1e-18;
    if(model->mbe_Given == FALSE) model->mbe=1.0;
    if(model->ires_Given == FALSE) model->ires=0.0;
    if(model->mre_Given == FALSE) model->mre=2.0;
    if(model->ibcs_Given == FALSE) model->ibcs=0.0;
    if(model->mbc_Given == FALSE) model->mbc=1.0;
    if(model->cje0_Given == FALSE) model->cje0=1.0e-20;
    if(model->vde_Given == FALSE) model->vde=0.9;
    if(model->ze_Given == FALSE) model->ze=0.5;
    if(model->aje_Given == FALSE) model->aje=2.5;
    if(model->t0_Given == FALSE) model->t0=0.0;
    if(model->dt0h_Given == FALSE) model->dt0h=0.0;
    if(model->tbvl_Given == FALSE) model->tbvl=0.0;
    if(model->tef0_Given == FALSE) model->tef0=0.0;
    if(model->gte_Given == FALSE) model->gte=1.0;
    if(model->thcs_Given == FALSE) model->thcs=0.0;
    if(model->ahc_Given == FALSE) model->ahc=0.1;
    if(model->tr_Given == FALSE) model->tr=0.0;
    if(model->rci0_Given == FALSE) model->rci0=150;
    if(model->vlim_Given == FALSE) model->vlim=0.5;
    if(model->vpt_Given == FALSE) model->vpt=100;
    if(model->vces_Given == FALSE) model->vces=0.1;
    if(model->cjci0_Given == FALSE) model->cjci0=1.0e-20;
    if(model->vdci_Given == FALSE) model->vdci=0.7;
    if(model->zci_Given == FALSE) model->zci=0.333;
    if(model->vptci_Given == FALSE) model->vptci=100;
    if(model->cjcx0_Given == FALSE) model->cjcx0=1.0e-20;
    if(model->vdcx_Given == FALSE) model->vdcx=0.7;
    if(model->zcx_Given == FALSE) model->zcx=0.333;
    if(model->vptcx_Given == FALSE) model->vptcx=100;
    if(model->fbc_Given == FALSE) model->fbc=1.0;
    if(model->rbi0_Given == FALSE) model->rbi0=0.0;
    if(model->vr0e_Given == FALSE) model->vr0e=2.5;
    if(model->vr0c_Given == FALSE) model->vr0c=1.0e6;
    if(model->fgeo_Given == FALSE) model->fgeo=0.656;
    if(model->rbx_Given == FALSE) model->rbx=0.0;
    if(model->rcx_Given == FALSE) model->rcx=0.0;
    if(model->re_Given == FALSE) model->re=0.0;
    if(model->itss_Given == FALSE) model->itss=0.0;
    if(model->msf_Given == FALSE) model->msf=1.0;
    if(model->iscs_Given == FALSE) model->iscs=0.0;
    if(model->msc_Given == FALSE) model->msc=1.0;
    if(model->cjs0_Given == FALSE) model->cjs0=1.0e-20;
    if(model->vds_Given == FALSE) model->vds=0.3;
    if(model->zs_Given == FALSE) model->zs=0.3;
    if(model->vpts_Given == FALSE) model->vpts=100;
    if(model->cbcpar_Given == FALSE) model->cbcpar=0.0;
    if(model->cbepar_Given == FALSE) model->cbepar=0.0;
    if(model->eavl_Given == FALSE) model->eavl=0.0;
    if(model->kavl_Given == FALSE) model->kavl=0.0;
    if(model->kf_Given == FALSE) model->kf=0.0;
    if(model->af_Given == FALSE) model->af=2.0;
    if(model->vgb_Given == FALSE) model->vgb=1.2;
    if(model->vge_Given == FALSE) model->vge=1.17;
    if(model->vgc_Given == FALSE) model->vgc=1.17;
    if(model->vgs_Given == FALSE) model->vgs=1.17;
    if(model->f1vg_Given == FALSE) model->f1vg=(-1.02377e-4);
    if(model->f2vg_Given == FALSE) model->f2vg=4.3215e-4;
    if(model->alt0_Given == FALSE) model->alt0=0.0;
    if(model->kt0_Given == FALSE) model->kt0=0.0;
    if(model->zetact_Given == FALSE) model->zetact=3.0;
    if(model->zetabet_Given == FALSE) model->zetabet=3.5;
    if(model->zetaci_Given == FALSE) model->zetaci=0.0;
    if(model->alvs_Given == FALSE) model->alvs=0.0;
    if(model->alces_Given == FALSE) model->alces=0.0;
    if(model->zetarbi_Given == FALSE) model->zetarbi=0.0;
    if(model->zetarbx_Given == FALSE) model->zetarbx=0.0;
    if(model->zetarcx_Given == FALSE) model->zetarcx=0.0;
    if(model->zetare_Given == FALSE) model->zetare=0.0;
    if(model->alkav_Given == FALSE) model->alkav=0.0;
    if(model->aleav_Given == FALSE) model->aleav=0.0;
    if(model->flsh_Given == FALSE) model->flsh=0;
    if(model->rth_Given == FALSE) model->rth=0.0;
    if(model->cth_Given == FALSE) model->cth=0.0;
    if(model->npn_Given == FALSE) model->npn=1;
    if(model->pnp_Given == FALSE) model->pnp=0;
    if(model->tnom_Given == FALSE) model->tnom=27;
    for ( here = model->hicum0instances ;here != NULL ; here = here->hicum0nextInstance )
    {
    if(here->dt_Given == FALSE) here->dt=0.0;
      /* Internal Nodes */
      {
        here->ciNode = -1;
        here->biNode = -1;
        here->eiNode = -1;
      }
      /* set states */
      here->state_tnode_GND = *states; *states += 2;
      here->state_b_e = *states; *states += 2;
      here->state_b_ci = *states; *states += 2;
      here->state_s_ci = *states; *states += 2;
      here->state_bi_ci = *states; *states += 2;
      here->state_bi_ei = *states; *states += 2;
      /* set Sparse Matrix Pointers */
      here->PTR_J_bi_bi_required=0;
      here->PTR_J_bi_ei_required=0;
      here->PTR_J_ei_bi_required=0;
      here->PTR_J_ei_ei_required=0;
      here->PTR_J_bi_ci_required=0;
      here->PTR_J_ci_bi_required=0;
      here->PTR_J_ci_ci_required=0;
      here->PTR_J_b_s_required=0;
      here->PTR_J_b_ci_required=0;
      here->PTR_J_s_s_required=0;
      here->PTR_J_s_ci_required=0;
      here->PTR_J_b_b_required=0;
      here->PTR_J_s_b_required=0;
      here->PTR_J_b_tnode_required=0;
      here->PTR_J_s_tnode_required=0;
      here->PTR_J_s_bi_required=0;
      here->PTR_J_s_ei_required=0;
      here->PTR_J_ci_ei_required=0;
      here->PTR_J_ci_s_required=0;
      here->PTR_J_ci_tnode_required=0;
      here->PTR_J_ci_b_required=0;
      here->PTR_J_b_bi_required=0;
      here->PTR_J_b_e_required=0;
      here->PTR_J_e_b_required=0;
      here->PTR_J_e_e_required=0;
      here->PTR_J_ei_tnode_required=0;
      here->PTR_J_e_tnode_required=0;
      here->PTR_J_ei_e_required=0;
      here->PTR_J_e_ei_required=0;
      here->PTR_J_c_tnode_required=0;
      here->PTR_J_ci_c_required=0;
      here->PTR_J_c_ci_required=0;
      here->PTR_J_c_c_required=0;
      here->PTR_J_b_ei_required=0;
      here->PTR_J_bi_b_required=0;
      here->PTR_J_bi_tnode_required=0;
      here->PTR_J_ei_ci_required=0;
      here->PTR_J_ei_b_required=0;
      here->PTR_J_tnode_tnode_required=0;
      here->PTR_J_tnode_ci_required=0;
      here->PTR_J_tnode_ei_required=0;
      here->PTR_J_tnode_b_required=0;
      here->PTR_J_tnode_bi_required=0;
      here->PTR_J_tnode_s_required=0;
      here->PTR_J_tnode_e_required=0;
      here->PTR_J_tnode_c_required=0;

      hicum0guesstopology(matrix,ckt,model,here);

      /* Internal Nodes */
      {
        int error;
        CKTnode *tmp;
        if(here->ciNode == -1)
        {
          error=CKTmkVolt(ckt,&tmp,here->hicum0name,"internal collector node");
          if(error) return(error);
          here->ciNode = tmp->number;
        }
        if(here->biNode == -1)
        {
          error=CKTmkVolt(ckt,&tmp,here->hicum0name,"internal base node");
          if(error) return(error);
          here->biNode = tmp->number;
        }
        if(here->eiNode == -1)
        {
          error=CKTmkVolt(ckt,&tmp,here->hicum0name,"internal emitter node");
          if(error) return(error);
          here->eiNode = tmp->number;
        }
      }
        if(here->PTR_J_bi_bi_required==1)
        {
          here->PTR_J_bi_bi=SMPmakeElt(matrix,here->biNode,here->biNode);
        }
        if(here->PTR_J_bi_ei_required==1)
        {
          here->PTR_J_bi_ei=SMPmakeElt(matrix,here->biNode,here->eiNode);
        }
        if(here->PTR_J_ei_bi_required==1)
        {
          here->PTR_J_ei_bi=SMPmakeElt(matrix,here->eiNode,here->biNode);
        }
        if(here->PTR_J_ei_ei_required==1)
        {
          here->PTR_J_ei_ei=SMPmakeElt(matrix,here->eiNode,here->eiNode);
        }
        if(here->PTR_J_bi_ci_required==1)
        {
          here->PTR_J_bi_ci=SMPmakeElt(matrix,here->biNode,here->ciNode);
        }
        if(here->PTR_J_ci_bi_required==1)
        {
          here->PTR_J_ci_bi=SMPmakeElt(matrix,here->ciNode,here->biNode);
        }
        if(here->PTR_J_ci_ci_required==1)
        {
          here->PTR_J_ci_ci=SMPmakeElt(matrix,here->ciNode,here->ciNode);
        }
        if(here->PTR_J_b_s_required==1)
        {
          here->PTR_J_b_s=SMPmakeElt(matrix,here->bNode,here->sNode);
        }
        if(here->PTR_J_b_ci_required==1)
        {
          here->PTR_J_b_ci=SMPmakeElt(matrix,here->bNode,here->ciNode);
        }
        if(here->PTR_J_s_s_required==1)
        {
          here->PTR_J_s_s=SMPmakeElt(matrix,here->sNode,here->sNode);
        }
        if(here->PTR_J_s_ci_required==1)
        {
          here->PTR_J_s_ci=SMPmakeElt(matrix,here->sNode,here->ciNode);
        }
        if(here->PTR_J_b_b_required==1)
        {
          here->PTR_J_b_b=SMPmakeElt(matrix,here->bNode,here->bNode);
        }
        if(here->PTR_J_s_b_required==1)
        {
          here->PTR_J_s_b=SMPmakeElt(matrix,here->sNode,here->bNode);
        }
        if(here->PTR_J_b_tnode_required==1)
        {
          here->PTR_J_b_tnode=SMPmakeElt(matrix,here->bNode,here->tnodeNode);
        }
        if(here->PTR_J_s_tnode_required==1)
        {
          here->PTR_J_s_tnode=SMPmakeElt(matrix,here->sNode,here->tnodeNode);
        }
        if(here->PTR_J_s_bi_required==1)
        {
          here->PTR_J_s_bi=SMPmakeElt(matrix,here->sNode,here->biNode);
        }
        if(here->PTR_J_s_ei_required==1)
        {
          here->PTR_J_s_ei=SMPmakeElt(matrix,here->sNode,here->eiNode);
        }
        if(here->PTR_J_ci_ei_required==1)
        {
          here->PTR_J_ci_ei=SMPmakeElt(matrix,here->ciNode,here->eiNode);
        }
        if(here->PTR_J_ci_s_required==1)
        {
          here->PTR_J_ci_s=SMPmakeElt(matrix,here->ciNode,here->sNode);
        }
        if(here->PTR_J_ci_tnode_required==1)
        {
          here->PTR_J_ci_tnode=SMPmakeElt(matrix,here->ciNode,here->tnodeNode);
        }
        if(here->PTR_J_ci_b_required==1)
        {
          here->PTR_J_ci_b=SMPmakeElt(matrix,here->ciNode,here->bNode);
        }
        if(here->PTR_J_b_bi_required==1)
        {
          here->PTR_J_b_bi=SMPmakeElt(matrix,here->bNode,here->biNode);
        }
        if(here->PTR_J_b_e_required==1)
        {
          here->PTR_J_b_e=SMPmakeElt(matrix,here->bNode,here->eNode);
        }
        if(here->PTR_J_e_b_required==1)
        {
          here->PTR_J_e_b=SMPmakeElt(matrix,here->eNode,here->bNode);
        }
        if(here->PTR_J_e_e_required==1)
        {
          here->PTR_J_e_e=SMPmakeElt(matrix,here->eNode,here->eNode);
        }
        if(here->PTR_J_ei_tnode_required==1)
        {
          here->PTR_J_ei_tnode=SMPmakeElt(matrix,here->eiNode,here->tnodeNode);
        }
        if(here->PTR_J_e_tnode_required==1)
        {
          here->PTR_J_e_tnode=SMPmakeElt(matrix,here->eNode,here->tnodeNode);
        }
        if(here->PTR_J_ei_e_required==1)
        {
          here->PTR_J_ei_e=SMPmakeElt(matrix,here->eiNode,here->eNode);
        }
        if(here->PTR_J_e_ei_required==1)
        {
          here->PTR_J_e_ei=SMPmakeElt(matrix,here->eNode,here->eiNode);
        }
        if(here->PTR_J_c_tnode_required==1)
        {
          here->PTR_J_c_tnode=SMPmakeElt(matrix,here->cNode,here->tnodeNode);
        }
        if(here->PTR_J_ci_c_required==1)
        {
          here->PTR_J_ci_c=SMPmakeElt(matrix,here->ciNode,here->cNode);
        }
        if(here->PTR_J_c_ci_required==1)
        {
          here->PTR_J_c_ci=SMPmakeElt(matrix,here->cNode,here->ciNode);
        }
        if(here->PTR_J_c_c_required==1)
        {
          here->PTR_J_c_c=SMPmakeElt(matrix,here->cNode,here->cNode);
        }
        if(here->PTR_J_b_ei_required==1)
        {
          here->PTR_J_b_ei=SMPmakeElt(matrix,here->bNode,here->eiNode);
        }
        if(here->PTR_J_bi_b_required==1)
        {
          here->PTR_J_bi_b=SMPmakeElt(matrix,here->biNode,here->bNode);
        }
        if(here->PTR_J_bi_tnode_required==1)
        {
          here->PTR_J_bi_tnode=SMPmakeElt(matrix,here->biNode,here->tnodeNode);
        }
        if(here->PTR_J_ei_ci_required==1)
        {
          here->PTR_J_ei_ci=SMPmakeElt(matrix,here->eiNode,here->ciNode);
        }
        if(here->PTR_J_ei_b_required==1)
        {
          here->PTR_J_ei_b=SMPmakeElt(matrix,here->eiNode,here->bNode);
        }
        if(here->PTR_J_tnode_tnode_required==1)
        {
          here->PTR_J_tnode_tnode=SMPmakeElt(matrix,here->tnodeNode,here->tnodeNode);
        }
        if(here->PTR_J_tnode_ci_required==1)
        {
          here->PTR_J_tnode_ci=SMPmakeElt(matrix,here->tnodeNode,here->ciNode);
        }
        if(here->PTR_J_tnode_ei_required==1)
        {
          here->PTR_J_tnode_ei=SMPmakeElt(matrix,here->tnodeNode,here->eiNode);
        }
        if(here->PTR_J_tnode_b_required==1)
        {
          here->PTR_J_tnode_b=SMPmakeElt(matrix,here->tnodeNode,here->bNode);
        }
        if(here->PTR_J_tnode_bi_required==1)
        {
          here->PTR_J_tnode_bi=SMPmakeElt(matrix,here->tnodeNode,here->biNode);
        }
        if(here->PTR_J_tnode_s_required==1)
        {
          here->PTR_J_tnode_s=SMPmakeElt(matrix,here->tnodeNode,here->sNode);
        }
        if(here->PTR_J_tnode_e_required==1)
        {
          here->PTR_J_tnode_e=SMPmakeElt(matrix,here->tnodeNode,here->eNode);
        }
        if(here->PTR_J_tnode_c_required==1)
        {
          here->PTR_J_tnode_c=SMPmakeElt(matrix,here->tnodeNode,here->cNode);
        }

    }
  }
  return(OK);
}
Ejemplo n.º 3
0
int
BSIM3v32setup (SMPmatrix *matrix, GENmodel *inModel, CKTcircuit *ckt,
            int *states)
{
BSIM3v32model *model = (BSIM3v32model*)inModel;
BSIM3v32instance *here;
int error;
CKTnode *tmp;

CKTnode *tmpNode;
IFuid tmpName;

#ifdef USE_OMP
int idx, InstCount;
BSIM3v32instance **InstArray;
#endif

    /*  loop through all the BSIM3v32 device models */
    for( ; model != NULL; model = model->BSIM3v32nextModel )
    {
/* Default value Processing for BSIM3v32 MOSFET Models */
        if (!model->BSIM3v32typeGiven)
            model->BSIM3v32type = NMOS;
        if (!model->BSIM3v32mobModGiven)
            model->BSIM3v32mobMod = 1;
        if (!model->BSIM3v32binUnitGiven)
            model->BSIM3v32binUnit = 1;
        if (!model->BSIM3v32paramChkGiven)
            model->BSIM3v32paramChk = 0;
        if (!model->BSIM3v32capModGiven)
            model->BSIM3v32capMod = 3;
        if (!model->BSIM3v32acmModGiven)
            model->BSIM3v32acmMod = 0;
        if (!model->BSIM3v32calcacmGiven)
            model->BSIM3v32calcacm = 0;
        if (!model->BSIM3v32noiModGiven)
            model->BSIM3v32noiMod = 1;
        if (!model->BSIM3v32nqsModGiven)
            model->BSIM3v32nqsMod = 0;
        else if ((model->BSIM3v32nqsMod != 0) && (model->BSIM3v32nqsMod != 1))
        {   model->BSIM3v32nqsMod = 0;
            printf("Warning: nqsMod has been set to its default value: 0.\n");
        }

        /* If the user does not provide the model revision,
         * we always choose the most recent.
         */
        if (!model->BSIM3v32versionGiven)
                model->BSIM3v32version = copy("3.2.4");

        /* I have added below the code that translate model string
         * into an integer. This trick is meant to speed up the
         * revision testing instruction, since comparing integer
         * is faster than comparing strings.
         * Paolo Nenzi 2002
         */
        if ((!strcmp(model->BSIM3v32version, "3.2.4"))||(!strcmp(model->BSIM3v32version, "3.24")))
                model->BSIM3v32intVersion = BSIM3v32V324;
        else if ((!strcmp(model->BSIM3v32version, "3.2.3"))||(!strcmp(model->BSIM3v32version, "3.23")))
                model->BSIM3v32intVersion = BSIM3v32V323;
        else if ((!strcmp(model->BSIM3v32version, "3.2.2"))||(!strcmp(model->BSIM3v32version, "3.22")))
                model->BSIM3v32intVersion = BSIM3v32V322;
        else if ((!strcmp(model->BSIM3v32version, "3.2"))||(!strcmp(model->BSIM3v32version, "3.20")))
                model->BSIM3v32intVersion = BSIM3v32V32;
        else
                model->BSIM3v32intVersion = BSIM3v32V3OLD;
        /* BSIM3v32V3OLD is a placeholder for pre 3.2 revision
         * This model should not be used for pre 3.2 models.
         */

        if (!model->BSIM3v32toxGiven)
            model->BSIM3v32tox = 150.0e-10;
        model->BSIM3v32cox = 3.453133e-11 / model->BSIM3v32tox;
        if (!model->BSIM3v32toxmGiven)
            model->BSIM3v32toxm = model->BSIM3v32tox;

        if (!model->BSIM3v32cdscGiven)
            model->BSIM3v32cdsc = 2.4e-4;   /* unit Q/V/m^2  */
        if (!model->BSIM3v32cdscbGiven)
            model->BSIM3v32cdscb = 0.0;   /* unit Q/V/m^2  */
        if (!model->BSIM3v32cdscdGiven)
            model->BSIM3v32cdscd = 0.0;   /* unit Q/V/m^2  */
        if (!model->BSIM3v32citGiven)
            model->BSIM3v32cit = 0.0;   /* unit Q/V/m^2  */
        if (!model->BSIM3v32nfactorGiven)
            model->BSIM3v32nfactor = 1;
        if (!model->BSIM3v32xjGiven)
            model->BSIM3v32xj = .15e-6;
        if (!model->BSIM3v32vsatGiven)
            model->BSIM3v32vsat = 8.0e4;    /* unit m/s */
        if (!model->BSIM3v32atGiven)
            model->BSIM3v32at = 3.3e4;    /* unit m/s */
        if (!model->BSIM3v32a0Given)
            model->BSIM3v32a0 = 1.0;
        if (!model->BSIM3v32agsGiven)
            model->BSIM3v32ags = 0.0;
        if (!model->BSIM3v32a1Given)
            model->BSIM3v32a1 = 0.0;
        if (!model->BSIM3v32a2Given)
            model->BSIM3v32a2 = 1.0;
        if (!model->BSIM3v32ketaGiven)
            model->BSIM3v32keta = -0.047;    /* unit  / V */
        if (!model->BSIM3v32nsubGiven)
            model->BSIM3v32nsub = 6.0e16;   /* unit 1/cm3 */
        if (!model->BSIM3v32npeakGiven)
            model->BSIM3v32npeak = 1.7e17;   /* unit 1/cm3 */
        if (!model->BSIM3v32ngateGiven)
            model->BSIM3v32ngate = 0;   /* unit 1/cm3 */
        if (!model->BSIM3v32vbmGiven)
            model->BSIM3v32vbm = -3.0;
        if (!model->BSIM3v32xtGiven)
            model->BSIM3v32xt = 1.55e-7;
        if (!model->BSIM3v32kt1Given)
            model->BSIM3v32kt1 = -0.11;      /* unit V */
        if (!model->BSIM3v32kt1lGiven)
            model->BSIM3v32kt1l = 0.0;      /* unit V*m */
        if (!model->BSIM3v32kt2Given)
            model->BSIM3v32kt2 = 0.022;      /* No unit */
        if (!model->BSIM3v32k3Given)
            model->BSIM3v32k3 = 80.0;
        if (!model->BSIM3v32k3bGiven)
            model->BSIM3v32k3b = 0.0;
        if (!model->BSIM3v32w0Given)
            model->BSIM3v32w0 = 2.5e-6;
        if (!model->BSIM3v32nlxGiven)
            model->BSIM3v32nlx = 1.74e-7;
        if (!model->BSIM3v32dvt0Given)
            model->BSIM3v32dvt0 = 2.2;
        if (!model->BSIM3v32dvt1Given)
            model->BSIM3v32dvt1 = 0.53;
        if (!model->BSIM3v32dvt2Given)
            model->BSIM3v32dvt2 = -0.032;   /* unit 1 / V */

        if (!model->BSIM3v32dvt0wGiven)
            model->BSIM3v32dvt0w = 0.0;
        if (!model->BSIM3v32dvt1wGiven)
            model->BSIM3v32dvt1w = 5.3e6;
        if (!model->BSIM3v32dvt2wGiven)
            model->BSIM3v32dvt2w = -0.032;

        if (!model->BSIM3v32droutGiven)
            model->BSIM3v32drout = 0.56;
        if (!model->BSIM3v32dsubGiven)
            model->BSIM3v32dsub = model->BSIM3v32drout;
        if (!model->BSIM3v32vth0Given)
            model->BSIM3v32vth0 = (model->BSIM3v32type == NMOS) ? 0.7 : -0.7;
        if (!model->BSIM3v32uaGiven)
            model->BSIM3v32ua = 2.25e-9;      /* unit m/V */
        if (!model->BSIM3v32ua1Given)
            model->BSIM3v32ua1 = 4.31e-9;      /* unit m/V */
        if (!model->BSIM3v32ubGiven)
            model->BSIM3v32ub = 5.87e-19;     /* unit (m/V)**2 */
        if (!model->BSIM3v32ub1Given)
            model->BSIM3v32ub1 = -7.61e-18;     /* unit (m/V)**2 */
        if (!model->BSIM3v32ucGiven)
            model->BSIM3v32uc = (model->BSIM3v32mobMod == 3) ? -0.0465 : -0.0465e-9;
        if (!model->BSIM3v32uc1Given)
            model->BSIM3v32uc1 = (model->BSIM3v32mobMod == 3) ? -0.056 : -0.056e-9;
        if (!model->BSIM3v32u0Given)
            model->BSIM3v32u0 = (model->BSIM3v32type == NMOS) ? 0.067 : 0.025;
        if (!model->BSIM3v32uteGiven)
            model->BSIM3v32ute = -1.5;
        if (!model->BSIM3v32voffGiven)
            model->BSIM3v32voff = -0.08;
        if (!model->BSIM3v32deltaGiven)
           model->BSIM3v32delta = 0.01;
        if (!model->BSIM3v32rdswGiven)
            model->BSIM3v32rdsw = 0;
        if (!model->BSIM3v32prwgGiven)
            model->BSIM3v32prwg = 0.0;      /* unit 1/V */
        if (!model->BSIM3v32prwbGiven)
            model->BSIM3v32prwb = 0.0;
        if (!model->BSIM3v32prtGiven)
            model->BSIM3v32prt = 0.0;
        if (!model->BSIM3v32eta0Given)
            model->BSIM3v32eta0 = 0.08;      /* no unit  */
        if (!model->BSIM3v32etabGiven)
            model->BSIM3v32etab = -0.07;      /* unit  1/V */
        if (!model->BSIM3v32pclmGiven)
            model->BSIM3v32pclm = 1.3;      /* no unit  */
        if (!model->BSIM3v32pdibl1Given)
            model->BSIM3v32pdibl1 = .39;    /* no unit  */
        if (!model->BSIM3v32pdibl2Given)
            model->BSIM3v32pdibl2 = 0.0086;    /* no unit  */
        if (!model->BSIM3v32pdiblbGiven)
            model->BSIM3v32pdiblb = 0.0;    /* 1/V  */
        if (!model->BSIM3v32pscbe1Given)
            model->BSIM3v32pscbe1 = 4.24e8;
        if (!model->BSIM3v32pscbe2Given)
            model->BSIM3v32pscbe2 = 1.0e-5;
        if (!model->BSIM3v32pvagGiven)
            model->BSIM3v32pvag = 0.0;
        if (!model->BSIM3v32wrGiven)
            model->BSIM3v32wr = 1.0;
        if (!model->BSIM3v32dwgGiven)
            model->BSIM3v32dwg = 0.0;
        if (!model->BSIM3v32dwbGiven)
            model->BSIM3v32dwb = 0.0;
        if (!model->BSIM3v32b0Given)
            model->BSIM3v32b0 = 0.0;
        if (!model->BSIM3v32b1Given)
            model->BSIM3v32b1 = 0.0;
        if (!model->BSIM3v32alpha0Given)
            model->BSIM3v32alpha0 = 0.0;
        if (!model->BSIM3v32alpha1Given)
            model->BSIM3v32alpha1 = 0.0;
        if (!model->BSIM3v32beta0Given)
            model->BSIM3v32beta0 = 30.0;
        if (!model->BSIM3v32ijthGiven)
            model->BSIM3v32ijth = 0.1; /* unit A */

        if (!model->BSIM3v32elmGiven)
            model->BSIM3v32elm = 5.0;
        if (!model->BSIM3v32cgslGiven)
            model->BSIM3v32cgsl = 0.0;
        if (!model->BSIM3v32cgdlGiven)
            model->BSIM3v32cgdl = 0.0;
        if (!model->BSIM3v32ckappaGiven)
            model->BSIM3v32ckappa = 0.6;
        if (!model->BSIM3v32clcGiven)
            model->BSIM3v32clc = 0.1e-6;
        if (!model->BSIM3v32cleGiven)
            model->BSIM3v32cle = 0.6;
        if (!model->BSIM3v32vfbcvGiven)
            model->BSIM3v32vfbcv = -1.0;
        if (!model->BSIM3v32acdeGiven)
            model->BSIM3v32acde = 1.0;
        if (!model->BSIM3v32moinGiven)
            model->BSIM3v32moin = 15.0;
        if (!model->BSIM3v32noffGiven)
            model->BSIM3v32noff = 1.0;
        if (!model->BSIM3v32voffcvGiven)
            model->BSIM3v32voffcv = 0.0;
        if (!model->BSIM3v32tcjGiven)
            model->BSIM3v32tcj = 0.0;
        if (!model->BSIM3v32tpbGiven)
            model->BSIM3v32tpb = 0.0;
        if (!model->BSIM3v32tcjswGiven)
            model->BSIM3v32tcjsw = 0.0;
        if (!model->BSIM3v32tpbswGiven)
            model->BSIM3v32tpbsw = 0.0;
        if (!model->BSIM3v32tcjswgGiven)
            model->BSIM3v32tcjswg = 0.0;
        if (!model->BSIM3v32tpbswgGiven)
            model->BSIM3v32tpbswg = 0.0;

        /* ACM model */
        if (!model->BSIM3v32hdifGiven)
          model->BSIM3v32hdif = 0.0;
        if (!model->BSIM3v32ldifGiven)
          model->BSIM3v32ldif = 0.0;
        if (!model->BSIM3v32ldGiven)
          model->BSIM3v32ld = 0.0;
        if (!model->BSIM3v32rdGiven)
          model->BSIM3v32rd = 0.0;
        if (!model->BSIM3v32rsGiven)
          model->BSIM3v32rs = 0.0;
        if (!model->BSIM3v32rdcGiven)
          model->BSIM3v32rdc = 0.0;
        if (!model->BSIM3v32rscGiven)
          model->BSIM3v32rsc = 0.0;
        if (!model->BSIM3v32wmltGiven)
          model->BSIM3v32wmlt = 1.0;

        if (!model->BSIM3v32lmltGiven)
          model->BSIM3v32lmlt = 1.0;

        /* Length dependence */
        if (!model->BSIM3v32lcdscGiven)
            model->BSIM3v32lcdsc = 0.0;
        if (!model->BSIM3v32lcdscbGiven)
            model->BSIM3v32lcdscb = 0.0;
        if (!model->BSIM3v32lcdscdGiven)
            model->BSIM3v32lcdscd = 0.0;
        if (!model->BSIM3v32lcitGiven)
            model->BSIM3v32lcit = 0.0;
        if (!model->BSIM3v32lnfactorGiven)
            model->BSIM3v32lnfactor = 0.0;
        if (!model->BSIM3v32lxjGiven)
            model->BSIM3v32lxj = 0.0;
        if (!model->BSIM3v32lvsatGiven)
            model->BSIM3v32lvsat = 0.0;
        if (!model->BSIM3v32latGiven)
            model->BSIM3v32lat = 0.0;
        if (!model->BSIM3v32la0Given)
            model->BSIM3v32la0 = 0.0;
        if (!model->BSIM3v32lagsGiven)
            model->BSIM3v32lags = 0.0;
        if (!model->BSIM3v32la1Given)
            model->BSIM3v32la1 = 0.0;
        if (!model->BSIM3v32la2Given)
            model->BSIM3v32la2 = 0.0;
        if (!model->BSIM3v32lketaGiven)
            model->BSIM3v32lketa = 0.0;
        if (!model->BSIM3v32lnsubGiven)
            model->BSIM3v32lnsub = 0.0;
        if (!model->BSIM3v32lnpeakGiven)
            model->BSIM3v32lnpeak = 0.0;
        if (!model->BSIM3v32lngateGiven)
            model->BSIM3v32lngate = 0.0;
        if (!model->BSIM3v32lvbmGiven)
            model->BSIM3v32lvbm = 0.0;
        if (!model->BSIM3v32lxtGiven)
            model->BSIM3v32lxt = 0.0;
        if (!model->BSIM3v32lkt1Given)
            model->BSIM3v32lkt1 = 0.0;
        if (!model->BSIM3v32lkt1lGiven)
            model->BSIM3v32lkt1l = 0.0;
        if (!model->BSIM3v32lkt2Given)
            model->BSIM3v32lkt2 = 0.0;
        if (!model->BSIM3v32lk3Given)
            model->BSIM3v32lk3 = 0.0;
        if (!model->BSIM3v32lk3bGiven)
            model->BSIM3v32lk3b = 0.0;
        if (!model->BSIM3v32lw0Given)
            model->BSIM3v32lw0 = 0.0;
        if (!model->BSIM3v32lnlxGiven)
            model->BSIM3v32lnlx = 0.0;
        if (!model->BSIM3v32ldvt0Given)
            model->BSIM3v32ldvt0 = 0.0;
        if (!model->BSIM3v32ldvt1Given)
            model->BSIM3v32ldvt1 = 0.0;
        if (!model->BSIM3v32ldvt2Given)
            model->BSIM3v32ldvt2 = 0.0;
        if (!model->BSIM3v32ldvt0wGiven)
            model->BSIM3v32ldvt0w = 0.0;
        if (!model->BSIM3v32ldvt1wGiven)
            model->BSIM3v32ldvt1w = 0.0;
        if (!model->BSIM3v32ldvt2wGiven)
            model->BSIM3v32ldvt2w = 0.0;
        if (!model->BSIM3v32ldroutGiven)
            model->BSIM3v32ldrout = 0.0;
        if (!model->BSIM3v32ldsubGiven)
            model->BSIM3v32ldsub = 0.0;
        if (!model->BSIM3v32lvth0Given)
           model->BSIM3v32lvth0 = 0.0;
        if (!model->BSIM3v32luaGiven)
            model->BSIM3v32lua = 0.0;
        if (!model->BSIM3v32lua1Given)
            model->BSIM3v32lua1 = 0.0;
        if (!model->BSIM3v32lubGiven)
            model->BSIM3v32lub = 0.0;
        if (!model->BSIM3v32lub1Given)
            model->BSIM3v32lub1 = 0.0;
        if (!model->BSIM3v32lucGiven)
            model->BSIM3v32luc = 0.0;
        if (!model->BSIM3v32luc1Given)
            model->BSIM3v32luc1 = 0.0;
        if (!model->BSIM3v32lu0Given)
            model->BSIM3v32lu0 = 0.0;
        if (!model->BSIM3v32luteGiven)
            model->BSIM3v32lute = 0.0;
        if (!model->BSIM3v32lvoffGiven)
            model->BSIM3v32lvoff = 0.0;
        if (!model->BSIM3v32ldeltaGiven)
            model->BSIM3v32ldelta = 0.0;
        if (!model->BSIM3v32lrdswGiven)
            model->BSIM3v32lrdsw = 0.0;
        if (!model->BSIM3v32lprwbGiven)
            model->BSIM3v32lprwb = 0.0;
        if (!model->BSIM3v32lprwgGiven)
            model->BSIM3v32lprwg = 0.0;
        if (!model->BSIM3v32lprtGiven)
            model->BSIM3v32lprt = 0.0;
        if (!model->BSIM3v32leta0Given)
            model->BSIM3v32leta0 = 0.0;
        if (!model->BSIM3v32letabGiven)
            model->BSIM3v32letab = -0.0;
        if (!model->BSIM3v32lpclmGiven)
            model->BSIM3v32lpclm = 0.0;
        if (!model->BSIM3v32lpdibl1Given)
            model->BSIM3v32lpdibl1 = 0.0;
        if (!model->BSIM3v32lpdibl2Given)
            model->BSIM3v32lpdibl2 = 0.0;
        if (!model->BSIM3v32lpdiblbGiven)
            model->BSIM3v32lpdiblb = 0.0;
        if (!model->BSIM3v32lpscbe1Given)
            model->BSIM3v32lpscbe1 = 0.0;
        if (!model->BSIM3v32lpscbe2Given)
            model->BSIM3v32lpscbe2 = 0.0;
        if (!model->BSIM3v32lpvagGiven)
            model->BSIM3v32lpvag = 0.0;
        if (!model->BSIM3v32lwrGiven)
            model->BSIM3v32lwr = 0.0;
        if (!model->BSIM3v32ldwgGiven)
            model->BSIM3v32ldwg = 0.0;
        if (!model->BSIM3v32ldwbGiven)
            model->BSIM3v32ldwb = 0.0;
        if (!model->BSIM3v32lb0Given)
            model->BSIM3v32lb0 = 0.0;
        if (!model->BSIM3v32lb1Given)
            model->BSIM3v32lb1 = 0.0;
        if (!model->BSIM3v32lalpha0Given)
            model->BSIM3v32lalpha0 = 0.0;
        if (!model->BSIM3v32lalpha1Given)
            model->BSIM3v32lalpha1 = 0.0;
        if (!model->BSIM3v32lbeta0Given)
            model->BSIM3v32lbeta0 = 0.0;
        if (!model->BSIM3v32lvfbGiven)
            model->BSIM3v32lvfb = 0.0;

        if (!model->BSIM3v32lelmGiven)
            model->BSIM3v32lelm = 0.0;
        if (!model->BSIM3v32lcgslGiven)
            model->BSIM3v32lcgsl = 0.0;
        if (!model->BSIM3v32lcgdlGiven)
            model->BSIM3v32lcgdl = 0.0;
        if (!model->BSIM3v32lckappaGiven)
            model->BSIM3v32lckappa = 0.0;
        if (!model->BSIM3v32lclcGiven)
            model->BSIM3v32lclc = 0.0;
        if (!model->BSIM3v32lcleGiven)
            model->BSIM3v32lcle = 0.0;
        if (!model->BSIM3v32lcfGiven)
            model->BSIM3v32lcf = 0.0;
        if (!model->BSIM3v32lvfbcvGiven)
            model->BSIM3v32lvfbcv = 0.0;
        if (!model->BSIM3v32lacdeGiven)
            model->BSIM3v32lacde = 0.0;
        if (!model->BSIM3v32lmoinGiven)
            model->BSIM3v32lmoin = 0.0;
        if (!model->BSIM3v32lnoffGiven)
            model->BSIM3v32lnoff = 0.0;
        if (!model->BSIM3v32lvoffcvGiven)
            model->BSIM3v32lvoffcv = 0.0;

        /* Width dependence */
        if (!model->BSIM3v32wcdscGiven)
            model->BSIM3v32wcdsc = 0.0;
        if (!model->BSIM3v32wcdscbGiven)
            model->BSIM3v32wcdscb = 0.0;
            if (!model->BSIM3v32wcdscdGiven)
            model->BSIM3v32wcdscd = 0.0;
        if (!model->BSIM3v32wcitGiven)
            model->BSIM3v32wcit = 0.0;
        if (!model->BSIM3v32wnfactorGiven)
            model->BSIM3v32wnfactor = 0.0;
        if (!model->BSIM3v32wxjGiven)
            model->BSIM3v32wxj = 0.0;
        if (!model->BSIM3v32wvsatGiven)
            model->BSIM3v32wvsat = 0.0;
        if (!model->BSIM3v32watGiven)
            model->BSIM3v32wat = 0.0;
        if (!model->BSIM3v32wa0Given)
            model->BSIM3v32wa0 = 0.0;
        if (!model->BSIM3v32wagsGiven)
            model->BSIM3v32wags = 0.0;
        if (!model->BSIM3v32wa1Given)
            model->BSIM3v32wa1 = 0.0;
        if (!model->BSIM3v32wa2Given)
            model->BSIM3v32wa2 = 0.0;
        if (!model->BSIM3v32wketaGiven)
            model->BSIM3v32wketa = 0.0;
        if (!model->BSIM3v32wnsubGiven)
            model->BSIM3v32wnsub = 0.0;
        if (!model->BSIM3v32wnpeakGiven)
            model->BSIM3v32wnpeak = 0.0;
        if (!model->BSIM3v32wngateGiven)
            model->BSIM3v32wngate = 0.0;
        if (!model->BSIM3v32wvbmGiven)
            model->BSIM3v32wvbm = 0.0;
        if (!model->BSIM3v32wxtGiven)
            model->BSIM3v32wxt = 0.0;
        if (!model->BSIM3v32wkt1Given)
            model->BSIM3v32wkt1 = 0.0;
        if (!model->BSIM3v32wkt1lGiven)
            model->BSIM3v32wkt1l = 0.0;
        if (!model->BSIM3v32wkt2Given)
            model->BSIM3v32wkt2 = 0.0;
        if (!model->BSIM3v32wk3Given)
            model->BSIM3v32wk3 = 0.0;
        if (!model->BSIM3v32wk3bGiven)
            model->BSIM3v32wk3b = 0.0;
        if (!model->BSIM3v32ww0Given)
            model->BSIM3v32ww0 = 0.0;
        if (!model->BSIM3v32wnlxGiven)
            model->BSIM3v32wnlx = 0.0;
        if (!model->BSIM3v32wdvt0Given)
            model->BSIM3v32wdvt0 = 0.0;
        if (!model->BSIM3v32wdvt1Given)
            model->BSIM3v32wdvt1 = 0.0;
        if (!model->BSIM3v32wdvt2Given)
            model->BSIM3v32wdvt2 = 0.0;
        if (!model->BSIM3v32wdvt0wGiven)
            model->BSIM3v32wdvt0w = 0.0;
        if (!model->BSIM3v32wdvt1wGiven)
            model->BSIM3v32wdvt1w = 0.0;
        if (!model->BSIM3v32wdvt2wGiven)
            model->BSIM3v32wdvt2w = 0.0;
        if (!model->BSIM3v32wdroutGiven)
            model->BSIM3v32wdrout = 0.0;
        if (!model->BSIM3v32wdsubGiven)
            model->BSIM3v32wdsub = 0.0;
        if (!model->BSIM3v32wvth0Given)
           model->BSIM3v32wvth0 = 0.0;
        if (!model->BSIM3v32wuaGiven)
            model->BSIM3v32wua = 0.0;
        if (!model->BSIM3v32wua1Given)
            model->BSIM3v32wua1 = 0.0;
        if (!model->BSIM3v32wubGiven)
            model->BSIM3v32wub = 0.0;
        if (!model->BSIM3v32wub1Given)
            model->BSIM3v32wub1 = 0.0;
        if (!model->BSIM3v32wucGiven)
            model->BSIM3v32wuc = 0.0;
        if (!model->BSIM3v32wuc1Given)
            model->BSIM3v32wuc1 = 0.0;
        if (!model->BSIM3v32wu0Given)
            model->BSIM3v32wu0 = 0.0;
        if (!model->BSIM3v32wuteGiven)
            model->BSIM3v32wute = 0.0;
        if (!model->BSIM3v32wvoffGiven)
            model->BSIM3v32wvoff = 0.0;
        if (!model->BSIM3v32wdeltaGiven)
            model->BSIM3v32wdelta = 0.0;
        if (!model->BSIM3v32wrdswGiven)
            model->BSIM3v32wrdsw = 0.0;
        if (!model->BSIM3v32wprwbGiven)
            model->BSIM3v32wprwb = 0.0;
        if (!model->BSIM3v32wprwgGiven)
            model->BSIM3v32wprwg = 0.0;
        if (!model->BSIM3v32wprtGiven)
            model->BSIM3v32wprt = 0.0;
        if (!model->BSIM3v32weta0Given)
            model->BSIM3v32weta0 = 0.0;
        if (!model->BSIM3v32wetabGiven)
            model->BSIM3v32wetab = 0.0;
        if (!model->BSIM3v32wpclmGiven)
            model->BSIM3v32wpclm = 0.0;
        if (!model->BSIM3v32wpdibl1Given)
            model->BSIM3v32wpdibl1 = 0.0;
        if (!model->BSIM3v32wpdibl2Given)
            model->BSIM3v32wpdibl2 = 0.0;
        if (!model->BSIM3v32wpdiblbGiven)
            model->BSIM3v32wpdiblb = 0.0;
        if (!model->BSIM3v32wpscbe1Given)
            model->BSIM3v32wpscbe1 = 0.0;
        if (!model->BSIM3v32wpscbe2Given)
            model->BSIM3v32wpscbe2 = 0.0;
        if (!model->BSIM3v32wpvagGiven)
            model->BSIM3v32wpvag = 0.0;
        if (!model->BSIM3v32wwrGiven)
            model->BSIM3v32wwr = 0.0;
        if (!model->BSIM3v32wdwgGiven)
            model->BSIM3v32wdwg = 0.0;
        if (!model->BSIM3v32wdwbGiven)
            model->BSIM3v32wdwb = 0.0;
        if (!model->BSIM3v32wb0Given)
            model->BSIM3v32wb0 = 0.0;
        if (!model->BSIM3v32wb1Given)
            model->BSIM3v32wb1 = 0.0;
        if (!model->BSIM3v32walpha0Given)
            model->BSIM3v32walpha0 = 0.0;
        if (!model->BSIM3v32walpha1Given)
            model->BSIM3v32walpha1 = 0.0;
        if (!model->BSIM3v32wbeta0Given)
            model->BSIM3v32wbeta0 = 0.0;
        if (!model->BSIM3v32wvfbGiven)
            model->BSIM3v32wvfb = 0.0;

        if (!model->BSIM3v32welmGiven)
            model->BSIM3v32welm = 0.0;
        if (!model->BSIM3v32wcgslGiven)
            model->BSIM3v32wcgsl = 0.0;
        if (!model->BSIM3v32wcgdlGiven)
            model->BSIM3v32wcgdl = 0.0;
        if (!model->BSIM3v32wckappaGiven)
            model->BSIM3v32wckappa = 0.0;
        if (!model->BSIM3v32wcfGiven)
            model->BSIM3v32wcf = 0.0;
        if (!model->BSIM3v32wclcGiven)
            model->BSIM3v32wclc = 0.0;
        if (!model->BSIM3v32wcleGiven)
            model->BSIM3v32wcle = 0.0;
        if (!model->BSIM3v32wvfbcvGiven)
            model->BSIM3v32wvfbcv = 0.0;
        if (!model->BSIM3v32wacdeGiven)
            model->BSIM3v32wacde = 0.0;
        if (!model->BSIM3v32wmoinGiven)
            model->BSIM3v32wmoin = 0.0;
        if (!model->BSIM3v32wnoffGiven)
            model->BSIM3v32wnoff = 0.0;
        if (!model->BSIM3v32wvoffcvGiven)
            model->BSIM3v32wvoffcv = 0.0;

        /* Cross-term dependence */
        if (!model->BSIM3v32pcdscGiven)
            model->BSIM3v32pcdsc = 0.0;
        if (!model->BSIM3v32pcdscbGiven)
            model->BSIM3v32pcdscb = 0.0;
        if (!model->BSIM3v32pcdscdGiven)
            model->BSIM3v32pcdscd = 0.0;
        if (!model->BSIM3v32pcitGiven)
            model->BSIM3v32pcit = 0.0;
        if (!model->BSIM3v32pnfactorGiven)
            model->BSIM3v32pnfactor = 0.0;
        if (!model->BSIM3v32pxjGiven)
            model->BSIM3v32pxj = 0.0;
        if (!model->BSIM3v32pvsatGiven)
            model->BSIM3v32pvsat = 0.0;
        if (!model->BSIM3v32patGiven)
            model->BSIM3v32pat = 0.0;
        if (!model->BSIM3v32pa0Given)
            model->BSIM3v32pa0 = 0.0;

        if (!model->BSIM3v32pagsGiven)
            model->BSIM3v32pags = 0.0;
        if (!model->BSIM3v32pa1Given)
            model->BSIM3v32pa1 = 0.0;
        if (!model->BSIM3v32pa2Given)
            model->BSIM3v32pa2 = 0.0;
        if (!model->BSIM3v32pketaGiven)
            model->BSIM3v32pketa = 0.0;
        if (!model->BSIM3v32pnsubGiven)
            model->BSIM3v32pnsub = 0.0;
        if (!model->BSIM3v32pnpeakGiven)
            model->BSIM3v32pnpeak = 0.0;
        if (!model->BSIM3v32pngateGiven)
            model->BSIM3v32pngate = 0.0;
        if (!model->BSIM3v32pvbmGiven)
            model->BSIM3v32pvbm = 0.0;
        if (!model->BSIM3v32pxtGiven)
            model->BSIM3v32pxt = 0.0;
        if (!model->BSIM3v32pkt1Given)
            model->BSIM3v32pkt1 = 0.0;
        if (!model->BSIM3v32pkt1lGiven)
            model->BSIM3v32pkt1l = 0.0;
        if (!model->BSIM3v32pkt2Given)
            model->BSIM3v32pkt2 = 0.0;
        if (!model->BSIM3v32pk3Given)
            model->BSIM3v32pk3 = 0.0;
        if (!model->BSIM3v32pk3bGiven)
            model->BSIM3v32pk3b = 0.0;
        if (!model->BSIM3v32pw0Given)
            model->BSIM3v32pw0 = 0.0;
        if (!model->BSIM3v32pnlxGiven)
            model->BSIM3v32pnlx = 0.0;
        if (!model->BSIM3v32pdvt0Given)
            model->BSIM3v32pdvt0 = 0.0;
        if (!model->BSIM3v32pdvt1Given)
            model->BSIM3v32pdvt1 = 0.0;
        if (!model->BSIM3v32pdvt2Given)
            model->BSIM3v32pdvt2 = 0.0;
        if (!model->BSIM3v32pdvt0wGiven)
            model->BSIM3v32pdvt0w = 0.0;
        if (!model->BSIM3v32pdvt1wGiven)
            model->BSIM3v32pdvt1w = 0.0;
        if (!model->BSIM3v32pdvt2wGiven)
            model->BSIM3v32pdvt2w = 0.0;
        if (!model->BSIM3v32pdroutGiven)
            model->BSIM3v32pdrout = 0.0;
        if (!model->BSIM3v32pdsubGiven)
            model->BSIM3v32pdsub = 0.0;
        if (!model->BSIM3v32pvth0Given)
           model->BSIM3v32pvth0 = 0.0;
        if (!model->BSIM3v32puaGiven)
            model->BSIM3v32pua = 0.0;
        if (!model->BSIM3v32pua1Given)
            model->BSIM3v32pua1 = 0.0;
        if (!model->BSIM3v32pubGiven)
            model->BSIM3v32pub = 0.0;
        if (!model->BSIM3v32pub1Given)
            model->BSIM3v32pub1 = 0.0;
        if (!model->BSIM3v32pucGiven)
            model->BSIM3v32puc = 0.0;
        if (!model->BSIM3v32puc1Given)
            model->BSIM3v32puc1 = 0.0;
        if (!model->BSIM3v32pu0Given)
            model->BSIM3v32pu0 = 0.0;
        if (!model->BSIM3v32puteGiven)
            model->BSIM3v32pute = 0.0;
        if (!model->BSIM3v32pvoffGiven)
            model->BSIM3v32pvoff = 0.0;
        if (!model->BSIM3v32pdeltaGiven)
            model->BSIM3v32pdelta = 0.0;
        if (!model->BSIM3v32prdswGiven)
            model->BSIM3v32prdsw = 0.0;
        if (!model->BSIM3v32pprwbGiven)
            model->BSIM3v32pprwb = 0.0;
        if (!model->BSIM3v32pprwgGiven)
            model->BSIM3v32pprwg = 0.0;
        if (!model->BSIM3v32pprtGiven)
            model->BSIM3v32pprt = 0.0;
        if (!model->BSIM3v32peta0Given)
            model->BSIM3v32peta0 = 0.0;
        if (!model->BSIM3v32petabGiven)
            model->BSIM3v32petab = 0.0;
        if (!model->BSIM3v32ppclmGiven)
            model->BSIM3v32ppclm = 0.0;
        if (!model->BSIM3v32ppdibl1Given)
            model->BSIM3v32ppdibl1 = 0.0;
        if (!model->BSIM3v32ppdibl2Given)
            model->BSIM3v32ppdibl2 = 0.0;
        if (!model->BSIM3v32ppdiblbGiven)
            model->BSIM3v32ppdiblb = 0.0;
        if (!model->BSIM3v32ppscbe1Given)
            model->BSIM3v32ppscbe1 = 0.0;
        if (!model->BSIM3v32ppscbe2Given)
            model->BSIM3v32ppscbe2 = 0.0;
        if (!model->BSIM3v32ppvagGiven)
            model->BSIM3v32ppvag = 0.0;
        if (!model->BSIM3v32pwrGiven)
            model->BSIM3v32pwr = 0.0;
        if (!model->BSIM3v32pdwgGiven)
            model->BSIM3v32pdwg = 0.0;
        if (!model->BSIM3v32pdwbGiven)
            model->BSIM3v32pdwb = 0.0;
        if (!model->BSIM3v32pb0Given)
            model->BSIM3v32pb0 = 0.0;
        if (!model->BSIM3v32pb1Given)
            model->BSIM3v32pb1 = 0.0;
        if (!model->BSIM3v32palpha0Given)
            model->BSIM3v32palpha0 = 0.0;
        if (!model->BSIM3v32palpha1Given)
            model->BSIM3v32palpha1 = 0.0;
        if (!model->BSIM3v32pbeta0Given)
            model->BSIM3v32pbeta0 = 0.0;
        if (!model->BSIM3v32pvfbGiven)
            model->BSIM3v32pvfb = 0.0;

        if (!model->BSIM3v32pelmGiven)
            model->BSIM3v32pelm = 0.0;
        if (!model->BSIM3v32pcgslGiven)
            model->BSIM3v32pcgsl = 0.0;
        if (!model->BSIM3v32pcgdlGiven)
            model->BSIM3v32pcgdl = 0.0;
        if (!model->BSIM3v32pckappaGiven)
            model->BSIM3v32pckappa = 0.0;
        if (!model->BSIM3v32pcfGiven)
            model->BSIM3v32pcf = 0.0;
        if (!model->BSIM3v32pclcGiven)
            model->BSIM3v32pclc = 0.0;
        if (!model->BSIM3v32pcleGiven)
            model->BSIM3v32pcle = 0.0;
        if (!model->BSIM3v32pvfbcvGiven)
            model->BSIM3v32pvfbcv = 0.0;
        if (!model->BSIM3v32pacdeGiven)
            model->BSIM3v32pacde = 0.0;
        if (!model->BSIM3v32pmoinGiven)
            model->BSIM3v32pmoin = 0.0;
        if (!model->BSIM3v32pnoffGiven)
            model->BSIM3v32pnoff = 0.0;
        if (!model->BSIM3v32pvoffcvGiven)
            model->BSIM3v32pvoffcv = 0.0;

        /* unit degree celcius */
        if (!model->BSIM3v32tnomGiven)
            model->BSIM3v32tnom = ckt->CKTnomTemp;
/*        else
            model->BSIM3v32tnom = model->BSIM3v32tnom + 273.15; we make this transform in b3v32mpar.c in the first run */
        if (!model->BSIM3v32LintGiven)
           model->BSIM3v32Lint = 0.0;
        if (!model->BSIM3v32LlGiven)
           model->BSIM3v32Ll = 0.0;
        if (!model->BSIM3v32LlcGiven)
           model->BSIM3v32Llc = model->BSIM3v32Ll;
        if (!model->BSIM3v32LlnGiven)
           model->BSIM3v32Lln = 1.0;
        if (!model->BSIM3v32LwGiven)
           model->BSIM3v32Lw = 0.0;
        if (!model->BSIM3v32LwcGiven)
           model->BSIM3v32Lwc = model->BSIM3v32Lw;
        if (!model->BSIM3v32LwnGiven)
           model->BSIM3v32Lwn = 1.0;
        if (!model->BSIM3v32LwlGiven)
           model->BSIM3v32Lwl = 0.0;
        if (!model->BSIM3v32LwlcGiven)
           model->BSIM3v32Lwlc = model->BSIM3v32Lwl;
        if (!model->BSIM3v32LminGiven)
           model->BSIM3v32Lmin = 0.0;
        if (!model->BSIM3v32LmaxGiven)
           model->BSIM3v32Lmax = 1.0;
        if (!model->BSIM3v32WintGiven)
           model->BSIM3v32Wint = 0.0;
        if (!model->BSIM3v32WlGiven)
           model->BSIM3v32Wl = 0.0;
        if (!model->BSIM3v32WlcGiven)
           model->BSIM3v32Wlc = model->BSIM3v32Wl;
        if (!model->BSIM3v32WlnGiven)
           model->BSIM3v32Wln = 1.0;
        if (!model->BSIM3v32WwGiven)
           model->BSIM3v32Ww = 0.0;
        if (!model->BSIM3v32WwcGiven)
           model->BSIM3v32Wwc = model->BSIM3v32Ww;
        if (!model->BSIM3v32WwnGiven)
           model->BSIM3v32Wwn = 1.0;
        if (!model->BSIM3v32WwlGiven)
           model->BSIM3v32Wwl = 0.0;
        if (!model->BSIM3v32WwlcGiven)
           model->BSIM3v32Wwlc = model->BSIM3v32Wwl;
        if (!model->BSIM3v32WminGiven)
           model->BSIM3v32Wmin = 0.0;
        if (!model->BSIM3v32WmaxGiven)
           model->BSIM3v32Wmax = 1.0;
        if (!model->BSIM3v32dwcGiven)
           model->BSIM3v32dwc = model->BSIM3v32Wint;
        if (!model->BSIM3v32dlcGiven)
           model->BSIM3v32dlc = model->BSIM3v32Lint;

        if (!model->BSIM3v32xlGiven)
           model->BSIM3v32xl = 0.0;
        if (!model->BSIM3v32xwGiven)
           model->BSIM3v32xw = 0.0;

        if (!model->BSIM3v32cfGiven)
            model->BSIM3v32cf = 2.0 * EPSOX / PI
                           * log(1.0 + 0.4e-6 / model->BSIM3v32tox);
        if (!model->BSIM3v32cgdoGiven)
        {   if (model->BSIM3v32dlcGiven && (model->BSIM3v32dlc > 0.0))
            {   model->BSIM3v32cgdo = model->BSIM3v32dlc * model->BSIM3v32cox
                                 - model->BSIM3v32cgdl ;
            }
            else
                model->BSIM3v32cgdo = 0.6 * model->BSIM3v32xj * model->BSIM3v32cox;
        }
        if (!model->BSIM3v32cgsoGiven)
        {   if (model->BSIM3v32dlcGiven && (model->BSIM3v32dlc > 0.0))
            {   model->BSIM3v32cgso = model->BSIM3v32dlc * model->BSIM3v32cox
                                 - model->BSIM3v32cgsl ;
            }
            else
                model->BSIM3v32cgso = 0.6 * model->BSIM3v32xj * model->BSIM3v32cox;
        }

        if (!model->BSIM3v32cgboGiven)
        {   model->BSIM3v32cgbo = 2.0 * model->BSIM3v32dwc * model->BSIM3v32cox;
        }
        if (!model->BSIM3v32xpartGiven)
            model->BSIM3v32xpart = 0.0;
        if (!model->BSIM3v32sheetResistanceGiven)
            model->BSIM3v32sheetResistance = 0.0;
        if (!model->BSIM3v32unitAreaJctCapGiven)
            model->BSIM3v32unitAreaJctCap = 5.0E-4;
        if (!model->BSIM3v32unitLengthSidewallJctCapGiven)
            model->BSIM3v32unitLengthSidewallJctCap = 5.0E-10;
        if (!model->BSIM3v32unitLengthGateSidewallJctCapGiven)
            model->BSIM3v32unitLengthGateSidewallJctCap = model->BSIM3v32unitLengthSidewallJctCap ;
        if (!model->BSIM3v32jctSatCurDensityGiven)
            model->BSIM3v32jctSatCurDensity = 1.0E-4;
        if (!model->BSIM3v32jctSidewallSatCurDensityGiven)
            model->BSIM3v32jctSidewallSatCurDensity = 0.0;
        if (!model->BSIM3v32bulkJctPotentialGiven)
            model->BSIM3v32bulkJctPotential = 1.0;
        if (!model->BSIM3v32sidewallJctPotentialGiven)
            model->BSIM3v32sidewallJctPotential = 1.0;
        if (!model->BSIM3v32GatesidewallJctPotentialGiven)
            model->BSIM3v32GatesidewallJctPotential = model->BSIM3v32sidewallJctPotential;
        if (!model->BSIM3v32bulkJctBotGradingCoeffGiven)
            model->BSIM3v32bulkJctBotGradingCoeff = 0.5;
        if (!model->BSIM3v32bulkJctSideGradingCoeffGiven)
            model->BSIM3v32bulkJctSideGradingCoeff = 0.33;
        if (!model->BSIM3v32bulkJctGateSideGradingCoeffGiven)
            model->BSIM3v32bulkJctGateSideGradingCoeff = model->BSIM3v32bulkJctSideGradingCoeff;
        if (!model->BSIM3v32jctEmissionCoeffGiven)
            model->BSIM3v32jctEmissionCoeff = 1.0;
        if (!model->BSIM3v32jctTempExponentGiven)
            model->BSIM3v32jctTempExponent = 3.0;
        if (!model->BSIM3v32oxideTrapDensityAGiven)
        {   if (model->BSIM3v32type == NMOS)
                model->BSIM3v32oxideTrapDensityA = 1e20;
            else
                model->BSIM3v32oxideTrapDensityA=9.9e18;
        }
        if (!model->BSIM3v32oxideTrapDensityBGiven)
        {   if (model->BSIM3v32type == NMOS)
                model->BSIM3v32oxideTrapDensityB = 5e4;
            else
                model->BSIM3v32oxideTrapDensityB = 2.4e3;
        }
        if (!model->BSIM3v32oxideTrapDensityCGiven)
        {   if (model->BSIM3v32type == NMOS)
                model->BSIM3v32oxideTrapDensityC = -1.4e-12;
            else
                model->BSIM3v32oxideTrapDensityC = 1.4e-12;

        }
        if (!model->BSIM3v32emGiven)
            model->BSIM3v32em = 4.1e7; /* V/m */
        if (!model->BSIM3v32efGiven)
            model->BSIM3v32ef = 1.0;
        if (!model->BSIM3v32afGiven)
            model->BSIM3v32af = 1.0;
        if (!model->BSIM3v32kfGiven)
            model->BSIM3v32kf = 0.0;

        if (!model->BSIM3v32vgsMaxGiven)
            model->BSIM3v32vgsMax = 1e99;
        if (!model->BSIM3v32vgdMaxGiven)
            model->BSIM3v32vgdMax = 1e99;
        if (!model->BSIM3v32vgbMaxGiven)
            model->BSIM3v32vgbMax = 1e99;
        if (!model->BSIM3v32vdsMaxGiven)
            model->BSIM3v32vdsMax = 1e99;
        if (!model->BSIM3v32vbsMaxGiven)
            model->BSIM3v32vbsMax = 1e99;
        if (!model->BSIM3v32vbdMaxGiven)
            model->BSIM3v32vbdMax = 1e99;

        /* loop through all the instances of the model */
        for (here = model->BSIM3v32instances; here != NULL ;
             here=here->BSIM3v32nextInstance)
        {
            /* allocate a chunk of the state vector */
            here->BSIM3v32states = *states;
            *states += BSIM3v32numStates;

            /* perform the parameter defaulting */
            if (!here->BSIM3v32drainAreaGiven)
                here->BSIM3v32drainArea = 0.0;
            if (!here->BSIM3v32drainPerimeterGiven)
                here->BSIM3v32drainPerimeter = 0.0;
            if (!here->BSIM3v32drainSquaresGiven)
            {
                if (model->BSIM3v32acmMod == 0)
                  here->BSIM3v32drainSquares = 1.0;
                else
                  here->BSIM3v32drainSquares = 0.0;
            }
            if (!here->BSIM3v32delvtoGiven)
                here->BSIM3v32delvto = 0.0;
            if (!here->BSIM3v32mulu0Given)
                here->BSIM3v32mulu0 = 1.0;
            if (!here->BSIM3v32icVBSGiven)
                here->BSIM3v32icVBS = 0.0;
            if (!here->BSIM3v32icVDSGiven)
                here->BSIM3v32icVDS = 0.0;
            if (!here->BSIM3v32icVGSGiven)
                here->BSIM3v32icVGS = 0.0;
            if (!here->BSIM3v32lGiven)
                here->BSIM3v32l = 5.0e-6;
            if (!here->BSIM3v32sourceAreaGiven)
                here->BSIM3v32sourceArea = 0.0;
            if (!here->BSIM3v32sourcePerimeterGiven)
                here->BSIM3v32sourcePerimeter = 0.0;
            if (!here->BSIM3v32sourceSquaresGiven)
            {
                if (model->BSIM3v32acmMod == 0)
                  here->BSIM3v32sourceSquares = 1.0;
                else
                  here->BSIM3v32sourceSquares = 0.0;
            }
            if (!here->BSIM3v32wGiven)
                here->BSIM3v32w = 5.0e-6;
            if (!here->BSIM3v32nqsModGiven)
                here->BSIM3v32nqsMod = model->BSIM3v32nqsMod;
            else if ((here->BSIM3v32nqsMod != 0) && (here->BSIM3v32nqsMod != 1))
            {   here->BSIM3v32nqsMod = model->BSIM3v32nqsMod;
                printf("Warning: nqsMod has been set to its global value %d.\n",
                model->BSIM3v32nqsMod);
            }
            if (!here->BSIM3v32geoGiven)
                here->BSIM3v32geo = 0;

            if (!here->BSIM3v32mGiven)
                here->BSIM3v32m = 1;

            /* process drain series resistance */
            if (  ((model->BSIM3v32sheetResistance > 0.0) && (here->BSIM3v32drainSquares > 0.0))
                ||((model->BSIM3v32sheetResistance > 0.0) && (model->BSIM3v32hdif > 0.0))
                ||((model->BSIM3v32rd > 0.0) && (model->BSIM3v32ldif > 0.0))
                ||((model->BSIM3v32rd > 0.0) && (model->BSIM3v32ld > 0.0))
                ||((model->BSIM3v32rdc > 0.0))
               )
            {
               if(here->BSIM3v32dNodePrime == 0) {
                 error = CKTmkVolt(ckt,&tmp,here->BSIM3v32name,"drain");
                 if(error) return(error);
                 here->BSIM3v32dNodePrime = tmp->number;
                 if (ckt->CKTcopyNodesets) {
                   if (CKTinst2Node(ckt,here,1,&tmpNode,&tmpName)==OK) {
                      if (tmpNode->nsGiven) {
                        tmp->nodeset=tmpNode->nodeset;
                        tmp->nsGiven=tmpNode->nsGiven;
                      }
                   }
                 }
               }
            }
            else
            {  here->BSIM3v32dNodePrime = here->BSIM3v32dNode;
            }

            /* process source series resistance */
            if (  ((model->BSIM3v32sheetResistance > 0.0) && (here->BSIM3v32sourceSquares > 0.0))
                ||((model->BSIM3v32sheetResistance > 0.0) && (model->BSIM3v32hdif > 0.0))
                ||((model->BSIM3v32rs > 0.0) && (model->BSIM3v32ldif > 0.0))
                ||((model->BSIM3v32rs > 0.0) && (model->BSIM3v32ld > 0.0))
                ||((model->BSIM3v32rsc > 0.0))
               )
            {
               if(here->BSIM3v32sNodePrime == 0) {
                 error = CKTmkVolt(ckt,&tmp,here->BSIM3v32name,"source");
                 if(error) return(error);
                 here->BSIM3v32sNodePrime = tmp->number;
                 if (ckt->CKTcopyNodesets) {
                   if (CKTinst2Node(ckt,here,3,&tmpNode,&tmpName)==OK) {
                      if (tmpNode->nsGiven) {
                        tmp->nodeset=tmpNode->nodeset;
                        tmp->nsGiven=tmpNode->nsGiven;
                      }
                   }
                 }
               }
            }
            else
            {  here->BSIM3v32sNodePrime = here->BSIM3v32sNode;
            }

            /* internal charge node */

            if (here->BSIM3v32nqsMod)
            {   if(here->BSIM3v32qNode == 0)
                {   error = CKTmkVolt(ckt,&tmp,here->BSIM3v32name,"charge");
                    if(error) return(error);
                    here->BSIM3v32qNode = tmp->number;
                }
            }
            else
            {   here->BSIM3v32qNode = 0;
            }

            /* Channel length scaling with lmlt model parameter */

            if (model->BSIM3v32lmltGiven)
                here->BSIM3v32l *= model->BSIM3v32lmlt;

        /* set Sparse Matrix Pointers */

/* macro to make elements with built in test for out of memory */
#define TSTALLOC(ptr,first,second) \
do { if((here->ptr = SMPmakeElt(matrix, here->first, here->second)) == NULL){\
    return(E_NOMEM);\
} } while(0)

            TSTALLOC(BSIM3v32DdPtr, BSIM3v32dNode, BSIM3v32dNode);
            TSTALLOC(BSIM3v32GgPtr, BSIM3v32gNode, BSIM3v32gNode);
            TSTALLOC(BSIM3v32SsPtr, BSIM3v32sNode, BSIM3v32sNode);
            TSTALLOC(BSIM3v32BbPtr, BSIM3v32bNode, BSIM3v32bNode);
            TSTALLOC(BSIM3v32DPdpPtr, BSIM3v32dNodePrime, BSIM3v32dNodePrime);
            TSTALLOC(BSIM3v32SPspPtr, BSIM3v32sNodePrime, BSIM3v32sNodePrime);
            TSTALLOC(BSIM3v32DdpPtr, BSIM3v32dNode, BSIM3v32dNodePrime);
            TSTALLOC(BSIM3v32GbPtr, BSIM3v32gNode, BSIM3v32bNode);
            TSTALLOC(BSIM3v32GdpPtr, BSIM3v32gNode, BSIM3v32dNodePrime);
            TSTALLOC(BSIM3v32GspPtr, BSIM3v32gNode, BSIM3v32sNodePrime);
            TSTALLOC(BSIM3v32SspPtr, BSIM3v32sNode, BSIM3v32sNodePrime);
            TSTALLOC(BSIM3v32BdpPtr, BSIM3v32bNode, BSIM3v32dNodePrime);
            TSTALLOC(BSIM3v32BspPtr, BSIM3v32bNode, BSIM3v32sNodePrime);
            TSTALLOC(BSIM3v32DPspPtr, BSIM3v32dNodePrime, BSIM3v32sNodePrime);
            TSTALLOC(BSIM3v32DPdPtr, BSIM3v32dNodePrime, BSIM3v32dNode);
            TSTALLOC(BSIM3v32BgPtr, BSIM3v32bNode, BSIM3v32gNode);
            TSTALLOC(BSIM3v32DPgPtr, BSIM3v32dNodePrime, BSIM3v32gNode);
            TSTALLOC(BSIM3v32SPgPtr, BSIM3v32sNodePrime, BSIM3v32gNode);
            TSTALLOC(BSIM3v32SPsPtr, BSIM3v32sNodePrime, BSIM3v32sNode);
            TSTALLOC(BSIM3v32DPbPtr, BSIM3v32dNodePrime, BSIM3v32bNode);
            TSTALLOC(BSIM3v32SPbPtr, BSIM3v32sNodePrime, BSIM3v32bNode);
            TSTALLOC(BSIM3v32SPdpPtr, BSIM3v32sNodePrime, BSIM3v32dNodePrime);

            TSTALLOC(BSIM3v32QqPtr, BSIM3v32qNode, BSIM3v32qNode);

            TSTALLOC(BSIM3v32QdpPtr, BSIM3v32qNode, BSIM3v32dNodePrime);
            TSTALLOC(BSIM3v32QspPtr, BSIM3v32qNode, BSIM3v32sNodePrime);
            TSTALLOC(BSIM3v32QgPtr, BSIM3v32qNode, BSIM3v32gNode);
            TSTALLOC(BSIM3v32QbPtr, BSIM3v32qNode, BSIM3v32bNode);
            TSTALLOC(BSIM3v32DPqPtr, BSIM3v32dNodePrime, BSIM3v32qNode);
            TSTALLOC(BSIM3v32SPqPtr, BSIM3v32sNodePrime, BSIM3v32qNode);
            TSTALLOC(BSIM3v32GqPtr, BSIM3v32gNode, BSIM3v32qNode);
            TSTALLOC(BSIM3v32BqPtr, BSIM3v32bNode, BSIM3v32qNode);

        }
    }
#ifdef USE_OMP
    InstCount = 0;
    model = (BSIM3v32model*)inModel;
    /* loop through all the BSIM3 device models
    to count the number of instances */

    for (; model != NULL; model = model->BSIM3v32nextModel)
    {
        /* loop through all the instances of the model */
        for (here = model->BSIM3v32instances; here != NULL;
            here = here->BSIM3v32nextInstance)
        {
            InstCount++;
        }
    }
    InstArray = TMALLOC(BSIM3v32instance*, InstCount);
    model = (BSIM3v32model*)inModel;
    idx = 0;
    for (; model != NULL; model = model->BSIM3v32nextModel)
    {
        /* loop through all the instances of the model */
        for (here = model->BSIM3v32instances; here != NULL;
            here = here->BSIM3v32nextInstance)
        {
            InstArray[idx] = here;
            idx++;
        }
        /* set the array pointer and instance count into each model */
        model->BSIM3v32InstCount = InstCount;
        model->BSIM3v32InstanceArray = InstArray;
    }

#endif
    return(OK);
}
Ejemplo n.º 4
0
int
JFET2setup(SMPmatrix *matrix, GENmodel *inModel, CKTcircuit *ckt, int *states)
        /* load the diode structure with those pointers needed later 
         * for fast matrix loading 
         */
{
    JFET2model *model = (JFET2model*)inModel;
    JFET2instance *here;
    int error;
    CKTnode *tmp;

    /*  loop through all the diode models */
    for( ; model != NULL; model = model->JFET2nextModel ) {

        if( (model->JFET2type != NJF) && (model->JFET2type != PJF) ) {
            model->JFET2type = NJF;
        }
#define  PARAM(code,id,flag,ref,default,descrip) \
              if(!model->flag) {model->ref = default;}
#include "jfet2parm.h"

        /* loop through all the instances of the model */
        for (here = model->JFET2instances; here != NULL ;
                here=here->JFET2nextInstance) {
            
            if(!here->JFET2areaGiven) {
                here->JFET2area = 1;
            }
  
            if(!here->JFET2mGiven) {
                here->JFET2m = 1;
            }

            here->JFET2state = *states;
            *states += JFET2_STATE_COUNT + 1;

            if(model->JFET2rs != 0) {
                if(here->JFET2sourcePrimeNode == 0) {
                error = CKTmkVolt(ckt,&tmp,here->JFET2name,"source");
                if(error) return(error);
                here->JFET2sourcePrimeNode = tmp->number;

                if (ckt->CKTcopyNodesets) {
		    CKTnode *tmpNode;
		    IFuid tmpName;
            
                  if (CKTinst2Node(ckt,here,3,&tmpNode,&tmpName)==OK) {
                     if (tmpNode->nsGiven) {
                       tmp->nodeset=tmpNode->nodeset; 
                       tmp->nsGiven=tmpNode->nsGiven; 
                     }
                  }
                }
                }

            } else {
                here->JFET2sourcePrimeNode = here->JFET2sourceNode;
            }
            if(model->JFET2rd != 0) {
                if(here->JFET2drainPrimeNode == 0) {
                error = CKTmkVolt(ckt,&tmp,here->JFET2name,"drain");
                if(error) return(error);
                here->JFET2drainPrimeNode = tmp->number;

                if (ckt->CKTcopyNodesets) {
		    CKTnode *tmpNode;
		    IFuid tmpName;

                  if (CKTinst2Node(ckt,here,1,&tmpNode,&tmpName)==OK) {
                     if (tmpNode->nsGiven) {
                       tmp->nodeset=tmpNode->nodeset; 
                       tmp->nsGiven=tmpNode->nsGiven; 
                     }
                  }
                }
                }
                
            } else {
                here->JFET2drainPrimeNode = here->JFET2drainNode;
            }

/* macro to make elements with built in test for out of memory */
#define TSTALLOC(ptr,first,second) \
do { if((here->ptr = SMPmakeElt(matrix, here->first, here->second)) == NULL){\
    return(E_NOMEM);\
} } while(0)

            TSTALLOC(JFET2drainDrainPrimePtr,JFET2drainNode,JFET2drainPrimeNode);
            TSTALLOC(JFET2gateDrainPrimePtr,JFET2gateNode,JFET2drainPrimeNode);
            TSTALLOC(JFET2gateSourcePrimePtr,JFET2gateNode,JFET2sourcePrimeNode);
            TSTALLOC(JFET2sourceSourcePrimePtr,JFET2sourceNode,JFET2sourcePrimeNode);
            TSTALLOC(JFET2drainPrimeDrainPtr,JFET2drainPrimeNode,JFET2drainNode);
            TSTALLOC(JFET2drainPrimeGatePtr,JFET2drainPrimeNode,JFET2gateNode);
            TSTALLOC(JFET2drainPrimeSourcePrimePtr,JFET2drainPrimeNode,JFET2sourcePrimeNode);
            TSTALLOC(JFET2sourcePrimeGatePtr,JFET2sourcePrimeNode,JFET2gateNode);
            TSTALLOC(JFET2sourcePrimeSourcePtr,JFET2sourcePrimeNode,JFET2sourceNode);
            TSTALLOC(JFET2sourcePrimeDrainPrimePtr,JFET2sourcePrimeNode,JFET2drainPrimeNode);
            TSTALLOC(JFET2drainDrainPtr,JFET2drainNode,JFET2drainNode);
            TSTALLOC(JFET2gateGatePtr,JFET2gateNode,JFET2gateNode);
            TSTALLOC(JFET2sourceSourcePtr,JFET2sourceNode,JFET2sourceNode);
            TSTALLOC(JFET2drainPrimeDrainPrimePtr,JFET2drainPrimeNode,JFET2drainPrimeNode);
            TSTALLOC(JFET2sourcePrimeSourcePrimePtr,JFET2sourcePrimeNode,JFET2sourcePrimeNode);
        }
    }
    return(OK);
}
Ejemplo n.º 5
0
int HFET2setup(SMPmatrix *matrix, GENmodel *inModel, CKTcircuit *ckt, int *states)
{
    
  HFET2model *model = (HFET2model*)inModel;
  HFET2instance *here;
  int error;
  CKTnode *tmp;

  for( ; model != NULL; model = model->HFET2nextModel ) {
    if((TYPE != NHFET) && (TYPE != PHFET) )
      TYPE = NHFET;
    if(!model->HFET2cfGiven)
      CF = 0;
    if(!model->HFET2d1Given)
      D1 = 0.03e-6;
    if(!model->HFET2d2Given)
      D2 = 0.2e-6;
    if(!model->HFET2delGiven)
      DEL = 0.04;
    if(!model->HFET2deltaGiven)
      DELTA = 3.0;
    if(!model->HFET2deltadGiven)
      DELTAD = 4.5e-9;
    if(!model->HFET2diGiven)
      DI = 0.04e-6;
    if(!model->HFET2epsiGiven)
      EPSI = 12.244*8.85418e-12;
    if(!model->HFET2etaGiven)
      {
        if(TYPE == NHFET)
          ETA = 1.28;
        else
          ETA = 1.4;
      }
    if(!model->HFET2eta1Given)
      ETA1 = 2;
    if(!model->HFET2eta2Given)
      ETA2 = 2;
    if(!model->HFET2gammaGiven)
      GAMMA = 3.0;
    if(!model->HFET2ggrGiven)
      GGR = 0;
    if(!model->HFET2jsGiven)
      JS = 0;
    if(!model->HFET2klambdaGiven)
      KLAMBDA = 0;
    if(!model->HFET2kmuGiven)
      KMU = 0;
    if(!model->HFET2knmaxGiven)
      KNMAX = 0;
    if(!model->HFET2kvtoGiven)
      KVTO = 0;
    if(!model->HFET2lambdaGiven)
      LAMBDA = 0.15;
    if(!model->HFET2mGiven)
      M = 3.0;
    if(!model->HFET2mcGiven)
      MC = 3.0;
    if(!model->HFET2muGiven)
      {
        if(TYPE == NHFET)
          MU = 0.4;
        else
          MU = 0.03;
      }
    if(!model->HFET2nGiven)
      N = 5.0;
    if(!model->HFET2nmaxGiven)
      NMAX = 2e16;
    if(!model->HFET2pGiven)
      PP = 1;
    if(!model->HFET2rdGiven)
      RD = 0;
    if(!model->HFET2rdiGiven)
      RDI = 0;
    if(!model->HFET2rsGiven)
      RS = 0;
    if(!model->HFET2rsiGiven)
      RSI = 0;
    if(!model->HFET2sigma0Given)
      SIGMA0 = 0.057;
    if(!model->HFET2vsGiven)
      {
        if(TYPE == NHFET)
          VS = 1.5e5;
        else
          VS = 0.8e5;
      }
    if(!model->HFET2vsigmaGiven)
      VSIGMA = 0.1;
    if(!model->HFET2vsigmatGiven)
      VSIGMAT = 0.3;
    if(!model->HFET2vt1Given)
      /* initialized in HFET2temp */
      HFET2_VT1 = 0;
    if(!model->HFET2vt2Given)
      /* initialized in HFET2temp */
      VT2 = 0;
    if(!model->HFET2vtoGiven) {
      if(model->HFET2type == NHFET)
        VTO = 0.15;
      else
        VTO = -0.15;
    }
   
    /* loop through all the instances of the model */
   
    
    for (here = model->HFET2instances; here != NULL; 
         here=here->HFET2nextInstance) {
      
      CKTnode *tmpNode;
      IFuid tmpName;
   
      here->HFET2state = *states;
      *states += 13;
      
      if(!here->HFET2lengthGiven)
        L = 1e-6;      
      if(!here->HFET2widthGiven)
        W = 20e-6;
      if(!here->HFET2mGiven)
        here->HFET2m = 1.0;

      if(model->HFET2rs != 0) {
          if(here->HFET2sourcePrimeNode == 0) {
        error = CKTmkVolt(ckt,&tmp,here->HFET2name,"source");
        if(error) return(error);
        here->HFET2sourcePrimeNode = tmp->number;
        
        if (ckt->CKTcopyNodesets) {
                  if (CKTinst2Node(ckt,here,3,&tmpNode,&tmpName)==OK) {
                     if (tmpNode->nsGiven) {
                       tmp->nodeset=tmpNode->nodeset; 
                       tmp->nsGiven=tmpNode->nsGiven; 
                     }
                  }
                }
          }
        
      } else {
        here->HFET2sourcePrimeNode = here->HFET2sourceNode;
      }
      if(model->HFET2rd != 0) {
          if(here->HFET2drainPrimeNode == 0) {
        error = CKTmkVolt(ckt,&tmp,here->HFET2name,"drain");
        if(error) return(error);
        here->HFET2drainPrimeNode = tmp->number;
        
        if (ckt->CKTcopyNodesets) {
                  if (CKTinst2Node(ckt,here,1,&tmpNode,&tmpName)==OK) {
                     if (tmpNode->nsGiven) {
                       tmp->nodeset=tmpNode->nodeset; 
                       tmp->nsGiven=tmpNode->nsGiven; 
                     }
                  }
                }
          }
        
      } else {
        here->HFET2drainPrimeNode = here->HFET2drainNode;
      }

/* macro to make elements with built in test for out of memory */
#define TSTALLOC(ptr,first,second) \
do { if((here->ptr = SMPmakeElt(matrix, here->first, here->second)) == NULL){\
    return(E_NOMEM);\
} } while(0)

    TSTALLOC(HFET2drainDrainPrimePtr,HFET2drainNode,HFET2drainPrimeNode);
    TSTALLOC(HFET2gateDrainPrimePtr,HFET2gateNode,HFET2drainPrimeNode);
    TSTALLOC(HFET2gateSourcePrimePtr,HFET2gateNode,HFET2sourcePrimeNode);
    TSTALLOC(HFET2sourceSourcePrimePtr,HFET2sourceNode,HFET2sourcePrimeNode);
    TSTALLOC(HFET2drainPrimeDrainPtr,HFET2drainPrimeNode,HFET2drainNode);
    TSTALLOC(HFET2drainPrimeGatePtr,HFET2drainPrimeNode,HFET2gateNode);
    TSTALLOC(HFET2drainPriHFET2ourcePrimePtr,HFET2drainPrimeNode,HFET2sourcePrimeNode);
    TSTALLOC(HFET2sourcePrimeGatePtr,HFET2sourcePrimeNode,HFET2gateNode);
    TSTALLOC(HFET2sourcePriHFET2ourcePtr,HFET2sourcePrimeNode,HFET2sourceNode);
    TSTALLOC(HFET2sourcePrimeDrainPrimePtr,HFET2sourcePrimeNode,HFET2drainPrimeNode);
    TSTALLOC(HFET2drainDrainPtr,HFET2drainNode,HFET2drainNode);
    TSTALLOC(HFET2gateGatePtr,HFET2gateNode,HFET2gateNode);
    TSTALLOC(HFET2sourceSourcePtr,HFET2sourceNode,HFET2sourceNode);
    TSTALLOC(HFET2drainPrimeDrainPrimePtr,HFET2drainPrimeNode,HFET2drainPrimeNode);
    TSTALLOC(HFET2sourcePriHFET2ourcePrimePtr,HFET2sourcePrimeNode,HFET2sourcePrimeNode);
  
    }
  }
  return(OK);
  
}
Ejemplo n.º 6
0
int hicum2setup (SMPmatrix *matrix, GENmodel *inModel, CKTcircuit *ckt, int *states)
     /* load the hicum2 device structure with those pointers needed later 
      * for fast matrix loading 
      */
{
  register hicum2model *model = (hicum2model*)inModel;
  register hicum2instance *here;
  
  /*  loop through all the hicum2 device models */
  for ( ;model != NULL ;model = model->hicum2nextModel )
  {
    if(model->pnp_Given == FALSE) model->pnp=0;
    if(model->npn_Given == FALSE) model->npn=1;
    if(model->c10_Given == FALSE) model->c10=2.0E-30;
    if(model->qp0_Given == FALSE) model->qp0=2.0E-14;
    if(model->ich_Given == FALSE) model->ich=0.0;
    if(model->hfe_Given == FALSE) model->hfe=1.0;
    if(model->hfc_Given == FALSE) model->hfc=1.0;
    if(model->hjei_Given == FALSE) model->hjei=1.0;
    if(model->hjci_Given == FALSE) model->hjci=1.0;
    if(model->ibeis_Given == FALSE) model->ibeis=1.0E-18;
    if(model->mbei_Given == FALSE) model->mbei=1.0;
    if(model->ireis_Given == FALSE) model->ireis=0.0;
    if(model->mrei_Given == FALSE) model->mrei=2.0;
    if(model->ibeps_Given == FALSE) model->ibeps=0.0;
    if(model->mbep_Given == FALSE) model->mbep=1.0;
    if(model->ireps_Given == FALSE) model->ireps=0.0;
    if(model->mrep_Given == FALSE) model->mrep=2.0;
    if(model->mcf_Given == FALSE) model->mcf=1.0;
    if(model->tbhrec_Given == FALSE) model->tbhrec=0.0;
    if(model->ibcis_Given == FALSE) model->ibcis=1.0E-16;
    if(model->mbci_Given == FALSE) model->mbci=1.0;
    if(model->ibcxs_Given == FALSE) model->ibcxs=0.0;
    if(model->mbcx_Given == FALSE) model->mbcx=1.0;
    if(model->ibets_Given == FALSE) model->ibets=0.0;
    if(model->abet_Given == FALSE) model->abet=40;
    if(model->tunode_Given == FALSE) model->tunode=1;
    if(model->favl_Given == FALSE) model->favl=0.0;
    if(model->qavl_Given == FALSE) model->qavl=0.0;
    if(model->alfav_Given == FALSE) model->alfav=0.0;
    if(model->alqav_Given == FALSE) model->alqav=0.0;
    if(model->rbi0_Given == FALSE) model->rbi0=0.0;
    if(model->rbx_Given == FALSE) model->rbx=0.0;
    if(model->fgeo_Given == FALSE) model->fgeo=0.6557;
    if(model->fdqr0_Given == FALSE) model->fdqr0=0.0;
    if(model->fcrbi_Given == FALSE) model->fcrbi=0.0;
    if(model->fqi_Given == FALSE) model->fqi=1.0;
    if(model->re_Given == FALSE) model->re=0.0;
    if(model->rcx_Given == FALSE) model->rcx=0.0;
    if(model->itss_Given == FALSE) model->itss=0.0;
    if(model->msf_Given == FALSE) model->msf=1.0;
    if(model->iscs_Given == FALSE) model->iscs=0.0;
    if(model->msc_Given == FALSE) model->msc=1.0;
    if(model->tsf_Given == FALSE) model->tsf=0.0;
    if(model->rsu_Given == FALSE) model->rsu=0.0;
    if(model->csu_Given == FALSE) model->csu=0.0;
    if(model->cjei0_Given == FALSE) model->cjei0=1.0E-20;
    if(model->vdei_Given == FALSE) model->vdei=0.9;
    if(model->zei_Given == FALSE) model->zei=0.5;
    if(model->ajei_Given == FALSE) model->ajei=2.5;
    if(model->cjep0_Given == FALSE) model->cjep0=1.0E-20;
    if(model->vdep_Given == FALSE) model->vdep=0.9;
    if(model->zep_Given == FALSE) model->zep=0.5;
    if(model->ajep_Given == FALSE) model->ajep=2.5;
    if(model->cjci0_Given == FALSE) model->cjci0=1.0E-20;
    if(model->vdci_Given == FALSE) model->vdci=0.7;
    if(model->zci_Given == FALSE) model->zci=0.4;
    if(model->vptci_Given == FALSE) model->vptci=100;
    if(model->cjcx0_Given == FALSE) model->cjcx0=1.0E-20;
    if(model->vdcx_Given == FALSE) model->vdcx=0.7;
    if(model->zcx_Given == FALSE) model->zcx=0.4;
    if(model->vptcx_Given == FALSE) model->vptcx=100;
    if(model->fbcpar_Given == FALSE) model->fbcpar=0.0;
    if(model->fbepar_Given == FALSE) model->fbepar=1.0;
    if(model->cjs0_Given == FALSE) model->cjs0=0.0;
    if(model->vds_Given == FALSE) model->vds=0.6;
    if(model->zs_Given == FALSE) model->zs=0.5;
    if(model->vpts_Given == FALSE) model->vpts=100;
    if(model->t0_Given == FALSE) model->t0=0.0;
    if(model->dt0h_Given == FALSE) model->dt0h=0.0;
    if(model->tbvl_Given == FALSE) model->tbvl=0.0;
    if(model->tef0_Given == FALSE) model->tef0=0.0;
    if(model->gtfe_Given == FALSE) model->gtfe=1.0;
    if(model->thcs_Given == FALSE) model->thcs=0.0;
    if(model->ahc_Given == FALSE) model->ahc=0.1;
    if(model->fthc_Given == FALSE) model->fthc=0.0;
    if(model->rci0_Given == FALSE) model->rci0=150;
    if(model->vlim_Given == FALSE) model->vlim=0.5;
    if(model->vces_Given == FALSE) model->vces=0.1;
    if(model->vpt_Given == FALSE) model->vpt=0.0;
    if(model->tr_Given == FALSE) model->tr=0.0;
    if(model->cbepar_Given == FALSE) model->cbepar=0.0;
    if(model->cbcpar_Given == FALSE) model->cbcpar=0.0;
    if(model->alqf_Given == FALSE) model->alqf=0.0;
    if(model->alit_Given == FALSE) model->alit=0.0;
    if(model->flnqs_Given == FALSE) model->flnqs=0;
    if(model->kf_Given == FALSE) model->kf=0.0;
    if(model->af_Given == FALSE) model->af=2.0;
    if(model->cfbe_Given == FALSE) model->cfbe=(-1);
    if(model->latb_Given == FALSE) model->latb=0.0;
    if(model->latl_Given == FALSE) model->latl=0.0;
    if(model->vgb_Given == FALSE) model->vgb=1.17;
    if(model->alt0_Given == FALSE) model->alt0=0.0;
    if(model->kt0_Given == FALSE) model->kt0=0.0;
    if(model->zetaci_Given == FALSE) model->zetaci=0.0;
    if(model->alvs_Given == FALSE) model->alvs=0.0;
    if(model->alces_Given == FALSE) model->alces=0.0;
    if(model->zetarbi_Given == FALSE) model->zetarbi=0.0;
    if(model->zetarbx_Given == FALSE) model->zetarbx=0.0;
    if(model->zetarcx_Given == FALSE) model->zetarcx=0.0;
    if(model->zetare_Given == FALSE) model->zetare=0.0;
    if(model->zetacx_Given == FALSE) model->zetacx=1.0;
    if(model->vge_Given == FALSE) model->vge=1.17;
    if(model->vgc_Given == FALSE) model->vgc=1.17;
    if(model->vgs_Given == FALSE) model->vgs=1.17;
    if(model->f1vg_Given == FALSE) model->f1vg=(-1.02377e-4);
    if(model->f2vg_Given == FALSE) model->f2vg=4.3215e-4;
    if(model->zetact_Given == FALSE) model->zetact=3.0;
    if(model->zetabet_Given == FALSE) model->zetabet=3.5;
    if(model->alb_Given == FALSE) model->alb=0.0;
    if(model->flsh_Given == FALSE) model->flsh=0;
    if(model->rth_Given == FALSE) model->rth=0.0;
    if(model->cth_Given == FALSE) model->cth=0.0;
    if(model->flcomp_Given == FALSE) model->flcomp=0.0;
    if(model->tnom_Given == FALSE) model->tnom=27.0;
    if(model->dt_Given == FALSE) model->dt=0.0;
    for ( here = model->hicum2instances ;here != NULL ; here = here->hicum2nextInstance )
    {
      /* Internal Nodes */
      {
        here->ciNode = -1;
        here->eiNode = -1;
        here->bpNode = -1;
        here->biNode = -1;
        here->siNode = -1;
        here->xf1Node = -1;
        here->xf2Node = -1;
        here->xf3Node = -1;
        here->xf4Node = -1;
        here->n1Node = -1;
        here->n2Node = -1;
      }
      /* set states */
      here->state_xf4_GND = *states; *states += 2;
      here->state_xf3_GND = *states; *states += 2;
      here->state_xf2_GND = *states; *states += 2;
      here->state_xf1_GND = *states; *states += 2;
      here->state_tnode_GND = *states; *states += 2;
      here->state_si_s = *states; *states += 2;
      here->state_si_ci = *states; *states += 2;
      here->state_bp_e = *states; *states += 2;
      here->state_b_e = *states; *states += 2;
      here->state_bp_ei = *states; *states += 2;
      here->state_bp_bi = *states; *states += 2;
      here->state_bp_ci = *states; *states += 2;
      here->state_b_ci = *states; *states += 2;
      here->state_bi_ci = *states; *states += 2;
      here->state_bi_ei = *states; *states += 2;
      /* set Sparse Matrix Pointers */
      here->PTR_J_bi_bi_required=0;
      here->PTR_J_bi_ei_required=0;
      here->PTR_J_ei_bi_required=0;
      here->PTR_J_ei_ei_required=0;
      here->PTR_J_bi_ci_required=0;
      here->PTR_J_ci_bi_required=0;
      here->PTR_J_ci_ci_required=0;
      here->PTR_J_b_bi_required=0;
      here->PTR_J_b_ci_required=0;
      here->PTR_J_b_tnode_required=0;
      here->PTR_J_ci_tnode_required=0;
      here->PTR_J_b_bp_required=0;
      here->PTR_J_ci_bp_required=0;
      here->PTR_J_b_b_required=0;
      here->PTR_J_ci_b_required=0;
      here->PTR_J_b_ei_required=0;
      here->PTR_J_ci_ei_required=0;
      here->PTR_J_bp_bp_required=0;
      here->PTR_J_bp_ci_required=0;
      here->PTR_J_bp_tnode_required=0;
      here->PTR_J_bp_b_required=0;
      here->PTR_J_bp_ei_required=0;
      here->PTR_J_bp_bi_required=0;
      here->PTR_J_bi_tnode_required=0;
      here->PTR_J_bi_bp_required=0;
      here->PTR_J_ei_tnode_required=0;
      here->PTR_J_ei_bp_required=0;
      here->PTR_J_ei_ci_required=0;
      here->PTR_J_bi_xf4_required=0;
      here->PTR_J_ei_xf4_required=0;
      here->PTR_J_bp_si_required=0;
      here->PTR_J_si_si_required=0;
      here->PTR_J_si_ci_required=0;
      here->PTR_J_si_bp_required=0;
      here->PTR_J_si_tnode_required=0;
      here->PTR_J_b_e_required=0;
      here->PTR_J_e_b_required=0;
      here->PTR_J_e_e_required=0;
      here->PTR_J_bp_e_required=0;
      here->PTR_J_e_bp_required=0;
      here->PTR_J_si_bi_required=0;
      here->PTR_J_si_ei_required=0;
      here->PTR_J_ci_si_required=0;
      here->PTR_J_si_b_required=0;
      here->PTR_J_ci_xf2_required=0;
      here->PTR_J_ei_xf2_required=0;
      here->PTR_J_c_tnode_required=0;
      here->PTR_J_ci_c_required=0;
      here->PTR_J_c_ci_required=0;
      here->PTR_J_c_c_required=0;
      here->PTR_J_e_tnode_required=0;
      here->PTR_J_ei_e_required=0;
      here->PTR_J_e_ei_required=0;
      here->PTR_J_si_s_required=0;
      here->PTR_J_s_si_required=0;
      here->PTR_J_s_s_required=0;
      here->PTR_J_tnode_tnode_required=0;
      here->PTR_J_tnode_bi_required=0;
      here->PTR_J_tnode_ci_required=0;
      here->PTR_J_tnode_ei_required=0;
      here->PTR_J_tnode_si_required=0;
      here->PTR_J_tnode_bp_required=0;
      here->PTR_J_tnode_e_required=0;
      here->PTR_J_tnode_c_required=0;
      here->PTR_J_tnode_b_required=0;
      here->PTR_J_xf1_bi_required=0;
      here->PTR_J_xf1_ci_required=0;
      here->PTR_J_xf1_ei_required=0;
      here->PTR_J_xf1_tnode_required=0;
      here->PTR_J_xf1_xf2_required=0;
      here->PTR_J_xf1_xf1_required=0;
      here->PTR_J_xf2_bi_required=0;
      here->PTR_J_xf2_ci_required=0;
      here->PTR_J_xf2_tnode_required=0;
      here->PTR_J_xf2_ei_required=0;
      here->PTR_J_xf2_xf1_required=0;
      here->PTR_J_xf2_xf2_required=0;
      here->PTR_J_xf3_bi_required=0;
      here->PTR_J_xf3_ei_required=0;
      here->PTR_J_xf3_ci_required=0;
      here->PTR_J_xf3_tnode_required=0;
      here->PTR_J_xf3_xf4_required=0;
      here->PTR_J_xf3_xf3_required=0;
      here->PTR_J_xf4_bi_required=0;
      here->PTR_J_xf4_ci_required=0;
      here->PTR_J_xf4_tnode_required=0;
      here->PTR_J_xf4_ei_required=0;
      here->PTR_J_xf4_xf3_required=0;
      here->PTR_J_xf4_xf4_required=0;
      here->PTR_J_s_tnode_required=0;

      hicum2guesstopology(matrix,ckt,model,here);

      /* Internal Nodes */
      {
        int error;
        CKTnode *tmp;
        if(here->ciNode == -1)
        {
          error=CKTmkVolt(ckt,&tmp,here->hicum2name,"no info");
          if(error) return(error);
          here->ciNode = tmp->number;
        }
        if(here->eiNode == -1)
        {
          error=CKTmkVolt(ckt,&tmp,here->hicum2name,"no info");
          if(error) return(error);
          here->eiNode = tmp->number;
        }
        if(here->bpNode == -1)
        {
          error=CKTmkVolt(ckt,&tmp,here->hicum2name,"no info");
          if(error) return(error);
          here->bpNode = tmp->number;
        }
        if(here->biNode == -1)
        {
          error=CKTmkVolt(ckt,&tmp,here->hicum2name,"no info");
          if(error) return(error);
          here->biNode = tmp->number;
        }
        if(here->siNode == -1)
        {
          error=CKTmkVolt(ckt,&tmp,here->hicum2name,"no info");
          if(error) return(error);
          here->siNode = tmp->number;
        }
        if(here->xf1Node == -1)
        {
          error=CKTmkVolt(ckt,&tmp,here->hicum2name,"no info");
          if(error) return(error);
          here->xf1Node = tmp->number;
        }
        if(here->xf2Node == -1)
        {
          error=CKTmkVolt(ckt,&tmp,here->hicum2name,"no info");
          if(error) return(error);
          here->xf2Node = tmp->number;
        }
        if(here->xf3Node == -1)
        {
          error=CKTmkVolt(ckt,&tmp,here->hicum2name,"no info");
          if(error) return(error);
          here->xf3Node = tmp->number;
        }
        if(here->xf4Node == -1)
        {
          error=CKTmkVolt(ckt,&tmp,here->hicum2name,"no info");
          if(error) return(error);
          here->xf4Node = tmp->number;
        }
        if(here->n1Node == -1)
        {
          error=CKTmkVolt(ckt,&tmp,here->hicum2name,"no info");
          if(error) return(error);
          here->n1Node = tmp->number;
        }
        if(here->n2Node == -1)
        {
          error=CKTmkVolt(ckt,&tmp,here->hicum2name,"no info");
          if(error) return(error);
          here->n2Node = tmp->number;
        }
      }
        if(here->PTR_J_bi_bi_required==1)
        {
          here->PTR_J_bi_bi=SMPmakeElt(matrix,here->biNode,here->biNode);
        }
        if(here->PTR_J_bi_ei_required==1)
        {
          here->PTR_J_bi_ei=SMPmakeElt(matrix,here->biNode,here->eiNode);
        }
        if(here->PTR_J_ei_bi_required==1)
        {
          here->PTR_J_ei_bi=SMPmakeElt(matrix,here->eiNode,here->biNode);
        }
        if(here->PTR_J_ei_ei_required==1)
        {
          here->PTR_J_ei_ei=SMPmakeElt(matrix,here->eiNode,here->eiNode);
        }
        if(here->PTR_J_bi_ci_required==1)
        {
          here->PTR_J_bi_ci=SMPmakeElt(matrix,here->biNode,here->ciNode);
        }
        if(here->PTR_J_ci_bi_required==1)
        {
          here->PTR_J_ci_bi=SMPmakeElt(matrix,here->ciNode,here->biNode);
        }
        if(here->PTR_J_ci_ci_required==1)
        {
          here->PTR_J_ci_ci=SMPmakeElt(matrix,here->ciNode,here->ciNode);
        }
        if(here->PTR_J_b_bi_required==1)
        {
          here->PTR_J_b_bi=SMPmakeElt(matrix,here->bNode,here->biNode);
        }
        if(here->PTR_J_b_ci_required==1)
        {
          here->PTR_J_b_ci=SMPmakeElt(matrix,here->bNode,here->ciNode);
        }
        if(here->PTR_J_b_tnode_required==1)
        {
          here->PTR_J_b_tnode=SMPmakeElt(matrix,here->bNode,here->tnodeNode);
        }
        if(here->PTR_J_ci_tnode_required==1)
        {
          here->PTR_J_ci_tnode=SMPmakeElt(matrix,here->ciNode,here->tnodeNode);
        }
        if(here->PTR_J_b_bp_required==1)
        {
          here->PTR_J_b_bp=SMPmakeElt(matrix,here->bNode,here->bpNode);
        }
        if(here->PTR_J_ci_bp_required==1)
        {
          here->PTR_J_ci_bp=SMPmakeElt(matrix,here->ciNode,here->bpNode);
        }
        if(here->PTR_J_b_b_required==1)
        {
          here->PTR_J_b_b=SMPmakeElt(matrix,here->bNode,here->bNode);
        }
        if(here->PTR_J_ci_b_required==1)
        {
          here->PTR_J_ci_b=SMPmakeElt(matrix,here->ciNode,here->bNode);
        }
        if(here->PTR_J_b_ei_required==1)
        {
          here->PTR_J_b_ei=SMPmakeElt(matrix,here->bNode,here->eiNode);
        }
        if(here->PTR_J_ci_ei_required==1)
        {
          here->PTR_J_ci_ei=SMPmakeElt(matrix,here->ciNode,here->eiNode);
        }
        if(here->PTR_J_bp_bp_required==1)
        {
          here->PTR_J_bp_bp=SMPmakeElt(matrix,here->bpNode,here->bpNode);
        }
        if(here->PTR_J_bp_ci_required==1)
        {
          here->PTR_J_bp_ci=SMPmakeElt(matrix,here->bpNode,here->ciNode);
        }
        if(here->PTR_J_bp_tnode_required==1)
        {
          here->PTR_J_bp_tnode=SMPmakeElt(matrix,here->bpNode,here->tnodeNode);
        }
        if(here->PTR_J_bp_b_required==1)
        {
          here->PTR_J_bp_b=SMPmakeElt(matrix,here->bpNode,here->bNode);
        }
        if(here->PTR_J_bp_ei_required==1)
        {
          here->PTR_J_bp_ei=SMPmakeElt(matrix,here->bpNode,here->eiNode);
        }
        if(here->PTR_J_bp_bi_required==1)
        {
          here->PTR_J_bp_bi=SMPmakeElt(matrix,here->bpNode,here->biNode);
        }
        if(here->PTR_J_bi_tnode_required==1)
        {
          here->PTR_J_bi_tnode=SMPmakeElt(matrix,here->biNode,here->tnodeNode);
        }
        if(here->PTR_J_bi_bp_required==1)
        {
          here->PTR_J_bi_bp=SMPmakeElt(matrix,here->biNode,here->bpNode);
        }
        if(here->PTR_J_ei_tnode_required==1)
        {
          here->PTR_J_ei_tnode=SMPmakeElt(matrix,here->eiNode,here->tnodeNode);
        }
        if(here->PTR_J_ei_bp_required==1)
        {
          here->PTR_J_ei_bp=SMPmakeElt(matrix,here->eiNode,here->bpNode);
        }
        if(here->PTR_J_ei_ci_required==1)
        {
          here->PTR_J_ei_ci=SMPmakeElt(matrix,here->eiNode,here->ciNode);
        }
        if(here->PTR_J_bi_xf4_required==1)
        {
          here->PTR_J_bi_xf4=SMPmakeElt(matrix,here->biNode,here->xf4Node);
        }
        if(here->PTR_J_ei_xf4_required==1)
        {
          here->PTR_J_ei_xf4=SMPmakeElt(matrix,here->eiNode,here->xf4Node);
        }
        if(here->PTR_J_bp_si_required==1)
        {
          here->PTR_J_bp_si=SMPmakeElt(matrix,here->bpNode,here->siNode);
        }
        if(here->PTR_J_si_si_required==1)
        {
          here->PTR_J_si_si=SMPmakeElt(matrix,here->siNode,here->siNode);
        }
        if(here->PTR_J_si_ci_required==1)
        {
          here->PTR_J_si_ci=SMPmakeElt(matrix,here->siNode,here->ciNode);
        }
        if(here->PTR_J_si_bp_required==1)
        {
          here->PTR_J_si_bp=SMPmakeElt(matrix,here->siNode,here->bpNode);
        }
        if(here->PTR_J_si_tnode_required==1)
        {
          here->PTR_J_si_tnode=SMPmakeElt(matrix,here->siNode,here->tnodeNode);
        }
        if(here->PTR_J_b_e_required==1)
        {
          here->PTR_J_b_e=SMPmakeElt(matrix,here->bNode,here->eNode);
        }
        if(here->PTR_J_e_b_required==1)
        {
          here->PTR_J_e_b=SMPmakeElt(matrix,here->eNode,here->bNode);
        }
        if(here->PTR_J_e_e_required==1)
        {
          here->PTR_J_e_e=SMPmakeElt(matrix,here->eNode,here->eNode);
        }
        if(here->PTR_J_bp_e_required==1)
        {
          here->PTR_J_bp_e=SMPmakeElt(matrix,here->bpNode,here->eNode);
        }
        if(here->PTR_J_e_bp_required==1)
        {
          here->PTR_J_e_bp=SMPmakeElt(matrix,here->eNode,here->bpNode);
        }
        if(here->PTR_J_si_bi_required==1)
        {
          here->PTR_J_si_bi=SMPmakeElt(matrix,here->siNode,here->biNode);
        }
        if(here->PTR_J_si_ei_required==1)
        {
          here->PTR_J_si_ei=SMPmakeElt(matrix,here->siNode,here->eiNode);
        }
        if(here->PTR_J_ci_si_required==1)
        {
          here->PTR_J_ci_si=SMPmakeElt(matrix,here->ciNode,here->siNode);
        }
        if(here->PTR_J_si_b_required==1)
        {
          here->PTR_J_si_b=SMPmakeElt(matrix,here->siNode,here->bNode);
        }
        if(here->PTR_J_ci_xf2_required==1)
        {
          here->PTR_J_ci_xf2=SMPmakeElt(matrix,here->ciNode,here->xf2Node);
        }
        if(here->PTR_J_ei_xf2_required==1)
        {
          here->PTR_J_ei_xf2=SMPmakeElt(matrix,here->eiNode,here->xf2Node);
        }
        if(here->PTR_J_c_tnode_required==1)
        {
          here->PTR_J_c_tnode=SMPmakeElt(matrix,here->cNode,here->tnodeNode);
        }
        if(here->PTR_J_ci_c_required==1)
        {
          here->PTR_J_ci_c=SMPmakeElt(matrix,here->ciNode,here->cNode);
        }
        if(here->PTR_J_c_ci_required==1)
        {
          here->PTR_J_c_ci=SMPmakeElt(matrix,here->cNode,here->ciNode);
        }
        if(here->PTR_J_c_c_required==1)
        {
          here->PTR_J_c_c=SMPmakeElt(matrix,here->cNode,here->cNode);
        }
        if(here->PTR_J_e_tnode_required==1)
        {
          here->PTR_J_e_tnode=SMPmakeElt(matrix,here->eNode,here->tnodeNode);
        }
        if(here->PTR_J_ei_e_required==1)
        {
          here->PTR_J_ei_e=SMPmakeElt(matrix,here->eiNode,here->eNode);
        }
        if(here->PTR_J_e_ei_required==1)
        {
          here->PTR_J_e_ei=SMPmakeElt(matrix,here->eNode,here->eiNode);
        }
        if(here->PTR_J_si_s_required==1)
        {
          here->PTR_J_si_s=SMPmakeElt(matrix,here->siNode,here->sNode);
        }
        if(here->PTR_J_s_si_required==1)
        {
          here->PTR_J_s_si=SMPmakeElt(matrix,here->sNode,here->siNode);
        }
        if(here->PTR_J_s_s_required==1)
        {
          here->PTR_J_s_s=SMPmakeElt(matrix,here->sNode,here->sNode);
        }
        if(here->PTR_J_tnode_tnode_required==1)
        {
          here->PTR_J_tnode_tnode=SMPmakeElt(matrix,here->tnodeNode,here->tnodeNode);
        }
        if(here->PTR_J_tnode_bi_required==1)
        {
          here->PTR_J_tnode_bi=SMPmakeElt(matrix,here->tnodeNode,here->biNode);
        }
        if(here->PTR_J_tnode_ci_required==1)
        {
          here->PTR_J_tnode_ci=SMPmakeElt(matrix,here->tnodeNode,here->ciNode);
        }
        if(here->PTR_J_tnode_ei_required==1)
        {
          here->PTR_J_tnode_ei=SMPmakeElt(matrix,here->tnodeNode,here->eiNode);
        }
        if(here->PTR_J_tnode_si_required==1)
        {
          here->PTR_J_tnode_si=SMPmakeElt(matrix,here->tnodeNode,here->siNode);
        }
        if(here->PTR_J_tnode_bp_required==1)
        {
          here->PTR_J_tnode_bp=SMPmakeElt(matrix,here->tnodeNode,here->bpNode);
        }
        if(here->PTR_J_tnode_e_required==1)
        {
          here->PTR_J_tnode_e=SMPmakeElt(matrix,here->tnodeNode,here->eNode);
        }
        if(here->PTR_J_tnode_c_required==1)
        {
          here->PTR_J_tnode_c=SMPmakeElt(matrix,here->tnodeNode,here->cNode);
        }
        if(here->PTR_J_tnode_b_required==1)
        {
          here->PTR_J_tnode_b=SMPmakeElt(matrix,here->tnodeNode,here->bNode);
        }
        if(here->PTR_J_xf1_bi_required==1)
        {
          here->PTR_J_xf1_bi=SMPmakeElt(matrix,here->xf1Node,here->biNode);
        }
        if(here->PTR_J_xf1_ci_required==1)
        {
          here->PTR_J_xf1_ci=SMPmakeElt(matrix,here->xf1Node,here->ciNode);
        }
        if(here->PTR_J_xf1_ei_required==1)
        {
          here->PTR_J_xf1_ei=SMPmakeElt(matrix,here->xf1Node,here->eiNode);
        }
        if(here->PTR_J_xf1_tnode_required==1)
        {
          here->PTR_J_xf1_tnode=SMPmakeElt(matrix,here->xf1Node,here->tnodeNode);
        }
        if(here->PTR_J_xf1_xf2_required==1)
        {
          here->PTR_J_xf1_xf2=SMPmakeElt(matrix,here->xf1Node,here->xf2Node);
        }
        if(here->PTR_J_xf1_xf1_required==1)
        {
          here->PTR_J_xf1_xf1=SMPmakeElt(matrix,here->xf1Node,here->xf1Node);
        }
        if(here->PTR_J_xf2_bi_required==1)
        {
          here->PTR_J_xf2_bi=SMPmakeElt(matrix,here->xf2Node,here->biNode);
        }
        if(here->PTR_J_xf2_ci_required==1)
        {
          here->PTR_J_xf2_ci=SMPmakeElt(matrix,here->xf2Node,here->ciNode);
        }
        if(here->PTR_J_xf2_tnode_required==1)
        {
          here->PTR_J_xf2_tnode=SMPmakeElt(matrix,here->xf2Node,here->tnodeNode);
        }
        if(here->PTR_J_xf2_ei_required==1)
        {
          here->PTR_J_xf2_ei=SMPmakeElt(matrix,here->xf2Node,here->eiNode);
        }
        if(here->PTR_J_xf2_xf1_required==1)
        {
          here->PTR_J_xf2_xf1=SMPmakeElt(matrix,here->xf2Node,here->xf1Node);
        }
        if(here->PTR_J_xf2_xf2_required==1)
        {
          here->PTR_J_xf2_xf2=SMPmakeElt(matrix,here->xf2Node,here->xf2Node);
        }
        if(here->PTR_J_xf3_bi_required==1)
        {
          here->PTR_J_xf3_bi=SMPmakeElt(matrix,here->xf3Node,here->biNode);
        }
        if(here->PTR_J_xf3_ei_required==1)
        {
          here->PTR_J_xf3_ei=SMPmakeElt(matrix,here->xf3Node,here->eiNode);
        }
        if(here->PTR_J_xf3_ci_required==1)
        {
          here->PTR_J_xf3_ci=SMPmakeElt(matrix,here->xf3Node,here->ciNode);
        }
        if(here->PTR_J_xf3_tnode_required==1)
        {
          here->PTR_J_xf3_tnode=SMPmakeElt(matrix,here->xf3Node,here->tnodeNode);
        }
        if(here->PTR_J_xf3_xf4_required==1)
        {
          here->PTR_J_xf3_xf4=SMPmakeElt(matrix,here->xf3Node,here->xf4Node);
        }
        if(here->PTR_J_xf3_xf3_required==1)
        {
          here->PTR_J_xf3_xf3=SMPmakeElt(matrix,here->xf3Node,here->xf3Node);
        }
        if(here->PTR_J_xf4_bi_required==1)
        {
          here->PTR_J_xf4_bi=SMPmakeElt(matrix,here->xf4Node,here->biNode);
        }
        if(here->PTR_J_xf4_ci_required==1)
        {
          here->PTR_J_xf4_ci=SMPmakeElt(matrix,here->xf4Node,here->ciNode);
        }
        if(here->PTR_J_xf4_tnode_required==1)
        {
          here->PTR_J_xf4_tnode=SMPmakeElt(matrix,here->xf4Node,here->tnodeNode);
        }
        if(here->PTR_J_xf4_ei_required==1)
        {
          here->PTR_J_xf4_ei=SMPmakeElt(matrix,here->xf4Node,here->eiNode);
        }
        if(here->PTR_J_xf4_xf3_required==1)
        {
          here->PTR_J_xf4_xf3=SMPmakeElt(matrix,here->xf4Node,here->xf3Node);
        }
        if(here->PTR_J_xf4_xf4_required==1)
        {
          here->PTR_J_xf4_xf4=SMPmakeElt(matrix,here->xf4Node,here->xf4Node);
        }
        if(here->PTR_J_s_tnode_required==1)
        {
          here->PTR_J_s_tnode=SMPmakeElt(matrix,here->sNode,here->tnodeNode);
        }

    }
  }
  return(OK);
}
Ejemplo n.º 7
0
/* ARGSUSED */
int
TRAsetup(SMPmatrix *matrix, GENmodel *inModel, CKTcircuit *ckt, int *state)
        /* load the transmission line structure with those pointers needed later
         * for fast matrix loading 
         */
{
    TRAmodel *model = (TRAmodel *)inModel;
    TRAinstance *here;
    int error;
    CKTnode *tmp;

    NG_IGNORE(state);

    /*  loop through all the transmission line models */
    for( ; model != NULL; model = model->TRAnextModel ) {

        /* loop through all the instances of the model */
        for (here = model->TRAinstances; here != NULL ;
                here=here->TRAnextInstance) {
            
            if(here->TRAbrEq1==0) {
                error = CKTmkVolt(ckt,&tmp,here->TRAname,"i1");
                if(error) return(error);
                here->TRAbrEq1 = tmp->number;
            }

            if(here->TRAbrEq2==0) {
                error = CKTmkVolt(ckt,&tmp,here->TRAname,"i2");
                if(error) return(error);
                here->TRAbrEq2 = tmp->number;
            }

            if(here->TRAintNode1==0) {
                error = CKTmkVolt(ckt,&tmp,here->TRAname,"int1");
                if(error) return(error);
                here->TRAintNode1 = tmp->number;
            }

            if(here->TRAintNode2==0) {
                error = CKTmkVolt(ckt,&tmp,here->TRAname,"int2");
                if(error) return(error);
                here->TRAintNode2 = tmp->number;
            }

            /* allocate the delay table */
            here->TRAdelays = TMALLOC(double, 15);
            here->TRAallocDelay = 4;

/* macro to make elements with built in test for out of memory */
#define TSTALLOC(ptr,first,second) \
do { if((here->ptr = SMPmakeElt(matrix, here->first, here->second)) == NULL){\
    return(E_NOMEM);\
} } while(0)

            TSTALLOC(TRAibr1Ibr2Ptr, TRAbrEq1, TRAbrEq2);
            TSTALLOC(TRAibr1Int1Ptr, TRAbrEq1, TRAintNode1);
            TSTALLOC(TRAibr1Neg1Ptr, TRAbrEq1, TRAnegNode1);
            TSTALLOC(TRAibr1Neg2Ptr, TRAbrEq1, TRAnegNode2);
            TSTALLOC(TRAibr1Pos2Ptr, TRAbrEq1, TRAposNode2);
            TSTALLOC(TRAibr2Ibr1Ptr, TRAbrEq2, TRAbrEq1);
            TSTALLOC(TRAibr2Int2Ptr, TRAbrEq2, TRAintNode2);
            TSTALLOC(TRAibr2Neg1Ptr, TRAbrEq2, TRAnegNode1);
            TSTALLOC(TRAibr2Neg2Ptr, TRAbrEq2, TRAnegNode2);
            TSTALLOC(TRAibr2Pos1Ptr, TRAbrEq2, TRAposNode1);
            TSTALLOC(TRAint1Ibr1Ptr, TRAintNode1, TRAbrEq1);
            TSTALLOC(TRAint1Int1Ptr, TRAintNode1, TRAintNode1);
            TSTALLOC(TRAint1Pos1Ptr, TRAintNode1, TRAposNode1);
            TSTALLOC(TRAint2Ibr2Ptr, TRAintNode2, TRAbrEq2);
            TSTALLOC(TRAint2Int2Ptr, TRAintNode2, TRAintNode2);
            TSTALLOC(TRAint2Pos2Ptr, TRAintNode2, TRAposNode2);
            TSTALLOC(TRAneg1Ibr1Ptr, TRAnegNode1, TRAbrEq1);
            TSTALLOC(TRAneg2Ibr2Ptr, TRAnegNode2, TRAbrEq2);
            TSTALLOC(TRApos1Int1Ptr, TRAposNode1, TRAintNode1);
            TSTALLOC(TRApos1Pos1Ptr, TRAposNode1, TRAposNode1);
            TSTALLOC(TRApos2Int2Ptr, TRAposNode2, TRAintNode2);
            TSTALLOC(TRApos2Pos2Ptr, TRAposNode2, TRAposNode2);

            if(!here->TRAnlGiven) {
                here->TRAnl = .25;
            }
            if(!here->TRAfGiven) {
                here->TRAf = 1e9;
            }
            if(!here->TRAreltolGiven) {
                here->TRAreltol = 1;
            }
            if(!here->TRAabstolGiven) {
                here->TRAabstol = 1;
            }
            if(!here->TRAimpedGiven) {
                SPfrontEnd->IFerror (ERR_FATAL,
                        "%s: transmission line z0 must be given",
                        &(here->TRAname));
                return(E_BADPARM);
            }
        }
    }
    return(OK);
}
Ejemplo n.º 8
0
int psp102setup (SMPmatrix *matrix, GENmodel *inModel, CKTcircuit *ckt, int *states)
     /* load the psp102 device structure with those pointers needed later 
      * for fast matrix loading 
      */
{
  register psp102model *model = (psp102model*)inModel;
  register psp102instance *here;
  
  /*  loop through all the psp102 device models */
  for ( ;model != NULL ;model = model->psp102nextModel )
  {
    if(model->pmos_Given == FALSE) model->pmos=1;
    if(model->nmos_Given == FALSE) model->nmos=1;
    if(model->LEVEL_Given == FALSE) model->LEVEL=1020;
    if(model->TYPE_Given == FALSE) model->TYPE=1.0;
    if(model->TR_Given == FALSE) model->TR=21.0;
    if(model->SWIGATE_Given == FALSE) model->SWIGATE=0.0;
    if(model->SWIMPACT_Given == FALSE) model->SWIMPACT=0.0;
    if(model->SWGIDL_Given == FALSE) model->SWGIDL=0.0;
    if(model->SWJUNCAP_Given == FALSE) model->SWJUNCAP=0.0;
    if(model->QMC_Given == FALSE) model->QMC=1.0;
    if(model->LVARO_Given == FALSE) model->LVARO=0.0;
    if(model->LVARL_Given == FALSE) model->LVARL=0.0;
    if(model->LVARW_Given == FALSE) model->LVARW=0.0;
    if(model->LAP_Given == FALSE) model->LAP=0.0;
    if(model->WVARO_Given == FALSE) model->WVARO=0.0;
    if(model->WVARL_Given == FALSE) model->WVARL=0.0;
    if(model->WVARW_Given == FALSE) model->WVARW=0.0;
    if(model->WOT_Given == FALSE) model->WOT=0.0;
    if(model->DLQ_Given == FALSE) model->DLQ=0.0;
    if(model->DWQ_Given == FALSE) model->DWQ=0.0;
    if(model->VFBO_Given == FALSE) model->VFBO=(-1.0);
    if(model->VFBL_Given == FALSE) model->VFBL=0.0;
    if(model->VFBW_Given == FALSE) model->VFBW=0.0;
    if(model->VFBLW_Given == FALSE) model->VFBLW=0.0;
    if(model->STVFBO_Given == FALSE) model->STVFBO=5e-4;
    if(model->STVFBL_Given == FALSE) model->STVFBL=0.0;
    if(model->STVFBW_Given == FALSE) model->STVFBW=0.0;
    if(model->STVFBLW_Given == FALSE) model->STVFBLW=0.0;
    if(model->TOXO_Given == FALSE) model->TOXO=2e-9;
    if(model->NSUBO_Given == FALSE) model->NSUBO=3e23;
    if(model->NSUBW_Given == FALSE) model->NSUBW=0.0;
    if(model->WSEG_Given == FALSE) model->WSEG=1e-8;
    if(model->NPCK_Given == FALSE) model->NPCK=1e24;
    if(model->NPCKW_Given == FALSE) model->NPCKW=0.0;
    if(model->WSEGP_Given == FALSE) model->WSEGP=1e-8;
    if(model->LPCK_Given == FALSE) model->LPCK=1e-8;
    if(model->LPCKW_Given == FALSE) model->LPCKW=0.0;
    if(model->FOL1_Given == FALSE) model->FOL1=0.0;
    if(model->FOL2_Given == FALSE) model->FOL2=0.0;
    if(model->VNSUBO_Given == FALSE) model->VNSUBO=0.0;
    if(model->NSLPO_Given == FALSE) model->NSLPO=0.05;
    if(model->DNSUBO_Given == FALSE) model->DNSUBO=0.0;
    if(model->DPHIBO_Given == FALSE) model->DPHIBO=0.0;
    if(model->DPHIBL_Given == FALSE) model->DPHIBL=0.0;
    if(model->DPHIBLEXP_Given == FALSE) model->DPHIBLEXP=1.0;
    if(model->DPHIBW_Given == FALSE) model->DPHIBW=0.0;
    if(model->DPHIBLW_Given == FALSE) model->DPHIBLW=0.0;
    if(model->NPO_Given == FALSE) model->NPO=1e26;
    if(model->NPL_Given == FALSE) model->NPL=0.0;
    if(model->CTO_Given == FALSE) model->CTO=0.0;
    if(model->CTL_Given == FALSE) model->CTL=0.0;
    if(model->CTLEXP_Given == FALSE) model->CTLEXP=1.0;
    if(model->CTW_Given == FALSE) model->CTW=0.0;
    if(model->CTLW_Given == FALSE) model->CTLW=0.0;
    if(model->TOXOVO_Given == FALSE) model->TOXOVO=2e-9;
    if(model->LOV_Given == FALSE) model->LOV=0;
    if(model->NOVO_Given == FALSE) model->NOVO=5e25;
    if(model->CFL_Given == FALSE) model->CFL=0.0;
    if(model->CFLEXP_Given == FALSE) model->CFLEXP=2.0;
    if(model->CFW_Given == FALSE) model->CFW=0.0;
    if(model->CFBO_Given == FALSE) model->CFBO=0.0;
    if(model->UO_Given == FALSE) model->UO=5e-2;
    if(model->FBET1_Given == FALSE) model->FBET1=0.0;
    if(model->FBET1W_Given == FALSE) model->FBET1W=0.0;
    if(model->LP1_Given == FALSE) model->LP1=1e-8;
    if(model->LP1W_Given == FALSE) model->LP1W=0.0;
    if(model->FBET2_Given == FALSE) model->FBET2=0.0;
    if(model->LP2_Given == FALSE) model->LP2=1e-8;
    if(model->BETW1_Given == FALSE) model->BETW1=0.0;
    if(model->BETW2_Given == FALSE) model->BETW2=0.0;
    if(model->WBET_Given == FALSE) model->WBET=1e-9;
    if(model->STBETO_Given == FALSE) model->STBETO=1.0;
    if(model->STBETL_Given == FALSE) model->STBETL=0.0;
    if(model->STBETW_Given == FALSE) model->STBETW=0.0;
    if(model->STBETLW_Given == FALSE) model->STBETLW=0.0;
    if(model->MUEO_Given == FALSE) model->MUEO=0.5;
    if(model->MUEW_Given == FALSE) model->MUEW=0.0;
    if(model->STMUEO_Given == FALSE) model->STMUEO=0.0;
    if(model->THEMUO_Given == FALSE) model->THEMUO=1.5;
    if(model->STTHEMUO_Given == FALSE) model->STTHEMUO=1.5;
    if(model->CSO_Given == FALSE) model->CSO=0.0;
    if(model->CSL_Given == FALSE) model->CSL=0.0;
    if(model->CSLEXP_Given == FALSE) model->CSLEXP=0.0;
    if(model->CSW_Given == FALSE) model->CSW=0.0;
    if(model->CSLW_Given == FALSE) model->CSLW=0.0;
    if(model->STCSO_Given == FALSE) model->STCSO=0.0;
    if(model->XCORO_Given == FALSE) model->XCORO=0.0;
    if(model->XCORL_Given == FALSE) model->XCORL=0.0;
    if(model->XCORW_Given == FALSE) model->XCORW=0.0;
    if(model->XCORLW_Given == FALSE) model->XCORLW=0.0;
    if(model->STXCORO_Given == FALSE) model->STXCORO=0.0;
    if(model->FETAO_Given == FALSE) model->FETAO=1.0;
    if(model->RSW1_Given == FALSE) model->RSW1=2.5e3;
    if(model->RSW2_Given == FALSE) model->RSW2=0.0;
    if(model->STRSO_Given == FALSE) model->STRSO=1.0;
    if(model->RSBO_Given == FALSE) model->RSBO=0.0;
    if(model->RSGO_Given == FALSE) model->RSGO=0.0;
    if(model->THESATO_Given == FALSE) model->THESATO=0.0;
    if(model->THESATL_Given == FALSE) model->THESATL=0.05;
    if(model->THESATLEXP_Given == FALSE) model->THESATLEXP=1.0;
    if(model->THESATW_Given == FALSE) model->THESATW=0.0;
    if(model->THESATLW_Given == FALSE) model->THESATLW=0.0;
    if(model->STTHESATO_Given == FALSE) model->STTHESATO=1.0;
    if(model->STTHESATL_Given == FALSE) model->STTHESATL=0.0;
    if(model->STTHESATW_Given == FALSE) model->STTHESATW=0.0;
    if(model->STTHESATLW_Given == FALSE) model->STTHESATLW=0.0;
    if(model->THESATBO_Given == FALSE) model->THESATBO=0.0;
    if(model->THESATGO_Given == FALSE) model->THESATGO=0.0;
    if(model->AXO_Given == FALSE) model->AXO=18;
    if(model->AXL_Given == FALSE) model->AXL=0.4;
    if(model->ALPL_Given == FALSE) model->ALPL=5e-4;
    if(model->ALPLEXP_Given == FALSE) model->ALPLEXP=1.0;
    if(model->ALPW_Given == FALSE) model->ALPW=0.0;
    if(model->ALP1L1_Given == FALSE) model->ALP1L1=0.0;
    if(model->ALP1LEXP_Given == FALSE) model->ALP1LEXP=0.5;
    if(model->ALP1L2_Given == FALSE) model->ALP1L2=0.0;
    if(model->ALP1W_Given == FALSE) model->ALP1W=0.0;
    if(model->ALP2L1_Given == FALSE) model->ALP2L1=0.0;
    if(model->ALP2LEXP_Given == FALSE) model->ALP2LEXP=0.5;
    if(model->ALP2L2_Given == FALSE) model->ALP2L2=0.0;
    if(model->ALP2W_Given == FALSE) model->ALP2W=0.0;
    if(model->VPO_Given == FALSE) model->VPO=0.05;
    if(model->A1O_Given == FALSE) model->A1O=1.0;
    if(model->A1L_Given == FALSE) model->A1L=0.0;
    if(model->A1W_Given == FALSE) model->A1W=0.0;
    if(model->A2O_Given == FALSE) model->A2O=10;
    if(model->STA2O_Given == FALSE) model->STA2O=0.0;
    if(model->A3O_Given == FALSE) model->A3O=1.0;
    if(model->A3L_Given == FALSE) model->A3L=0.0;
    if(model->A3W_Given == FALSE) model->A3W=0.0;
    if(model->A4O_Given == FALSE) model->A4O=0.0;
    if(model->A4L_Given == FALSE) model->A4L=0.0;
    if(model->A4W_Given == FALSE) model->A4W=0.0;
    if(model->GCOO_Given == FALSE) model->GCOO=0.0;
    if(model->IGINVLW_Given == FALSE) model->IGINVLW=0.0;
    if(model->IGOVW_Given == FALSE) model->IGOVW=0.0;
    if(model->STIGO_Given == FALSE) model->STIGO=2.0;
    if(model->GC2O_Given == FALSE) model->GC2O=0.375;
    if(model->GC3O_Given == FALSE) model->GC3O=0.063;
    if(model->CHIBO_Given == FALSE) model->CHIBO=3.1;
    if(model->AGIDLW_Given == FALSE) model->AGIDLW=0.0;
    if(model->BGIDLO_Given == FALSE) model->BGIDLO=41;
    if(model->STBGIDLO_Given == FALSE) model->STBGIDLO=0.0;
    if(model->CGIDLO_Given == FALSE) model->CGIDLO=0.0;
    if(model->CGBOVL_Given == FALSE) model->CGBOVL=0.0;
    if(model->CFRW_Given == FALSE) model->CFRW=0.0;
    if(model->FNTO_Given == FALSE) model->FNTO=1.0;
    if(model->NFALW_Given == FALSE) model->NFALW=8e22;
    if(model->NFBLW_Given == FALSE) model->NFBLW=3e7;
    if(model->NFCLW_Given == FALSE) model->NFCLW=0.0;
    if(model->DTA_Given == FALSE) model->DTA=0;
    if(model->SAREF_Given == FALSE) model->SAREF=1.0e-6;
    if(model->SBREF_Given == FALSE) model->SBREF=1.0e-6;
    if(model->WLOD_Given == FALSE) model->WLOD=0;
    if(model->KUO_Given == FALSE) model->KUO=0;
    if(model->KVSAT_Given == FALSE) model->KVSAT=0;
    if(model->TKUO_Given == FALSE) model->TKUO=0;
    if(model->LKUO_Given == FALSE) model->LKUO=0;
    if(model->WKUO_Given == FALSE) model->WKUO=0;
    if(model->PKUO_Given == FALSE) model->PKUO=0;
    if(model->LLODKUO_Given == FALSE) model->LLODKUO=0;
    if(model->WLODKUO_Given == FALSE) model->WLODKUO=0;
    if(model->KVTHO_Given == FALSE) model->KVTHO=0;
    if(model->LKVTHO_Given == FALSE) model->LKVTHO=0;
    if(model->WKVTHO_Given == FALSE) model->WKVTHO=0;
    if(model->PKVTHO_Given == FALSE) model->PKVTHO=0;
    if(model->LLODVTH_Given == FALSE) model->LLODVTH=0;
    if(model->WLODVTH_Given == FALSE) model->WLODVTH=0;
    if(model->STETAO_Given == FALSE) model->STETAO=0;
    if(model->LODETAO_Given == FALSE) model->LODETAO=1.0;
    if(model->TRJ_Given == FALSE) model->TRJ=21;
    if(model->IMAX_Given == FALSE) model->IMAX=1000;
    if(model->CJORBOT_Given == FALSE) model->CJORBOT=1E-3;
    if(model->CJORSTI_Given == FALSE) model->CJORSTI=1E-9;
    if(model->CJORGAT_Given == FALSE) model->CJORGAT=1E-9;
    if(model->VBIRBOT_Given == FALSE) model->VBIRBOT=1;
    if(model->VBIRSTI_Given == FALSE) model->VBIRSTI=1;
    if(model->VBIRGAT_Given == FALSE) model->VBIRGAT=1;
    if(model->PBOT_Given == FALSE) model->PBOT=0.5;
    if(model->PSTI_Given == FALSE) model->PSTI=0.5;
    if(model->PGAT_Given == FALSE) model->PGAT=0.5;
    if(model->PHIGBOT_Given == FALSE) model->PHIGBOT=1.16;
    if(model->PHIGSTI_Given == FALSE) model->PHIGSTI=1.16;
    if(model->PHIGGAT_Given == FALSE) model->PHIGGAT=1.16;
    if(model->IDSATRBOT_Given == FALSE) model->IDSATRBOT=1E-12;
    if(model->IDSATRSTI_Given == FALSE) model->IDSATRSTI=1E-18;
    if(model->IDSATRGAT_Given == FALSE) model->IDSATRGAT=1E-18;
    if(model->CSRHBOT_Given == FALSE) model->CSRHBOT=1E2;
    if(model->CSRHSTI_Given == FALSE) model->CSRHSTI=1E-4;
    if(model->CSRHGAT_Given == FALSE) model->CSRHGAT=1E-4;
    if(model->XJUNSTI_Given == FALSE) model->XJUNSTI=100E-9;
    if(model->XJUNGAT_Given == FALSE) model->XJUNGAT=100E-9;
    if(model->CTATBOT_Given == FALSE) model->CTATBOT=1E2;
    if(model->CTATSTI_Given == FALSE) model->CTATSTI=1E-4;
    if(model->CTATGAT_Given == FALSE) model->CTATGAT=1E-4;
    if(model->MEFFTATBOT_Given == FALSE) model->MEFFTATBOT=0.25;
    if(model->MEFFTATSTI_Given == FALSE) model->MEFFTATSTI=0.25;
    if(model->MEFFTATGAT_Given == FALSE) model->MEFFTATGAT=0.25;
    if(model->CBBTBOT_Given == FALSE) model->CBBTBOT=1E-12;
    if(model->CBBTSTI_Given == FALSE) model->CBBTSTI=1E-18;
    if(model->CBBTGAT_Given == FALSE) model->CBBTGAT=1E-18;
    if(model->FBBTRBOT_Given == FALSE) model->FBBTRBOT=1E9;
    if(model->FBBTRSTI_Given == FALSE) model->FBBTRSTI=1E9;
    if(model->FBBTRGAT_Given == FALSE) model->FBBTRGAT=1E9;
    if(model->STFBBTBOT_Given == FALSE) model->STFBBTBOT=(-1E-3);
    if(model->STFBBTSTI_Given == FALSE) model->STFBBTSTI=(-1E-3);
    if(model->STFBBTGAT_Given == FALSE) model->STFBBTGAT=(-1E-3);
    if(model->VBRBOT_Given == FALSE) model->VBRBOT=10;
    if(model->VBRSTI_Given == FALSE) model->VBRSTI=10;
    if(model->VBRGAT_Given == FALSE) model->VBRGAT=10;
    if(model->PBRBOT_Given == FALSE) model->PBRBOT=4;
    if(model->PBRSTI_Given == FALSE) model->PBRSTI=4;
    if(model->PBRGAT_Given == FALSE) model->PBRGAT=4;
    for ( here = model->psp102instances ;here != NULL ; here = here->psp102nextInstance )
    {
    if(here->L_Given == FALSE) here->L=10e-6;
    if(here->W_Given == FALSE) here->W=10e-6;
    if(here->SA_Given == FALSE) here->SA=0.0;
    if(here->SB_Given == FALSE) here->SB=0.0;
    if(here->ABSOURCE_Given == FALSE) here->ABSOURCE=1E-12;
    if(here->LSSOURCE_Given == FALSE) here->LSSOURCE=1E-6;
    if(here->LGSOURCE_Given == FALSE) here->LGSOURCE=1E-6;
    if(here->ABDRAIN_Given == FALSE) here->ABDRAIN=1E-12;
    if(here->LSDRAIN_Given == FALSE) here->LSDRAIN=1E-6;
    if(here->LGDRAIN_Given == FALSE) here->LGDRAIN=1E-6;
    if(here->AS_Given == FALSE) here->AS=1E-12;
    if(here->PS_Given == FALSE) here->PS=1E-6;
    if(here->AD_Given == FALSE) here->AD=1E-12;
    if(here->PD_Given == FALSE) here->PD=1E-6;
    if(here->MULT_Given == FALSE) here->MULT=1.0;
      /* Internal Nodes */
      {
        here->NOINode = -1;
        here->NOI2Node = -1;
      }
      /* set states */
      here->state_NOI_GND = *states; *states += 2;
      here->state_B_D = *states; *states += 2;
      here->state_B_S = *states; *states += 2;
      here->state_G_B = *states; *states += 2;
      here->state_G_D = *states; *states += 2;
      here->state_G_S = *states; *states += 2;
      here->state_D_S = *states; *states += 2;
      /* set Sparse Matrix Pointers */
      here->PTR_J_D_S_required=0;
      here->PTR_J_D_B_required=0;
      here->PTR_J_B_S_required=0;
      here->PTR_J_B_B_required=0;
      here->PTR_J_D_D_required=0;
      here->PTR_J_B_D_required=0;
      here->PTR_J_D_G_required=0;
      here->PTR_J_B_G_required=0;
      here->PTR_J_S_S_required=0;
      here->PTR_J_S_B_required=0;
      here->PTR_J_S_D_required=0;
      here->PTR_J_S_G_required=0;
      here->PTR_J_G_S_required=0;
      here->PTR_J_G_B_required=0;
      here->PTR_J_G_D_required=0;
      here->PTR_J_G_G_required=0;
      here->PTR_J_NOI2_NOI2_required=0;
      here->PTR_J_NOI2_S_required=0;
      here->PTR_J_NOI2_B_required=0;
      here->PTR_J_NOI2_G_required=0;
      here->PTR_J_NOI2_D_required=0;
      here->PTR_J_NOI_S_required=0;
      here->PTR_J_NOI_B_required=0;
      here->PTR_J_NOI_G_required=0;
      here->PTR_J_NOI_D_required=0;
      here->PTR_J_NOI_NOI2_required=0;
      here->PTR_J_NOI_NOI_required=0;
      here->PTR_J_D_NOI2_required=0;
      here->PTR_J_S_NOI2_required=0;
      here->PTR_J_G_NOI_required=0;
      here->PTR_J_S_NOI_required=0;
      here->PTR_J_D_NOI_required=0;

      psp102guesstopology(matrix,ckt,model,here);

      /* Internal Nodes */
      {
        int error;
        CKTnode *tmp;
        if(here->NOINode == -1)
        {
          error=CKTmkVolt(ckt,&tmp,here->psp102name,"no info");
          if(error) return(error);
          here->NOINode = tmp->number;
        }
        if(here->NOI2Node == -1)
        {
          error=CKTmkVolt(ckt,&tmp,here->psp102name,"no info");
          if(error) return(error);
          here->NOI2Node = tmp->number;
        }
      }
        if(here->PTR_J_D_S_required==1)
        {
          here->PTR_J_D_S=SMPmakeElt(matrix,here->DNode,here->SNode);
        }
        if(here->PTR_J_D_B_required==1)
        {
          here->PTR_J_D_B=SMPmakeElt(matrix,here->DNode,here->BNode);
        }
        if(here->PTR_J_B_S_required==1)
        {
          here->PTR_J_B_S=SMPmakeElt(matrix,here->BNode,here->SNode);
        }
        if(here->PTR_J_B_B_required==1)
        {
          here->PTR_J_B_B=SMPmakeElt(matrix,here->BNode,here->BNode);
        }
        if(here->PTR_J_D_D_required==1)
        {
          here->PTR_J_D_D=SMPmakeElt(matrix,here->DNode,here->DNode);
        }
        if(here->PTR_J_B_D_required==1)
        {
          here->PTR_J_B_D=SMPmakeElt(matrix,here->BNode,here->DNode);
        }
        if(here->PTR_J_D_G_required==1)
        {
          here->PTR_J_D_G=SMPmakeElt(matrix,here->DNode,here->GNode);
        }
        if(here->PTR_J_B_G_required==1)
        {
          here->PTR_J_B_G=SMPmakeElt(matrix,here->BNode,here->GNode);
        }
        if(here->PTR_J_S_S_required==1)
        {
          here->PTR_J_S_S=SMPmakeElt(matrix,here->SNode,here->SNode);
        }
        if(here->PTR_J_S_B_required==1)
        {
          here->PTR_J_S_B=SMPmakeElt(matrix,here->SNode,here->BNode);
        }
        if(here->PTR_J_S_D_required==1)
        {
          here->PTR_J_S_D=SMPmakeElt(matrix,here->SNode,here->DNode);
        }
        if(here->PTR_J_S_G_required==1)
        {
          here->PTR_J_S_G=SMPmakeElt(matrix,here->SNode,here->GNode);
        }
        if(here->PTR_J_G_S_required==1)
        {
          here->PTR_J_G_S=SMPmakeElt(matrix,here->GNode,here->SNode);
        }
        if(here->PTR_J_G_B_required==1)
        {
          here->PTR_J_G_B=SMPmakeElt(matrix,here->GNode,here->BNode);
        }
        if(here->PTR_J_G_D_required==1)
        {
          here->PTR_J_G_D=SMPmakeElt(matrix,here->GNode,here->DNode);
        }
        if(here->PTR_J_G_G_required==1)
        {
          here->PTR_J_G_G=SMPmakeElt(matrix,here->GNode,here->GNode);
        }
        if(here->PTR_J_NOI2_NOI2_required==1)
        {
          here->PTR_J_NOI2_NOI2=SMPmakeElt(matrix,here->NOI2Node,here->NOI2Node);
        }
        if(here->PTR_J_NOI2_S_required==1)
        {
          here->PTR_J_NOI2_S=SMPmakeElt(matrix,here->NOI2Node,here->SNode);
        }
        if(here->PTR_J_NOI2_B_required==1)
        {
          here->PTR_J_NOI2_B=SMPmakeElt(matrix,here->NOI2Node,here->BNode);
        }
        if(here->PTR_J_NOI2_G_required==1)
        {
          here->PTR_J_NOI2_G=SMPmakeElt(matrix,here->NOI2Node,here->GNode);
        }
        if(here->PTR_J_NOI2_D_required==1)
        {
          here->PTR_J_NOI2_D=SMPmakeElt(matrix,here->NOI2Node,here->DNode);
        }
        if(here->PTR_J_NOI_S_required==1)
        {
          here->PTR_J_NOI_S=SMPmakeElt(matrix,here->NOINode,here->SNode);
        }
        if(here->PTR_J_NOI_B_required==1)
        {
          here->PTR_J_NOI_B=SMPmakeElt(matrix,here->NOINode,here->BNode);
        }
        if(here->PTR_J_NOI_G_required==1)
        {
          here->PTR_J_NOI_G=SMPmakeElt(matrix,here->NOINode,here->GNode);
        }
        if(here->PTR_J_NOI_D_required==1)
        {
          here->PTR_J_NOI_D=SMPmakeElt(matrix,here->NOINode,here->DNode);
        }
        if(here->PTR_J_NOI_NOI2_required==1)
        {
          here->PTR_J_NOI_NOI2=SMPmakeElt(matrix,here->NOINode,here->NOI2Node);
        }
        if(here->PTR_J_NOI_NOI_required==1)
        {
          here->PTR_J_NOI_NOI=SMPmakeElt(matrix,here->NOINode,here->NOINode);
        }
        if(here->PTR_J_D_NOI2_required==1)
        {
          here->PTR_J_D_NOI2=SMPmakeElt(matrix,here->DNode,here->NOI2Node);
        }
        if(here->PTR_J_S_NOI2_required==1)
        {
          here->PTR_J_S_NOI2=SMPmakeElt(matrix,here->SNode,here->NOI2Node);
        }
        if(here->PTR_J_G_NOI_required==1)
        {
          here->PTR_J_G_NOI=SMPmakeElt(matrix,here->GNode,here->NOINode);
        }
        if(here->PTR_J_S_NOI_required==1)
        {
          here->PTR_J_S_NOI=SMPmakeElt(matrix,here->SNode,here->NOINode);
        }
        if(here->PTR_J_D_NOI_required==1)
        {
          here->PTR_J_D_NOI=SMPmakeElt(matrix,here->DNode,here->NOINode);
        }

    }
  }
  return(OK);
}
Ejemplo n.º 9
0
int
HFETAsetup(SMPmatrix *matrix, GENmodel *inModel, CKTcircuit *ckt, int *states)
        /* load the diode structure with those pointers needed later 
         * for fast matrix loading 
         */
{
    HFETAmodel *model = (HFETAmodel*)inModel;
    HFETAinstance *here;
    int error;
    CKTnode *tmp;


    /*  loop through all the diode models */
    for( ; model != NULL; model = model->HFETAnextModel ) {
        if( (model->HFETAtype != NHFET) && (model->HFETAtype != PHFET) ) {
            model->HFETAtype = NHFET;
        }
        if(!model->HFETAthresholdGiven) {
            if(model->HFETAtype == NHFET)
              model->HFETAthreshold = 0.15;
            else
              model->HFETAthreshold = -0.15;
        }
        if(!model->HFETAdiGiven) {
            model->HFETAdi = 0.04e-6;
        }
        if(!model->HFETAlambdaGiven) {
            model->HFETAlambda = 0.15;
        }
        if(!model->HFETAetaGiven) {
            if(model->HFETAtype == NHFET)
              model->HFETAeta = 1.28;
            else
              model->HFETAeta = 1.4;
        }
        if(!model->HFETAmGiven) {
            model->HFETAm = 3.0;
        }
        if(!model->HFETAmcGiven) {
            model->HFETAmc = 3.0;
        }
        if(!model->HFETAgammaGiven) {
            model->HFETAgamma = 3.0;
        }
        if(!model->HFETAsigma0Given) {
            model->HFETAsigma0 = 0.057;
        }
        if(!model->HFETAvsigmatGiven) {
            model->HFETAvsigmat = 0.3;
        }
        if(!model->HFETAvsigmaGiven) {
            model->HFETAvsigma = 0.1;
        }
        if(!model->HFETAmuGiven) {
            if(model->HFETAtype == NHFET)
              model->HFETAmu = 0.4;
            else
              model->HFETAmu = 0.03;
        }
        if(!model->HFETAdeltaGiven) {
            model->HFETAdelta = 3.0;
        }
        if(!model->HFETAvsGiven) {
            if(model->HFETAtype == NHFET)
              model->HFETAvs = 1.5e5;
            else
              model->HFETAvs = 0.8e5;
        }
        if(!model->HFETAnmaxGiven) {
            model->HFETAnmax = 2e16;
        }
        if(!model->HFETAdeltadGiven) {
            model->HFETAdeltad = 4.5e-9;
        }
        if(!model->HFETAjs1dGiven) {
            model->HFETAjs1d = 1.0;
        }
        if(!model->HFETAjs2dGiven) {
            model->HFETAjs2d = 1.15e6;
        }
        if(!model->HFETAjs1sGiven) {
            model->HFETAjs1s = 1.0;
        }
        if(!model->HFETAjs2sGiven) {
            model->HFETAjs2s = 1.15e6;
        }
        if(!model->HFETAm1dGiven) {
            model->HFETAm1d = 1.32;
        }
        if(!model->HFETAm2dGiven) {
            model->HFETAm2d = 6.9;
        }
        if(!model->HFETAm1sGiven) {
            model->HFETAm1s = 1.32;
        }
        if(!model->HFETAm2sGiven) {
            model->HFETAm2s = 6.9;
        }
        if(!model->HFETArdGiven) {
            model->HFETArd = 0;
        }
        if(!model->HFETArsGiven) {
            model->HFETArs = 0;
        }
        if(!model->HFETArdiGiven) {
            model->HFETArdi = 0;
        }
        if(!model->HFETArsiGiven) {
            model->HFETArsi = 0;
        }
        if(!model->HFETArgsGiven) {
            model->HFETArgs = 90;
        }
        if(!model->HFETArgdGiven) {
            model->HFETArgd = 90;
        }
        if(!model->HFETAriGiven) {
            model->HFETAri = 0;
        }
        if(!model->HFETArfGiven) {
            model->HFETArf = 0;
        }
        if(!model->HFETAepsiGiven) {
            model->HFETAepsi = 12.244*8.85418e-12;
        }
        if(!model->HFETAa1Given) {
            model->HFETAa1 = 0;
        }
        if(!model->HFETAa2Given) {
            model->HFETAa2 = 0;
        }
        if(!model->HFETAmv1Given) {
            model->HFETAmv1 = 3;
        }
        if(!model->HFETApGiven) {
            model->HFETAp = 1;
        }
        if(!model->HFETAkappaGiven) {
            model->HFETAkappa = 0;
        }
        if(!model->HFETAdelfGiven) {
            model->HFETAdelf = 0;
        }
        if(!model->HFETAfgdsGiven) {
            model->HFETAfgds = 0;
        }
        if(!model->HFETAtfGiven) {
            model->HFETAtf = ckt->CKTtemp;
        }
        if(!model->HFETAcdsGiven) {
            model->HFETAcds = 0;
        }
        if(!model->HFETAphibGiven) {
            model->HFETAphib = 0.5*CHARGE;
        }
        if(!model->HFETAtalphaGiven) {
            model->HFETAtalpha = 1200;
        }
        if(!model->HFETAmt1Given) {
            model->HFETAmt1 = 3.5;
        }
        if(!model->HFETAmt2Given) {
            model->HFETAmt2 = 9.9;
        }
        if(!model->HFETAck1Given) {
            model->HFETAck1 = 1;
        }
        if(!model->HFETAck2Given) {
            model->HFETAck2 = 0;
        }
        if(!model->HFETAcm1Given) {
            model->HFETAcm1 = 3;
        }
        if(!model->HFETAcm2Given) {
            model->HFETAcm2 = 0;
        }
        if(!model->HFETAcm3Given) {
            model->HFETAcm3 = 0.17;
        }
        if(!model->HFETAastarGiven) {
            model->HFETAastar = 4.0e4;
        }
        if(!model->HFETAeta1Given) {
            model->HFETAeta1 = 2;
        }
        if(!model->HFETAd1Given) {
            model->HFETAd1 = 0.03e-6;
        }
        if(!model->HFETAeta2Given) {
            model->HFETAeta2 = 2;
        }
        if(!model->HFETAd2Given) {
            model->HFETAd2 = 0.2e-6;
        }
        if(!model->HFETAvt2Given) {
            /* initialized in HFETAtemp */
            model->HFETAvt2 = 0;
        }
        
        if(!model->HFETAggrGiven) {
            model->HFETAggr = 40;
        }
        if(!model->HFETAdelGiven) {
            model->HFETAdel = 0.04;
        }             
       if(!model->HFETAklambdaGiven)
         KLAMBDA = 0;
       if(!model->HFETAkmuGiven)
         KMU = 0;
       if(!model->HFETAkvtoGiven)
         KVTO = 0;
        
        /* loop through all the instances of the model */
        for (here = model->HFETAinstances; here != NULL ;
                here=here->HFETAnextInstance) {
           
            if(!here->HFETAlengthGiven) {
                here->HFETAlength = 1e-6;
            }
            if(!here->HFETAwidthGiven) {
                here->HFETAwidth = 20e-6;
            }
            if(!here->HFETAmGiven) {
                here->HFETAm = 1.0;
            }
 
            here->HFETAstate = *states;
            /* *states += 24; */
            *states += HFETAnumStates;

            if(model->HFETArs != 0) {
                if(here->HFETAsourcePrimeNode == 0) {
                error = CKTmkVolt(ckt,&tmp,here->HFETAname,"source");
                if(error) return(error);
                here->HFETAsourcePrimeNode = tmp->number;

/* XXX: Applied AlansFixes  */               
                 if (ckt->CKTcopyNodesets) {
		    CKTnode *tmpNode;
		    IFuid tmpName;

                  if (CKTinst2Node(ckt,here,3,&tmpNode,&tmpName)==OK) {
                     if (tmpNode->nsGiven) {
                       tmp->nodeset=tmpNode->nodeset; 
                       tmp->nsGiven=tmpNode->nsGiven; 
                     }
                  }
                }
                }

            } else {
                here->HFETAsourcePrimeNode = here->HFETAsourceNode;
            }
            
            if(model->HFETArd != 0) {
                if(here->HFETAdrainPrimeNode == 0) {
                error = CKTmkVolt(ckt,&tmp,here->HFETAname,"drain");
                if(error) return(error);
                here->HFETAdrainPrimeNode = tmp->number;
       
                if (ckt->CKTcopyNodesets) {
		    CKTnode *tmpNode;
		    IFuid tmpName;

                  if (CKTinst2Node(ckt,here,1,&tmpNode,&tmpName)==OK) {
                     if (tmpNode->nsGiven) {
                       tmp->nodeset=tmpNode->nodeset; 
                       tmp->nsGiven=tmpNode->nsGiven; 
                     }
                  }
                }
                }
                
            } else {
                here->HFETAdrainPrimeNode = here->HFETAdrainNode;
            }
            
            if(model->HFETArg != 0) {
                if(here->HFETAgatePrimeNode == 0) {
                error = CKTmkVolt(ckt,&tmp,here->HFETAname,"gate");
                if(error) return(error);
                here->HFETAgatePrimeNode = tmp->number;
         
                if (ckt->CKTcopyNodesets) {
		    CKTnode *tmpNode;
		    IFuid tmpName;

                  if (CKTinst2Node(ckt,here,2,&tmpNode,&tmpName)==OK) {
                     if (tmpNode->nsGiven) {
                       tmp->nodeset=tmpNode->nodeset; 
                       tmp->nsGiven=tmpNode->nsGiven; 
                     }
                  }
                }     
                }     
                
            } else {
                here->HFETAgatePrimeNode = here->HFETAgateNode;
            }
            if(model->HFETArf != 0) {
                if(here->HFETAdrainPrmPrmNode == 0) {
                error = CKTmkVolt(ckt,&tmp,here->HFETAname,"gd");
                if(error) return(error);
                here->HFETAdrainPrmPrmNode = tmp->number;
         
                if (ckt->CKTcopyNodesets) {
		    CKTnode *tmpNode;
		    IFuid tmpName;

                  if (CKTinst2Node(ckt,here,1,&tmpNode,&tmpName)==OK) {
                     if (tmpNode->nsGiven) {
                       tmp->nodeset=tmpNode->nodeset; 
                       tmp->nsGiven=tmpNode->nsGiven; 
                     }
                  }
                }
                }

            } else {
                here->HFETAdrainPrmPrmNode = here->HFETAdrainPrimeNode;
            }
            
            if(model->HFETAri != 0) {
                if(here->HFETAsourcePrmPrmNode == 0) {
                error = CKTmkVolt(ckt,&tmp,here->HFETAname,"gs");
                if(error) return(error);
                here->HFETAsourcePrmPrmNode = tmp->number;
      
                if (ckt->CKTcopyNodesets) {
		    CKTnode *tmpNode;
		    IFuid tmpName;

                  if (CKTinst2Node(ckt,here,3,&tmpNode,&tmpName)==OK) {
                     if (tmpNode->nsGiven) {
                       tmp->nodeset=tmpNode->nodeset; 
                       tmp->nsGiven=tmpNode->nsGiven; 
                     }
                  }
                }
                }
                
            } else {
                here->HFETAsourcePrmPrmNode = here->HFETAsourcePrimeNode;
            }

/* macro to make elements with built in test for out of memory */
#define TSTALLOC(ptr,first,second) \
do { if((here->ptr = SMPmakeElt(matrix, here->first, here->second)) == NULL){\
    return(E_NOMEM);\
} } while(0)

  TSTALLOC(HFETAdrainDrainPrimePtr,HFETAdrainNode,HFETAdrainPrimeNode);
  TSTALLOC(HFETAgatePrimeDrainPrimePtr,HFETAgatePrimeNode,HFETAdrainPrimeNode);
  TSTALLOC(HFETAgatePrimeSourcePrimePtr,HFETAgatePrimeNode,HFETAsourcePrimeNode);
  TSTALLOC(HFETAsourceSourcePrimePtr,HFETAsourceNode,HFETAsourcePrimeNode);
  TSTALLOC(HFETAdrainPrimeDrainPtr,HFETAdrainPrimeNode,HFETAdrainNode);
  TSTALLOC(HFETAdrainPrimeGatePrimePtr,HFETAdrainPrimeNode,HFETAgatePrimeNode);
  TSTALLOC(HFETAdrainPrimeSourcePrimePtr,HFETAdrainPrimeNode,HFETAsourcePrimeNode);
  TSTALLOC(HFETAsourcePrimeGatePrimePtr,HFETAsourcePrimeNode,HFETAgatePrimeNode);
  TSTALLOC(HFETAsourcePrimeSourcePtr,HFETAsourcePrimeNode,HFETAsourceNode);
  TSTALLOC(HFETAsourcePrimeDrainPrimePtr,HFETAsourcePrimeNode,HFETAdrainPrimeNode);
  TSTALLOC(HFETAdrainDrainPtr,HFETAdrainNode,HFETAdrainNode);
  TSTALLOC(HFETAgatePrimeGatePrimePtr,HFETAgatePrimeNode,HFETAgatePrimeNode);
  TSTALLOC(HFETAsourceSourcePtr,HFETAsourceNode,HFETAsourceNode);
  TSTALLOC(HFETAdrainPrimeDrainPrimePtr,HFETAdrainPrimeNode,HFETAdrainPrimeNode);
  TSTALLOC(HFETAsourcePrimeSourcePrimePtr,HFETAsourcePrimeNode,HFETAsourcePrimeNode);
  TSTALLOC(HFETAdrainPrimeDrainPrmPrmPtr,HFETAdrainPrimeNode,HFETAdrainPrmPrmNode);
  TSTALLOC(HFETAdrainPrmPrmDrainPrimePtr,HFETAdrainPrmPrmNode,HFETAdrainPrimeNode);
  TSTALLOC(HFETAdrainPrmPrmGatePrimePtr,HFETAdrainPrmPrmNode,HFETAgatePrimeNode);
  TSTALLOC(HFETAgatePrimeDrainPrmPrmPtr,HFETAgatePrimeNode,HFETAdrainPrmPrmNode);
  TSTALLOC(HFETAdrainPrmPrmDrainPrmPrmPtr,HFETAdrainPrmPrmNode,HFETAdrainPrmPrmNode);
  TSTALLOC(HFETAsourcePrimeSourcePrmPrmPtr,HFETAsourcePrimeNode,HFETAsourcePrmPrmNode);
  TSTALLOC(HFETAsourcePrmPrmSourcePrimePtr,HFETAsourcePrmPrmNode,HFETAsourcePrimeNode);
  TSTALLOC(HFETAsourcePrmPrmGatePrimePtr,HFETAsourcePrmPrmNode,HFETAgatePrimeNode);
  TSTALLOC(HFETAgatePrimeSourcePrmPrmPtr,HFETAgatePrimeNode,HFETAsourcePrmPrmNode);
  TSTALLOC(HFETAsourcePrmPrmSourcePrmPrmPtr,HFETAsourcePrmPrmNode,HFETAsourcePrmPrmNode);
  TSTALLOC(HFETAgateGatePtr,HFETAgateNode,HFETAgateNode);
  TSTALLOC(HFETAgateGatePrimePtr,HFETAgateNode,HFETAgatePrimeNode);
  TSTALLOC(HFETAgatePrimeGatePtr,HFETAgatePrimeNode,HFETAgateNode);
 }
 }
 return(OK);
}
Ejemplo n.º 10
0
int
MESAsetup(SMPmatrix *matrix, GENmodel *inModel, CKTcircuit *ckt, int *states)

        /* load the diode structure with those pointers needed later 
         * for fast matrix loading 
         */
{
    MESAmodel *model = (MESAmodel*)inModel;
    MESAinstance *here;
    int error;
    CKTnode *tmp;


    /*  loop through all the diode models */
    for( ; model != NULL; model = model->MESAnextModel ) {
        if( (model->MESAtype != NMF) ) {
            fprintf(stderr, "Only nmf model type supported, set to nmf\n");
            model->MESAtype = NMF;
        }
        if(!model->MESAthresholdGiven) {
            model->MESAthreshold = -1.26;
        }
        if(!model->MESAdGiven) {
            model->MESAd = 0.12e-6;
        }
        if(!model->MESAduGiven) {
            model->MESAdu = 0.035e-6;
        }
        if(!model->MESAlambdaGiven) {
            model->MESAlambda = 0.045;
        }
        if(!model->MESAvsGiven) {
            model->MESAvs = 1.5e5;
        }
        if(!model->MESAbetaGiven) {
            model->MESAbeta = 0.0085;
        }
        if(!model->MESAetaGiven) {
            model->MESAeta = 1.73;
        }
        if(!model->MESAmGiven) {
            model->MESAm = 2.5;
        }
        if(!model->MESAmcGiven) {
            model->MESAmc = 3.0;
        }
        if(!model->MESAalphaGiven) {
            model->MESAalpha = 0.0;
        }
        if(!model->MESAsigma0Given) {
            model->MESAsigma0 = 0.081;
        }
        if(!model->MESAvsigmatGiven) {
            model->MESAvsigmat = 1.01;
        }
        if(!model->MESAvsigmaGiven) {
            model->MESAvsigma = 0.1;
        }
        if(!model->MESAmuGiven) {
            model->MESAmu = 0.23;
        }
        if(!model->MESAthetaGiven) {
            model->MESAtheta = 0;
        }
        if(!model->MESAmu1Given) {
            model->MESAmu1 = 0;
        }
        if(!model->MESAmu2Given) {
            model->MESAmu2 = 0;
        }
        if(!model->MESAndGiven) {
            model->MESAnd = 2.0e23;
        }
        if(!model->MESAnduGiven) {
            model->MESAndu = 1e22;
        }
        if(!model->MESAndeltaGiven) {
            model->MESAndelta = 6e24;
        }
        if(!model->MESAthGiven) {
            model->MESAth = 0.01e-6;
        }
        if(!model->MESAdeltaGiven) {
            model->MESAdelta = 5.0;
        }
        if(!model->MESAtcGiven) {
            model->MESAtc = 0.0;
        }
        if(!model->MESAdrainResistGiven) {
            model->MESAdrainResist = 0;
        }
        if(!model->MESAsourceResistGiven) {
            model->MESAsourceResist = 0;
        }
        if(!model->MESAgateResistGiven) {
            model->MESAgateResist = 0;
        }
        if(!model->MESAriGiven) {
            model->MESAri = 0;
        }
        if(!model->MESArfGiven) {
            model->MESArf = 0;
        }
        if(!model->MESArdiGiven) {
            model->MESArdi = 0;
        }
        if(!model->MESArsiGiven) {
            model->MESArsi = 0;
        }
        if(!model->MESAphibGiven) {
            model->MESAphib = 0.5*CHARGE;
        }
        if(!model->MESAphib1Given) {
            model->MESAphib1 = 0;
        }
        if(!model->MESAastarGiven) {
            model->MESAastar = 4.0e4;
        }
        if(!model->MESAggrGiven) {
            model->MESAggr = 40;
        }
        if(!model->MESAdelGiven) {
            model->MESAdel = 0.04;
        }
        if(!model->MESAxchiGiven) {
            model->MESAxchi = 0.033;
        }
        if(!model->MESAnGiven) {
            model->MESAn = 1;
        }
        if(!model->MESAtvtoGiven) {
            model->MESAtvto = 0;
        }
        if(!model->MESAtlambdaGiven) {
            model->MESAtlambda = DBL_MAX;
        }
        if(!model->MESAteta0Given) {
            model->MESAteta0 = DBL_MAX;
        }
        if(!model->MESAteta1Given) {
            model->MESAteta1 = 0;
        }
        if(!model->MESAtmuGiven) {
            model->MESAtmu = 300.15;
        }
        if(!model->MESAxtm0Given) {
            model->MESAxtm0 = 0;
        }
        if(!model->MESAxtm1Given) {
            model->MESAxtm1 = 0;
        }
        if(!model->MESAxtm2Given) {
            model->MESAxtm2 = 0;
        }
        if(!model->MESAksGiven) {
            model->MESAks = 0;
        }
        if(!model->MESAvsgGiven) {
            model->MESAvsg = 0;
        }
        if(!model->MESAtfGiven) {
            model->MESAtf = ckt->CKTtemp;
        }
        if(!model->MESAfloGiven) {
            model->MESAflo = 0;
        }
        if(!model->MESAdelfoGiven) {
            model->MESAdelfo = 0;
        }
        if(!model->MESAagGiven) {
            model->MESAag = 0;
        }
        if(!model->MESAtc1Given) {
            model->MESAtc1 = 0;
        }
        if(!model->MESAtc2Given) {
            model->MESAtc2 = 0;
        }
        if(!model->MESAzetaGiven) {
            model->MESAzeta = 1;
        }
        if(!model->MESAlevelGiven) {
            model->MESAlevel = 2;
        }
        if(!model->MESAnmaxGiven) {
            model->MESAnmax = 2e16;
        }
        if(!model->MESAgammaGiven) {
            model->MESAgamma = 3.0;
        }
        if(!model->MESAepsiGiven) {
            model->MESAepsi = 12.244*8.85418e-12;
        }
        if(!model->MESAcasGiven) {
            model->MESAcas = 1;
        }                      
        if(!model->MESAcbsGiven) {
            model->MESAcbs = 1;
        }

        if(model->MESAdrainResist != 0) {
            model->MESAdrainConduct = 1./model->MESAdrainResist;
        } else {
            model->MESAdrainConduct = DBL_MAX;
        }
        if(model->MESAsourceResist != 0) {
            model->MESAsourceConduct = 1./model->MESAsourceResist;
        } else {
            model->MESAsourceConduct = DBL_MAX;
        }

        model->MESAvcrit = 0.; /* until model has changed */

        /* loop through all the instances of the model */
        for (here = model->MESAinstances; here != NULL ;
                here=here->MESAnextInstance) {
         
            if(!here->MESAlengthGiven) {
                here->MESAlength = 1e-6;
            }
            if(!here->MESAwidthGiven) {
                here->MESAwidth = 20e-6;
            }
            if(!here->MESAmGiven) {
                here->MESAm = 1.0;
            }
            if(!here->MESAdtempGiven) {
                here->MESAdtemp = 0.0;
            }
            if(!here->MESAtdGiven) {
                here->MESAtd = ckt->CKTtemp + here->MESAdtemp;
            }
            if(!here->MESAtsGiven) {
                here->MESAts = ckt->CKTtemp + here->MESAdtemp;
            }


            here->MESAstate = *states;
            *states += 20;

            if(model->MESAsourceResist != 0) {
                if(here->MESAsourcePrimeNode == 0) {
                error = CKTmkVolt(ckt,&tmp,here->MESAname,"source");
                if(error) return(error);
                here->MESAsourcePrimeNode = tmp->number;
                
                if (ckt->CKTcopyNodesets) {
		    CKTnode *tmpNode;
		    IFuid tmpName;

                  if (CKTinst2Node(ckt,here,3,&tmpNode,&tmpName)==OK) {
                     if (tmpNode->nsGiven) {
                       tmp->nodeset=tmpNode->nodeset; 
                       tmp->nsGiven=tmpNode->nsGiven; 
                     }
                  }
                }
                }
                
            } else {
                here->MESAsourcePrimeNode = here->MESAsourceNode;
            }
            
            if(model->MESAdrainResist != 0) {
                if(here->MESAdrainPrimeNode == 0) {
                error = CKTmkVolt(ckt,&tmp,here->MESAname,"drain");
                if(error) return(error);
                here->MESAdrainPrimeNode = tmp->number;
                
                if (ckt->CKTcopyNodesets) {
		    CKTnode *tmpNode;
		    IFuid tmpName;

                  if (CKTinst2Node(ckt,here,1,&tmpNode,&tmpName)==OK) {
                     if (tmpNode->nsGiven) {
                       tmp->nodeset=tmpNode->nodeset; 
                       tmp->nsGiven=tmpNode->nsGiven; 
                     }
                  }
                }
                }
                
            } else {
                here->MESAdrainPrimeNode = here->MESAdrainNode;
            }
            if(model->MESAgateResist != 0) {
                if(here->MESAgatePrimeNode == 0) {
                error = CKTmkVolt(ckt,&tmp,here->MESAname,"gate");
                if(error) return(error);
                here->MESAgatePrimeNode = tmp->number;
                
                if (ckt->CKTcopyNodesets) {
		    CKTnode *tmpNode;
		    IFuid tmpName;

                  if (CKTinst2Node(ckt,here,2,&tmpNode,&tmpName)==OK) {
                     if (tmpNode->nsGiven) {
                       tmp->nodeset=tmpNode->nodeset; 
                       tmp->nsGiven=tmpNode->nsGiven; 
                     }
                  }
                }
                }
                
                
            } else {
                here->MESAgatePrimeNode = here->MESAgateNode;
            }
            
            
            if(model->MESAri != 0) {
                if(here->MESAsourcePrmPrmNode == 0) {
                error = CKTmkVolt(ckt,&tmp,here->MESAname,"gs");
                if(error) return(error);
                here->MESAsourcePrmPrmNode = tmp->number;
                
                if (ckt->CKTcopyNodesets) {
		    CKTnode *tmpNode;
		    IFuid tmpName;

                  if (CKTinst2Node(ckt,here,3,&tmpNode,&tmpName)==OK) {
                     if (tmpNode->nsGiven) {
                       tmp->nodeset=tmpNode->nodeset; 
                       tmp->nsGiven=tmpNode->nsGiven; 
                     }
                  }
                }
                }
                
            } else {
                here->MESAsourcePrmPrmNode = here->MESAsourcePrimeNode;
            }
            if(model->MESArf != 0) {
                if(here->MESAdrainPrmPrmNode == 0) {
                error = CKTmkVolt(ckt,&tmp,here->MESAname,"gd");
                if(error) return(error);
                here->MESAdrainPrmPrmNode = tmp->number;
                
                if (ckt->CKTcopyNodesets) {
		    CKTnode *tmpNode;
		    IFuid tmpName;

                  if (CKTinst2Node(ckt,here,1,&tmpNode,&tmpName)==OK) {
                     if (tmpNode->nsGiven) {
                       tmp->nodeset=tmpNode->nodeset; 
                       tmp->nsGiven=tmpNode->nsGiven; 
                     }
                  }
                }
                }
                
            } else {
                here->MESAdrainPrmPrmNode = here->MESAdrainPrimeNode;
            }

#define TSTALLOC(ptr,first,second) \
do { if((here->ptr = SMPmakeElt(matrix, here->first, here->second)) == NULL){\
    return(E_NOMEM);\
} } while(0)

  TSTALLOC(MESAdrainDrainPtr,MESAdrainNode,MESAdrainNode);
  TSTALLOC(MESAdrainPrimeDrainPrimePtr,MESAdrainPrimeNode,MESAdrainPrimeNode);
  TSTALLOC(MESAdrainPrmPrmDrainPrmPrmPtr,MESAdrainPrmPrmNode,MESAdrainPrmPrmNode);
  TSTALLOC(MESAgateGatePtr,MESAgateNode,MESAgateNode);
  TSTALLOC(MESAgatePrimeGatePrimePtr,MESAgatePrimeNode,MESAgatePrimeNode);
  TSTALLOC(MESAsourceSourcePtr,MESAsourceNode,MESAsourceNode);
  TSTALLOC(MESAsourcePrimeSourcePrimePtr,MESAsourcePrimeNode,MESAsourcePrimeNode);
  TSTALLOC(MESAsourcePrmPrmSourcePrmPrmPtr,MESAsourcePrmPrmNode,MESAsourcePrmPrmNode);
  TSTALLOC(MESAdrainDrainPrimePtr,MESAdrainNode,MESAdrainPrimeNode);
  TSTALLOC(MESAdrainPrimeDrainPtr,MESAdrainPrimeNode,MESAdrainNode);
  TSTALLOC(MESAgatePrimeDrainPrimePtr,MESAgatePrimeNode,MESAdrainPrimeNode);
  TSTALLOC(MESAdrainPrimeGatePrimePtr,MESAdrainPrimeNode,MESAgatePrimeNode);
  TSTALLOC(MESAgatePrimeSourcePrimePtr,MESAgatePrimeNode,MESAsourcePrimeNode);
  TSTALLOC(MESAsourcePrimeGatePrimePtr,MESAsourcePrimeNode,MESAgatePrimeNode)  ;
  TSTALLOC(MESAsourceSourcePrimePtr,MESAsourceNode,MESAsourcePrimeNode);
  TSTALLOC(MESAsourcePrimeSourcePtr,MESAsourcePrimeNode,MESAsourceNode);
  TSTALLOC(MESAdrainPrimeSourcePrimePtr,MESAdrainPrimeNode,MESAsourcePrimeNode);
  TSTALLOC(MESAsourcePrimeDrainPrimePtr,MESAsourcePrimeNode,MESAdrainPrimeNode);
  TSTALLOC(MESAgatePrimeGatePtr,MESAgatePrimeNode,MESAgateNode);
  TSTALLOC(MESAgateGatePrimePtr,MESAgateNode,MESAgatePrimeNode);
  TSTALLOC(MESAsourcePrmPrmSourcePrimePtr,MESAsourcePrmPrmNode,MESAsourcePrimeNode);
  TSTALLOC(MESAsourcePrimeSourcePrmPrmPtr,MESAsourcePrimeNode,MESAsourcePrmPrmNode);
  TSTALLOC(MESAsourcePrmPrmGatePrimePtr,MESAsourcePrmPrmNode,MESAgatePrimeNode);
  TSTALLOC(MESAgatePrimeSourcePrmPrmPtr,MESAgatePrimeNode,MESAsourcePrmPrmNode);
  TSTALLOC(MESAdrainPrmPrmDrainPrimePtr,MESAdrainPrmPrmNode,MESAdrainPrimeNode);
  TSTALLOC(MESAdrainPrimeDrainPrmPrmPtr,MESAdrainPrimeNode,MESAdrainPrmPrmNode);
  TSTALLOC(MESAdrainPrmPrmGatePrimePtr,MESAdrainPrmPrmNode,MESAgatePrimeNode);
  TSTALLOC(MESAgatePrimeDrainPrmPrmPtr,MESAgatePrimeNode,MESAdrainPrmPrmNode);
  }
  }
    return(OK);
}
Ejemplo n.º 11
0
int
DIOsetup(SMPmatrix *matrix, GENmodel *inModel, CKTcircuit *ckt, int *states)
{
    DIOmodel *model = (DIOmodel*)inModel;
    DIOinstance *here;
    int error;
    CKTnode *tmp;

    /*  loop through all the diode models */
    for( ; model != NULL; model = model->DIOnextModel ) {

        if(!model->DIOlevelGiven) {
            model->DIOlevel = 1;
        }
        if(!model->DIOemissionCoeffGiven) {
            model->DIOemissionCoeff = 1;
        }
        if(!model->DIOsatCurGiven) {
            model->DIOsatCur = 1e-14;
        }
        if(!model->DIOsatSWCurGiven) {
            model->DIOsatSWCur = 0.0;
        }
        if(!model->DIOswEmissionCoeffGiven) {
            model->DIOswEmissionCoeff = 1;
        }
        if(!model->DIObreakdownCurrentGiven) {
            model->DIObreakdownCurrent = 1e-3;
        }
        if(!model->DIOjunctionPotGiven){
            model->DIOjunctionPot = 1;
        }
        if(!model->DIOgradingCoeffGiven) {
            model->DIOgradingCoeff = .5;
        }
        if(!model->DIOgradCoeffTemp1Given) {
            model->DIOgradCoeffTemp1 = 0.0;
        }
        if(!model->DIOgradCoeffTemp2Given) {
            model->DIOgradCoeffTemp2 = 0.0;
        }
        if(!model->DIOdepletionCapCoeffGiven) {
            model->DIOdepletionCapCoeff = .5;
        }
        if(!model->DIOdepletionSWcapCoeffGiven) {
            model->DIOdepletionSWcapCoeff = .5;
        }
        if(!model->DIOtransitTimeGiven) {
            model->DIOtransitTime = 0;
        }
        if(!model->DIOtranTimeTemp1Given) {
            model->DIOtranTimeTemp1 = 0.0;
        }
        if(!model->DIOtranTimeTemp2Given) {
            model->DIOtranTimeTemp2 = 0.0;
        }
        if(!model->DIOjunctionCapGiven) {
            model->DIOjunctionCap = 0;
        }
        if(!model->DIOjunctionSWCapGiven) {
            model->DIOjunctionSWCap = 0;
        }
        if(!model->DIOjunctionSWPotGiven){
            model->DIOjunctionSWPot = 1;
        }
        if(!model->DIOgradingSWCoeffGiven) {
            model->DIOgradingSWCoeff = .33;
        }
        if(!model->DIOforwardKneeCurrentGiven) {
            model->DIOforwardKneeCurrent = 0.0;
        }
        if(!model->DIOreverseKneeCurrentGiven) {
            model->DIOreverseKneeCurrent = 0.0;
        }
        if(!model->DIObrkdEmissionCoeffGiven) {
            model->DIObrkdEmissionCoeff = model->DIOemissionCoeff;
        }
        if(!model->DIOtlevGiven) {
            model->DIOtlev = 0;
        }
        if(!model->DIOtlevcGiven) {
            model->DIOtlevc = 0;
        }
        if(!model->DIOactivationEnergyGiven) {
            model->DIOactivationEnergy = 1.11;
        }
        if(!model->DIOsaturationCurrentExpGiven) {
            model->DIOsaturationCurrentExp = 3;
        }
        if(!model->DIOctaGiven) {
            model->DIOcta = 0.0;
        }
        if(!model->DIOctpGiven) {
            model->DIOctp = 0.0;
        }
        if(!model->DIOtpbGiven) {
            model->DIOtpb = 0.0;
        }
        if(!model->DIOtphpGiven) {
            model->DIOtphp = 0.0;
        }
        if(!model->DIOfNcoefGiven) {
            model->DIOfNcoef = 0.0;
        }
        if(!model->DIOfNexpGiven) {
            model->DIOfNexp = 1.0;
        }
        if(!model->DIOresistTemp1Given) {
            model->DIOresistTemp1 = 0.0;
        }
        if(!model->DIOresistTemp2Given) {
            model->DIOresistTemp2 = 0.0;
        }
        if(!model->DIOtcvGiven) {
            model->DIOtcv = 0.0;
        }
        if(!model->DIOareaGiven) {
            model->DIOarea = 1.0;
        }
        if(!model->DIOpjGiven) {
            model->DIOpj = 0.0;
        }
        if(!model->DIOtunSatCurGiven) {
            model->DIOtunSatCur = 0.0;
        }
        if(!model->DIOtunSatSWCurGiven) {
            model->DIOtunSatSWCur = 0.0;
        }
        if(!model->DIOtunEmissionCoeffGiven) {
            model->DIOtunEmissionCoeff = 30.0;
        }
        if(!model->DIOtunSaturationCurrentExpGiven) {
            model->DIOtunSaturationCurrentExp = 3.0;
        }
        if(!model->DIOtunEGcorrectionFactorGiven) {
            model->DIOtunEGcorrectionFactor = 1.0;
        }
        if(!model->DIOfv_maxGiven) {
            model->DIOfv_max = 1e99;
        }
        if(!model->DIObv_maxGiven) {
            model->DIObv_max = 1e99;
        }

        /* loop through all the instances of the model */
        for (here = model->DIOinstances; here != NULL ;
                here=here->DIOnextInstance) {

            if(!here->DIOareaGiven) {
                if((!here->DIOwGiven) && (!here->DIOlGiven))  {
                    here->DIOarea = model->DIOarea;
                } else {
                    here->DIOarea = 1;
                }
            }
            if(!here->DIOpjGiven) {
                if((!here->DIOwGiven) && (!here->DIOlGiven))  {
                    here->DIOpj = model->DIOpj;
                } else {
                    here->DIOpj = 0;
                }
            }
            if(!here->DIOmGiven) {
                here->DIOm = 1;
            }

            here->DIOarea = here->DIOarea * here->DIOm;
            here->DIOpj = here->DIOpj * here->DIOm;
            if (model->DIOlevel == 3) {
                if((here->DIOwGiven) && (here->DIOlGiven))  {
                    here->DIOarea = here->DIOw * here->DIOl * here->DIOm;
                    here->DIOpj = (2 * here->DIOw + 2 * here->DIOl) * here->DIOm;
                }
            }
            here->DIOforwardKneeCurrent = model->DIOforwardKneeCurrent * here->DIOarea;
            here->DIOreverseKneeCurrent = model->DIOreverseKneeCurrent * here->DIOarea;
            here->DIOjunctionCap = model->DIOjunctionCap * here->DIOarea;
            here->DIOjunctionSWCap = model->DIOjunctionSWCap * here->DIOpj;

            here->DIOstate = *states;
            *states += 5;
            if(ckt->CKTsenInfo && (ckt->CKTsenInfo->SENmode & TRANSEN) ){
                *states += 2 * (ckt->CKTsenInfo->SENparms);
            }

            if(model->DIOresist == 0) {

                here->DIOposPrimeNode = here->DIOposNode;

            } else if(here->DIOposPrimeNode == 0) {

               CKTnode *tmpNode;
               IFuid tmpName;

                error = CKTmkVolt(ckt,&tmp,here->DIOname,"internal");
                if(error) return(error);
                here->DIOposPrimeNode = tmp->number;
                if (ckt->CKTcopyNodesets) {
                  if (CKTinst2Node(ckt,here,1,&tmpNode,&tmpName)==OK) {
                     if (tmpNode->nsGiven) {
                       tmp->nodeset=tmpNode->nodeset;
                       tmp->nsGiven=tmpNode->nsGiven;
                     }
                  }
                }
            }

/* macro to make elements with built in test for out of memory */
#define TSTALLOC(ptr,first,second) \
do { if((here->ptr = SMPmakeElt(matrix, here->first, here->second)) == NULL){\
    return(E_NOMEM);\
} } while(0)

            TSTALLOC(DIOposPosPrimePtr,DIOposNode,DIOposPrimeNode);
            TSTALLOC(DIOnegPosPrimePtr,DIOnegNode,DIOposPrimeNode);
            TSTALLOC(DIOposPrimePosPtr,DIOposPrimeNode,DIOposNode);
            TSTALLOC(DIOposPrimeNegPtr,DIOposPrimeNode,DIOnegNode);
            TSTALLOC(DIOposPosPtr,DIOposNode,DIOposNode);
            TSTALLOC(DIOnegNegPtr,DIOnegNode,DIOnegNode);
            TSTALLOC(DIOposPrimePosPrimePtr,DIOposPrimeNode,DIOposPrimeNode);
        }
    }
    return(OK);
}
Ejemplo n.º 12
0
Archivo: b3set.c Proyecto: imr/ngspice
int
BSIM3setup(
SMPmatrix *matrix,
GENmodel *inModel,
CKTcircuit *ckt,
int *states)
{
BSIM3model *model = (BSIM3model*)inModel;
BSIM3instance *here;
int error;
CKTnode *tmp;
CKTnode *tmpNode;
IFuid tmpName;

#ifdef USE_OMP
int idx, InstCount;
BSIM3instance **InstArray;
#endif


    /*  loop through all the BSIM3 device models */
    for( ; model != NULL; model = BSIM3nextModel(model))
    {
/* Default value Processing for BSIM3 MOSFET Models */
        if (!model->BSIM3typeGiven)
            model->BSIM3type = NMOS;
        if (!model->BSIM3mobModGiven)
            model->BSIM3mobMod = 1;
        if (!model->BSIM3binUnitGiven)
            model->BSIM3binUnit = 1;
        if (!model->BSIM3paramChkGiven)
            model->BSIM3paramChk = 0;
        if (!model->BSIM3capModGiven)
            model->BSIM3capMod = 3;
        if (!model->BSIM3acmModGiven)
            model->BSIM3acmMod = 0;
        if (!model->BSIM3calcacmGiven)
            model->BSIM3calcacm = 0;
        if (!model->BSIM3noiModGiven)
            model->BSIM3noiMod = 1;
        if (!model->BSIM3nqsModGiven)
            model->BSIM3nqsMod = 0;
        else if ((model->BSIM3nqsMod != 0) && (model->BSIM3nqsMod != 1))
        {   model->BSIM3nqsMod = 0;
            printf("Warning: nqsMod has been set to its default value: 0.\n");
        }
        if (!model->BSIM3acnqsModGiven)
            model->BSIM3acnqsMod = 0;
        else if ((model->BSIM3acnqsMod != 0) && (model->BSIM3acnqsMod != 1))
        {   model->BSIM3acnqsMod = 0;
            printf("Warning: acnqsMod has been set to its default value: 0.\n");
        }
        if (!model->BSIM3versionGiven)
            model->BSIM3version = copy("3.3.0");
        if (!model->BSIM3toxGiven)
            model->BSIM3tox = 150.0e-10;
        model->BSIM3cox = 3.453133e-11 / model->BSIM3tox;
        if (!model->BSIM3toxmGiven)
            model->BSIM3toxm = model->BSIM3tox;

        if (!model->BSIM3cdscGiven)
            model->BSIM3cdsc = 2.4e-4;   /* unit Q/V/m^2  */
        if (!model->BSIM3cdscbGiven)
            model->BSIM3cdscb = 0.0;   /* unit Q/V/m^2  */
        if (!model->BSIM3cdscdGiven)
            model->BSIM3cdscd = 0.0;   /* unit Q/V/m^2  */
        if (!model->BSIM3citGiven)
            model->BSIM3cit = 0.0;   /* unit Q/V/m^2  */
        if (!model->BSIM3nfactorGiven)
            model->BSIM3nfactor = 1;
        if (!model->BSIM3xjGiven)
            model->BSIM3xj = .15e-6;
        if (!model->BSIM3vsatGiven)
            model->BSIM3vsat = 8.0e4;    /* unit m/s */
        if (!model->BSIM3atGiven)
            model->BSIM3at = 3.3e4;    /* unit m/s */
        if (!model->BSIM3a0Given)
            model->BSIM3a0 = 1.0;
        if (!model->BSIM3agsGiven)
            model->BSIM3ags = 0.0;
        if (!model->BSIM3a1Given)
            model->BSIM3a1 = 0.0;
        if (!model->BSIM3a2Given)
            model->BSIM3a2 = 1.0;
        if (!model->BSIM3ketaGiven)
            model->BSIM3keta = -0.047;    /* unit  / V */
        if (!model->BSIM3nsubGiven)
            model->BSIM3nsub = 6.0e16;   /* unit 1/cm3 */
        if (!model->BSIM3npeakGiven)
            model->BSIM3npeak = 1.7e17;   /* unit 1/cm3 */
        if (!model->BSIM3ngateGiven)
            model->BSIM3ngate = 0;   /* unit 1/cm3 */
        if (!model->BSIM3vbmGiven)
            model->BSIM3vbm = -3.0;
        if (!model->BSIM3xtGiven)
            model->BSIM3xt = 1.55e-7;
        if (!model->BSIM3kt1Given)
            model->BSIM3kt1 = -0.11;      /* unit V */
        if (!model->BSIM3kt1lGiven)
            model->BSIM3kt1l = 0.0;      /* unit V*m */
        if (!model->BSIM3kt2Given)
            model->BSIM3kt2 = 0.022;      /* No unit */
        if (!model->BSIM3k3Given)
            model->BSIM3k3 = 80.0;
        if (!model->BSIM3k3bGiven)
            model->BSIM3k3b = 0.0;
        if (!model->BSIM3w0Given)
            model->BSIM3w0 = 2.5e-6;
        if (!model->BSIM3nlxGiven)
            model->BSIM3nlx = 1.74e-7;
        if (!model->BSIM3dvt0Given)
            model->BSIM3dvt0 = 2.2;
        if (!model->BSIM3dvt1Given)
            model->BSIM3dvt1 = 0.53;
        if (!model->BSIM3dvt2Given)
            model->BSIM3dvt2 = -0.032;   /* unit 1 / V */

        if (!model->BSIM3dvt0wGiven)
            model->BSIM3dvt0w = 0.0;
        if (!model->BSIM3dvt1wGiven)
            model->BSIM3dvt1w = 5.3e6;
        if (!model->BSIM3dvt2wGiven)
            model->BSIM3dvt2w = -0.032;

        if (!model->BSIM3droutGiven)
            model->BSIM3drout = 0.56;
        if (!model->BSIM3dsubGiven)
            model->BSIM3dsub = model->BSIM3drout;
        if (!model->BSIM3vth0Given)
            model->BSIM3vth0 = (model->BSIM3type == NMOS) ? 0.7 : -0.7;
        if (!model->BSIM3uaGiven)
            model->BSIM3ua = 2.25e-9;      /* unit m/V */
        if (!model->BSIM3ua1Given)
            model->BSIM3ua1 = 4.31e-9;      /* unit m/V */
        if (!model->BSIM3ubGiven)
            model->BSIM3ub = 5.87e-19;     /* unit (m/V)**2 */
        if (!model->BSIM3ub1Given)
            model->BSIM3ub1 = -7.61e-18;     /* unit (m/V)**2 */
        if (!model->BSIM3ucGiven)
            model->BSIM3uc = (model->BSIM3mobMod == 3) ? -0.0465 : -0.0465e-9;
        if (!model->BSIM3uc1Given)
            model->BSIM3uc1 = (model->BSIM3mobMod == 3) ? -0.056 : -0.056e-9;
        if (!model->BSIM3u0Given)
            model->BSIM3u0 = (model->BSIM3type == NMOS) ? 0.067 : 0.025;
        if (!model->BSIM3uteGiven)
            model->BSIM3ute = -1.5;
        if (!model->BSIM3voffGiven)
            model->BSIM3voff = -0.08;
        if (!model->BSIM3deltaGiven)
           model->BSIM3delta = 0.01;
        if (!model->BSIM3rdswGiven)
            model->BSIM3rdsw = 0;
        if (!model->BSIM3prwgGiven)
            model->BSIM3prwg = 0.0;      /* unit 1/V */
        if (!model->BSIM3prwbGiven)
            model->BSIM3prwb = 0.0;
        if (!model->BSIM3prtGiven)
            model->BSIM3prt = 0.0;
        if (!model->BSIM3eta0Given)
            model->BSIM3eta0 = 0.08;      /* no unit  */
        if (!model->BSIM3etabGiven)
            model->BSIM3etab = -0.07;      /* unit  1/V */
        if (!model->BSIM3pclmGiven)
            model->BSIM3pclm = 1.3;      /* no unit  */
        if (!model->BSIM3pdibl1Given)
            model->BSIM3pdibl1 = .39;    /* no unit  */
        if (!model->BSIM3pdibl2Given)
            model->BSIM3pdibl2 = 0.0086;    /* no unit  */
        if (!model->BSIM3pdiblbGiven)
            model->BSIM3pdiblb = 0.0;    /* 1/V  */
        if (!model->BSIM3pscbe1Given)
            model->BSIM3pscbe1 = 4.24e8;
        if (!model->BSIM3pscbe2Given)
            model->BSIM3pscbe2 = 1.0e-5;
        if (!model->BSIM3pvagGiven)
            model->BSIM3pvag = 0.0;
        if (!model->BSIM3wrGiven)
            model->BSIM3wr = 1.0;
        if (!model->BSIM3dwgGiven)
            model->BSIM3dwg = 0.0;
        if (!model->BSIM3dwbGiven)
            model->BSIM3dwb = 0.0;
        if (!model->BSIM3b0Given)
            model->BSIM3b0 = 0.0;
        if (!model->BSIM3b1Given)
            model->BSIM3b1 = 0.0;
        if (!model->BSIM3alpha0Given)
            model->BSIM3alpha0 = 0.0;
        if (!model->BSIM3alpha1Given)
            model->BSIM3alpha1 = 0.0;
        if (!model->BSIM3beta0Given)
            model->BSIM3beta0 = 30.0;
        if (!model->BSIM3ijthGiven)
            model->BSIM3ijth = 0.1; /* unit A */

        if (!model->BSIM3elmGiven)
            model->BSIM3elm = 5.0;
        if (!model->BSIM3cgslGiven)
            model->BSIM3cgsl = 0.0;
        if (!model->BSIM3cgdlGiven)
            model->BSIM3cgdl = 0.0;
        if (!model->BSIM3ckappaGiven)
            model->BSIM3ckappa = 0.6;
        if (!model->BSIM3clcGiven)
            model->BSIM3clc = 0.1e-6;
        if (!model->BSIM3cleGiven)
            model->BSIM3cle = 0.6;
        if (!model->BSIM3vfbcvGiven)
            model->BSIM3vfbcv = -1.0;
        if (!model->BSIM3acdeGiven)
            model->BSIM3acde = 1.0;
        if (!model->BSIM3moinGiven)
            model->BSIM3moin = 15.0;
        if (!model->BSIM3noffGiven)
            model->BSIM3noff = 1.0;
        if (!model->BSIM3voffcvGiven)
            model->BSIM3voffcv = 0.0;
        if (!model->BSIM3tcjGiven)
            model->BSIM3tcj = 0.0;
        if (!model->BSIM3tpbGiven)
            model->BSIM3tpb = 0.0;
        if (!model->BSIM3tcjswGiven)
            model->BSIM3tcjsw = 0.0;
        if (!model->BSIM3tpbswGiven)
            model->BSIM3tpbsw = 0.0;
        if (!model->BSIM3tcjswgGiven)
            model->BSIM3tcjswg = 0.0;
        if (!model->BSIM3tpbswgGiven)
            model->BSIM3tpbswg = 0.0;

        /* ACM model */
        if (!model->BSIM3hdifGiven)
          model->BSIM3hdif = 0.0;
        if (!model->BSIM3ldifGiven)
          model->BSIM3ldif = 0.0;
        if (!model->BSIM3ldGiven)
          model->BSIM3ld = 0.0;
        if (!model->BSIM3rdGiven)
          model->BSIM3rd = 0.0;
        if (!model->BSIM3rsGiven)
          model->BSIM3rs = 0.0;
        if (!model->BSIM3rdcGiven)
          model->BSIM3rdc = 0.0;
        if (!model->BSIM3rscGiven)
          model->BSIM3rsc = 0.0;
        if (!model->BSIM3wmltGiven)
          model->BSIM3wmlt = 1.0;

        /* Length dependence */
        if (!model->BSIM3lcdscGiven)
            model->BSIM3lcdsc = 0.0;
        if (!model->BSIM3lcdscbGiven)
            model->BSIM3lcdscb = 0.0;
        if (!model->BSIM3lcdscdGiven)
            model->BSIM3lcdscd = 0.0;
        if (!model->BSIM3lcitGiven)
            model->BSIM3lcit = 0.0;
        if (!model->BSIM3lnfactorGiven)
            model->BSIM3lnfactor = 0.0;
        if (!model->BSIM3lxjGiven)
            model->BSIM3lxj = 0.0;
        if (!model->BSIM3lvsatGiven)
            model->BSIM3lvsat = 0.0;
        if (!model->BSIM3latGiven)
            model->BSIM3lat = 0.0;
        if (!model->BSIM3la0Given)
            model->BSIM3la0 = 0.0;
        if (!model->BSIM3lagsGiven)
            model->BSIM3lags = 0.0;
        if (!model->BSIM3la1Given)
            model->BSIM3la1 = 0.0;
        if (!model->BSIM3la2Given)
            model->BSIM3la2 = 0.0;
        if (!model->BSIM3lketaGiven)
            model->BSIM3lketa = 0.0;
        if (!model->BSIM3lnsubGiven)
            model->BSIM3lnsub = 0.0;
        if (!model->BSIM3lnpeakGiven)
            model->BSIM3lnpeak = 0.0;
        if (!model->BSIM3lngateGiven)
            model->BSIM3lngate = 0.0;
        if (!model->BSIM3lvbmGiven)
            model->BSIM3lvbm = 0.0;
        if (!model->BSIM3lxtGiven)
            model->BSIM3lxt = 0.0;
        if (!model->BSIM3lkt1Given)
            model->BSIM3lkt1 = 0.0;
        if (!model->BSIM3lkt1lGiven)
            model->BSIM3lkt1l = 0.0;
        if (!model->BSIM3lkt2Given)
            model->BSIM3lkt2 = 0.0;
        if (!model->BSIM3lk3Given)
            model->BSIM3lk3 = 0.0;
        if (!model->BSIM3lk3bGiven)
            model->BSIM3lk3b = 0.0;
        if (!model->BSIM3lw0Given)
            model->BSIM3lw0 = 0.0;
        if (!model->BSIM3lnlxGiven)
            model->BSIM3lnlx = 0.0;
        if (!model->BSIM3ldvt0Given)
            model->BSIM3ldvt0 = 0.0;
        if (!model->BSIM3ldvt1Given)
            model->BSIM3ldvt1 = 0.0;
        if (!model->BSIM3ldvt2Given)
            model->BSIM3ldvt2 = 0.0;
        if (!model->BSIM3ldvt0wGiven)
            model->BSIM3ldvt0w = 0.0;
        if (!model->BSIM3ldvt1wGiven)
            model->BSIM3ldvt1w = 0.0;
        if (!model->BSIM3ldvt2wGiven)
            model->BSIM3ldvt2w = 0.0;
        if (!model->BSIM3ldroutGiven)
            model->BSIM3ldrout = 0.0;
        if (!model->BSIM3ldsubGiven)
            model->BSIM3ldsub = 0.0;
        if (!model->BSIM3lvth0Given)
           model->BSIM3lvth0 = 0.0;
        if (!model->BSIM3luaGiven)
            model->BSIM3lua = 0.0;
        if (!model->BSIM3lua1Given)
            model->BSIM3lua1 = 0.0;
        if (!model->BSIM3lubGiven)
            model->BSIM3lub = 0.0;
        if (!model->BSIM3lub1Given)
            model->BSIM3lub1 = 0.0;
        if (!model->BSIM3lucGiven)
            model->BSIM3luc = 0.0;
        if (!model->BSIM3luc1Given)
            model->BSIM3luc1 = 0.0;
        if (!model->BSIM3lu0Given)
            model->BSIM3lu0 = 0.0;
        if (!model->BSIM3luteGiven)
            model->BSIM3lute = 0.0;
        if (!model->BSIM3lvoffGiven)
            model->BSIM3lvoff = 0.0;
        if (!model->BSIM3ldeltaGiven)
            model->BSIM3ldelta = 0.0;
        if (!model->BSIM3lrdswGiven)
            model->BSIM3lrdsw = 0.0;
        if (!model->BSIM3lprwbGiven)
            model->BSIM3lprwb = 0.0;
        if (!model->BSIM3lprwgGiven)
            model->BSIM3lprwg = 0.0;
        if (!model->BSIM3lprtGiven)
            model->BSIM3lprt = 0.0;
        if (!model->BSIM3leta0Given)
            model->BSIM3leta0 = 0.0;
        if (!model->BSIM3letabGiven)
            model->BSIM3letab = -0.0;
        if (!model->BSIM3lpclmGiven)
            model->BSIM3lpclm = 0.0;
        if (!model->BSIM3lpdibl1Given)
            model->BSIM3lpdibl1 = 0.0;
        if (!model->BSIM3lpdibl2Given)
            model->BSIM3lpdibl2 = 0.0;
        if (!model->BSIM3lpdiblbGiven)
            model->BSIM3lpdiblb = 0.0;
        if (!model->BSIM3lpscbe1Given)
            model->BSIM3lpscbe1 = 0.0;
        if (!model->BSIM3lpscbe2Given)
            model->BSIM3lpscbe2 = 0.0;
        if (!model->BSIM3lpvagGiven)
            model->BSIM3lpvag = 0.0;
        if (!model->BSIM3lwrGiven)
            model->BSIM3lwr = 0.0;
        if (!model->BSIM3ldwgGiven)
            model->BSIM3ldwg = 0.0;
        if (!model->BSIM3ldwbGiven)
            model->BSIM3ldwb = 0.0;
        if (!model->BSIM3lb0Given)
            model->BSIM3lb0 = 0.0;
        if (!model->BSIM3lb1Given)
            model->BSIM3lb1 = 0.0;
        if (!model->BSIM3lalpha0Given)
            model->BSIM3lalpha0 = 0.0;
        if (!model->BSIM3lalpha1Given)
            model->BSIM3lalpha1 = 0.0;
        if (!model->BSIM3lbeta0Given)
            model->BSIM3lbeta0 = 0.0;
        if (!model->BSIM3lvfbGiven)
            model->BSIM3lvfb = 0.0;

        if (!model->BSIM3lelmGiven)
            model->BSIM3lelm = 0.0;
        if (!model->BSIM3lcgslGiven)
            model->BSIM3lcgsl = 0.0;
        if (!model->BSIM3lcgdlGiven)
            model->BSIM3lcgdl = 0.0;
        if (!model->BSIM3lckappaGiven)
            model->BSIM3lckappa = 0.0;
        if (!model->BSIM3lclcGiven)
            model->BSIM3lclc = 0.0;
        if (!model->BSIM3lcleGiven)
            model->BSIM3lcle = 0.0;
        if (!model->BSIM3lcfGiven)
            model->BSIM3lcf = 0.0;
        if (!model->BSIM3lvfbcvGiven)
            model->BSIM3lvfbcv = 0.0;
        if (!model->BSIM3lacdeGiven)
            model->BSIM3lacde = 0.0;
        if (!model->BSIM3lmoinGiven)
            model->BSIM3lmoin = 0.0;
        if (!model->BSIM3lnoffGiven)
            model->BSIM3lnoff = 0.0;
        if (!model->BSIM3lvoffcvGiven)
            model->BSIM3lvoffcv = 0.0;

        /* Width dependence */
        if (!model->BSIM3wcdscGiven)
            model->BSIM3wcdsc = 0.0;
        if (!model->BSIM3wcdscbGiven)
            model->BSIM3wcdscb = 0.0;
        if (!model->BSIM3wcdscdGiven)
            model->BSIM3wcdscd = 0.0;
        if (!model->BSIM3wcitGiven)
            model->BSIM3wcit = 0.0;
        if (!model->BSIM3wnfactorGiven)
            model->BSIM3wnfactor = 0.0;
        if (!model->BSIM3wxjGiven)
            model->BSIM3wxj = 0.0;
        if (!model->BSIM3wvsatGiven)
            model->BSIM3wvsat = 0.0;
        if (!model->BSIM3watGiven)
            model->BSIM3wat = 0.0;
        if (!model->BSIM3wa0Given)
            model->BSIM3wa0 = 0.0;
        if (!model->BSIM3wagsGiven)
            model->BSIM3wags = 0.0;
        if (!model->BSIM3wa1Given)
            model->BSIM3wa1 = 0.0;
        if (!model->BSIM3wa2Given)
            model->BSIM3wa2 = 0.0;
        if (!model->BSIM3wketaGiven)
            model->BSIM3wketa = 0.0;
        if (!model->BSIM3wnsubGiven)
            model->BSIM3wnsub = 0.0;
        if (!model->BSIM3wnpeakGiven)
            model->BSIM3wnpeak = 0.0;
        if (!model->BSIM3wngateGiven)
            model->BSIM3wngate = 0.0;
        if (!model->BSIM3wvbmGiven)
            model->BSIM3wvbm = 0.0;
        if (!model->BSIM3wxtGiven)
            model->BSIM3wxt = 0.0;
        if (!model->BSIM3wkt1Given)
            model->BSIM3wkt1 = 0.0;
        if (!model->BSIM3wkt1lGiven)
            model->BSIM3wkt1l = 0.0;
        if (!model->BSIM3wkt2Given)
            model->BSIM3wkt2 = 0.0;
        if (!model->BSIM3wk3Given)
            model->BSIM3wk3 = 0.0;
        if (!model->BSIM3wk3bGiven)
            model->BSIM3wk3b = 0.0;
        if (!model->BSIM3ww0Given)
            model->BSIM3ww0 = 0.0;
        if (!model->BSIM3wnlxGiven)
            model->BSIM3wnlx = 0.0;
        if (!model->BSIM3wdvt0Given)
            model->BSIM3wdvt0 = 0.0;
        if (!model->BSIM3wdvt1Given)
            model->BSIM3wdvt1 = 0.0;
        if (!model->BSIM3wdvt2Given)
            model->BSIM3wdvt2 = 0.0;
        if (!model->BSIM3wdvt0wGiven)
            model->BSIM3wdvt0w = 0.0;
        if (!model->BSIM3wdvt1wGiven)
            model->BSIM3wdvt1w = 0.0;
        if (!model->BSIM3wdvt2wGiven)
            model->BSIM3wdvt2w = 0.0;
        if (!model->BSIM3wdroutGiven)
            model->BSIM3wdrout = 0.0;
        if (!model->BSIM3wdsubGiven)
            model->BSIM3wdsub = 0.0;
        if (!model->BSIM3wvth0Given)
           model->BSIM3wvth0 = 0.0;
        if (!model->BSIM3wuaGiven)
            model->BSIM3wua = 0.0;
        if (!model->BSIM3wua1Given)
            model->BSIM3wua1 = 0.0;
        if (!model->BSIM3wubGiven)
            model->BSIM3wub = 0.0;
        if (!model->BSIM3wub1Given)
            model->BSIM3wub1 = 0.0;
        if (!model->BSIM3wucGiven)
            model->BSIM3wuc = 0.0;
        if (!model->BSIM3wuc1Given)
            model->BSIM3wuc1 = 0.0;
        if (!model->BSIM3wu0Given)
            model->BSIM3wu0 = 0.0;
        if (!model->BSIM3wuteGiven)
            model->BSIM3wute = 0.0;
        if (!model->BSIM3wvoffGiven)
            model->BSIM3wvoff = 0.0;
        if (!model->BSIM3wdeltaGiven)
            model->BSIM3wdelta = 0.0;
        if (!model->BSIM3wrdswGiven)
            model->BSIM3wrdsw = 0.0;
        if (!model->BSIM3wprwbGiven)
            model->BSIM3wprwb = 0.0;
        if (!model->BSIM3wprwgGiven)
            model->BSIM3wprwg = 0.0;
        if (!model->BSIM3wprtGiven)
            model->BSIM3wprt = 0.0;
        if (!model->BSIM3weta0Given)
            model->BSIM3weta0 = 0.0;
        if (!model->BSIM3wetabGiven)
            model->BSIM3wetab = 0.0;
        if (!model->BSIM3wpclmGiven)
            model->BSIM3wpclm = 0.0;
        if (!model->BSIM3wpdibl1Given)
            model->BSIM3wpdibl1 = 0.0;
        if (!model->BSIM3wpdibl2Given)
            model->BSIM3wpdibl2 = 0.0;
        if (!model->BSIM3wpdiblbGiven)
            model->BSIM3wpdiblb = 0.0;
        if (!model->BSIM3wpscbe1Given)
            model->BSIM3wpscbe1 = 0.0;
        if (!model->BSIM3wpscbe2Given)
            model->BSIM3wpscbe2 = 0.0;
        if (!model->BSIM3wpvagGiven)
            model->BSIM3wpvag = 0.0;
        if (!model->BSIM3wwrGiven)
            model->BSIM3wwr = 0.0;
        if (!model->BSIM3wdwgGiven)
            model->BSIM3wdwg = 0.0;
        if (!model->BSIM3wdwbGiven)
            model->BSIM3wdwb = 0.0;
        if (!model->BSIM3wb0Given)
            model->BSIM3wb0 = 0.0;
        if (!model->BSIM3wb1Given)
            model->BSIM3wb1 = 0.0;
        if (!model->BSIM3walpha0Given)
            model->BSIM3walpha0 = 0.0;
        if (!model->BSIM3walpha1Given)
            model->BSIM3walpha1 = 0.0;
        if (!model->BSIM3wbeta0Given)
            model->BSIM3wbeta0 = 0.0;
        if (!model->BSIM3wvfbGiven)
            model->BSIM3wvfb = 0.0;

        if (!model->BSIM3welmGiven)
            model->BSIM3welm = 0.0;
        if (!model->BSIM3wcgslGiven)
            model->BSIM3wcgsl = 0.0;
        if (!model->BSIM3wcgdlGiven)
            model->BSIM3wcgdl = 0.0;
        if (!model->BSIM3wckappaGiven)
            model->BSIM3wckappa = 0.0;
        if (!model->BSIM3wcfGiven)
            model->BSIM3wcf = 0.0;
        if (!model->BSIM3wclcGiven)
            model->BSIM3wclc = 0.0;
        if (!model->BSIM3wcleGiven)
            model->BSIM3wcle = 0.0;
        if (!model->BSIM3wvfbcvGiven)
            model->BSIM3wvfbcv = 0.0;
        if (!model->BSIM3wacdeGiven)
            model->BSIM3wacde = 0.0;
        if (!model->BSIM3wmoinGiven)
            model->BSIM3wmoin = 0.0;
        if (!model->BSIM3wnoffGiven)
            model->BSIM3wnoff = 0.0;
        if (!model->BSIM3wvoffcvGiven)
            model->BSIM3wvoffcv = 0.0;

        /* Cross-term dependence */
        if (!model->BSIM3pcdscGiven)
            model->BSIM3pcdsc = 0.0;
        if (!model->BSIM3pcdscbGiven)
            model->BSIM3pcdscb = 0.0;
        if (!model->BSIM3pcdscdGiven)
            model->BSIM3pcdscd = 0.0;
        if (!model->BSIM3pcitGiven)
            model->BSIM3pcit = 0.0;
        if (!model->BSIM3pnfactorGiven)
            model->BSIM3pnfactor = 0.0;
        if (!model->BSIM3pxjGiven)
            model->BSIM3pxj = 0.0;
        if (!model->BSIM3pvsatGiven)
            model->BSIM3pvsat = 0.0;
        if (!model->BSIM3patGiven)
            model->BSIM3pat = 0.0;
        if (!model->BSIM3pa0Given)
            model->BSIM3pa0 = 0.0;

        if (!model->BSIM3pagsGiven)
            model->BSIM3pags = 0.0;
        if (!model->BSIM3pa1Given)
            model->BSIM3pa1 = 0.0;
        if (!model->BSIM3pa2Given)
            model->BSIM3pa2 = 0.0;
        if (!model->BSIM3pketaGiven)
            model->BSIM3pketa = 0.0;
        if (!model->BSIM3pnsubGiven)
            model->BSIM3pnsub = 0.0;
        if (!model->BSIM3pnpeakGiven)
            model->BSIM3pnpeak = 0.0;
        if (!model->BSIM3pngateGiven)
            model->BSIM3pngate = 0.0;
        if (!model->BSIM3pvbmGiven)
            model->BSIM3pvbm = 0.0;
        if (!model->BSIM3pxtGiven)
            model->BSIM3pxt = 0.0;
        if (!model->BSIM3pkt1Given)
            model->BSIM3pkt1 = 0.0;
        if (!model->BSIM3pkt1lGiven)
            model->BSIM3pkt1l = 0.0;
        if (!model->BSIM3pkt2Given)
            model->BSIM3pkt2 = 0.0;
        if (!model->BSIM3pk3Given)
            model->BSIM3pk3 = 0.0;
        if (!model->BSIM3pk3bGiven)
            model->BSIM3pk3b = 0.0;
        if (!model->BSIM3pw0Given)
            model->BSIM3pw0 = 0.0;
        if (!model->BSIM3pnlxGiven)
            model->BSIM3pnlx = 0.0;
        if (!model->BSIM3pdvt0Given)
            model->BSIM3pdvt0 = 0.0;
        if (!model->BSIM3pdvt1Given)
            model->BSIM3pdvt1 = 0.0;
        if (!model->BSIM3pdvt2Given)
            model->BSIM3pdvt2 = 0.0;
        if (!model->BSIM3pdvt0wGiven)
            model->BSIM3pdvt0w = 0.0;
        if (!model->BSIM3pdvt1wGiven)
            model->BSIM3pdvt1w = 0.0;
        if (!model->BSIM3pdvt2wGiven)
            model->BSIM3pdvt2w = 0.0;
        if (!model->BSIM3pdroutGiven)
            model->BSIM3pdrout = 0.0;
        if (!model->BSIM3pdsubGiven)
            model->BSIM3pdsub = 0.0;
        if (!model->BSIM3pvth0Given)
           model->BSIM3pvth0 = 0.0;
        if (!model->BSIM3puaGiven)
            model->BSIM3pua = 0.0;
        if (!model->BSIM3pua1Given)
            model->BSIM3pua1 = 0.0;
        if (!model->BSIM3pubGiven)
            model->BSIM3pub = 0.0;
        if (!model->BSIM3pub1Given)
            model->BSIM3pub1 = 0.0;
        if (!model->BSIM3pucGiven)
            model->BSIM3puc = 0.0;
        if (!model->BSIM3puc1Given)
            model->BSIM3puc1 = 0.0;
        if (!model->BSIM3pu0Given)
            model->BSIM3pu0 = 0.0;
        if (!model->BSIM3puteGiven)
            model->BSIM3pute = 0.0;
        if (!model->BSIM3pvoffGiven)
            model->BSIM3pvoff = 0.0;
        if (!model->BSIM3pdeltaGiven)
            model->BSIM3pdelta = 0.0;
        if (!model->BSIM3prdswGiven)
            model->BSIM3prdsw = 0.0;
        if (!model->BSIM3pprwbGiven)
            model->BSIM3pprwb = 0.0;
        if (!model->BSIM3pprwgGiven)
            model->BSIM3pprwg = 0.0;
        if (!model->BSIM3pprtGiven)
            model->BSIM3pprt = 0.0;
        if (!model->BSIM3peta0Given)
            model->BSIM3peta0 = 0.0;
        if (!model->BSIM3petabGiven)
            model->BSIM3petab = 0.0;
        if (!model->BSIM3ppclmGiven)
            model->BSIM3ppclm = 0.0;
        if (!model->BSIM3ppdibl1Given)
            model->BSIM3ppdibl1 = 0.0;
        if (!model->BSIM3ppdibl2Given)
            model->BSIM3ppdibl2 = 0.0;
        if (!model->BSIM3ppdiblbGiven)
            model->BSIM3ppdiblb = 0.0;
        if (!model->BSIM3ppscbe1Given)
            model->BSIM3ppscbe1 = 0.0;
        if (!model->BSIM3ppscbe2Given)
            model->BSIM3ppscbe2 = 0.0;
        if (!model->BSIM3ppvagGiven)
            model->BSIM3ppvag = 0.0;
        if (!model->BSIM3pwrGiven)
            model->BSIM3pwr = 0.0;
        if (!model->BSIM3pdwgGiven)
            model->BSIM3pdwg = 0.0;
        if (!model->BSIM3pdwbGiven)
            model->BSIM3pdwb = 0.0;
        if (!model->BSIM3pb0Given)
            model->BSIM3pb0 = 0.0;
        if (!model->BSIM3pb1Given)
            model->BSIM3pb1 = 0.0;
        if (!model->BSIM3palpha0Given)
            model->BSIM3palpha0 = 0.0;
        if (!model->BSIM3palpha1Given)
            model->BSIM3palpha1 = 0.0;
        if (!model->BSIM3pbeta0Given)
            model->BSIM3pbeta0 = 0.0;
        if (!model->BSIM3pvfbGiven)
            model->BSIM3pvfb = 0.0;

        if (!model->BSIM3pelmGiven)
            model->BSIM3pelm = 0.0;
        if (!model->BSIM3pcgslGiven)
            model->BSIM3pcgsl = 0.0;
        if (!model->BSIM3pcgdlGiven)
            model->BSIM3pcgdl = 0.0;
        if (!model->BSIM3pckappaGiven)
            model->BSIM3pckappa = 0.0;
        if (!model->BSIM3pcfGiven)
            model->BSIM3pcf = 0.0;
        if (!model->BSIM3pclcGiven)
            model->BSIM3pclc = 0.0;
        if (!model->BSIM3pcleGiven)
            model->BSIM3pcle = 0.0;
        if (!model->BSIM3pvfbcvGiven)
            model->BSIM3pvfbcv = 0.0;
        if (!model->BSIM3pacdeGiven)
            model->BSIM3pacde = 0.0;
        if (!model->BSIM3pmoinGiven)
            model->BSIM3pmoin = 0.0;
        if (!model->BSIM3pnoffGiven)
            model->BSIM3pnoff = 0.0;
        if (!model->BSIM3pvoffcvGiven)
            model->BSIM3pvoffcv = 0.0;

        /* unit degree celcius */
        if (!model->BSIM3tnomGiven)
            model->BSIM3tnom = ckt->CKTnomTemp;
/*        else
            model->BSIM3tnom = model->BSIM3tnom + 273.15; we make this transform in b3mpar.c in the first run */
        if (!model->BSIM3lintnoiGiven)
            model->BSIM3lintnoi = 0.0;  /* unit m */
        if (!model->BSIM3LintGiven)
           model->BSIM3Lint = 0.0;
        if (!model->BSIM3LlGiven)
           model->BSIM3Ll = 0.0;
        if (!model->BSIM3LlcGiven)
           model->BSIM3Llc = model->BSIM3Ll;
        if (!model->BSIM3LlnGiven)
           model->BSIM3Lln = 1.0;
        if (!model->BSIM3LwGiven)
           model->BSIM3Lw = 0.0;
        if (!model->BSIM3LwcGiven)
           model->BSIM3Lwc = model->BSIM3Lw;
        if (!model->BSIM3LwnGiven)
           model->BSIM3Lwn = 1.0;
        if (!model->BSIM3LwlGiven)
           model->BSIM3Lwl = 0.0;
        if (!model->BSIM3LwlcGiven)
           model->BSIM3Lwlc = model->BSIM3Lwl;
        if (!model->BSIM3LminGiven)
           model->BSIM3Lmin = 0.0;
        if (!model->BSIM3LmaxGiven)
           model->BSIM3Lmax = 1.0;
        if (!model->BSIM3WintGiven)
           model->BSIM3Wint = 0.0;
        if (!model->BSIM3WlGiven)
           model->BSIM3Wl = 0.0;
        if (!model->BSIM3WlcGiven)
           model->BSIM3Wlc = model->BSIM3Wl;
        if (!model->BSIM3WlnGiven)
           model->BSIM3Wln = 1.0;
        if (!model->BSIM3WwGiven)
           model->BSIM3Ww = 0.0;
        if (!model->BSIM3WwcGiven)
           model->BSIM3Wwc = model->BSIM3Ww;
        if (!model->BSIM3WwnGiven)
           model->BSIM3Wwn = 1.0;
        if (!model->BSIM3WwlGiven)
           model->BSIM3Wwl = 0.0;
        if (!model->BSIM3WwlcGiven)
           model->BSIM3Wwlc = model->BSIM3Wwl;
        if (!model->BSIM3WminGiven)
           model->BSIM3Wmin = 0.0;
        if (!model->BSIM3WmaxGiven)
           model->BSIM3Wmax = 1.0;
        if (!model->BSIM3dwcGiven)
           model->BSIM3dwc = model->BSIM3Wint;
        if (!model->BSIM3dlcGiven)
           model->BSIM3dlc = model->BSIM3Lint;

        if (!model->BSIM3xlGiven)
           model->BSIM3xl = 0.0;
        if (!model->BSIM3xwGiven)
           model->BSIM3xw = 0.0;

        if (!model->BSIM3cfGiven)
            model->BSIM3cf = 2.0 * EPSOX / PI
                           * log(1.0 + 0.4e-6 / model->BSIM3tox);
        if (!model->BSIM3cgdoGiven)
        {   if (model->BSIM3dlcGiven && (model->BSIM3dlc > 0.0))
            {   model->BSIM3cgdo = model->BSIM3dlc * model->BSIM3cox
                                 - model->BSIM3cgdl ;
            }
            else
                model->BSIM3cgdo = 0.6 * model->BSIM3xj * model->BSIM3cox;
        }
        if (!model->BSIM3cgsoGiven)
        {   if (model->BSIM3dlcGiven && (model->BSIM3dlc > 0.0))
            {   model->BSIM3cgso = model->BSIM3dlc * model->BSIM3cox
                                 - model->BSIM3cgsl ;
            }
            else
                model->BSIM3cgso = 0.6 * model->BSIM3xj * model->BSIM3cox;
        }

        if (!model->BSIM3cgboGiven)
        {   model->BSIM3cgbo = 2.0 * model->BSIM3dwc * model->BSIM3cox;
        }
        if (!model->BSIM3xpartGiven)
            model->BSIM3xpart = 0.0;
        if (!model->BSIM3sheetResistanceGiven)
            model->BSIM3sheetResistance = 0.0;
        if (!model->BSIM3unitAreaJctCapGiven)
            model->BSIM3unitAreaJctCap = 5.0E-4;
        if (!model->BSIM3unitLengthSidewallJctCapGiven)
            model->BSIM3unitLengthSidewallJctCap = 5.0E-10;
        if (!model->BSIM3unitLengthGateSidewallJctCapGiven)
            model->BSIM3unitLengthGateSidewallJctCap = model->BSIM3unitLengthSidewallJctCap ;
        if (!model->BSIM3jctSatCurDensityGiven)
            model->BSIM3jctSatCurDensity = 1.0E-4;
        if (!model->BSIM3jctSidewallSatCurDensityGiven)
            model->BSIM3jctSidewallSatCurDensity = 0.0;
        if (!model->BSIM3bulkJctPotentialGiven)
            model->BSIM3bulkJctPotential = 1.0;
        if (!model->BSIM3sidewallJctPotentialGiven)
            model->BSIM3sidewallJctPotential = 1.0;
        if (!model->BSIM3GatesidewallJctPotentialGiven)
            model->BSIM3GatesidewallJctPotential = model->BSIM3sidewallJctPotential;
        if (!model->BSIM3bulkJctBotGradingCoeffGiven)
            model->BSIM3bulkJctBotGradingCoeff = 0.5;
        if (!model->BSIM3bulkJctSideGradingCoeffGiven)
            model->BSIM3bulkJctSideGradingCoeff = 0.33;
        if (!model->BSIM3bulkJctGateSideGradingCoeffGiven)
            model->BSIM3bulkJctGateSideGradingCoeff = model->BSIM3bulkJctSideGradingCoeff;
        if (!model->BSIM3jctEmissionCoeffGiven)
            model->BSIM3jctEmissionCoeff = 1.0;
        if (!model->BSIM3jctTempExponentGiven)
            model->BSIM3jctTempExponent = 3.0;
        if (!model->BSIM3oxideTrapDensityAGiven)
        {   if (model->BSIM3type == NMOS)
                model->BSIM3oxideTrapDensityA = 1e20;
            else
                model->BSIM3oxideTrapDensityA=9.9e18;
        }
        if (!model->BSIM3oxideTrapDensityBGiven)
        {   if (model->BSIM3type == NMOS)
                model->BSIM3oxideTrapDensityB = 5e4;
            else
                model->BSIM3oxideTrapDensityB = 2.4e3;
        }
        if (!model->BSIM3oxideTrapDensityCGiven)
        {   if (model->BSIM3type == NMOS)
                model->BSIM3oxideTrapDensityC = -1.4e-12;
            else
                model->BSIM3oxideTrapDensityC = 1.4e-12;

        }
        if (!model->BSIM3emGiven)
            model->BSIM3em = 4.1e7; /* V/m */
        if (!model->BSIM3efGiven)
            model->BSIM3ef = 1.0;
        if (!model->BSIM3afGiven)
            model->BSIM3af = 1.0;
        if (!model->BSIM3kfGiven)
            model->BSIM3kf = 0.0;

        if (!model->BSIM3vgsMaxGiven)
            model->BSIM3vgsMax = 1e99;
        if (!model->BSIM3vgdMaxGiven)
            model->BSIM3vgdMax = 1e99;
        if (!model->BSIM3vgbMaxGiven)
            model->BSIM3vgbMax = 1e99;
        if (!model->BSIM3vdsMaxGiven)
            model->BSIM3vdsMax = 1e99;
        if (!model->BSIM3vbsMaxGiven)
            model->BSIM3vbsMax = 1e99;
        if (!model->BSIM3vbdMaxGiven)
            model->BSIM3vbdMax = 1e99;
        if (!model->BSIM3vgsrMaxGiven)
            model->BSIM3vgsrMax = 1e99;
        if (!model->BSIM3vgdrMaxGiven)
            model->BSIM3vgdrMax = 1e99;
        if (!model->BSIM3vgbrMaxGiven)
            model->BSIM3vgbrMax = 1e99;
        if (!model->BSIM3vbsrMaxGiven)
            model->BSIM3vbsrMax = 1e99;
        if (!model->BSIM3vbdrMaxGiven)
            model->BSIM3vbdrMax = 1e99;

        /* loop through all the instances of the model */
        for (here = BSIM3instances(model); here != NULL ;
             here=BSIM3nextInstance(here))
        {
            /* allocate a chunk of the state vector */
            here->BSIM3states = *states;
            *states += BSIM3numStates;

            /* perform the parameter defaulting */
            if (!here->BSIM3drainAreaGiven)
                here->BSIM3drainArea = 0.0;
            if (!here->BSIM3drainPerimeterGiven)
                here->BSIM3drainPerimeter = 0.0;
            if (!here->BSIM3drainSquaresGiven)
            {
                if (model->BSIM3acmMod == 0)
                  here->BSIM3drainSquares = 1.0;
                else
                  here->BSIM3drainSquares = 0.0;
            }
            if (!here->BSIM3delvtoGiven)
                here->BSIM3delvto = 0.0;
            if (!here->BSIM3mulu0Given)
                here->BSIM3mulu0 = 1.0;
            if (!here->BSIM3icVBSGiven)
                here->BSIM3icVBS = 0.0;
            if (!here->BSIM3icVDSGiven)
                here->BSIM3icVDS = 0.0;
            if (!here->BSIM3icVGSGiven)
                here->BSIM3icVGS = 0.0;
            if (!here->BSIM3lGiven)
                here->BSIM3l = 5.0e-6;
            if (!here->BSIM3sourceAreaGiven)
                here->BSIM3sourceArea = 0.0;
            if (!here->BSIM3sourcePerimeterGiven)
                here->BSIM3sourcePerimeter = 0.0;
            if (!here->BSIM3sourceSquaresGiven)
            {
                if (model->BSIM3acmMod == 0)
                  here->BSIM3sourceSquares = 1.0;
                else
                  here->BSIM3sourceSquares = 0.0;
            }
            if (!here->BSIM3wGiven)
                here->BSIM3w = 5.0e-6;
            if (!here->BSIM3nqsModGiven)
                here->BSIM3nqsMod = model->BSIM3nqsMod;
            else if ((here->BSIM3nqsMod != 0) && (here->BSIM3nqsMod != 1))
            {   here->BSIM3nqsMod = model->BSIM3nqsMod;
                printf("Warning: nqsMod has been set to its global value %d.\n",
                model->BSIM3nqsMod);
            }
            if (!here->BSIM3acnqsModGiven)
                here->BSIM3acnqsMod = model->BSIM3acnqsMod;
            else if ((here->BSIM3acnqsMod != 0) && (here->BSIM3acnqsMod != 1))
            {   here->BSIM3acnqsMod = model->BSIM3acnqsMod;
                printf("Warning: acnqsMod has been set to its global value %d.\n",
                model->BSIM3acnqsMod);
            }

            if (!here->BSIM3geoGiven)
                here->BSIM3geo = 0;

            if (!here->BSIM3mGiven)
                here->BSIM3m = 1;

            /* process source/drain series resistance */
            /* ACM model */

            double drainResistance, sourceResistance;

            if (model->BSIM3acmMod == 0)
            {
                drainResistance = model->BSIM3sheetResistance
                    * here->BSIM3drainSquares;
                sourceResistance = model->BSIM3sheetResistance
                    * here->BSIM3sourceSquares;
            }
            else /* ACM > 0 */
            {
                error = ACM_SourceDrainResistances(
                    model->BSIM3acmMod,
                    model->BSIM3ld,
                    model->BSIM3ldif,
                    model->BSIM3hdif,
                    model->BSIM3wmlt,
                    here->BSIM3w,
                    model->BSIM3xw,
                    model->BSIM3sheetResistance,
                    here->BSIM3drainSquaresGiven,
                    model->BSIM3rd,
                    model->BSIM3rdc,
                    here->BSIM3drainSquares,
                    here->BSIM3sourceSquaresGiven,
                    model->BSIM3rs,
                    model->BSIM3rsc,
                    here->BSIM3sourceSquares,
                    &drainResistance,
                    &sourceResistance
                    );
                if (error)
                    return(error);
            }

            /* process drain series resistance */
            if (drainResistance != 0.0)
            {
              if(here->BSIM3dNodePrime == 0) {
                error = CKTmkVolt(ckt,&tmp,here->BSIM3name,"drain");
                if(error) return(error);
                here->BSIM3dNodePrime = tmp->number;
                if (ckt->CKTcopyNodesets) {
                  if (CKTinst2Node(ckt,here,1,&tmpNode,&tmpName)==OK) {
                    if (tmpNode->nsGiven) {
                      tmp->nodeset=tmpNode->nodeset;
                      tmp->nsGiven=tmpNode->nsGiven;
                    }
                  }
                }
              }
            }
            else
            {   here->BSIM3dNodePrime = here->BSIM3dNode;
            }

            /* process source series resistance */
            if (sourceResistance != 0.0)
            {
              if(here->BSIM3sNodePrime == 0) {
                error = CKTmkVolt(ckt,&tmp,here->BSIM3name,"source");
                if(error) return(error);
                here->BSIM3sNodePrime = tmp->number;
                if (ckt->CKTcopyNodesets) {
                  if (CKTinst2Node(ckt,here,3,&tmpNode,&tmpName)==OK) {
                    if (tmpNode->nsGiven) {
                      tmp->nodeset=tmpNode->nodeset;
                      tmp->nsGiven=tmpNode->nsGiven;
                    }
                  }
                }
              }
            }
            else
            {   here->BSIM3sNodePrime = here->BSIM3sNode;
            }

 /* internal charge node */

            if (here->BSIM3nqsMod)
            {   if (here->BSIM3qNode == 0)
            {   error = CKTmkVolt(ckt,&tmp,here->BSIM3name,"charge");
                if(error) return(error);
                here->BSIM3qNode = tmp->number;
            }
            }
            else
            {   here->BSIM3qNode = 0;
            }

        /* set Sparse Matrix Pointers */


/* macro to make elements with built in test for out of memory */
#define TSTALLOC(ptr,first,second) \
do { if((here->ptr = SMPmakeElt(matrix, here->first, here->second)) == NULL){\
    return(E_NOMEM);\
} } while(0)

            TSTALLOC(BSIM3DdPtr, BSIM3dNode, BSIM3dNode);
            TSTALLOC(BSIM3GgPtr, BSIM3gNode, BSIM3gNode);
            TSTALLOC(BSIM3SsPtr, BSIM3sNode, BSIM3sNode);
            TSTALLOC(BSIM3BbPtr, BSIM3bNode, BSIM3bNode);
            TSTALLOC(BSIM3DPdpPtr, BSIM3dNodePrime, BSIM3dNodePrime);
            TSTALLOC(BSIM3SPspPtr, BSIM3sNodePrime, BSIM3sNodePrime);
            TSTALLOC(BSIM3DdpPtr, BSIM3dNode, BSIM3dNodePrime);
            TSTALLOC(BSIM3GbPtr, BSIM3gNode, BSIM3bNode);
            TSTALLOC(BSIM3GdpPtr, BSIM3gNode, BSIM3dNodePrime);
            TSTALLOC(BSIM3GspPtr, BSIM3gNode, BSIM3sNodePrime);
            TSTALLOC(BSIM3SspPtr, BSIM3sNode, BSIM3sNodePrime);
            TSTALLOC(BSIM3BdpPtr, BSIM3bNode, BSIM3dNodePrime);
            TSTALLOC(BSIM3BspPtr, BSIM3bNode, BSIM3sNodePrime);
            TSTALLOC(BSIM3DPspPtr, BSIM3dNodePrime, BSIM3sNodePrime);
            TSTALLOC(BSIM3DPdPtr, BSIM3dNodePrime, BSIM3dNode);
            TSTALLOC(BSIM3BgPtr, BSIM3bNode, BSIM3gNode);
            TSTALLOC(BSIM3DPgPtr, BSIM3dNodePrime, BSIM3gNode);
            TSTALLOC(BSIM3SPgPtr, BSIM3sNodePrime, BSIM3gNode);
            TSTALLOC(BSIM3SPsPtr, BSIM3sNodePrime, BSIM3sNode);
            TSTALLOC(BSIM3DPbPtr, BSIM3dNodePrime, BSIM3bNode);
            TSTALLOC(BSIM3SPbPtr, BSIM3sNodePrime, BSIM3bNode);
            TSTALLOC(BSIM3SPdpPtr, BSIM3sNodePrime, BSIM3dNodePrime);

            TSTALLOC(BSIM3QqPtr, BSIM3qNode, BSIM3qNode);

            TSTALLOC(BSIM3QdpPtr, BSIM3qNode, BSIM3dNodePrime);
            TSTALLOC(BSIM3QspPtr, BSIM3qNode, BSIM3sNodePrime);
            TSTALLOC(BSIM3QgPtr, BSIM3qNode, BSIM3gNode);
            TSTALLOC(BSIM3QbPtr, BSIM3qNode, BSIM3bNode);
            TSTALLOC(BSIM3DPqPtr, BSIM3dNodePrime, BSIM3qNode);
            TSTALLOC(BSIM3SPqPtr, BSIM3sNodePrime, BSIM3qNode);
            TSTALLOC(BSIM3GqPtr, BSIM3gNode, BSIM3qNode);
            TSTALLOC(BSIM3BqPtr, BSIM3bNode, BSIM3qNode);

        }
    }
#ifdef USE_OMP
    InstCount = 0;
    model = (BSIM3model*)inModel;
    /* loop through all the BSIM3 device models
       to count the number of instances */

    for( ; model != NULL; model = BSIM3nextModel(model))
    {
        /* loop through all the instances of the model */
        for (here = BSIM3instances(model); here != NULL ;
             here=BSIM3nextInstance(here))
        {
            InstCount++;
        }
        model->BSIM3InstCount = 0;
        model->BSIM3InstanceArray = NULL;
    }
    InstArray = TMALLOC(BSIM3instance*, InstCount);
    model = (BSIM3model*)inModel;
    /* store this in the first model only */
    model->BSIM3InstCount = InstCount;
    model->BSIM3InstanceArray = InstArray;
    idx = 0;
    for( ; model != NULL; model = BSIM3nextModel(model))
    {
        /* loop through all the instances of the model */
        for (here = BSIM3instances(model); here != NULL ;
             here=BSIM3nextInstance(here))
        {
            InstArray[idx] = here;
            idx++;
        }
    }

#endif
    return(OK);
}