You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
146 lines
4.8 KiB
146 lines
4.8 KiB
/**********
|
|
Copyright 1990 Regents of the University of California. All rights reserved.
|
|
Author: 1985 Thomas L. Quarles
|
|
Modified by Dietmar Warning 2003 and Paolo Nenzi 2003
|
|
**********/
|
|
|
|
#include "ngspice.h"
|
|
#include "const.h"
|
|
#include "devdefs.h"
|
|
#include "ifsim.h"
|
|
#include "cktdefs.h"
|
|
#include "diodefs.h"
|
|
#include "sperror.h"
|
|
#include "suffix.h"
|
|
|
|
/* ARGSUSED */
|
|
int
|
|
DIOask (CKTcircuit *ckt, GENinstance *inst, int which, IFvalue *value,
|
|
IFvalue *select)
|
|
{
|
|
DIOinstance *here = (DIOinstance*)inst;
|
|
double vr;
|
|
double vi;
|
|
double sr;
|
|
double si;
|
|
double vm;
|
|
static char *msg = "Current and power not available for ac analysis";
|
|
|
|
switch (which) {
|
|
case DIO_OFF:
|
|
value->iValue = here->DIOoff;
|
|
return(OK);
|
|
case DIO_IC:
|
|
value->rValue = here->DIOinitCond;
|
|
return(OK);
|
|
case DIO_AREA:
|
|
value->rValue = here->DIOarea;
|
|
return(OK);
|
|
case DIO_PJ:
|
|
value->rValue = here->DIOpj;
|
|
return(OK);
|
|
case DIO_M:
|
|
value->rValue = here->DIOm;
|
|
return(OK);
|
|
|
|
case DIO_TEMP:
|
|
value->rValue = here->DIOtemp-CONSTCtoK;
|
|
return(OK);
|
|
case DIO_DTEMP:
|
|
value->rValue = here->DIOdtemp;
|
|
return(OK);
|
|
case DIO_VOLTAGE:
|
|
value->rValue = *(ckt->CKTstate0+here->DIOvoltage);
|
|
return(OK);
|
|
case DIO_CURRENT:
|
|
value->rValue = *(ckt->CKTstate0+here->DIOcurrent);
|
|
return(OK);
|
|
case DIO_CAP:
|
|
value->rValue = here->DIOcap;
|
|
return(OK);
|
|
case DIO_CHARGE:
|
|
value->rValue = *(ckt->CKTstate0+here->DIOcapCharge);
|
|
return(OK);
|
|
case DIO_CAPCUR:
|
|
value->rValue = *(ckt->CKTstate0+here->DIOcapCurrent);
|
|
return(OK);
|
|
case DIO_CONDUCT:
|
|
value->rValue = *(ckt->CKTstate0+here->DIOconduct);
|
|
return(OK);
|
|
case DIO_POWER :
|
|
if (ckt->CKTcurrentAnalysis & DOING_AC) {
|
|
errMsg = MALLOC(strlen(msg)+1);
|
|
errRtn = "DIOask";
|
|
strcpy(errMsg,msg);
|
|
return(E_ASKPOWER);
|
|
} else {
|
|
value->rValue = *(ckt->CKTstate0 + here->DIOcurrent) *
|
|
*(ckt->CKTstate0 + here->DIOvoltage);
|
|
}
|
|
return(OK);
|
|
case DIO_QUEST_SENS_DC:
|
|
if(ckt->CKTsenInfo){
|
|
value->rValue = *(ckt->CKTsenInfo->SEN_Sap[select->iValue + 1]+
|
|
here->DIOsenParmNo);
|
|
}
|
|
return(OK);
|
|
case DIO_QUEST_SENS_REAL:
|
|
if(ckt->CKTsenInfo){
|
|
value->rValue = *(ckt->CKTsenInfo->SEN_RHS[select->iValue + 1]+
|
|
here->DIOsenParmNo);
|
|
}
|
|
return(OK);
|
|
case DIO_QUEST_SENS_IMAG:
|
|
if(ckt->CKTsenInfo){
|
|
value->rValue = *(ckt->CKTsenInfo->SEN_iRHS[select->iValue + 1]+
|
|
here->DIOsenParmNo);
|
|
}
|
|
return(OK);
|
|
case DIO_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->DIOsenParmNo);
|
|
si = *(ckt->CKTsenInfo->SEN_iRHS[select->iValue + 1]+
|
|
here->DIOsenParmNo);
|
|
value->rValue = (vr * sr + vi * si)/vm;
|
|
}
|
|
return(OK);
|
|
case DIO_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->DIOsenParmNo);
|
|
si = *(ckt->CKTsenInfo->SEN_iRHS[select->iValue + 1]+
|
|
here->DIOsenParmNo);
|
|
|
|
value->rValue = (vr * si - vi * sr)/vm;
|
|
}
|
|
return(OK);
|
|
case DIO_QUEST_SENS_CPLX:
|
|
if(ckt->CKTsenInfo){
|
|
value->cValue.real=
|
|
*(ckt->CKTsenInfo->SEN_RHS[select->iValue + 1]+
|
|
here->DIOsenParmNo);
|
|
value->cValue.imag=
|
|
*(ckt->CKTsenInfo->SEN_iRHS[select->iValue + 1]+
|
|
here->DIOsenParmNo);
|
|
}
|
|
return(OK);
|
|
default:
|
|
return(E_BADPARM);
|
|
}
|
|
}
|
|
|