mos2ask.c

来自「ngspice又一个电子CAD仿真软件代码.功能更全」· C语言 代码 · 共 438 行 · 第 1/2 页

            value->rValue = *(ckt->CKTstate0 + here->MOS2cqgb);
            return(OK);
        case MOS2_QBD:
            value->rValue = *(ckt->CKTstate0 + here->MOS2qbd);
            return(OK);
        case MOS2_CQBD:
            value->rValue = *(ckt->CKTstate0 + here->MOS2cqbd);
            return(OK);
        case MOS2_QBS:
            value->rValue = *(ckt->CKTstate0 + here->MOS2qbs);
            return(OK);
        case MOS2_CQBS:
            value->rValue = *(ckt->CKTstate0 + here->MOS2cqbs);
            return(OK);
        case MOS2_L_SENS_DC:
            if(ckt->CKTsenInfo){
               value->rValue = *(ckt->CKTsenInfo->SEN_Sap[select->iValue + 1]+
                here->MOS2senParmNo);
            }
            return(OK);
        case MOS2_L_SENS_REAL:
            if(ckt->CKTsenInfo){
               value->rValue = *(ckt->CKTsenInfo->SEN_RHS[select->iValue + 1]+
                here->MOS2senParmNo);
            }
            return(OK);
        case MOS2_L_SENS_IMAG:
            if(ckt->CKTsenInfo){
               value->rValue = *(ckt->CKTsenInfo->SEN_iRHS[select->iValue + 1]+
                here->MOS2senParmNo);
            }
            return(OK);
        case MOS2_L_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->MOS2senParmNo);
                si = *(ckt->CKTsenInfo->SEN_iRHS[select->iValue + 1]+
                        here->MOS2senParmNo);
                value->rValue = (vr * sr + vi * si)/vm;
            }
            return(OK);
        case MOS2_L_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->MOS2senParmNo);
                si = *(ckt->CKTsenInfo->SEN_iRHS[select->iValue + 1]+
                    here->MOS2senParmNo);
                value->rValue =  (vr * si - vi * sr)/vm;
            }
            return(OK);
        case MOS2_L_SENS_CPLX:
            if(ckt->CKTsenInfo){
                value->cValue.real= 
                        *(ckt->CKTsenInfo->SEN_RHS[select->iValue + 1]+
                        here->MOS2senParmNo);
                value->cValue.imag= 
                        *(ckt->CKTsenInfo->SEN_iRHS[select->iValue + 1]+
                        here->MOS2senParmNo);
            }
            return(OK);
        case MOS2_W_SENS_DC:
            if(ckt->CKTsenInfo){
                value->rValue = *(ckt->CKTsenInfo->SEN_Sap[select->iValue + 1]+
                        here->MOS2senParmNo + here->MOS2sens_l);
            }
            return(OK);
        case MOS2_W_SENS_REAL:
            if(ckt->CKTsenInfo){
                value->rValue = *(ckt->CKTsenInfo->SEN_RHS[select->iValue + 1]+
                        here->MOS2senParmNo + here->MOS2sens_l);
            }
            return(OK);
        case MOS2_W_SENS_IMAG:
            if(ckt->CKTsenInfo){
                value->rValue = *(ckt->CKTsenInfo->SEN_iRHS[select->iValue + 1]+
                        here->MOS2senParmNo + here->MOS2sens_l);
            }
            return(OK);
        case MOS2_W_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->MOS2senParmNo + here->MOS2sens_l);
                si = *(ckt->CKTsenInfo->SEN_iRHS[select->iValue + 1]+
                        here->MOS2senParmNo + here->MOS2sens_l);
                        value->rValue = (vr * sr + vi * si)/vm;
            }
            return(OK);
        case MOS2_W_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->MOS2senParmNo + here->MOS2sens_l);
                si = *(ckt->CKTsenInfo->SEN_iRHS[select->iValue + 1]+
                        here->MOS2senParmNo + here->MOS2sens_l);
                value->rValue =  (vr * si - vi * sr)/vm;
            }
            return(OK);
        case MOS2_W_SENS_CPLX:
            if(ckt->CKTsenInfo){
                value->cValue.real= 
                        *(ckt->CKTsenInfo->SEN_RHS[select->iValue + 1]+
                        here->MOS2senParmNo + here->MOS2sens_l);
                value->cValue.imag= 
                        *(ckt->CKTsenInfo->SEN_iRHS[select->iValue + 1]+
                        here->MOS2senParmNo + here->MOS2sens_l);
             }
             return(OK);
        case MOS2_CB :
            if (ckt->CKTcurrentAnalysis & DOING_AC) {
                errMsg = MALLOC(strlen(msg)+1);
                errRtn = "MOS2ask.c";
                strcpy(errMsg,msg);
                return(E_ASKCURRENT);
            } else {
                value->rValue = here->MOS2cbd + here->MOS2cbs - *(ckt->CKTstate0
                        + here->MOS2cqgb);
            }
            return(OK);
        case MOS2_CG :
            if (ckt->CKTcurrentAnalysis & DOING_AC) {
                errMsg = MALLOC(strlen(msg)+1);
                errRtn = "MOS2ask.c";
                strcpy(errMsg,msg);
                return(E_ASKCURRENT);
            } else if (ckt->CKTcurrentAnalysis & (DOING_DCOP | DOING_TRCV)) {
                value->rValue = 0;
            } else if ((ckt->CKTcurrentAnalysis & DOING_TRAN) && 
                   (ckt->CKTmode & MODETRANOP)) {
                value->rValue = 0;
            } else {
                value->rValue = *(ckt->CKTstate0 + here->MOS2cqgb) +
                        *(ckt->CKTstate0 + here->MOS2cqgd) + *(ckt->CKTstate0 + 
                        here->MOS2cqgs);
            }
            return(OK);
        case MOS2_CS :
            if (ckt->CKTcurrentAnalysis & DOING_AC) {
                errMsg = MALLOC(strlen(msg)+1);
                errRtn = "MOS2ask.c";
                strcpy(errMsg,msg);
                return(E_ASKCURRENT);
            } else {
                value->rValue = -here->MOS2cd;
                value->rValue -= here->MOS2cbd + here->MOS2cbs -
                        *(ckt->CKTstate0 + here->MOS2cqgb);
                if ((ckt->CKTcurrentAnalysis & DOING_TRAN) && 
                        !(ckt->CKTmode & MODETRANOP)) {
                    value->rValue -= *(ckt->CKTstate0 + here->MOS2cqgb) + 
                            *(ckt->CKTstate0 + here->MOS2cqgd) + 
                            *(ckt->CKTstate0 + here->MOS2cqgs);
                }
            }
            return(OK);
        case MOS2_POWER :
            if (ckt->CKTcurrentAnalysis & DOING_AC) {
                errMsg = MALLOC(strlen(msg)+1); 
                errRtn = "MOS2ask.c";
                strcpy(errMsg,msg);
                return(E_ASKPOWER);
            } else {
                double temp;

                value->rValue = here->MOS2cd * 
                        *(ckt->CKTrhsOld + here->MOS2dNode);
                value->rValue += (here->MOS2cbd + here->MOS2cbs -
                        *(ckt->CKTstate0 + here->MOS2cqgb)) * 
                        *(ckt->CKTrhsOld + here->MOS2bNode);
                if ((ckt->CKTcurrentAnalysis & DOING_TRAN) && 
                        !(ckt->CKTmode & MODETRANOP)) {
                    value->rValue += (*(ckt->CKTstate0 + here->MOS2cqgb) + 
                            *(ckt->CKTstate0 + here->MOS2cqgd) + 
                            *(ckt->CKTstate0 + here->MOS2cqgs)) *
                            *(ckt->CKTrhsOld + here->MOS2gNode);
                }
                temp = -here->MOS2cd;
                temp -= here->MOS2cbd + here->MOS2cbs ;
                if ((ckt->CKTcurrentAnalysis & DOING_TRAN) && 
                        !(ckt->CKTmode & MODETRANOP)) {
                    temp -= *(ckt->CKTstate0 + here->MOS2cqgb) + 
                            *(ckt->CKTstate0 + here->MOS2cqgd) + 
                            *(ckt->CKTstate0 + here->MOS2cqgs);
                }
                value->rValue += temp * *(ckt->CKTrhsOld + here->MOS2sNode);
            }
            return(OK);
        default:
            return(E_BADPARM);
    }
    /* NOTREACHED */
}