cswtrunc.c

ngspice又一个电子CAD仿真软件代码.功能更全

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

#include "ngspice.h"
#include "cktdefs.h"
#include "sperror.h"
#include "suffix.h"

#include "cswdefs.h"

int
CSWtrunc(GENmodel *inModel, CKTcircuit *ckt, double *timeStep)
{
    CSWmodel *model = (CSWmodel*)inModel;
    CSWinstance *here;

    double   lastChange, maxChange, maxStep, ref;

    for( ; model!= NULL; model = model->CSWnextModel) {
        for(here = model->CSWinstances ; here != NULL ;
                here = here->CSWnextInstance) {
            lastChange = *(ckt->CKTstate0+(here->CSWstate+1)) -
                          *(ckt->CKTstate1+(here->CSWstate+1));
            if (*(ckt->CKTstate0+(here->CSWstate))==0) {
              ref = (model->CSWiThreshold + model->CSWiHysteresis);
              if ((*(ckt->CKTstate0+(here->CSWstate+1))<ref) && (lastChange>0)) {
                   maxChange = (ref - *(ckt->CKTstate0+(here->CSWstate+1))) *
                                0.75 + 0.00005;
                   maxStep = maxChange/lastChange * ckt->CKTdeltaOld[0];
                   if (*timeStep > maxStep) { *timeStep = maxStep; }
              }
            } else {
              ref = (model->CSWiThreshold - model->CSWiHysteresis);
              if ((*(ckt->CKTstate0+(here->CSWstate+1))>ref) && (lastChange<0)) {
                   maxChange = (ref - *(ckt->CKTstate0+(here->CSWstate+1))) *
                                0.75 - 0.00005;
                   maxStep = maxChange/lastChange * ckt->CKTdeltaOld[0];
                   if (*timeStep > maxStep) { *timeStep = maxStep; }
              }
            }
        }
    }
    return(OK);
}