nbjtask.c

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

/**********
Copyright 1990 Regents of the University of California.  All rights reserved.
Author: 1987 Thomas L. Quarles
**********/

#include "spice.h"
#include <stdio.h>
#include "const.h"
#include "util.h"
#include "ifsim.h"
#include "cktdefs.h"
#include "devdefs.h"
#include "nbjtdefs.h"
#include "sperror.h"
#include "suffix.h"

/* ARGSUSED */
int
NBJTask(ckt, inInst, which, value, select)
  CKTcircuit *ckt;
  GENinstance *inInst;
  int which;
  IFvalue *value;
  IFvalue *select;
{
  NBJTinstance *inst = (NBJTinstance *) inInst;
  switch (which) {
  case NBJT_AREA:
    value->rValue = inst->NBJTarea;
    return (OK);
  case NBJT_TEMP:
    value->rValue = inst->NBJTtemp - CONSTCtoK;
    return (OK);
  case NBJT_G11:
    value->rValue = *(ckt->CKTstate0 + inst->NBJTdIcDVce);
    return (OK);
  case NBJT_G12:
    value->rValue = *(ckt->CKTstate0 + inst->NBJTdIcDVbe);
    return (OK);
  case NBJT_G13:
    value->rValue = -*(ckt->CKTstate0 + inst->NBJTdIcDVce)
	- *(ckt->CKTstate0 + inst->NBJTdIcDVbe);
    return (OK);
  case NBJT_G21:
    value->rValue = *(ckt->CKTstate0 + inst->NBJTdIeDVce)
	- *(ckt->CKTstate0 + inst->NBJTdIcDVce);
    return (OK);
  case NBJT_G22:
    value->rValue = *(ckt->CKTstate0 + inst->NBJTdIeDVbe)
	- *(ckt->CKTstate0 + inst->NBJTdIcDVbe);
    return (OK);
  case NBJT_G23:
    value->rValue = -*(ckt->CKTstate0 + inst->NBJTdIeDVce)
	+ *(ckt->CKTstate0 + inst->NBJTdIcDVce);
    -*(ckt->CKTstate0 + inst->NBJTdIeDVbe)
	+ *(ckt->CKTstate0 + inst->NBJTdIcDVbe);
    return (OK);
  case NBJT_G31:
    value->rValue = -*(ckt->CKTstate0 + inst->NBJTdIeDVce);
    return (OK);
  case NBJT_G32:
    value->rValue = -*(ckt->CKTstate0 + inst->NBJTdIeDVbe);
    return (OK);
  case NBJT_G33:
    value->rValue = *(ckt->CKTstate0 + inst->NBJTdIeDVce)
	+ *(ckt->CKTstate0 + inst->NBJTdIeDVbe);
    return (OK);
  case NBJT_C11:
    if (!inst->NBJTsmSigAvail
	&& ckt->CKTcurrentAnalysis != DOING_TRAN) {
      NBJTinitSmSig(inst);
    }
    value->rValue = inst->NBJTc11;
    return (OK);
  case NBJT_C12:
    if (!inst->NBJTsmSigAvail
	&& ckt->CKTcurrentAnalysis != DOING_TRAN) {
      NBJTinitSmSig(inst);
    }
    value->rValue = inst->NBJTc12;
    return (OK);
  case NBJT_C13:
    if (!inst->NBJTsmSigAvail
	&& ckt->CKTcurrentAnalysis != DOING_TRAN) {
      NBJTinitSmSig(inst);
    }
    value->rValue = -inst->NBJTc11 - inst->NBJTc12;
    return (OK);
  case NBJT_C21:
    if (!inst->NBJTsmSigAvail
	&& ckt->CKTcurrentAnalysis != DOING_TRAN) {
      NBJTinitSmSig(inst);
    }
    value->rValue = inst->NBJTc21;
    return (OK);
  case NBJT_C22:
    if (!inst->NBJTsmSigAvail
	&& ckt->CKTcurrentAnalysis != DOING_TRAN) {
      NBJTinitSmSig(inst);
    }
    value->rValue = inst->NBJTc22;
    return (OK);
  case NBJT_C23:
    if (!inst->NBJTsmSigAvail
	&& ckt->CKTcurrentAnalysis != DOING_TRAN) {
      NBJTinitSmSig(inst);
    }
    value->rValue = -inst->NBJTc21 - inst->NBJTc22;
    return (OK);
  case NBJT_C31:
    if (!inst->NBJTsmSigAvail
	&& ckt->CKTcurrentAnalysis != DOING_TRAN) {
      NBJTinitSmSig(inst);
    }
    value->rValue = -inst->NBJTc11 - inst->NBJTc21;
    return (OK);
  case NBJT_C32:
    if (!inst->NBJTsmSigAvail
	&& ckt->CKTcurrentAnalysis != DOING_TRAN) {
      NBJTinitSmSig(inst);
    }
    value->rValue = -inst->NBJTc12 - inst->NBJTc22;
    return (OK);
  case NBJT_C33:
    if (!inst->NBJTsmSigAvail
	&& ckt->CKTcurrentAnalysis != DOING_TRAN) {
      NBJTinitSmSig(inst);
    }
    value->rValue = inst->NBJTc11 + inst->NBJTc21
	+ inst->NBJTc12 + inst->NBJTc22;
    return (OK);
  case NBJT_Y11:
    if (!inst->NBJTsmSigAvail
	&& ckt->CKTcurrentAnalysis != DOING_TRAN) {
      NBJTinitSmSig(inst);
    }
    value->cValue.real = inst->NBJTy11r;
    value->cValue.imag = inst->NBJTy11i;
    return (OK);
  case NBJT_Y12:
    if (!inst->NBJTsmSigAvail
	&& ckt->CKTcurrentAnalysis != DOING_TRAN) {
      NBJTinitSmSig(inst);
    }
    value->cValue.real = inst->NBJTy12r;
    value->cValue.imag = inst->NBJTy12i;
    return (OK);
  case NBJT_Y13:
    if (!inst->NBJTsmSigAvail
	&& ckt->CKTcurrentAnalysis != DOING_TRAN) {
      NBJTinitSmSig(inst);
    }
    value->cValue.real = -inst->NBJTy11r - inst->NBJTy12r;
    value->cValue.imag = -inst->NBJTy11i - inst->NBJTy12i;
    return (OK);
  case NBJT_Y21:
    if (!inst->NBJTsmSigAvail
	&& ckt->CKTcurrentAnalysis != DOING_TRAN) {
      NBJTinitSmSig(inst);
    }
    value->cValue.real = inst->NBJTy21r;
    value->cValue.imag = inst->NBJTy21i;
    return (OK);
  case NBJT_Y22:
    if (!inst->NBJTsmSigAvail
	&& ckt->CKTcurrentAnalysis != DOING_TRAN) {
      NBJTinitSmSig(inst);
    }
    value->cValue.real = inst->NBJTy22r;
    value->cValue.imag = inst->NBJTy22i;
    return (OK);
  case NBJT_Y23:
    if (!inst->NBJTsmSigAvail
	&& ckt->CKTcurrentAnalysis != DOING_TRAN) {
      NBJTinitSmSig(inst);
    }
    value->cValue.real = -inst->NBJTy21r - inst->NBJTy22r;
    value->cValue.imag = -inst->NBJTy21i - inst->NBJTy22i;
    return (OK);
  case NBJT_Y31:
    if (!inst->NBJTsmSigAvail
	&& ckt->CKTcurrentAnalysis != DOING_TRAN) {
      NBJTinitSmSig(inst);
    }
    value->cValue.real = -inst->NBJTy11r - inst->NBJTy21r;
    value->cValue.imag = -inst->NBJTy11i - inst->NBJTy21i;
    return (OK);
  case NBJT_Y32:
    if (!inst->NBJTsmSigAvail
	&& ckt->CKTcurrentAnalysis != DOING_TRAN) {
      NBJTinitSmSig(inst);
    }
    value->cValue.real = -inst->NBJTy12r - inst->NBJTy22r;
    value->cValue.imag = -inst->NBJTy12i - inst->NBJTy22i;
    return (OK);
  case NBJT_Y33:
    if (!inst->NBJTsmSigAvail
	&& ckt->CKTcurrentAnalysis != DOING_TRAN) {
      NBJTinitSmSig(inst);
    }
    value->cValue.real = inst->NBJTy11r + inst->NBJTy21r
	+ inst->NBJTy12r + inst->NBJTy22r;
    value->cValue.imag = inst->NBJTy11i + inst->NBJTy21i
	+ inst->NBJTy12i + inst->NBJTy22i;
    return (OK);
  default:
    return (E_BADPARM);
  }
  /* NOTREACHED */
}