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

#include "ngspice/ngspice.h"
#include "ngspice/ifsim.h"
#include "urcdefs.h"
#include "ngspice/cktdefs.h"
#include "ngspice/gendefs.h"
#include "ngspice/sperror.h"
#include "ngspice/suffix.h"

/* ARGSUSED */
int
URCsetup(SMPmatrix *matrix, GENmodel *inModel, CKTcircuit *ckt, int *state)
        /* create the resistors/capacitors used to model the URC
         */
{
    URCmodel *model = (URCmodel *)inModel;
    URCinstance *here;
    int rtype;
    int ctype;
    int dtype;
    CKTnode * lowl;
    CKTnode * lowr;
    CKTnode * hil;
    CKTnode * hir;
    char *nameelt;
    char *namehi;
    CKTnode *nodehi;
    CKTnode *nodelo;
    char *namelo;
    double r;
    double c;
    IFvalue ptemp;
    double p;
    double r0;
    double c0;
    double i0;
    double r1;
    double c1;
    double i1;
    double rd;
    double wnorm;
    double prop;
    int i;
    GENinstance  *fast;
    GENmodel *modfast; /* capacitor or diode model */
    GENmodel *rmodfast; /* resistor model */
    int error;
    IFuid dioUid;
    IFuid resUid;
    IFuid capUid;
    IFuid eltUid;

    NG_IGNORE(state);
    NG_IGNORE(matrix);

    rtype = CKTtypelook("Resistor");
    ctype = CKTtypelook("Capacitor");
    dtype = CKTtypelook("Diode");
    /*  loop through all the URC models */
    for( ; model != NULL; model = model->URCnextModel ) {
	if(!model->URCkGiven)
	    model->URCk = 1.5;
	if(!model->URCfmaxGiven)
	    model->URCfmax = 1e9;
	if(!model->URCrPerLGiven)
	    model->URCrPerL = 1000;
	if(!model->URCcPerLGiven)
	    model->URCcPerL = 1e-12;

/* may need to put in limits:  k>=1.1, freq <=1e9, rperl >=.1 */

        /* loop through all the instances of the model */
        for (here = model->URCinstances; here != NULL ;
                here=here->URCnextInstance)
	{
            p = model->URCk;
            r0 = here->URClength * model->URCrPerL;
            c0 = here->URClength * model->URCcPerL;
            i0 = here->URClength * model->URCisPerL;
            if(!here->URClumpsGiven) {
                wnorm = model->URCfmax * r0 * c0 * 2.0 * M_PI;
                here->URClumps=(int)MAX(3.0,log(wnorm*(((p-1)/p)*((p-1)/p)))/log(p));
                if(wnorm <35) here->URClumps=3;
                /* may want to limit lumps to <= 100 or something like that */
            }
            r1 = (r0*(p-1))/((2*(pow(p,(double)here->URClumps)))-2);
            c1 = (c0 * (p-1))/((pow(p,(double)(here->URClumps-1)))*(p+1)-2);
            i1 = (i0 * (p-1))/((pow(p,(double)(here->URClumps-1)))*(p+1)-2);
            rd = here->URClength * here->URClumps * model->URCrsPerL;
            /* may want to check that c1 > 0 */
            prop=1;

            if(model->URCisPerLGiven) {
                error = SPfrontEnd->IFnewUid (ckt, &dioUid, here->URCname,
                        "diodemod", UID_MODEL, NULL);
                if(error) return(error);
                modfast = NULL;
                error = CKTmodCrt(ckt,dtype,&modfast,
                        dioUid);
                if(error) return(error);
                ptemp.rValue = c1;
                error= CKTpModName("cjo",&ptemp,ckt,dtype,dioUid,&modfast);
                if(error) return(error);
                ptemp.rValue = rd;
                error = CKTpModName("rs",&ptemp,ckt,dtype,dioUid,&modfast);
                if(error) return(error);
                ptemp.rValue = i1;
                error = CKTpModName("is",&ptemp,ckt,dtype,dioUid,&modfast);
                if(error) return(error);
            } else {
                error = SPfrontEnd->IFnewUid (ckt, &capUid,
                        here->URCname, "capmod", UID_MODEL, NULL);
                if(error) return(error);
                modfast = NULL;
                error = CKTmodCrt(ckt,ctype,&modfast,
                        capUid);
                if(error) return(error);
            }

            error = SPfrontEnd->IFnewUid (ckt, &resUid, here->URCname,
                    "resmod", UID_MODEL, NULL);
            if(error) return(error);
            rmodfast = NULL;
            error = CKTmodCrt(ckt,rtype,&rmodfast,resUid);
            if(error) return(error);
            lowl = CKTnum2nod(ckt,here->URCposNode);
            hir = CKTnum2nod(ckt,here->URCnegNode);
            for(i=1;i<=here->URClumps;i++) {
                namehi = tprintf("hi%d", i);
                error = CKTmkVolt(ckt, &nodehi, here->URCname, namehi);
                if(error) return(error);
                hil = nodehi;
                if(i==here->URClumps) {
                    lowr = hil;
                } else {
                    namelo = tprintf("lo%d", i);
                    error = CKTmkVolt(ckt, &nodelo, here->URCname,
                            namelo);
                    if(error) return(error);
                    lowr = nodelo;
                }
                r = prop*r1;
                c = prop*c1;

                nameelt = tprintf("rlo%d", i);
                error = SPfrontEnd->IFnewUid (ckt, &eltUid, here->URCname,
                        nameelt, UID_INSTANCE, NULL);
                if(error) return(error);
                error = CKTcrtElt(ckt,rmodfast,
                        &fast,eltUid);
                if(error) return(error);
                error = CKTbindNode(ckt,fast,1,lowl);
                if(error) return(error);
                error = CKTbindNode(ckt,fast,2,lowr);
                if(error) return(error);
                ptemp.rValue = r;
                error = CKTpName("resistance",&ptemp,ckt,rtype,nameelt,&fast);
                if(error) return(error);

                nameelt = tprintf("rhi%d", i);
                error = SPfrontEnd->IFnewUid (ckt, &eltUid, here->URCname,
                        nameelt, UID_INSTANCE, NULL);
                if(error) return(error);
                error = CKTcrtElt(ckt,rmodfast,
                        &fast,eltUid);
                if(error) return(error);
                error = CKTbindNode(ckt,fast,1,hil);
                if(error) return(error);
                error = CKTbindNode(ckt,fast,2,hir);
                if(error) return(error);
                ptemp.rValue = r;
                error = CKTpName("resistance",&ptemp,ckt,rtype,nameelt,&fast);
                if(error) return(error);

                if(model->URCisPerLGiven) {
                    /* use diode */
                    nameelt = tprintf("dlo%d", i);
                    error = SPfrontEnd->IFnewUid (ckt, &eltUid,
                            here->URCname,nameelt,UID_INSTANCE, 
                            NULL);
                    if(error) return(error);
                    error = CKTcrtElt(ckt,modfast,
                            &fast, eltUid);
                    if(error) return(error);
                    error = CKTbindNode(ckt,fast,1,lowr);
                    if(error) return(error);
                    error = CKTbindNode(ckt,fast,2,
                            CKTnum2nod(ckt, here->URCgndNode));
                    if(error) return(error);
                    ptemp.rValue = prop;
                    error = CKTpName("area",&ptemp,ckt,dtype,nameelt,&fast);
                    if(error) return(error);
                } else {
                    /* use simple capacitor */
                    nameelt = tprintf("clo%d", i);
                    error = SPfrontEnd->IFnewUid (ckt, &eltUid, here->URCname
                            ,nameelt, UID_INSTANCE, NULL);
                    if(error) return(error);
                    error = CKTcrtElt(ckt,modfast,
                            &fast, eltUid);
                    if(error) return(error);
                    error = CKTbindNode(ckt,fast,1,lowr);
                    if(error) return(error);
                    error = CKTbindNode(ckt,fast,2,
                            CKTnum2nod(ckt, here->URCgndNode));
                    if(error) return(error);
                    ptemp.rValue = c;
                    error = CKTpName("capacitance",&ptemp,ckt,ctype,nameelt,
                            &fast);
                    if(error) return(error);
                }

                if(i!=here->URClumps){
                    if(model->URCisPerLGiven) {
                        /* use diode */
                        nameelt = tprintf("dhi%d", i);
                        error = SPfrontEnd->IFnewUid (ckt, &eltUid,
                                here->URCname,nameelt,UID_INSTANCE,
                                NULL);
                        if(error) return(error);
                        error = CKTcrtElt(ckt,modfast,
                                &fast,eltUid);
                        if(error) return(error);
                        error = CKTbindNode(ckt,fast,1,hil);
                        if(error) return(error);
                        error = CKTbindNode(ckt,fast,2,
                                CKTnum2nod(ckt, here->URCgndNode));
                        if(error) return(error);
                        ptemp.rValue = prop;
                        error=CKTpName("area",&ptemp,ckt,dtype,nameelt,&fast);
                        if(error) return(error);
                    } else {
                        /* use simple capacitor */
                        nameelt = tprintf("chi%d", i);
                        error = SPfrontEnd->IFnewUid (ckt, &eltUid,
                                here->URCname,nameelt,UID_INSTANCE,
                                NULL);
                        if(error) return(error);
                        error = CKTcrtElt(ckt,modfast,
                                &fast,eltUid);
                        if(error) return(error);
                        error = CKTbindNode(ckt,fast,1,hil);
                        if(error) return(error);
                        error = CKTbindNode(ckt,fast,2,
                                CKTnum2nod(ckt, here->URCgndNode));
                        if(error) return(error);
                        ptemp.rValue = c;
                        error =CKTpName("capacitance",&ptemp,ckt,ctype,nameelt,
                                &fast);
                        if(error) return(error);
                    }
                }
                prop *= p;
                lowl = lowr;
                hir = hil;
            }
        }
    }
    return(OK);
}

