restemp.c File Reference

#include "spice.h"
#include <stdio.h>
#include "resdefs.h"
#include "sperror.h"
#include "util.h"

Include dependency graph for restemp.c:

Go to the source code of this file.

Functions
int	REStemp (GENmodel *inModel, CKTcircuit *ckt)

Function Documentation

int REStemp	(	GENmodel *	inModel,
		CKTcircuit *	ckt
	)

Definition at line 16 of file restemp.c.

{
    register RESmodel *model =  (RESmodel *)inModel;
    register RESinstance *here;
    double factor;
    double difference;


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

        /* Default Value Processing for Resistor Models */
        if(!model->REStnomGiven) model->REStnom = ckt->CKTnomTemp;
        if(!model->RESsheetResGiven) model->RESsheetRes = 0;
        if(!model->RESdefWidthGiven) model->RESdefWidth = 10.e-6; /*M*/
        if(!model->REStc1Given) model->REStempCoeff1 = 0;
        if(!model->REStc2Given) model->REStempCoeff2 = 0;
        if(!model->RESnarrowGiven) model->RESnarrow = 0;

        /* loop through all the instances of the model */
        for (here = model->RESinstances; here != NULL ;
                here=here->RESnextInstance) {
            
            /* Default Value Processing for Resistor Instance */
            if(!here->REStempGiven) here->REStemp = ckt->CKTtemp;
            if(!here->RESwidthGiven) here->RESwidth = model->RESdefWidth;
            if(!here->RESlengthGiven)  here->RESlength = 0 ;
            if(!here->RESresGiven)  {
                if(model->RESsheetResGiven && (model->RESsheetRes != 0) &&
                        (here->RESlength != 0)) {
                    here->RESresist = model->RESsheetRes * (here->RESlength -
                        model->RESnarrow) / (here->RESwidth - model->RESnarrow);
                } else {
                    (*(SPfrontEnd->IFerror))(ERR_WARNING,
                            "%s: resistance=0, set to 1000",&(here->RESname));
                    here->RESresist=1000;
                }
            }

            difference = here->REStemp - model->REStnom;
            factor = 1.0 + (model->REStempCoeff1)*difference + 
                    (model->REStempCoeff2)*difference*difference;

            here->RESconduct = 1.0/(here->RESresist * factor);

            if (ckt->CKTpreload) {
                /* preload constants */

                if (!here->RESnegNode)
                    *(here->RESposPosptr) += here->RESconduct;
                else if (!here->RESposNode)
                    *(here->RESnegNegptr) += here->RESconduct;
                else {
                    *(here->RESposPosptr) += here->RESconduct;
                    *(here->RESnegNegptr) += here->RESconduct;
                    *(here->RESposNegptr) -= here->RESconduct;
                    *(here->RESnegPosptr) -= here->RESconduct;
                }
            }
        }
    }
    return (OK);
}