/**********
Copyright 1990 Regents of the University of California.  All rights reserved.
Author: 1985 Thomas L. Quarles
**********/

#include "ngspice/ngspice.h"
#include "ngspice/smpdefs.h"
#include "ngspice/cktdefs.h"
#include "inddefs.h"
#include "ngspice/sperror.h"
#include "ngspice/suffix.h"

int
INDsetup(SMPmatrix *matrix, GENmodel *inModel, CKTcircuit *ckt, int *states)
   /* load the inductor structure with those pointers needed later 
   * for fast matrix loading 
   */
{
    INDmodel *model = (INDmodel*)inModel;
    INDinstance *here;
    int error;
    CKTnode *tmp;

    /*  loop through all the inductor models */
    for( ; model != NULL; model = INDnextModel(model)) {
 
   /* Default Value Processing for Model Parameters */
        if (!model->INDmIndGiven) {
             model->INDmInd = 0.0;
        }
	if (!model->INDtnomGiven) {
             model->INDtnom = ckt->CKTnomTemp;
        }
	if (!model->INDtc1Given) {
             model->INDtempCoeff1 = 0.0;
        }
        if (!model->INDtc2Given) {
             model->INDtempCoeff2 = 0.0;
        }
        if (!model->INDcsectGiven){
             model->INDcsect = 0.0;
        }
        if (!model->INDlengthGiven) {
             model->INDlength = 0.0;
        }
        if (!model->INDmodNtGiven) {
             model->INDmodNt = 0.0;
        }
        if (!model->INDmuGiven) {
             model->INDmu = 0.0;
        }
          
	/* precompute specific inductance (one turn) */
	if((model->INDlengthGiven) 
              && (model->INDlength > 0.0)) {
                
		if (model->INDmuGiven)
                    model->INDspecInd = (model->INDmu * CONSTmuZero 
		     * model->INDcsect * model->INDcsect) / model->INDlength;   
		else
                   model->INDspecInd = (CONSTmuZero * model->INDcsect
		   * model->INDcsect ) / model->INDlength; 
	
	} else  {
	        model->INDspecInd = 0.0;
	}	
        
	if (!model->INDmIndGiven) 
            model->INDmInd = model->INDmodNt * model->INDmodNt * model->INDspecInd;
		   
            
        
	
        /* loop through all the instances of the model */
        for (here = INDinstances(model); here != NULL ;
                here=INDnextInstance(here)) {

            here->INDflux = *states;
            *states += INDnumStates ;
            if(ckt->CKTsenInfo && (ckt->CKTsenInfo->SENmode & TRANSEN) ){
                *states += INDnumSenStates * (ckt->CKTsenInfo->SENparms);
            }

            if(here->INDbrEq == 0) {
                error = CKTmkCur(ckt,&tmp,here->INDname,"branch");
                if(error) return(error);
                here->INDbrEq = tmp->number;
            }

            here->system = NULL;
            here->system_next_ind = NULL;

/* 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(INDposIbrPtr,INDposNode,INDbrEq);
            TSTALLOC(INDnegIbrPtr,INDnegNode,INDbrEq);
            TSTALLOC(INDibrNegPtr,INDbrEq,INDnegNode);
            TSTALLOC(INDibrPosPtr,INDbrEq,INDposNode);
            TSTALLOC(INDibrIbrPtr,INDbrEq,INDbrEq);
        }
    }
    return(OK);
}

int
INDunsetup(GENmodel *inModel, CKTcircuit *ckt)
{
    INDmodel *model;
    INDinstance *here;

    for (model = (INDmodel *)inModel; model != NULL;
	    model = INDnextModel(model))
    {
        for (here = INDinstances(model); here != NULL;
                here=INDnextInstance(here))
	{
	    if (here->INDbrEq > 0)
		CKTdltNNum(ckt, here->INDbrEq);
            here->INDbrEq = 0;
	}
    }
    return OK;
}
