asrcacld.c

linux平台下类似著名的电路板作图软件 Spice的源代码

/**********
Copyright 1990 Regents of the University of California.  All rights reserved.
Author: 1988 Kanwar Jit Singh
**********/
/*
 * [e-mail]
 */

#include "spice.h"
#include <stdio.h>
#include "cktdefs.h"
#include "asrcdefs.h"
#include "sperror.h"
#include "util.h"
#include "suffix.h"


/*ARGSUSED*/
int
ASRCacLoad(inModel,ckt)
GENmodel *inModel;
CKTcircuit *ckt;
{

    /*
     * Actually load the current voltage value into the 
     * sparse matrix previously provided. The values have
     * been precomputed and stored with the instance model.
     */

    register ASRCmodel *model = (ASRCmodel*)inModel;
    register ASRCinstance *here;
    int i, v_first, j;
    double *derivs;
    double rhs;

    /*  loop through all the Arbitrary source models */
    for( ; model != NULL; model = model->ASRCnextModel ) {

        /* loop through all the instances of the model */
        for (here = model->ASRCinstances; here != NULL ;
                here = here->ASRCnextInstance) {
            
        /*
         * Get the function and its derivatives from the
         * field in the instance structure. The field is 
         * an array of doubles holding the rhs, and the
         * entries of the jacobian.
         */

        v_first = 1;
        j=0;
        derivs = here->ASRCacValues;
        rhs = (here->ASRCacValues)[here->ASRCtree->numVars];

        for(i=0; i < here->ASRCtree->numVars; i++){
        switch(here->ASRCtree->varTypes[i]){
        case IF_INSTANCE:
            if( here->ASRCtype == ASRC_VOLTAGE){
            /* CCVS */
            if(v_first){
                *(here->ASRCposptr[j++]) += 1.0;
                *(here->ASRCposptr[j++]) -= 1.0;
                *(here->ASRCposptr[j++]) -= 1.0;
                *(here->ASRCposptr[j++]) += 1.0;
                v_first = 0;
               }
            *(here->ASRCposptr[j++]) -= derivs[i];
            } else{
            /* CCCS */
            *(here->ASRCposptr[j++]) += derivs[i];
            *(here->ASRCposptr[j++]) -= derivs[i];
            }
            break;
        case IF_NODE:
            if(here->ASRCtype == ASRC_VOLTAGE){
            /* VCVS */
            if( v_first){
                *(here->ASRCposptr[j++]) += 1.0;
                *(here->ASRCposptr[j++]) -= 1.0;
                *(here->ASRCposptr[j++]) -= 1.0;
                *(here->ASRCposptr[j++]) += 1.0;
                v_first = 0;
            }
            *(here->ASRCposptr[j++]) -= derivs[i];
            } else {
            /*VCCS*/
            *(here->ASRCposptr[j++]) += derivs[i];
            *(here->ASRCposptr[j++]) -= derivs[i];
            }
            break;
        default:
            return(E_BADPARM);
        }
        }
#ifdef notdef
        /* Insert the RHS */
        if( here->ASRCtype == ASRC_VOLTAGE){
        *(ckt->CKTrhs+(here->ASRCbranch)) += rhs;
        } else {
        *(ckt->CKTrhs+(here->ASRCposNode)) -= rhs;
        *(ckt->CKTrhs+(here->ASRCnegNode)) += rhs;
        }
#endif
    }
    }
    return(OK);
}