/**** BSIM3v3.2.4, Released by Xuemei Xi 12/21/2001 ****/

/**********
 * Copyright 2001 Regents of the University of California. All rights reserved.
 * File: b3pzld.c of BSIM3v3.2.4
 * Author: 1995 Min-Chie Jeng and Mansun Chan.
 * Author: 1997-1999 Weidong Liu.
 * Author: 2001  Xuemei Xi
 * Modified by Paolo Nenzi 2002
 **********/

#include "ngspice/ngspice.h"
#include "ngspice/cktdefs.h"
#include "ngspice/complex.h"
#include "ngspice/sperror.h"
#include "bsim3v32def.h"
#include "ngspice/suffix.h"

int
BSIM3v32pzLoad (GENmodel *inModel, CKTcircuit *ckt, SPcomplex *s)
{
BSIM3v32model *model = (BSIM3v32model*)inModel;
BSIM3v32instance *here;
double xcggb, xcgdb, xcgsb, xcgbb, xcbgb, xcbdb, xcbsb, xcbbb;
double xcdgb, xcddb, xcdsb, xcdbb, xcsgb, xcsdb, xcssb, xcsbb;
double gdpr, gspr, gds, gbd, gbs, capbd, capbs, FwdSum, RevSum, Gm, Gmbs;
double cggb, cgdb, cgsb, cbgb, cbdb, cbsb, cddb, cdgb, cdsb;
double GSoverlapCap, GDoverlapCap, GBoverlapCap;
double dxpart, sxpart, xgtg, xgtd, xgts, xgtb;
double xcqgb = 0.0, xcqdb = 0.0, xcqsb = 0.0, xcqbb = 0.0;
double gbspsp, gbbdp, gbbsp, gbspg, gbspb;
double gbspdp, gbdpdp, gbdpg, gbdpb, gbdpsp;
double ddxpart_dVd, ddxpart_dVg, ddxpart_dVb, ddxpart_dVs;
double dsxpart_dVd, dsxpart_dVg, dsxpart_dVb, dsxpart_dVs;
double T1, CoxWL, qcheq, Cdg, Cdd, Cds, Csg, Csd, Css;
double ScalingFactor = 1.0e-9;
double m;

    for (; model != NULL; model = model->BSIM3v32nextModel)
    {    for (here = model->BSIM3v32instances; here!= NULL;
              here = here->BSIM3v32nextInstance)
         {
              if (here->BSIM3v32mode >= 0)
              {   Gm = here->BSIM3v32gm;
                  Gmbs = here->BSIM3v32gmbs;
                  FwdSum = Gm + Gmbs;
                  RevSum = 0.0;

                  gbbdp = -here->BSIM3v32gbds;
                  gbbsp = here->BSIM3v32gbds + here->BSIM3v32gbgs + here->BSIM3v32gbbs;

                  gbdpg = here->BSIM3v32gbgs;
                  gbdpdp = here->BSIM3v32gbds;
                  gbdpb = here->BSIM3v32gbbs;
                  gbdpsp = -(gbdpg + gbdpdp + gbdpb);

                  gbspg = 0.0;
                  gbspdp = 0.0;
                  gbspb = 0.0;
                  gbspsp = 0.0;

                  if (here->BSIM3v32nqsMod == 0)
                  {   cggb = here->BSIM3v32cggb;
                      cgsb = here->BSIM3v32cgsb;
                      cgdb = here->BSIM3v32cgdb;

                      cbgb = here->BSIM3v32cbgb;
                      cbsb = here->BSIM3v32cbsb;
                      cbdb = here->BSIM3v32cbdb;

                      cdgb = here->BSIM3v32cdgb;
                      cdsb = here->BSIM3v32cdsb;
                      cddb = here->BSIM3v32cddb;

                      xgtg = xgtd = xgts = xgtb = 0.0;
                      sxpart = 0.6;
                      dxpart = 0.4;
                      ddxpart_dVd = ddxpart_dVg = ddxpart_dVb
                                  = ddxpart_dVs = 0.0;
                      dsxpart_dVd = dsxpart_dVg = dsxpart_dVb
                                  = dsxpart_dVs = 0.0;
                  }
                  else
                  {   cggb = cgdb = cgsb = 0.0;
                      cbgb = cbdb = cbsb = 0.0;
                      cdgb = cddb = cdsb = 0.0;

                      xgtg = here->BSIM3v32gtg;
                      xgtd = here->BSIM3v32gtd;
                      xgts = here->BSIM3v32gts;
                      xgtb = here->BSIM3v32gtb;

                      xcqgb = here->BSIM3v32cqgb;
                      xcqdb = here->BSIM3v32cqdb;
                      xcqsb = here->BSIM3v32cqsb;
                      xcqbb = here->BSIM3v32cqbb;

                      CoxWL = model->BSIM3v32cox * here->pParam->BSIM3v32weffCV
                            * here->pParam->BSIM3v32leffCV;
                      qcheq = -(here->BSIM3v32qgate + here->BSIM3v32qbulk);
                      if (fabs(qcheq) <= 1.0e-5 * CoxWL)
                      {   if (model->BSIM3v32xpart < 0.5)
                          {   dxpart = 0.4;
                          }
                          else if (model->BSIM3v32xpart > 0.5)
                          {   dxpart = 0.0;
                          }
                          else
                          {   dxpart = 0.5;
                          }
                          ddxpart_dVd = ddxpart_dVg = ddxpart_dVb
                                      = ddxpart_dVs = 0.0;
                      }
                      else
                      {   dxpart = here->BSIM3v32qdrn / qcheq;
                          Cdd = here->BSIM3v32cddb;
                          Csd = -(here->BSIM3v32cgdb + here->BSIM3v32cddb
                              + here->BSIM3v32cbdb);
                          ddxpart_dVd = (Cdd - dxpart * (Cdd + Csd)) / qcheq;
                          Cdg = here->BSIM3v32cdgb;
                          Csg = -(here->BSIM3v32cggb + here->BSIM3v32cdgb
                              + here->BSIM3v32cbgb);
                          ddxpart_dVg = (Cdg - dxpart * (Cdg + Csg)) / qcheq;

                          Cds = here->BSIM3v32cdsb;
                          Css = -(here->BSIM3v32cgsb + here->BSIM3v32cdsb
                              + here->BSIM3v32cbsb);
                          ddxpart_dVs = (Cds - dxpart * (Cds + Css)) / qcheq;

                          ddxpart_dVb = -(ddxpart_dVd + ddxpart_dVg
                                      + ddxpart_dVs);
                      }
                      sxpart = 1.0 - dxpart;
                      dsxpart_dVd = -ddxpart_dVd;
                      dsxpart_dVg = -ddxpart_dVg;
                      dsxpart_dVs = -ddxpart_dVs;
                      dsxpart_dVb = -(dsxpart_dVd + dsxpart_dVg + dsxpart_dVs);
                  }
              }
              else
              {   Gm = -here->BSIM3v32gm;
                  Gmbs = -here->BSIM3v32gmbs;
                  FwdSum = 0.0;
                  RevSum = -(Gm + Gmbs);

                  gbbsp = -here->BSIM3v32gbds;
                  gbbdp = here->BSIM3v32gbds + here->BSIM3v32gbgs + here->BSIM3v32gbbs;

                  gbdpg = 0.0;
                  gbdpsp = 0.0;
                  gbdpb = 0.0;
                  gbdpdp = 0.0;

                  gbspg = here->BSIM3v32gbgs;
                  gbspsp = here->BSIM3v32gbds;
                  gbspb = here->BSIM3v32gbbs;
                  gbspdp = -(gbspg + gbspsp + gbspb);

                  if (here->BSIM3v32nqsMod == 0)
                  {   cggb = here->BSIM3v32cggb;
                      cgsb = here->BSIM3v32cgdb;
                      cgdb = here->BSIM3v32cgsb;

                      cbgb = here->BSIM3v32cbgb;
                      cbsb = here->BSIM3v32cbdb;
                      cbdb = here->BSIM3v32cbsb;

                      cdgb = -(here->BSIM3v32cdgb + cggb + cbgb);
                      cdsb = -(here->BSIM3v32cddb + cgsb + cbsb);
                      cddb = -(here->BSIM3v32cdsb + cgdb + cbdb);

                      xgtg = xgtd = xgts = xgtb = 0.0;
                      sxpart = 0.4;
                      dxpart = 0.6;
                      ddxpart_dVd = ddxpart_dVg = ddxpart_dVb
                                  = ddxpart_dVs = 0.0;
                      dsxpart_dVd = dsxpart_dVg = dsxpart_dVb
                                  = dsxpart_dVs = 0.0;
                  }
                  else
                  {   cggb = cgdb = cgsb = 0.0;
                      cbgb = cbdb = cbsb = 0.0;
                      cdgb = cddb = cdsb = 0.0;

                      xgtg = here->BSIM3v32gtg;
                      xgtd = here->BSIM3v32gts;
                      xgts = here->BSIM3v32gtd;
                      xgtb = here->BSIM3v32gtb;

                      xcqgb = here->BSIM3v32cqgb;
                      xcqdb = here->BSIM3v32cqsb;
                      xcqsb = here->BSIM3v32cqdb;
                      xcqbb = here->BSIM3v32cqbb;

                      CoxWL = model->BSIM3v32cox * here->pParam->BSIM3v32weffCV
                            * here->pParam->BSIM3v32leffCV;
                      qcheq = -(here->BSIM3v32qgate + here->BSIM3v32qbulk);
                      if (fabs(qcheq) <= 1.0e-5 * CoxWL)
                      {   if (model->BSIM3v32xpart < 0.5)
                          {   sxpart = 0.4;
                          }
                          else if (model->BSIM3v32xpart > 0.5)
                          {   sxpart = 0.0;
                          }
                          else
                          {   sxpart = 0.5;
                          }
                          dsxpart_dVd = dsxpart_dVg = dsxpart_dVb
                                      = dsxpart_dVs = 0.0;
                      }
                      else
                      {   sxpart = here->BSIM3v32qdrn / qcheq;
                          Css = here->BSIM3v32cddb;
                          Cds = -(here->BSIM3v32cgdb + here->BSIM3v32cddb
                              + here->BSIM3v32cbdb);
                          dsxpart_dVs = (Css - sxpart * (Css + Cds)) / qcheq;
                          Csg = here->BSIM3v32cdgb;
                          Cdg = -(here->BSIM3v32cggb + here->BSIM3v32cdgb
                              + here->BSIM3v32cbgb);
                          dsxpart_dVg = (Csg - sxpart * (Csg + Cdg)) / qcheq;

                          Csd = here->BSIM3v32cdsb;
                          Cdd = -(here->BSIM3v32cgsb + here->BSIM3v32cdsb
                              + here->BSIM3v32cbsb);
                          dsxpart_dVd = (Csd - sxpart * (Csd + Cdd)) / qcheq;

                          dsxpart_dVb = -(dsxpart_dVd + dsxpart_dVg
                                      + dsxpart_dVs);
                      }
                      dxpart = 1.0 - sxpart;
                      ddxpart_dVd = -dsxpart_dVd;
                      ddxpart_dVg = -dsxpart_dVg;
                      ddxpart_dVs = -dsxpart_dVs;
                      ddxpart_dVb = -(ddxpart_dVd + ddxpart_dVg + ddxpart_dVs);
                  }
              }


              T1 = *(ckt->CKTstate0 + here->BSIM3v32qdef) * here->BSIM3v32gtau;
              gdpr = here->BSIM3v32drainConductance;
              gspr = here->BSIM3v32sourceConductance;
              gds = here->BSIM3v32gds;
              gbd = here->BSIM3v32gbd;
              gbs = here->BSIM3v32gbs;
              capbd = here->BSIM3v32capbd;
              capbs = here->BSIM3v32capbs;

              GSoverlapCap = here->BSIM3v32cgso;
              GDoverlapCap = here->BSIM3v32cgdo;
              GBoverlapCap = here->pParam->BSIM3v32cgbo;

              xcdgb = (cdgb - GDoverlapCap);
              xcddb = (cddb + capbd + GDoverlapCap);
              xcdsb = cdsb;
              xcdbb = -(xcdgb + xcddb + xcdsb);
              xcsgb = -(cggb + cbgb + cdgb + GSoverlapCap);
              xcsdb = -(cgdb + cbdb + cddb);
              xcssb = (capbs + GSoverlapCap - (cgsb + cbsb + cdsb));
              xcsbb = -(xcsgb + xcsdb + xcssb);
              xcggb = (cggb + GDoverlapCap + GSoverlapCap + GBoverlapCap);
              xcgdb = (cgdb - GDoverlapCap);
              xcgsb = (cgsb - GSoverlapCap);
              xcgbb = -(xcggb + xcgdb + xcgsb);
              xcbgb = (cbgb - GBoverlapCap);
              xcbdb = (cbdb - capbd);
              xcbsb = (cbsb - capbs);
              xcbbb = -(xcbgb + xcbdb + xcbsb);

              m = here->BSIM3v32m;

              *(here->BSIM3v32GgPtr ) += m * (xcggb * s->real);
              *(here->BSIM3v32GgPtr +1) += m * (xcggb * s->imag);
              *(here->BSIM3v32BbPtr ) += m * (xcbbb * s->real);
              *(here->BSIM3v32BbPtr +1) += m * (xcbbb * s->imag);
              *(here->BSIM3v32DPdpPtr ) += m * (xcddb * s->real);
              *(here->BSIM3v32DPdpPtr +1) += m * (xcddb * s->imag);
              *(here->BSIM3v32SPspPtr ) += m * (xcssb * s->real);
              *(here->BSIM3v32SPspPtr +1) += m * (xcssb * s->imag);

              *(here->BSIM3v32GbPtr ) += m * (xcgbb * s->real);
              *(here->BSIM3v32GbPtr +1) += m * (xcgbb * s->imag);
              *(here->BSIM3v32GdpPtr ) += m * (xcgdb * s->real);
              *(here->BSIM3v32GdpPtr +1) += m * (xcgdb * s->imag);
              *(here->BSIM3v32GspPtr ) += m * (xcgsb * s->real);
              *(here->BSIM3v32GspPtr +1) += m * (xcgsb * s->imag);

              *(here->BSIM3v32BgPtr ) += m * (xcbgb * s->real);
              *(here->BSIM3v32BgPtr +1) += m * (xcbgb * s->imag);
              *(here->BSIM3v32BdpPtr ) += m * (xcbdb * s->real);
              *(here->BSIM3v32BdpPtr +1) += m * (xcbdb * s->imag);
              *(here->BSIM3v32BspPtr ) += m * (xcbsb * s->real);
              *(here->BSIM3v32BspPtr +1) += m * (xcbsb * s->imag);

              *(here->BSIM3v32DPgPtr ) += m * (xcdgb * s->real);
              *(here->BSIM3v32DPgPtr +1) += m * (xcdgb * s->imag);
              *(here->BSIM3v32DPbPtr ) += m * (xcdbb * s->real);
              *(here->BSIM3v32DPbPtr +1) += m * (xcdbb * s->imag);
              *(here->BSIM3v32DPspPtr ) += m * (xcdsb * s->real);
              *(here->BSIM3v32DPspPtr +1) += m * (xcdsb * s->imag);

              *(here->BSIM3v32SPgPtr ) += m * (xcsgb * s->real);
              *(here->BSIM3v32SPgPtr +1) += m * (xcsgb * s->imag);
              *(here->BSIM3v32SPbPtr ) += m * (xcsbb * s->real);
              *(here->BSIM3v32SPbPtr +1) += m * (xcsbb * s->imag);
              *(here->BSIM3v32SPdpPtr ) += m * (xcsdb * s->real);
              *(here->BSIM3v32SPdpPtr +1) += m * (xcsdb * s->imag);

              *(here->BSIM3v32DdPtr) += m * gdpr;
              *(here->BSIM3v32DdpPtr) -= m * gdpr;
              *(here->BSIM3v32DPdPtr) -= m * gdpr;

              *(here->BSIM3v32SsPtr) += m * gspr;
              *(here->BSIM3v32SspPtr) -= m * gspr;
              *(here->BSIM3v32SPsPtr) -= m * gspr;

              *(here->BSIM3v32BgPtr) -= m * here->BSIM3v32gbgs;
              *(here->BSIM3v32BbPtr) += m * (gbd + gbs - here->BSIM3v32gbbs);
              *(here->BSIM3v32BdpPtr) -= m * (gbd - gbbdp);
              *(here->BSIM3v32BspPtr) -= m * (gbs - gbbsp);

              *(here->BSIM3v32DPgPtr) += Gm + dxpart * xgtg
                                   + T1 * ddxpart_dVg + gbdpg;
              *(here->BSIM3v32DPdpPtr) += gdpr + gds + gbd + RevSum
                                    + dxpart * xgtd + T1 * ddxpart_dVd + gbdpdp;
              *(here->BSIM3v32DPspPtr) -= gds + FwdSum - dxpart * xgts
                                    - T1 * ddxpart_dVs - gbdpsp;
              *(here->BSIM3v32DPbPtr) -= gbd - Gmbs - dxpart * xgtb
                                   - T1 * ddxpart_dVb - gbdpb;

              *(here->BSIM3v32SPgPtr) -= Gm - sxpart * xgtg
                                   - T1 * dsxpart_dVg - gbspg;
              *(here->BSIM3v32SPspPtr) += gspr + gds + gbs + FwdSum
                                   + sxpart * xgts + T1 * dsxpart_dVs + gbspsp;
              *(here->BSIM3v32SPbPtr) -= gbs + Gmbs - sxpart * xgtb
                                   - T1 * dsxpart_dVb - gbspb;
              *(here->BSIM3v32SPdpPtr) -= gds + RevSum - sxpart * xgtd
                                    - T1 * dsxpart_dVd - gbspdp;

              *(here->BSIM3v32GgPtr) -= xgtg;
              *(here->BSIM3v32GbPtr) -= xgtb;
              *(here->BSIM3v32GdpPtr) -= xgtd;
              *(here->BSIM3v32GspPtr) -= xgts;

              if (here->BSIM3v32nqsMod)
              {   *(here->BSIM3v32QqPtr ) += m * (s->real * ScalingFactor);
                  *(here->BSIM3v32QqPtr +1) += m * (s->imag * ScalingFactor);
                  *(here->BSIM3v32QgPtr ) -= m * (xcqgb * s->real);
                  *(here->BSIM3v32QgPtr +1) -= m * (xcqgb * s->imag);
                  *(here->BSIM3v32QdpPtr ) -= m * (xcqdb * s->real);
                  *(here->BSIM3v32QdpPtr +1) -= m * (xcqdb * s->imag);
                  *(here->BSIM3v32QbPtr ) -= m * (xcqbb * s->real);
                  *(here->BSIM3v32QbPtr +1) -= m * (xcqbb * s->imag);
                  *(here->BSIM3v32QspPtr ) -= m * (xcqsb * s->real);
                  *(here->BSIM3v32QspPtr +1) -= m * (xcqsb * s->imag);

                  *(here->BSIM3v32GqPtr) -= m * (here->BSIM3v32gtau);
                  *(here->BSIM3v32DPqPtr) += m * (dxpart * here->BSIM3v32gtau);
                  *(here->BSIM3v32SPqPtr) += m * (sxpart * here->BSIM3v32gtau);

                  *(here->BSIM3v32QqPtr) += m * (here->BSIM3v32gtau);
                  *(here->BSIM3v32QgPtr) += m * xgtg;
                  *(here->BSIM3v32QdpPtr) += m * xgtd;
                  *(here->BSIM3v32QbPtr) += m * xgtb;
                  *(here->BSIM3v32QspPtr) += m * xgts;
              }
         }
    }
    return(OK);
}