mesdisto.c

spice中支持多层次元件模型仿真的可单独运行的插件源码

/**********
Copyright 1990 Regents of the University of California.  All rights reserved.
Author: 1988 Jaijeet S Roychowdhury
**********/

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

int
MESdisto(mode,genmodel,ckt)
 GENmodel *genmodel;
 register CKTcircuit *ckt;
 int mode;

/* assuming here that ckt->CKTomega has been initialised to 
 * the correct value
 */
{
 MESmodel *model = (MESmodel *) genmodel;
 DISTOAN* job = (DISTOAN*) ckt->CKTcurJob;
 DpassStr pass;
 double r1h1x,i1h1x;
 double r1h1y,i1h1y;
 double r1h1z, i1h1z;
 double r1h2x, i1h2x;
 double r1h2y, i1h2y;
 double r1h2z, i1h2z;
 double r1hm2x,i1hm2x;
 double r1hm2y,i1hm2y;
 double r1hm2z, i1hm2z;
 double r2h11x,i2h11x;
 double r2h11y,i2h11y;
 double r2h11z, i2h11z;
 double r2h1m2x,i2h1m2x;
 double r2h1m2y,i2h1m2y;
 double r2h1m2z, i2h1m2z;
 double temp, itemp;
 register MESinstance *here;

if (mode == D_SETUP)
 return(MESdSetup(model,ckt));

if ((mode == D_TWOF1) || (mode == D_THRF1) || 
 (mode == D_F1PF2) || (mode == D_F1MF2) ||
 (mode == D_2F1MF2)) {

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

  /* loop through all the instances of the model */
  for (here = model->MESinstances; here != NULL ;
      here=here->MESnextInstance) {
    if (here->MESowner != ARCHme) continue;

    /* loading starts here */

    switch (mode) {
    case D_TWOF1:
      /* x = vgs, y = vds */

      /* getting first order (linear) Volterra kernel */
      r1h1x = *(job->r1H1ptr + (here->MESgateNode)) -
          *(job->r1H1ptr + (here->MESsourcePrimeNode));
      i1h1x = *(job->i1H1ptr + (here->MESgateNode)) -
          *(job->i1H1ptr + (here->MESsourcePrimeNode));

      r1h1y = *(job->r1H1ptr + (here->MESdrainPrimeNode)) -
          *(job->r1H1ptr + (here->MESsourcePrimeNode));
      i1h1y = *(job->i1H1ptr + (here->MESdrainPrimeNode)) -
          *(job->i1H1ptr + (here->MESsourcePrimeNode));


      /* loading starts here */
      /* loading cdrain term */

      temp = DFn2F1(here->cdr_x2,
      here->cdr_z2,
      0.0,
      here->cdr_xz,
      0.0,
      0.0,
      r1h1x,
      i1h1x,
      r1h1y,
      i1h1y,
      0.0,
      0.0);

      itemp = DFi2F1(here->cdr_x2,
      here->cdr_z2,
      0.0,
      here->cdr_xz,
      0.0,
      0.0,
      r1h1x,
      i1h1x,
      r1h1y,
      i1h1y,
      0.0,
      0.0);

      *(ckt->CKTrhs + (here->MESdrainPrimeNode)) -= temp;
      *(ckt->CKTirhs + (here->MESdrainPrimeNode)) -= itemp;
      *(ckt->CKTrhs + (here->MESsourcePrimeNode)) += temp;
      *(ckt->CKTirhs + (here->MESsourcePrimeNode)) += itemp;

      /* cdrain term over */

      /* loading ggs term */

      temp = D1n2F1(here->ggs2,
      r1h1x,
      i1h1x);

      itemp = D1i2F1(here->ggs2,
      r1h1x,
      i1h1x);


      *(ckt->CKTrhs + (here->MESgateNode)) -= temp;
      *(ckt->CKTirhs + (here->MESgateNode)) -= itemp;
      *(ckt->CKTrhs + (here->MESsourcePrimeNode)) += temp;
      *(ckt->CKTirhs + (here->MESsourcePrimeNode)) += itemp;

      /* ggs over */

      /* loading ggd term */

      temp = D1n2F1(here->ggd2,
      r1h1x - r1h1y,
      i1h1x - i1h1y);

      itemp = D1i2F1(here->ggd2,
      r1h1x - r1h1y,
      i1h1x - i1h1y);


      *(ckt->CKTrhs + (here->MESgateNode)) -= temp;
      *(ckt->CKTirhs + (here->MESgateNode)) -= itemp;
      *(ckt->CKTrhs + (here->MESdrainPrimeNode)) += temp;
      *(ckt->CKTirhs + (here->MESdrainPrimeNode)) += itemp;

      /* ggd over */

      /* loading qgs term */

      temp = -ckt->CKTomega *
          DFi2F1(here->qgs_x2,
      here->qgs_y2,
      0.0,
      here->qgs_xy,
      0.0,
      0.0,
      r1h1x,
      i1h1x,
      r1h1x - r1h1y,
      i1h1x - i1h1y,
      0.0,
      0.0);

      itemp = ckt->CKTomega*
          DFn2F1(here->qgs_x2,
      here->qgs_y2,
      0.0,
      here->qgs_xy,
      0.0,
      0.0,
      r1h1x,
      i1h1x,
      r1h1x - r1h1y,
      i1h1x - i1h1y,
      0.0,
      0.0);

      *(ckt->CKTrhs + (here->MESgateNode)) -= temp;
      *(ckt->CKTirhs + (here->MESgateNode)) -= itemp;
      *(ckt->CKTrhs + (here->MESsourcePrimeNode)) += temp;
      *(ckt->CKTirhs + (here->MESsourcePrimeNode)) += itemp;

      /* qgs term over */

      /* loading qgd term */

      temp = -ckt->CKTomega *
          DFi2F1(here->qgd_x2,
      here->qgd_y2,
      0.0,
      here->qgd_xy,
      0.0,
      0.0,
      r1h1x,
      i1h1x,
      r1h1x - r1h1y,
      i1h1x - i1h1y,
      0.0,
      0.0);

      itemp = ckt->CKTomega*
          DFn2F1(here->qgd_x2,
      here->qgd_y2,
      0.0,
      here->qgd_xy,
      0.0,
      0.0,
      r1h1x,
      i1h1x,
      r1h1x - r1h1y,
      i1h1x - i1h1y,
      0.0,
      0.0);

      *(ckt->CKTrhs + (here->MESgateNode)) -= temp;
      *(ckt->CKTirhs + (here->MESgateNode)) -= itemp;
      *(ckt->CKTrhs + (here->MESdrainPrimeNode)) += temp;
      *(ckt->CKTirhs + (here->MESdrainPrimeNode)) += itemp;

      /* qgd term over */
      /* all done */

      break;

    case D_THRF1:
      /* x = vgs, y = vds */

      /* getting first order (linear) Volterra kernel */
      r1h1x = *(job->r1H1ptr + (here->MESgateNode)) -
          *(job->r1H1ptr + (here->MESsourcePrimeNode));
      i1h1x = *(job->i1H1ptr + (here->MESgateNode)) -
          *(job->i1H1ptr + (here->MESsourcePrimeNode));

      r1h1y = *(job->r1H1ptr + (here->MESdrainPrimeNode)) -
          *(job->r1H1ptr + (here->MESsourcePrimeNode));
      i1h1y = *(job->i1H1ptr + (here->MESdrainPrimeNode)) -
          *(job->i1H1ptr + (here->MESsourcePrimeNode));

      r2h11x = *(job->r2H11ptr + (here->MESgateNode)) -
          *(job->r2H11ptr + (here->MESsourcePrimeNode));
      i2h11x = *(job->i2H11ptr + (here->MESgateNode)) -
          *(job->i2H11ptr + (here->MESsourcePrimeNode));

      r2h11y = *(job->r2H11ptr + (here->MESdrainPrimeNode)) -
          *(job->r2H11ptr + (here->MESsourcePrimeNode));
      i2h11y = *(job->i2H11ptr + (here->MESdrainPrimeNode)) -
          *(job->i2H11ptr + (here->MESsourcePrimeNode));

      /* loading starts here */
      /* loading cdrain term */

      temp = DFn3F1(here->cdr_x2,
      here->cdr_z2,
      0.0,
      here->cdr_xz,
      0.0,
      0.0,
      here->cdr_x3,
      here->cdr_z3,
      0.0,
      here->cdr_x2z,
      0.0,
      here->cdr_xz2,
      0.0,
      0.0,
      0.0,
      0.0,
      r1h1x,
      i1h1x,
      r1h1y,
      i1h1y,
      0.0,
      0.0,
      r2h11x,
      i2h11x,
      r2h11y,
      i2h11y,
      0.0,
      0.0);

      itemp = DFi3F1(here->cdr_x2,
      here->cdr_z2,
      0.0,
      here->cdr_xz,
      0.0,
      0.0,
      here->cdr_x3,
      here->cdr_z3,
      0.0,
      here->cdr_x2z,
      0.0,
      here->cdr_xz2,
      0.0,
      0.0,
      0.0,
      0.0,
      r1h1x,
      i1h1x,
      r1h1y,
      i1h1y,
      0.0,
      0.0,
      r2h11x,
      i2h11x,
      r2h11y,
      i2h11y,
      0.0,
      0.0);

      *(ckt->CKTrhs + (here->MESdrainPrimeNode)) -= temp;
      *(ckt->CKTirhs + (here->MESdrainPrimeNode)) -= itemp;
      *(ckt->CKTrhs + (here->MESsourcePrimeNode)) += temp;
      *(ckt->CKTirhs + (here->MESsourcePrimeNode)) += itemp;

      /* cdrain term over */

      /* loading ggs term */

      temp = D1n3F1(here->ggs2,
      here->ggs3,
      r1h1x,
      i1h1x,
      r2h11x,
      i2h11x);

      itemp = D1i3F1(here->ggs2,
      here->ggs3,
      r1h1x,
      i1h1x,
      r2h11x,
      i2h11x);


      *(ckt->CKTrhs + (here->MESgateNode)) -= temp;
      *(ckt->CKTirhs + (here->MESgateNode)) -= itemp;
      *(ckt->CKTrhs + (here->MESsourcePrimeNode)) += temp;
      *(ckt->CKTirhs + (here->MESsourcePrimeNode)) += itemp;

      /* ggs over */

      /* loading ggd term */

      temp = D1n3F1(here->ggd2,
      here->ggd3,
      r1h1x - r1h1y,
      i1h1x - i1h1y,
      r2h11x - r2h11y,
      i2h11x - i2h11y);

      itemp = D1i3F1(here->ggd2,
      here->ggd3,
      r1h1x - r1h1y,
      i1h1x - i1h1y,
      r2h11x - r2h11y,
      i2h11x - i2h11y);



      *(ckt->CKTrhs + (here->MESgateNode)) -= temp;
      *(ckt->CKTirhs + (here->MESgateNode)) -= itemp;
      *(ckt->CKTrhs + (here->MESdrainPrimeNode)) += temp;
      *(ckt->CKTirhs + (here->MESdrainPrimeNode)) += itemp;

      /* ggd over */

      /* loading qgs term */

      temp = -ckt->CKTomega*
          DFi3F1(here->qgs_x2,
      here->qgs_y2,
      0.0,
      here->qgs_xy,
      0.0,
      0.0,
      here->qgs_x3,
      here->qgs_y3,
      0.0,
      here->qgs_x2y,
      0.0,
      here->qgs_xy2,
      0.0,
      0.0,
      0.0,
      0.0,
      r1h1x,
      i1h1x,
      r1h1x - r1h1y,
      i1h1x - i1h1y,
      0.0,
      0.0,
      r2h11x,
      i2h11x,
      r2h11x - r2h11y,
      i2h11x - i2h11y,
      0.0,
      0.0);

      itemp =ckt->CKTomega*
          DFn3F1(here->qgs_x2,
      here->qgs_y2,
      0.0,
      here->qgs_xy,
      0.0,
      0.0,
      here->qgs_x3,
      here->qgs_y3,
      0.0,
      here->qgs_x2y,
      0.0,
      here->qgs_xy2,
      0.0,
      0.0,
      0.0,
      0.0,
      r1h1x,
      i1h1x,
      r1h1x - r1h1y,
      i1h1x - i1h1y,
      0.0,
      0.0,
      r2h11x,
      i2h11x,
      r2h11x - r2h11y,
      i2h11x - i2h11y,
      0.0,
      0.0);

      *(ckt->CKTrhs + (here->MESgateNode)) -= temp;
      *(ckt->CKTirhs + (here->MESgateNode)) -= itemp;
      *(ckt->CKTrhs + (here->MESsourcePrimeNode)) += temp;
      *(ckt->CKTirhs + (here->MESsourcePrimeNode)) += itemp;

      /* qgs term over */

      /* loading qgd term */

      temp = -ckt->CKTomega*
          DFi3F1(here->qgd_x2,
      here->qgd_y2,
      0.0,
      here->qgd_xy,
      0.0,
      0.0,
      here->qgd_x3,
      here->qgd_y3,
      0.0,
      here->qgd_x2y,
      0.0,
      here->qgd_xy2,
      0.0,
      0.0,
      0.0,
      0.0,
      r1h1x,
      i1h1x,
      r1h1x - r1h1y,
      i1h1x - i1h1y,
      0.0,
      0.0,
      r2h11x,
      i2h11x,
      r2h11x - r2h11y,
      i2h11x - i2h11y,
      0.0,
      0.0);

      itemp =ckt->CKTomega*
          DFn3F1(here->qgd_x2,
      here->qgd_y2,
      0.0,
      here->qgd_xy,
      0.0,
      0.0,
      here->qgd_x3,
      here->qgd_y3,
      0.0,
      here->qgd_x2y,
      0.0,
      here->qgd_xy2,
      0.0,
      0.0,
      0.0,
      0.0,
      r1h1x,
      i1h1x,
      r1h1x - r1h1y,
      i1h1x - i1h1y,
      0.0,
      0.0,
      r2h11x,
      i2h11x,
      r2h11x - r2h11y,
      i2h11x - i2h11y,
      0.0,
      0.0);

      *(ckt->CKTrhs + (here->MESgateNode)) -= temp;
      *(ckt->CKTirhs + (here->MESgateNode)) -= itemp;
      *(ckt->CKTrhs + (here->MESdrainPrimeNode)) += temp;
      *(ckt->CKTirhs + (here->MESdrainPrimeNode)) += itemp;

      /* qgd term over */
      /* all done */

      break;
    case D_F1PF2:
      /* x = vgs, y = vds */

      /* getting first order (linear) Volterra transform */
      r1h1x = *(job->r1H1ptr + (here->MESgateNode)) -
          *(job->r1H1ptr + (here->MESsourcePrimeNode));
      i1h1x = *(job->i1H1ptr + (here->MESgateNode)) -
          *(job->i1H1ptr + (here->MESsourcePrimeNode));

      r1h1y = *(job->r1H1ptr + (here->MESdrainPrimeNode)) -
          *(job->r1H1ptr + (here->MESsourcePrimeNode));
      i1h1y = *(job->i1H1ptr + (here->MESdrainPrimeNode)) -
          *(job->i1H1ptr + (here->MESsourcePrimeNode));

      r1h2x = *(job->r1H2ptr + (here->MESgateNode)) -
          *(job->r1H2ptr + (here->MESsourcePrimeNode));
      i1h2x = *(job->i1H2ptr + (here->MESgateNode)) -
          *(job->i1H2ptr + (here->MESsourcePrimeNode));

      r1h2y = *(job->r1H2ptr + (here->MESdrainPrimeNode)) -
          *(job->r1H2ptr + (here->MESsourcePrimeNode));
      i1h2y = *(job->i1H2ptr + (here->MESdrainPrimeNode)) -
          *(job->i1H2ptr + (here->MESsourcePrimeNode));

      /* loading starts here */
      /* loading cdrain term */

      temp = DFnF12(here->cdr_x2,
      here->cdr_z2,
      0.0,
      here->cdr_xz,
      0.0,
      0.0,
      r1h1x,
      i1h1x,
      r1h1y,
      i1h1y,
      0.0,
      0.0,
      r1h2x,
      i1h2x,
      r1h2y,
      i1h2y,
      0.0,
      0.0);

      itemp = DFiF12(here->cdr_x2,
      here->cdr_z2,
      0.0,
      here->cdr_xz,
      0.0,
      0.0,
      r1h1x,
      i1h1x,
      r1h1y,
      i1h1y,
      0.0,
      0.0,
      r1h2x,
      i1h2x,
      r1h2y,
      i1h2y,
      0.0,
      0.0);

      *(ckt->CKTrhs + (here->MESdrainPrimeNode)) -= temp;
      *(ckt->CKTirhs + (here->MESdrainPrimeNode)) -= itemp;
      *(ckt->CKTrhs + (here->MESsourcePrimeNode)) += temp;
      *(ckt->CKTirhs + (here->MESsourcePrimeNode)) += itemp;

      /* cdrain term over */

      /* loading ggs term */

      temp = D1nF12(here->ggs2,
      r1h1x,
      i1h1x,
      r1h2x,
      i1h2x);

      itemp = D1iF12(here->ggs2,
      r1h1x,
      i1h1x,
      r1h2x,
      i1h2x);


      *(ckt->CKTrhs + (here->MESgateNode)) -= temp;
      *(ckt->CKTirhs + (here->MESgateNode)) -= itemp;
      *(ckt->CKTrhs + (here->MESsourcePrimeNode)) += temp;
      *(ckt->CKTirhs + (here->MESsourcePrimeNode)) += itemp;

      /* ggs over */

      /* loading ggd term */

      temp = D1nF12(here->ggd2,
      r1h1x - r1h1y,
      i1h1x - i1h1y,
      r1h2x - r1h2y,
      i1h2x - i1h2y);

      itemp = D1iF12(here->ggd2,