/********** Copyright 1990 Regents of the University of California. All rights reserved. Author: 1985 Mathew Lew and Thomas L. Quarles Model Author: 1995 Colin McAndrew Motorola Spice3 Implementation: 2003 Dietmar Warning DAnalyse GmbH **********/ /* * This routine gives access to the internal device * parameters for VBICs */ #include "ngspice/ngspice.h" #include "ngspice/const.h" #include "ngspice/cktdefs.h" #include "vbicdefs.h" #include "ngspice/ifsim.h" #include "ngspice/sperror.h" #include "ngspice/suffix.h" /*ARGSUSED*/ int VBICask(CKTcircuit *ckt, GENinstance *instPtr, int which, IFvalue *value, IFvalue *select) { VBICinstance *here = (VBICinstance*)instPtr; int itmp; double vr; double vi; double sr; double si; double vm; switch(which) { case VBIC_AREA: value->rValue = here->VBICarea; return(OK); case VBIC_OFF: value->iValue = here->VBICoff; return(OK); case VBIC_IC_VBE: value->rValue = here->VBICicVBE; return(OK); case VBIC_IC_VCE: value->rValue = here->VBICicVCE; return(OK); case VBIC_TEMP: value->rValue = here->VBICtemp - CONSTCtoK; return(OK); case VBIC_M: value->rValue = here->VBICm; return(OK); case VBIC_QUEST_COLLNODE: value->iValue = here->VBICcollNode; return(OK); case VBIC_QUEST_BASENODE: value->iValue = here->VBICbaseNode; return(OK); case VBIC_QUEST_EMITNODE: value->iValue = here->VBICemitNode; return(OK); case VBIC_QUEST_SUBSNODE: value->iValue = here->VBICsubsNode; return(OK); case VBIC_QUEST_COLLCXNODE: value->iValue = here->VBICcollCXNode; return(OK); case VBIC_QUEST_BASEBXNODE: value->iValue = here->VBICbaseBXNode; return(OK); case VBIC_QUEST_EMITEINODE: value->iValue = here->VBICemitEINode; return(OK); case VBIC_QUEST_SUBSSINODE: value->iValue = here->VBICsubsSINode; return(OK); case VBIC_QUEST_VBE: value->rValue = *(ckt->CKTstate0 + here->VBICvbei); return(OK); case VBIC_QUEST_VBC: value->rValue = *(ckt->CKTstate0 + here->VBICvbci); return(OK); case VBIC_QUEST_CC: value->rValue = *(ckt->CKTstate0 + here->VBICitzf) - *(ckt->CKTstate0 + here->VBICitzr) - *(ckt->CKTstate0 + here->VBICibc); return(OK); case VBIC_QUEST_CB: value->rValue = *(ckt->CKTstate0 + here->VBICibe) + *(ckt->CKTstate0 + here->VBICibc) + *(ckt->CKTstate0 + here->VBICibex) + *(ckt->CKTstate0 + here->VBICibep) + *(ckt->CKTstate0 + here->VBICiccp); return(OK); case VBIC_QUEST_CE: value->rValue = - *(ckt->CKTstate0 + here->VBICibe) - *(ckt->CKTstate0 + here->VBICibex) - *(ckt->CKTstate0 + here->VBICitzf) + *(ckt->CKTstate0 + here->VBICitzr); return(OK); case VBIC_QUEST_CS: value->rValue = *(ckt->CKTstate0 + here->VBICiccp) - *(ckt->CKTstate0 + here->VBICibcp); return(OK); case VBIC_QUEST_POWER: value->rValue = fabs(*(ckt->CKTstate0 + here->VBICitzf) - *(ckt->CKTstate0 + here->VBICitzr)) * fabs(*(ckt->CKTstate0 + here->VBICvbei) - *(ckt->CKTstate0 + here->VBICvbci)) + fabs(*(ckt->CKTstate0 + here->VBICibe) * *(ckt->CKTstate0 + here->VBICvbei)) + fabs(*(ckt->CKTstate0 + here->VBICibex) * *(ckt->CKTstate0 + here->VBICvbex)) + fabs(*(ckt->CKTstate0 + here->VBICibc) * *(ckt->CKTstate0 + here->VBICvbci)) + fabs(*(ckt->CKTstate0 + here->VBICibcp) * *(ckt->CKTstate0 + here->VBICvbcp)) + fabs(*(ckt->CKTstate0 + here->VBICiccp)) * fabs(*(ckt->CKTstate0 + here->VBICvbep) - *(ckt->CKTstate0 + here->VBICvbcp)); return(OK); case VBIC_QUEST_GM: value->rValue = *(ckt->CKTstate0 + here->VBICitzf_Vbei); return(OK); case VBIC_QUEST_GO: value->rValue = *(ckt->CKTstate0 + here->VBICitzf_Vbci); return(OK); case VBIC_QUEST_GPI: value->rValue = *(ckt->CKTstate0 + here->VBICibe_Vbei); return(OK); case VBIC_QUEST_GMU: value->rValue = *(ckt->CKTstate0 + here->VBICibc_Vbci); return(OK); case VBIC_QUEST_GX: value->rValue = *(ckt->CKTstate0 + here->VBICirbi_Vrbi); return(OK); case VBIC_QUEST_CBE: value->rValue = here->VBICcapbe; return(OK); case VBIC_QUEST_CBEX: value->rValue = here->VBICcapbex; return(OK); case VBIC_QUEST_CBC: value->rValue = here->VBICcapbc; return(OK); case VBIC_QUEST_CBCX: value->rValue = here->VBICcapbcx; return(OK); case VBIC_QUEST_CBEP: value->rValue = here->VBICcapbep; return(OK); case VBIC_QUEST_CBCP: value->rValue = here->VBICcapbcp; return(OK); case VBIC_QUEST_QBE: value->rValue = *(ckt->CKTstate0 + here->VBICqbe); return(OK); case VBIC_QUEST_QBC: value->rValue = *(ckt->CKTstate0 + here->VBICqbc); return(OK); case VBIC_QUEST_SENS_DC: if(ckt->CKTsenInfo){ value->rValue = *(ckt->CKTsenInfo->SEN_Sap[select->iValue + 1]+ here->VBICsenParmNo); } return(OK); case VBIC_QUEST_SENS_REAL: if(ckt->CKTsenInfo){ value->rValue = *(ckt->CKTsenInfo->SEN_RHS[select->iValue + 1]+ here->VBICsenParmNo); } return(OK); case VBIC_QUEST_SENS_IMAG: if(ckt->CKTsenInfo){ value->rValue = *(ckt->CKTsenInfo->SEN_iRHS[select->iValue + 1]+ here->VBICsenParmNo); } return(OK); case VBIC_QUEST_SENS_MAG: if(ckt->CKTsenInfo){ vr = *(ckt->CKTrhsOld + select->iValue + 1); vi = *(ckt->CKTirhsOld + select->iValue + 1); vm = sqrt(vr*vr + vi*vi); if(vm == 0){ value->rValue = 0; return(OK); } sr = *(ckt->CKTsenInfo->SEN_RHS[select->iValue + 1]+ here->VBICsenParmNo); si = *(ckt->CKTsenInfo->SEN_iRHS[select->iValue + 1]+ here->VBICsenParmNo); value->rValue = (vr * sr + vi * si)/vm; } return(OK); case VBIC_QUEST_SENS_PH: if(ckt->CKTsenInfo){ vr = *(ckt->CKTrhsOld + select->iValue + 1); vi = *(ckt->CKTirhsOld + select->iValue + 1); vm = vr*vr + vi*vi; if(vm == 0){ value->rValue = 0; return(OK); } sr = *(ckt->CKTsenInfo->SEN_RHS[select->iValue + 1]+ here->VBICsenParmNo); si = *(ckt->CKTsenInfo->SEN_iRHS[select->iValue + 1]+ here->VBICsenParmNo); value->rValue = (vr * si - vi * sr)/vm; } return(OK); case VBIC_QUEST_SENS_CPLX: if(ckt->CKTsenInfo){ itmp = select->iValue + 1; value->cValue.real= *(ckt->CKTsenInfo->SEN_RHS[itmp]+ here->VBICsenParmNo); value->cValue.imag= *(ckt->CKTsenInfo->SEN_iRHS[itmp]+ here->VBICsenParmNo); } return(OK); default: return(E_BADPARM); } /* NOTREACHED */ }