Switch code
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: