mesdisto.c

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

      r1h1x - r1h1y,
      i1h1x - i1h1y,
      r1h2x - r1h2y,
      i1h2x - i1h2y);


      *(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*
          DFiF12(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,
      r1h2x,
      i1h2x,
      r1h2x - r1h2y,
      i1h2x - i1h2y,
      0.0,
      0.0);

      itemp =ckt->CKTomega*
          DFnF12(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,
      r1h2x,
      i1h2x,
      r1h2x - r1h2y,
      i1h2x - i1h2y,
      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*
          DFiF12(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,
      r1h2x,
      i1h2x,
      r1h2x - r1h2y,
      i1h2x - i1h2y,
      0.0,
      0.0);

      itemp =ckt->CKTomega*
          DFnF12(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,
      r1h2x,
      i1h2x,
      r1h2x - r1h2y,
      i1h2x - i1h2y,
      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_F1MF2:
      /* 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));

      r1hm2x = *(job->r1H2ptr + (here->MESgateNode)) -
          *(job->r1H2ptr + (here->MESsourcePrimeNode));
      i1hm2x = -(*(job->i1H2ptr + (here->MESgateNode)) -
          *(job->i1H2ptr + (here->MESsourcePrimeNode)));

      r1hm2y = *(job->r1H2ptr + (here->MESdrainPrimeNode)) -
          *(job->r1H2ptr + (here->MESsourcePrimeNode));
      i1hm2y = -(*(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,
      r1hm2x,
      i1hm2x,
      r1hm2y,
      i1hm2y,
      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,
      r1hm2x,
      i1hm2x,
      r1hm2y,
      i1hm2y,
      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,
      r1hm2x,
      i1hm2x);

      itemp = D1iF12(here->ggs2,
      r1h1x,
      i1h1x,
      r1hm2x,
      i1hm2x);


      *(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,
      r1hm2x - r1hm2y,
      i1hm2x - i1hm2y);

      itemp = D1iF12(here->ggd2,
      r1h1x - r1h1y,
      i1h1x - i1h1y,
      r1hm2x - r1hm2y,
      i1hm2x - i1hm2y);


      *(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*
          DFiF12(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,
      r1hm2x,
      i1hm2x,
      r1hm2x - r1hm2y,
      i1hm2x - i1hm2y,
      0.0,
      0.0);

      itemp = ckt->CKTomega*
          DFnF12(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,
      r1hm2x,
      i1hm2x,
      r1hm2x - r1hm2y,
      i1hm2x - i1hm2y,
      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*
          DFiF12(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,
      r1hm2x,
      i1hm2x,
      r1hm2x - r1hm2y,
      i1hm2x - i1hm2y,
      0.0,
      0.0);

      itemp = ckt->CKTomega*
          DFnF12(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,
      r1hm2x,
      i1hm2x,
      r1hm2x - r1hm2y,
      i1hm2x - i1hm2y,
      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_2F1MF2:
      /* 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));

      r1hm2x = *(job->r1H2ptr + (here->MESgateNode)) -
          *(job->r1H2ptr + (here->MESsourcePrimeNode));
      i1hm2x = -(*(job->i1H2ptr + (here->MESgateNode)) -
          *(job->i1H2ptr + (here->MESsourcePrimeNode)));

      r1hm2y = *(job->r1H2ptr + (here->MESdrainPrimeNode)) -
          *(job->r1H2ptr + (here->MESsourcePrimeNode));
      i1hm2y = -(*(job->i1H2ptr + (here->MESdrainPrimeNode)) -
          *(job->i1H2ptr + (here->MESsourcePrimeNode)));

      r2h1m2x = *(job->r2H1m2ptr + (here->MESgateNode)) -
          *(job->r2H1m2ptr + (here->MESsourcePrimeNode));
      i2h1m2x = *(job->i2H1m2ptr + (here->MESgateNode)) -
          *(job->i2H1m2ptr + (here->MESsourcePrimeNode));

      r2h1m2y = *(job->r2H1m2ptr + (here->MESdrainPrimeNode)) -
          *(job->r2H1m2ptr + (here->MESsourcePrimeNode));
      i2h1m2y = *(job->i2H1m2ptr + (here->MESdrainPrimeNode))
        - *(job->i2H1m2ptr + (here->MESsourcePrimeNode));

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

      pass.cxx = here->cdr_x2;
      pass.cyy = here->cdr_z2;
      pass.czz = 0.0;
      pass.cxy = here->cdr_xz;
      pass.cyz = 0.0;
      pass.cxz = 0.0;
      pass.cxxx = here->cdr_x3;
      pass.cyyy = here->cdr_z3;
      pass.czzz = 0.0;
      pass.cxxy = here->cdr_x2z;
      pass.cxxz = 0.0;
      pass.cxyy = here->cdr_xz2;
      pass.cyyz = 0.0;
      pass.cxzz = 0.0;
      pass.cyzz = 0.0;
      pass.cxyz = 0.0;
      pass.r1h1x = r1h1x;
      pass.i1h1x = i1h1x;
      pass.r1h1y = r1h1y;
      pass.i1h1y = i1h1y;
      pass.r1h1z = 0.0;
      pass.i1h1z = 0.0;
      pass.r1h2x = r1hm2x;
      pass.i1h2x = i1hm2x;
      pass.r1h2y = r1hm2y;
      pass.i1h2y = i1hm2y;
      pass.r1h2z = 0.0;
      pass.i1h2z = 0.0;
      pass.r2h11x = r2h11x;
      pass.i2h11x = i2h11x;
      pass.r2h11y = r2h11y;
      pass.i2h11y = i2h11y;
      pass.r2h11z = 0.0;
      pass.i2h11z = 0.0;
      pass.h2f1f2x = r2h1m2x;
      pass.ih2f1f2x = i2h1m2x;
      pass.h2f1f2y = r2h1m2y;
      pass.ih2f1f2y = i2h1m2y;
      pass.h2f1f2z = 0.0;
      pass.ih2f1f2z = 0.0;
      temp = DFn2F12(&pass);

      itemp = DFi2F12(&pass);

      *(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 = D1n2F12(here->ggs2,
      here->ggs3,
      r1h1x,
      i1h1x,
      r1hm2x,
      i1hm2x,
      r2h11x,
      i2h11x,
      r2h1m2x,
      i2h1m2x);

      itemp = D1i2F12(here->ggs2,
      here->ggs3,
      r1h1x,
      i1h1x,
      r1hm2x,
      i1hm2x,
      r2h11x,
      i2h11x,
      r2h1m2x,
      i2h1m2x);


      *(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 = D1n2F12(here->ggd2,
      here->ggd3,
      r1h1x - r1h1y,
      i1h1x - i1h1y,
      r1hm2x - r1hm2y,
      i1hm2x - i1hm2y,
      r2h11x - r2h11y,
      i2h11x - i2h11y,
      r2h1m2x - r2h1m2y,
      i2h1m2x - i2h1m2y);

      itemp = D1i2F12(here->ggd2,
      here->ggd3,
      r1h1x - r1h1y,
      i1h1x - i1h1y,
      r1hm2x - r1hm2y,
      i1hm2x - i1hm2y,
      r2h11x - r2h11y,
      i2h11x - i2h11y,
      r2h1m2x - r2h1m2y,
      i2h1m2x - i2h1m2y);



      *(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 */

      pass.cxx = here->qgs_x2;
      pass.cyy = here->qgs_y2;
      pass.czz = 0.0;
      pass.cxy = here->qgs_xy;
      pass.cyz = 0.0;
      pass.cxz = 0.0;
      pass.cxxx = here->qgs_x3;
      pass.cyyy = here->qgs_y3;
      pass.czzz = 0.0;
      pass.cxxy = here->qgs_x2y;
      pass.cxxz = 0.0;
      pass.cxyy = here->qgs_xy2;
      pass.cyyz = 0.0;
      pass.cxzz = 0.0;
      pass.cyzz = 0.0;
      pass.cxyz = 0.0;
      pass.r1h1x = r1h1x;
      pass.i1h1x = i1h1x;
      pass.r1h1y = r1h1x - r1h1y;
      pass.i1h1y = i1h1x - i1h1y;
      pass.r1h1z = 0.0;
      pass.i1h1z = 0.0;
      pass.r1h2x = r1hm2x;
      pass.i1h2x = i1hm2x;
      pass.r1h2y = r1hm2x - r1hm2y;
      pass.i1h2y = i1hm2x - i1hm2y;
      pass.r1h2z = 0.0;
      pass.i1h2z = 0.0;
      pass.r2h11x = r2h11x;
      pass.i2h11x = i2h11x;
      pass.r2h11y = r2h11x - r2h11y;
      pass.i2h11y = i2h11x - i2h11y;
      pass.r2h11z = 0.0;
      pass.i2h11z = 0.0;
      pass.h2f1f2x = r2h1m2x;
      pass.ih2f1f2x = i2h1m2x;
      pass.h2f1f2y = r2h1m2x - r2h1m2y;
      pass.ih2f1f2y = i2h1m2x - i2h1m2y;
      pass.h2f1f2z = 0.0;
      pass.ih2f1f2z = 0.0;
      temp = -ckt->CKTomega*DFi2F12(&pass);

      itemp = ckt->CKTomega*DFn2F12(&pass);

      *(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 */

      pass.cxx = here->qgd_x2;
      pass.cyy = here->qgd_y2;
      pass.czz = 0.0;
      pass.cxy = here->qgd_xy;
      pass.cyz = 0.0;
      pass.cxz = 0.0;
      pass.cxxx = here->qgd_x3;
      pass.cyyy = here->qgd_y3;
      pass.czzz = 0.0;
      pass.cxxy = here->qgd_x2y;
      pass.cxxz = 0.0;
      pass.cxyy = here->qgd_xy2;
      pass.cyyz = 0.0;
      pass.cxzz = 0.0;
      pass.cyzz = 0.0;
      pass.cxyz = 0.0;
      pass.r1h1x = r1h1x;
      pass.i1h1x = i1h1x;
      pass.r1h1y = r1h1x - r1h1y;
      pass.i1h1y = i1h1x - i1h1y;
      pass.r1h1z = 0.0;
      pass.i1h1z = 0.0;
      pass.r1h2x = r1hm2x;
      pass.i1h2x = i1hm2x;
      pass.r1h2y = r1hm2x - r1hm2y;
      pass.i1h2y = i1hm2x - i1hm2y;
      pass.r1h2z = 0.0;
      pass.i1h2z = 0.0;
      pass.r2h11x = r2h11x;
      pass.i2h11x = i2h11x;
      pass.r2h11y = r2h11x - r2h11y;
      pass.i2h11y = i2h11x - i2h11y;
      pass.r2h11z = 0.0;
      pass.i2h11z = 0.0;
      pass.h2f1f2x = r2h1m2x;
      pass.ih2f1f2x = i2h1m2x;
      pass.h2f1f2y = r2h1m2x - r2h1m2y;
      pass.ih2f1f2y = i2h1m2x - i2h1m2y;
      pass.h2f1f2z = 0.0;
      pass.ih2f1f2z = 0.0;
      temp = -ckt->CKTomega*DFi2F12(&pass);

      itemp = ckt->CKTomega*DFn2F12(&pass);

      *(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;
    default:
;
    }
  }
}
return(OK);
}
else
  return(E_BADPARM);
}