int
URCunsetup(GENmodel *inModel, CKTcircuit *ckt)
{
    IFuid	varUid;
    int	error;
    URCinstance *here;
    URCmodel *model = (URCmodel *) inModel;
    GENinstance *in;
    GENmodel *modfast;

    /* Delete models, devices, and intermediate nodes; */

    for ( ; model; model = model->URCnextModel) {
	for (here = model->URCinstances; here;
	     here = here->URCnextInstance)
	{
	    if(model->URCisPerLGiven) {
		/* Diodes */
		error = SPfrontEnd->IFnewUid (ckt, &varUid,
						  here->URCname,
						  "diodemod",
						  UID_MODEL,
						  NULL);
	    } else {
		/* Capacitors */
		error = SPfrontEnd->IFnewUid (ckt, &varUid,
						  here->URCname,
						  "capmod",
						  UID_MODEL,
						  NULL);
	    }
	    if (error && error != E_EXISTS)
		return error;

	    modfast = CKTfndMod(ckt, varUid);
	    if (!modfast)
		return E_NOMOD;

	    for (in = modfast->GENinstances; in; in = in->GENnextInstance)
		CKTdltNNum(ckt, in->GENnode1);

	    CKTdltMod(ckt, modfast);	/* Does the elements too */

	    /* Resistors */
	    error = SPfrontEnd->IFnewUid (ckt, &varUid, here->URCname,
					      "resmod", UID_MODEL, NULL);
	    if (error && error != E_EXISTS)
		return error;

	    modfast = CKTfndMod(ckt, varUid);
	    if (!modfast)
		return E_NOMOD;
	    CKTdltMod(ckt, modfast);
	}
    }
    return OK;
}
