/**********
Copyright 1991 Regents of the University of California.  All rights reserved.
Author:	1987 Kartikeya Mayaram, U. C. Berkeley CAD Group
**********/

/*
 * Function to load the COMPLEX circuit matrix using the small signal
 * parameters saved during a previous DC operating point analysis.
 */

#include "ngspice/ngspice.h"
#include "ngspice/cktdefs.h"
#include "numosdef.h"
#include "ngspice/sperror.h"
#include "ngspice/suffix.h"
#include "ngspice/complex.h"
#include "../../../ciderlib/twod/twoddefs.h"
#include "../../../ciderlib/twod/twodext.h"
#include "ngspice/cidersupt.h"


int
NUMOSacLoad(GENmodel *inModel, CKTcircuit *ckt)
{
  register NUMOSmodel *model = (NUMOSmodel *) inModel;
  register NUMOSinstance *inst;
  struct mosAdmittances yAc;
  double startTime;

  for (; model != NULL; model = model->NUMOSnextModel) {
    FieldDepMobility = model->NUMOSmodels->MODLfieldDepMobility;
    TransDepMobility = model->NUMOSmodels->MODLtransDepMobility;
    SurfaceMobility = model->NUMOSmodels->MODLsurfaceMobility;
    Srh = model->NUMOSmodels->MODLsrh;
    Auger = model->NUMOSmodels->MODLauger;
    AvalancheGen = model->NUMOSmodels->MODLavalancheGen;
    OneCarrier = model->NUMOSmethods->METHoneCarrier;
    AcAnalysisMethod = model->NUMOSmethods->METHacAnalysisMethod;
    MobDeriv = model->NUMOSmethods->METHmobDeriv;
    TWOacDebug = model->NUMOSoutputs->OUTPacDebug;

    for (inst = model->NUMOSinstances; inst != NULL;
	inst = inst->NUMOSnextInstance) {

      startTime = SPfrontEnd->IFseconds();
      /* Get Temp.-Dep. Global Parameters */
      GLOBgetGlobals(&(inst->NUMOSglobals));

      model->NUMOSmethods->METHacAnalysisMethod =
	  NUMOSadmittance(inst->NUMOSpDevice,
	  ckt->CKTomega, &yAc);

      *(inst->NUMOSdrainDrainPtr) += yAc.yIdVdb.real;
      *(inst->NUMOSdrainDrainPtr + 1) += yAc.yIdVdb.imag;
      *(inst->NUMOSdrainSourcePtr) += yAc.yIdVsb.real;
      *(inst->NUMOSdrainSourcePtr + 1) += yAc.yIdVsb.imag;
      *(inst->NUMOSdrainGatePtr) += yAc.yIdVgb.real;
      *(inst->NUMOSdrainGatePtr + 1) += yAc.yIdVgb.imag;
      *(inst->NUMOSdrainBulkPtr) -=
	  yAc.yIdVdb.real + yAc.yIdVsb.real + yAc.yIdVgb.real;
      *(inst->NUMOSdrainBulkPtr + 1) -=
	  yAc.yIdVdb.imag + yAc.yIdVsb.imag + yAc.yIdVgb.imag;

      *(inst->NUMOSsourceDrainPtr) += yAc.yIsVdb.real;
      *(inst->NUMOSsourceDrainPtr + 1) += yAc.yIsVdb.imag;
      *(inst->NUMOSsourceSourcePtr) += yAc.yIsVsb.real;
      *(inst->NUMOSsourceSourcePtr + 1) += yAc.yIsVsb.imag;
      *(inst->NUMOSsourceGatePtr) += yAc.yIsVgb.real;
      *(inst->NUMOSsourceGatePtr + 1) += yAc.yIsVgb.imag;
      *(inst->NUMOSsourceBulkPtr) -=
	  yAc.yIsVdb.real + yAc.yIsVsb.real + yAc.yIsVgb.real;
      *(inst->NUMOSsourceBulkPtr + 1) -=
	  yAc.yIsVdb.imag + yAc.yIsVsb.imag + yAc.yIsVgb.imag;

      *(inst->NUMOSgateDrainPtr) += yAc.yIgVdb.real;
      *(inst->NUMOSgateDrainPtr + 1) += yAc.yIgVdb.imag;
      *(inst->NUMOSgateSourcePtr) += yAc.yIgVsb.real;
      *(inst->NUMOSgateSourcePtr + 1) += yAc.yIgVsb.imag;
      *(inst->NUMOSgateGatePtr) += yAc.yIgVgb.real;
      *(inst->NUMOSgateGatePtr + 1) += yAc.yIgVgb.imag;
      *(inst->NUMOSgateBulkPtr) -=
	  yAc.yIgVdb.real + yAc.yIgVsb.real + yAc.yIgVgb.real;
      *(inst->NUMOSgateBulkPtr + 1) -=
	  yAc.yIgVdb.imag + yAc.yIgVsb.imag + yAc.yIgVgb.imag;

      *(inst->NUMOSbulkDrainPtr) -=
	  yAc.yIdVdb.real + yAc.yIsVdb.real + yAc.yIgVdb.real;
      *(inst->NUMOSbulkDrainPtr + 1) -=
	  yAc.yIdVdb.imag + yAc.yIsVdb.imag + yAc.yIgVdb.imag;
      *(inst->NUMOSbulkSourcePtr) -=
	  yAc.yIdVsb.real + yAc.yIsVsb.real + yAc.yIgVsb.real;
      *(inst->NUMOSbulkSourcePtr + 1) -=
	  yAc.yIdVsb.imag + yAc.yIsVsb.imag + yAc.yIgVsb.imag;
      *(inst->NUMOSbulkGatePtr) -=
	  yAc.yIdVgb.real + yAc.yIsVgb.real + yAc.yIgVgb.real;
      *(inst->NUMOSbulkGatePtr + 1) -=
	  yAc.yIdVgb.imag + yAc.yIsVgb.imag + yAc.yIgVgb.imag;
      *(inst->NUMOSbulkBulkPtr) += yAc.yIdVdb.real + yAc.yIdVsb.real +
	  yAc.yIdVgb.real + yAc.yIsVdb.real +
	  yAc.yIsVsb.real + yAc.yIsVgb.real +
	  yAc.yIgVdb.real + yAc.yIgVsb.real +
	  yAc.yIgVgb.real;
      *(inst->NUMOSbulkBulkPtr + 1) -= yAc.yIdVdb.imag + yAc.yIdVsb.imag +
	  yAc.yIdVgb.imag + yAc.yIsVdb.imag +
	  yAc.yIsVsb.imag + yAc.yIsVgb.imag +
	  yAc.yIgVdb.imag + yAc.yIgVsb.imag +
	  yAc.yIgVgb.imag;
      if (ckt->CKTomega != 0.0) {
	inst->NUMOSc11 = yAc.yIdVdb.imag / ckt->CKTomega;
	inst->NUMOSc12 = yAc.yIdVgb.imag / ckt->CKTomega;
	inst->NUMOSc13 = yAc.yIdVsb.imag / ckt->CKTomega;
	inst->NUMOSc21 = yAc.yIgVdb.imag / ckt->CKTomega;
	inst->NUMOSc22 = yAc.yIgVgb.imag / ckt->CKTomega;
	inst->NUMOSc23 = yAc.yIgVsb.imag / ckt->CKTomega;
	inst->NUMOSc31 = yAc.yIsVdb.imag / ckt->CKTomega;
	inst->NUMOSc32 = yAc.yIsVgb.imag / ckt->CKTomega;
	inst->NUMOSc33 = yAc.yIsVsb.imag / ckt->CKTomega;
      } else {
	inst->NUMOSc11 = 0.0;	/* XXX What else can be done?! */
	inst->NUMOSc12 = 0.0;	/* XXX What else can be done?! */
	inst->NUMOSc13 = 0.0;	/* XXX What else can be done?! */
	inst->NUMOSc21 = 0.0;	/* XXX What else can be done?! */
	inst->NUMOSc22 = 0.0;	/* XXX What else can be done?! */
	inst->NUMOSc23 = 0.0;	/* XXX What else can be done?! */
	inst->NUMOSc31 = 0.0;	/* XXX What else can be done?! */
	inst->NUMOSc32 = 0.0;	/* XXX What else can be done?! */
	inst->NUMOSc33 = 0.0;	/* XXX What else can be done?! */
      }
      inst->NUMOSy11r = yAc.yIdVdb.real;
      inst->NUMOSy11i = yAc.yIdVdb.imag;
      inst->NUMOSy12r = yAc.yIdVgb.real;
      inst->NUMOSy12i = yAc.yIdVgb.imag;
      inst->NUMOSy13r = yAc.yIdVsb.real;
      inst->NUMOSy13i = yAc.yIdVsb.imag;
      inst->NUMOSy21r = yAc.yIgVdb.real;
      inst->NUMOSy21i = yAc.yIgVdb.imag;
      inst->NUMOSy22r = yAc.yIgVgb.real;
      inst->NUMOSy22i = yAc.yIgVgb.imag;
      inst->NUMOSy23r = yAc.yIgVsb.real;
      inst->NUMOSy23i = yAc.yIgVsb.imag;
      inst->NUMOSy31r = yAc.yIsVdb.real;
      inst->NUMOSy31i = yAc.yIsVdb.imag;
      inst->NUMOSy32r = yAc.yIsVgb.real;
      inst->NUMOSy32i = yAc.yIsVgb.imag;
      inst->NUMOSy33r = yAc.yIsVsb.real;
      inst->NUMOSy33i = yAc.yIsVsb.imag;
      inst->NUMOSsmSigAvail = TRUE;
      inst->NUMOSpDevice->pStats->totalTime[STAT_AC] +=
	  SPfrontEnd->IFseconds() - startTime;
    }
  }
  return (OK);
}