Switch code


To ng-spice@ieee.ing.uniroma1.it
From Herbert Widlok <widlok@uci.agh.edu.pl>
Date Mon, 30 Oct 2000 11:04:11 +0100
Delivered-To mailing list ng-spice@ieee.ing.uniroma1.it
Mailing-List contact ng-spice-help@ieee.ing.uniroma1.it; run by ezmlm
References <JOJMALEAILMIEAAA@hotbot.com > <00102720300901.01174@spiff.al.dynip.com >
Reply-To ng-spice@ieee.ing.uniroma1.it

Hi,

    These are my switches, please tell me what do You deside to do with
them. If You accept this, I can update manuals too.

We have some bug in "let" command I think. I spend last night lookin for
it without any results. I __suppose__ that something
overides kewords variable and then (form time to time) clookup function
gives a memory fault. If anyone knows anythink please tell me.

Michael
/**********
Copyright 1990 Regents of the University of California.  All rights reserved.
Author: 1985 Gordon Jacobs
**********/
/* MW. Current switch changed like voltage switch */

#include "ngspice.h"
#include <stdio.h>
#include "cktdefs.h"
#include "cswdefs.h"
#include "trandefs.h"
#include "sperror.h"
#include "suffix.h"

int
CSWload(inModel,ckt)
    GENmodel *inModel;
    CKTcircuit *ckt;

        /* actually load the current values into the 
         * sparse matrix previously provided 
         */
{
    CSWmodel *model = (CSWmodel*)inModel;
    CSWinstance *here;
    double g_now;
    double i_ctrl;
    double previous_state; 
    double current_state;

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

        /* loop through all the instances of the model */
        for (here = model->CSWinstances; here != NULL ;
                here=here->CSWnextInstance) {
            if (here->CSWowner != ARCHme) continue;

            /* decide the state of the switch */

            if(ckt->CKTmode & (MODEINITFIX|MODEINITJCT)) {

                if(here->CSWzero_stateGiven) {
                        /* switch specified "on" */
                    *(ckt->CKTstate0 + here->CSWstate) = 1.0;
                    current_state = 1.0;
                } else {
                    *(ckt->CKTstate0 + here->CSWstate) = 0.0;
                    current_state = 0.0;
                }
            } else if (ckt->CKTmode & (MODEINITSMSIG)) {

                previous_state = *(ckt->CKTstate0 + here->CSWstate);
                current_state = previous_state;

            } else if (ckt->CKTmode & (MODEINITFLOAT)) {
                /* use state0 since INITTRAN or INITPRED already called */
                previous_state = *(ckt->CKTstate0 + here->CSWstate);
                i_ctrl = *(ckt->CKTrhsOld + 
                        here->CSWcontBranch);
                if(i_ctrl > (model->CSWiThreshold+model->CSWiHysteresis)) {
                    *(ckt->CKTstate0 + here->CSWstate) = 1.0;
                    current_state = 1.0;
                }
                else if(i_ctrl < (model->CSWiThreshold - 
                        model->CSWiHysteresis)) {
                    *(ckt->CKTstate0 + here->CSWstate) = 0.0;
                    current_state = 0.0;
                }
                else {

/* MW. New Switch state 'middle' 0-Off, 1-On, and lineary between */
                        current_state = ((i_ctrl-model->CSWiThreshold) / 
                                                        
(2*model->CSWiHysteresis)) + 0.5;
                        *(ckt->CKTstate0 + here->CSWstate) = current_state;
                }

/* MW. Convergence check - use reltol only */
                if (fabs(current_state-previous_state) > ckt->CKTreltol) {
                    ckt->CKTnoncon++;       /* ensure one more iteration */
                                ckt->CKTtroubleElt = (GENinstance *) here;
                        }

            } else if (ckt->CKTmode & (MODEINITTRAN|MODEINITPRED)) {

                previous_state = *(ckt->CKTstate1 + here->CSWstate);
                i_ctrl = *(ckt->CKTrhsOld + 
                        here->CSWcontBranch);

                if(i_ctrl > (model->CSWiThreshold+model->CSWiHysteresis)) {
                    current_state = 1;
                } else if(i_ctrl < (model->CSWiThreshold - 
                        model->CSWiHysteresis))  {
                    current_state = 0;
                } else {

/* MW. New Switch state 'middle' 0-Off, 1-On, and lineary between */
                        current_state = ((i_ctrl-model->CSWiThreshold) / 
                                                        
(2*model->CSWiHysteresis)) + 0.5;
                        *(ckt->CKTstate0 + here->CSWstate) = current_state;
                }

                *(ckt->CKTstate0 + here->CSWstate) = current_state;

            }

/* MW. Calculate conductance */
            g_now = model->CSWoffConduct + (current_state *
                                (model->CSWonConduct - model->CSWoffConduct));
            here->CSWcond = g_now;

            *(here->CSWposPosptr) += g_now;
            *(here->CSWposNegptr) -= g_now;
            *(here->CSWnegPosptr) -= g_now;
            *(here->CSWnegNegptr) += g_now;
        }
    }
    return(OK);
}
/**********
Copyright 1990 Regents of the University of California.  All rights reserved.
Author: 1985 Gordon Jacobs
**********/
/* MW. Switch changed - it switches lineary */

#include "ngspice.h"
#include <stdio.h>
#include "cktdefs.h"
#include "swdefs.h"
#include "trandefs.h"
#include "sperror.h"
#include "suffix.h"

int
SWload(inModel,ckt)
    GENmodel *inModel;
    CKTcircuit *ckt;
        /* actually load the current values into the 
         * sparse matrix previously provided 
         */
{
    SWmodel *model = (SWmodel *) inModel;
    SWinstance *here;
    double g_now;
    double v_ctrl;
    double previous_state; 
    double current_state;

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

        /* loop through all the instances of the model */
        for (here = model->SWinstances; here != NULL ;
                here=here->SWnextInstance) {
            if (here->SWowner != ARCHme) continue;

            /* decide the state of the switch */

            if(ckt->CKTmode & (MODEINITFIX|MODEINITJCT)) {

                if(here->SWzero_stateGiven) {
                        /* switch specified "on" */
                    *(ckt->CKTstate0 + here->SWstate) = 1.0;
                    current_state = 1.0;
                } else {
                    *(ckt->CKTstate0 + here->SWstate) = 0.0;
                    current_state = 0.0;
                }

            } else if (ckt->CKTmode & (MODEINITSMSIG)) {

                previous_state = *(ckt->CKTstate0 + here->SWstate);
                current_state = previous_state;

            } else if (ckt->CKTmode & (MODEINITFLOAT)) {

                /* use state0 since INITTRAN or INITPRED already called */
                previous_state = *(ckt->CKTstate0 + here->SWstate);
                v_ctrl = *(ckt->CKTrhsOld + here->SWposCntrlNode)
                         - *(ckt->CKTrhsOld + here->SWnegCntrlNode);
                if(v_ctrl > (model->SWvThreshold + model->SWvHysteresis)) {
                    *(ckt->CKTstate0 + here->SWstate) = 1.0;
                    current_state = 1.0;
                } else if(v_ctrl < (model->SWvThreshold - 
                        model->SWvHysteresis)) {
                    *(ckt->CKTstate0 + here->SWstate) = 0.0;
                    current_state = 0.0;
                } else {
                   
/* MW. New Switch state 'middle' 0-Off, 1-On, and lineary between */
                        current_state = ((v_ctrl-model->SWvThreshold) / 
                                                        
(2*model->SWvHysteresis)) + 0.5;
                        *(ckt->CKTstate0 + here->SWstate) = current_state;
                }
                
/* MW. Convergence check - use reltol only */
                if (fabs(current_state-previous_state) > ckt->CKTreltol) {
                    ckt->CKTnoncon++;       /* ensure one more iteration */
                                ckt->CKTtroubleElt = (GENinstance *) here;
                }


            } else if(ckt->CKTmode & (MODEINITTRAN|MODEINITPRED) ) {

                previous_state = *(ckt->CKTstate1 + here->SWstate);
                v_ctrl = *(ckt->CKTrhsOld + here->SWposCntrlNode)
                         - *(ckt->CKTrhsOld + here->SWnegCntrlNode);

                if(v_ctrl > (model->SWvThreshold + model->SWvHysteresis))  {
                    current_state = 1.0;
                } else if(v_ctrl < (model->SWvThreshold - 
                        model->SWvHysteresis)) {
                    current_state = 0.0;
                } else {

/* MW. New Switch state 'middle' 0-Off, 1-On, and lineary between */
                        current_state = ((v_ctrl-model->SWvThreshold) / 
                                                        
(2*model->SWvHysteresis)) + 0.5;
                }

                    *(ckt->CKTstate0 + here->SWstate) = current_state;

            }


/* MW. Calculate conductance */
            g_now = model->SWoffConduct + (current_state *
                                (model->SWonConduct - model->SWoffConduct));
            here->SWcond = g_now;

            *(here->SWposPosptr) += g_now;
            *(here->SWposNegptr) -= g_now;
            *(here->SWnegPosptr) -= g_now;
            *(here->SWnegNegptr) += g_now;
        }
    }
    return(OK);
}

Partial thread listing: