示例#1
0
        public CALCFC(CommonBlock Default)
        {
            #region Common varaible Initialization

            #region Common Block: Default Initialization

            this._default = Default;
            NPROB = Default.intData[0];
            #endregion
            #endregion
        }
示例#2
0
        public CALCFC()
        {
            #region Initialization Common Blocks

            CommonBlock Default = new CommonBlock(0, 1, 0, 0);
            #endregion
            #region Common varaible Initialization

            #region Common Block: Default Initialization

            this._default = Default;
            NPROB = Default.intData[0];
            #endregion
            #endregion
        }
示例#3
0
        public CONTR5(CommonBlock CONRA5)
        {
            #region Common varaible Initialization

            #region Common Block: CONRA5 Initialization

            this._conra5 = CONRA5;
            NN = CONRA5.intData[0];
            NN2 = CONRA5.intData[1];
            NN3 = CONRA5.intData[2];
            NN4 = CONRA5.intData[3];
            XSOL = CONRA5.doubleData[0];
            HSOL = CONRA5.doubleData[1];
            C2M1 = CONRA5.doubleData[2];
            C1M1 = CONRA5.doubleData[3];
            #endregion
            #endregion
        }
示例#4
0
        public DECOMC(DECC decc, DECBC decbc, DECHC dechc, CommonBlock LINAL)
        {
            #region Set Dependencies

            this._decc = decc; this._decbc = decbc; this._dechc = dechc;
            #endregion
            #region Common varaible Initialization

            #region Common Block: LINAL Initialization

            this._linal = LINAL;
            MLE = LINAL.intData[0];
            MUE = LINAL.intData[1];
            MBJAC = LINAL.intData[2];
            MBB = LINAL.intData[3];
            MDIAG = LINAL.intData[4];
            MDIFF = LINAL.intData[5];
            MBDIAG = LINAL.intData[6];
            #endregion
            #endregion
        }
示例#5
0
        public CommonBlock(int DimensionDouble, int DimensionInt, int DimensionSingle, int DimensionComplex)
        {
            //this.MeName = TheName;
            if (DimensionDouble > 0)
            {
                this.MeDoubleData = new Odouble[DimensionDouble];

                for (int i = 0; i < DimensionDouble; i++)
                {
                    MeDoubleData[i] = new Odouble();
                }
            }
            if (DimensionInt > 0)
            {
                this.MeIntData = new Oint[DimensionInt];

                for (int i = 0; i < DimensionInt; i++)
                {
                    MeIntData[i] = new Oint();
                }
            }
        }
示例#6
0
        public CommonBlock(int DimensionDouble, int DimensionInt, int DimensionSingle, int DimensionComplex)
        {
            //this.MeName = TheName;
            if (DimensionDouble > 0)
            {
                this.MeDoubleData = new Odouble[DimensionDouble];

                for (int i = 0; i < DimensionDouble; i++)
                {
                    MeDoubleData[i] = new Odouble();
                }
            }
            if (DimensionInt > 0)
            {
                this.MeIntData = new Oint[DimensionInt];

                for (int i = 0; i < DimensionInt; i++)
                {
                    MeIntData[i] = new Oint();
                }
            }
        }
示例#7
0
文件: dvode.cs 项目: davidsiaw/neuron
        public DVODE(DVNLSD dvnlsd, DUMACH dumach, DVNORM dvnorm, DCOPY dcopy, DEWSET dewset, DVHIN dvhin, DSCAL dscal, DVINDY dvindy, XERRWD xerrwd, DVSTEP dvstep
            , CommonBlock DVOD01, CommonBlock DVOD02)
        {
            #region Set Dependencies

            this._dvnlsd = dvnlsd; this._dumach = dumach; this._dvnorm = dvnorm; this._dcopy = dcopy; this._dewset = dewset;
            this._dvhin = dvhin;this._dscal = dscal; this._dvindy = dvindy; this._xerrwd = xerrwd; this._dvstep = dvstep;
            #endregion
            #region Data Initialization

            //MORD(1)/12
            MORD[1 + o_mord] = 12;
            //MORD(2)/5
            MORD[2 + o_mord] = 5;
            //MXSTP0/500
            MXSTP0 = 500;
            //MXHNL0/10
            MXHNL0 = 10;
            //ZERO/0.0D0
            ZERO = 0.0E0;
            //ONE/1.0D0
            ONE = 1.0E0;
            //TWO/2.0D0
            TWO = 2.0E0;
            //FOUR/4.0D0
            FOUR = 4.0E0;
            //PT2/0.2D0
            PT2 = 0.2E0;
            //HUN/100.0D0
            HUN = 100.0E0;
            #endregion
            #region Common varaible Initialization

            #region Common Block: DVOD01 Initialization

            this._dvod01 = DVOD01;
            ACNRM = DVOD01.doubleData[0];
            CCMXJ = DVOD01.doubleData[1];
            CONP = DVOD01.doubleData[2];
            CRATE = DVOD01.doubleData[3];
            DRC = DVOD01.doubleData[4];
            //Start Array: EL  StartIndex: 5  LastIndex: 17
            EL = DVOD01.doubleData;
            o_el = 4;  //o_ = StartIndex-1
            //End Array: EL
            ETA = DVOD01.doubleData[18];
            ETAMAX = DVOD01.doubleData[19];
            H = DVOD01.doubleData[20];
            HMIN = DVOD01.doubleData[21];
            HMXI = DVOD01.doubleData[22];
            HNEW = DVOD01.doubleData[23];
            HSCAL = DVOD01.doubleData[24];
            PRL1 = DVOD01.doubleData[25];
            RC = DVOD01.doubleData[26];
            RL1 = DVOD01.doubleData[27];
            //Start Array: TAU  StartIndex: 28  LastIndex: 40
            TAU = DVOD01.doubleData;
            o_tau = 27;  //o_ = StartIndex-1
            //End Array: TAU
            //Start Array: TQ  StartIndex: 41  LastIndex: 45
            TQ = DVOD01.doubleData;
            o_tq = 40;  //o_ = StartIndex-1
            //End Array: TQ
            TN = DVOD01.doubleData[46];
            UROUND = DVOD01.doubleData[47];
            ICF = DVOD01.intData[0];
            INIT = DVOD01.intData[1];
            IPUP = DVOD01.intData[2];
            JCUR = DVOD01.intData[3];
            JSTART = DVOD01.intData[4];
            JSV = DVOD01.intData[5];
            KFLAG = DVOD01.intData[6];
            KUTH = DVOD01.intData[7];
            L = DVOD01.intData[8];
            LMAX = DVOD01.intData[9];
            LYH = DVOD01.intData[10];
            LEWT = DVOD01.intData[11];
            LACOR = DVOD01.intData[12];
            LSAVF = DVOD01.intData[13];
            LWM = DVOD01.intData[14];
            LIWM = DVOD01.intData[15];
            LOCJS = DVOD01.intData[16];
            MAXORD = DVOD01.intData[17];
            METH = DVOD01.intData[18];
            MITER = DVOD01.intData[19];
            MSBJ = DVOD01.intData[20];
            MXHNIL = DVOD01.intData[21];
            MXSTEP = DVOD01.intData[22];
            N = DVOD01.intData[23];
            NEWH = DVOD01.intData[24];
            NEWQ = DVOD01.intData[25];
            NHNIL = DVOD01.intData[26];
            NQ = DVOD01.intData[27];
            NQNYH = DVOD01.intData[28];
            NQWAIT = DVOD01.intData[29];
            NSLJ = DVOD01.intData[30];
            NSLP = DVOD01.intData[31];
            NYH = DVOD01.intData[32];
            #endregion
            #region Common Block: DVOD02 Initialization

            this._dvod02 = DVOD02;
            HU = DVOD02.doubleData[0];
            NCFN = DVOD02.intData[0];
            NETF = DVOD02.intData[1];
            NFE = DVOD02.intData[2];
            NJE = DVOD02.intData[3];
            NLU = DVOD02.intData[4];
            NNI = DVOD02.intData[5];
            NQU = DVOD02.intData[6];
            NST = DVOD02.intData[7];
            #endregion
            #endregion
        }
示例#8
0
文件: dvode.cs 项目: davidsiaw/neuron
        public DVJAC(DVNORM dvnorm, DCOPY dcopy, DSCAL dscal, DGEFA dgefa, DACOPY dacopy, DGBFA dgbfa, CommonBlock DVOD01, CommonBlock DVOD02)
        {
            #region Set Dependencies

            this._dvnorm = dvnorm; this._dcopy = dcopy; this._dscal = dscal; this._dgefa = dgefa; this._dacopy = dacopy;
            this._dgbfa = dgbfa;
            #endregion
            #region Data Initialization

            //ONE/1.0D0
            ONE = 1.0E0;
            //THOU/1000.0D0
            THOU = 1000.0E0;
            //ZERO/0.0D0
            ZERO = 0.0E0;
            //PT1/0.1D0
            PT1 = 0.1E0;
            #endregion
            #region Common varaible Initialization

            #region Common Block: DVOD01 Initialization

            this._dvod01 = DVOD01;
            ACNRM = DVOD01.doubleData[0];
            CCMXJ = DVOD01.doubleData[1];
            CONP = DVOD01.doubleData[2];
            CRATE = DVOD01.doubleData[3];
            DRC = DVOD01.doubleData[4];
            //Start Array: EL  StartIndex: 5  LastIndex: 17
            EL = DVOD01.doubleData;
            o_el = 4;  //o_ = StartIndex-1
            //End Array: EL
            ETA = DVOD01.doubleData[18];
            ETAMAX = DVOD01.doubleData[19];
            H = DVOD01.doubleData[20];
            HMIN = DVOD01.doubleData[21];
            HMXI = DVOD01.doubleData[22];
            HNEW = DVOD01.doubleData[23];
            HSCAL = DVOD01.doubleData[24];
            PRL1 = DVOD01.doubleData[25];
            RC = DVOD01.doubleData[26];
            RL1 = DVOD01.doubleData[27];
            //Start Array: TAU  StartIndex: 28  LastIndex: 40
            TAU = DVOD01.doubleData;
            o_tau = 27;  //o_ = StartIndex-1
            //End Array: TAU
            //Start Array: TQ  StartIndex: 41  LastIndex: 45
            TQ = DVOD01.doubleData;
            o_tq = 40;  //o_ = StartIndex-1
            //End Array: TQ
            TN = DVOD01.doubleData[46];
            UROUND = DVOD01.doubleData[47];
            ICF = DVOD01.intData[0];
            INIT = DVOD01.intData[1];
            IPUP = DVOD01.intData[2];
            JCUR = DVOD01.intData[3];
            JSTART = DVOD01.intData[4];
            JSV = DVOD01.intData[5];
            KFLAG = DVOD01.intData[6];
            KUTH = DVOD01.intData[7];
            L = DVOD01.intData[8];
            LMAX = DVOD01.intData[9];
            LYH = DVOD01.intData[10];
            LEWT = DVOD01.intData[11];
            LACOR = DVOD01.intData[12];
            LSAVF = DVOD01.intData[13];
            LWM = DVOD01.intData[14];
            LIWM = DVOD01.intData[15];
            LOCJS = DVOD01.intData[16];
            MAXORD = DVOD01.intData[17];
            METH = DVOD01.intData[18];
            MITER = DVOD01.intData[19];
            MSBJ = DVOD01.intData[20];
            MXHNIL = DVOD01.intData[21];
            MXSTEP = DVOD01.intData[22];
            N = DVOD01.intData[23];
            NEWH = DVOD01.intData[24];
            NEWQ = DVOD01.intData[25];
            NHNIL = DVOD01.intData[26];
            NQ = DVOD01.intData[27];
            NQNYH = DVOD01.intData[28];
            NQWAIT = DVOD01.intData[29];
            NSLJ = DVOD01.intData[30];
            NSLP = DVOD01.intData[31];
            NYH = DVOD01.intData[32];
            #endregion
            #region Common Block: DVOD02 Initialization

            this._dvod02 = DVOD02;
            HU = DVOD02.doubleData[0];
            NCFN = DVOD02.intData[0];
            NETF = DVOD02.intData[1];
            NFE = DVOD02.intData[2];
            NJE = DVOD02.intData[3];
            NLU = DVOD02.intData[4];
            NNI = DVOD02.intData[5];
            NQU = DVOD02.intData[6];
            NST = DVOD02.intData[7];
            #endregion
            #endregion
        }
示例#9
0
        public COBYLATEST(COBYLA cobyla, CommonBlock Default)
        {
            #region Set Dependencies

            this._cobyla = cobyla;
            #endregion
            #region Common varaible Initialization

            #region Common Block: Default Initialization

            this._default = Default;
            NPROB = Default.intData[0];
            #endregion
            #endregion
        }
示例#10
0
文件: TN.cs 项目: davidsiaw/neuron
        public MSOLVE()
        {
            #region Initialization Common Blocks

            CommonBlock SUBSCR = new CommonBlock(0, 15, 0, 0);
            #endregion
            #region Dependencies (Initialization)

            DDOT ddot = new DDOT();
            SSBFGS ssbfgs = new SSBFGS();
            MSLV mslv = new MSLV(ddot, ssbfgs);
            #endregion
            #region Set Dependencies

            this._mslv = mslv;
            #endregion
            #region Common varaible Initialization

            #region Common Block: SUBSCR Initialization

            this._subscr = SUBSCR;
            LGV = SUBSCR.intData[0];
            LZ1 = SUBSCR.intData[1];
            LZK = SUBSCR.intData[2];
            LV = SUBSCR.intData[3];
            LSK = SUBSCR.intData[4];
            LYK = SUBSCR.intData[5];
            LDIAGB = SUBSCR.intData[6];
            LSR = SUBSCR.intData[7];
            LYR = SUBSCR.intData[8];
            LHYR = SUBSCR.intData[9];
            LHG = SUBSCR.intData[10];
            LHYK = SUBSCR.intData[11];
            LPK = SUBSCR.intData[12];
            LEMAT = SUBSCR.intData[13];
            LWTEST = SUBSCR.intData[14];
            #endregion
            #endregion
        }
示例#11
0
文件: TN.cs 项目: davidsiaw/neuron
        public LMQNBC()
        {
            #region Initialization Common Blocks

            CommonBlock SUBSCR = new CommonBlock(0, 15, 0, 0);
            #endregion
            #region Dependencies (Initialization)

            DDOT ddot = new DDOT();
            DNRM2 dnrm2 = new DNRM2();
            MCHPR1 mchpr1 = new MCHPR1();
            CRASH crash = new CRASH();
            ZTIME ztime = new ZTIME();
            MONIT monit = new MONIT();
            DCOPY dcopy = new DCOPY();
            SSBFGS ssbfgs = new SSBFGS();
            DAXPY daxpy = new DAXPY();
            NEGVEC negvec = new NEGVEC();
            STPMAX stpmax = new STPMAX();
            GETPTC getptc = new GETPTC();
            LSOUT lsout = new LSOUT();
            MODZ modz = new MODZ();
            CNVTST cnvtst = new CNVTST();
            DXPY dxpy = new DXPY();
            SETPAR setpar = new SETPAR(SUBSCR);
            GTIMS gtims = new GTIMS(SUBSCR);
            STEP1 step1 = new STEP1(mchpr1);
            CHKUCP chkucp = new CHKUCP(mchpr1, dnrm2);
            SETUCR setucr = new SETUCR(ddot);
            INITP3 initp3 = new INITP3(ddot, dcopy);
            INITPC initpc = new INITPC(initp3, SUBSCR);
            MSLV mslv = new MSLV(ddot, ssbfgs);
            MSOLVE msolve = new MSOLVE(mslv, SUBSCR);
            NDIA3 ndia3 = new NDIA3(ddot);
            MODLNP modlnp = new MODLNP(ddot, initpc, ztime, msolve, gtims, ndia3, daxpy, negvec, dcopy);
            LINDER linder = new LINDER(ddot, getptc, lsout, dcopy);
            #endregion
            #region Set Dependencies

            this._ddot = ddot; this._dnrm2 = dnrm2; this._step1 = step1; this._crash = crash; this._setpar = setpar;
            this._chkucp = chkucp;this._setucr = setucr; this._ztime = ztime; this._monit = monit; this._modlnp = modlnp;
            this._dcopy = dcopy;this._stpmax = stpmax; this._linder = linder; this._modz = modz; this._cnvtst = cnvtst;
            this._dxpy = dxpy;
            #endregion
            #region Common varaible Initialization

            #region Common Block: SUBSCR Initialization

            this._subscr = SUBSCR;
            LGV = SUBSCR.intData[0];
            LZ1 = SUBSCR.intData[1];
            LZK = SUBSCR.intData[2];
            LV = SUBSCR.intData[3];
            LSK = SUBSCR.intData[4];
            LYK = SUBSCR.intData[5];
            LDIAGB = SUBSCR.intData[6];
            LSR = SUBSCR.intData[7];
            LYR = SUBSCR.intData[8];
            LOLDG = SUBSCR.intData[9];
            LHG = SUBSCR.intData[10];
            LHYK = SUBSCR.intData[11];
            LPK = SUBSCR.intData[12];
            LEMAT = SUBSCR.intData[13];
            LWTEST = SUBSCR.intData[14];
            #endregion
            #endregion
        }
示例#12
0
文件: TN.cs 项目: davidsiaw/neuron
        public LMQN(DNRM2 dnrm2, STEP1 step1, DDOT ddot, SETPAR setpar, CHKUCP chkucp, SETUCR setucr, MODLNP modlnp, DCOPY dcopy, LINDER linder, DXPY dxpy
            , CommonBlock SUBSCR)
        {
            #region Set Dependencies

            this._dnrm2 = dnrm2; this._step1 = step1; this._ddot = ddot; this._setpar = setpar; this._chkucp = chkucp;
            this._setucr = setucr;this._modlnp = modlnp; this._dcopy = dcopy; this._linder = linder; this._dxpy = dxpy;
            #endregion
            #region Common varaible Initialization

            #region Common Block: SUBSCR Initialization

            this._subscr = SUBSCR;
            LGV = SUBSCR.intData[0];
            LZ1 = SUBSCR.intData[1];
            LZK = SUBSCR.intData[2];
            LV = SUBSCR.intData[3];
            LSK = SUBSCR.intData[4];
            LYK = SUBSCR.intData[5];
            LDIAGB = SUBSCR.intData[6];
            LSR = SUBSCR.intData[7];
            LYR = SUBSCR.intData[8];
            LOLDG = SUBSCR.intData[9];
            LHG = SUBSCR.intData[10];
            LHYK = SUBSCR.intData[11];
            LPK = SUBSCR.intData[12];
            LEMAT = SUBSCR.intData[13];
            LWTEST = SUBSCR.intData[14];
            #endregion
            #endregion
        }
示例#13
0
        public SLVRAI()
        {
            #region Initialization Common Blocks

            CommonBlock LINAL = new CommonBlock(0, 7, 0, 0);
            #endregion
            #region Dependencies (Initialization)

            SOLC solc = new SOLC();
            SOLBC solbc = new SOLBC();
            SOLHC solhc = new SOLHC();
            #endregion
            #region Set Dependencies

            this._solc = solc; this._solbc = solbc; this._solhc = solhc;
            #endregion
            #region Common varaible Initialization

            #region Common Block: LINAL Initialization

            this._linal = LINAL;
            MLE = LINAL.intData[0];
            MUE = LINAL.intData[1];
            MBJAC = LINAL.intData[2];
            MBB = LINAL.intData[3];
            MDIAG = LINAL.intData[4];
            MDIFF = LINAL.intData[5];
            MBDIAG = LINAL.intData[6];
            #endregion
            #endregion
        }
示例#14
0
        public SLVRAI(SOLC solc, SOLBC solbc, SOLHC solhc, CommonBlock LINAL)
        {
            #region Set Dependencies

            this._solc = solc; this._solbc = solbc; this._solhc = solhc;
            #endregion
            #region Common varaible Initialization

            #region Common Block: LINAL Initialization

            this._linal = LINAL;
            MLE = LINAL.intData[0];
            MUE = LINAL.intData[1];
            MBJAC = LINAL.intData[2];
            MBB = LINAL.intData[3];
            MDIAG = LINAL.intData[4];
            MDIFF = LINAL.intData[5];
            MBDIAG = LINAL.intData[6];
            #endregion
            #endregion
        }
示例#15
0
文件: dvode.cs 项目: davidsiaw/neuron
        public DVSTEP()
        {
            #region Initialization Common Blocks

            CommonBlock DVOD01 = new CommonBlock(48, 33, 0, 0);
            CommonBlock DVOD02 = new CommonBlock(1, 8, 0, 0);
            #endregion
            #region Dependencies (Initialization)

            DVNORM dvnorm = new DVNORM();
            DAXPY daxpy = new DAXPY();
            DSCAL dscal = new DSCAL();
            DCOPY dcopy = new DCOPY();
            DVSET dvset = new DVSET(DVOD01);
            DVJUST dvjust = new DVJUST(daxpy, DVOD01);
            #endregion
            #region Set Dependencies

            this._dvnorm = dvnorm; this._dvjust = dvjust; this._dscal = dscal; this._dvset = dvset; this._daxpy = daxpy;
            this._dcopy = dcopy;
            #endregion
            #region Data Initialization

            //KFC/-3
            KFC =  - 3;
            //KFH/-7
            KFH =  - 7;
            //MXNCF/10
            MXNCF = 10;
            //ADDON/1.0D-6
            ADDON = 1.0E-6;
            //BIAS1/6.0D0
            BIAS1 = 6.0E0;
            //BIAS2/6.0D0
            BIAS2 = 6.0E0;
            //BIAS3/10.0D0
            BIAS3 = 10.0E0;
            //ETACF/0.25D0
            ETACF = 0.25E0;
            //ETAMIN/0.1D0
            ETAMIN = 0.1E0;
            //ETAMXF/0.2D0
            ETAMXF = 0.2E0;
            //ETAMX1/1.0D4
            ETAMX1 = 1.0E4;
            //ETAMX2/10.0D0
            ETAMX2 = 10.0E0;
            //ETAMX3/10.0D0
            ETAMX3 = 10.0E0;
            //ONEPSM/1.00001D0
            ONEPSM = 1.00001E0;
            //THRESH/1.5D0
            THRESH = 1.5E0;
            //ONE/1.0D0
            ONE = 1.0E0;
            //ZERO/0.0D0
            ZERO = 0.0E0;
            #endregion
            #region Common varaible Initialization

            #region Common Block: DVOD01 Initialization

            this._dvod01 = DVOD01;
            ACNRM = DVOD01.doubleData[0];
            CCMXJ = DVOD01.doubleData[1];
            CONP = DVOD01.doubleData[2];
            CRATE = DVOD01.doubleData[3];
            DRC = DVOD01.doubleData[4];
            //Start Array: EL  StartIndex: 5  LastIndex: 17
            EL = DVOD01.doubleData;
            o_el = 4;  //o_ = StartIndex-1
            //End Array: EL
            ETA = DVOD01.doubleData[18];
            ETAMAX = DVOD01.doubleData[19];
            H = DVOD01.doubleData[20];
            HMIN = DVOD01.doubleData[21];
            HMXI = DVOD01.doubleData[22];
            HNEW = DVOD01.doubleData[23];
            HSCAL = DVOD01.doubleData[24];
            PRL1 = DVOD01.doubleData[25];
            RC = DVOD01.doubleData[26];
            RL1 = DVOD01.doubleData[27];
            //Start Array: TAU  StartIndex: 28  LastIndex: 40
            TAU = DVOD01.doubleData;
            o_tau = 27;  //o_ = StartIndex-1
            //End Array: TAU
            //Start Array: TQ  StartIndex: 41  LastIndex: 45
            TQ = DVOD01.doubleData;
            o_tq = 40;  //o_ = StartIndex-1
            //End Array: TQ
            TN = DVOD01.doubleData[46];
            UROUND = DVOD01.doubleData[47];
            ICF = DVOD01.intData[0];
            INIT = DVOD01.intData[1];
            IPUP = DVOD01.intData[2];
            JCUR = DVOD01.intData[3];
            JSTART = DVOD01.intData[4];
            JSV = DVOD01.intData[5];
            KFLAG = DVOD01.intData[6];
            KUTH = DVOD01.intData[7];
            L = DVOD01.intData[8];
            LMAX = DVOD01.intData[9];
            LYH = DVOD01.intData[10];
            LEWT = DVOD01.intData[11];
            LACOR = DVOD01.intData[12];
            LSAVF = DVOD01.intData[13];
            LWM = DVOD01.intData[14];
            LIWM = DVOD01.intData[15];
            LOCJS = DVOD01.intData[16];
            MAXORD = DVOD01.intData[17];
            METH = DVOD01.intData[18];
            MITER = DVOD01.intData[19];
            MSBJ = DVOD01.intData[20];
            MXHNIL = DVOD01.intData[21];
            MXSTEP = DVOD01.intData[22];
            N = DVOD01.intData[23];
            NEWH = DVOD01.intData[24];
            NEWQ = DVOD01.intData[25];
            NHNIL = DVOD01.intData[26];
            NQ = DVOD01.intData[27];
            NQNYH = DVOD01.intData[28];
            NQWAIT = DVOD01.intData[29];
            NSLJ = DVOD01.intData[30];
            NSLP = DVOD01.intData[31];
            NYH = DVOD01.intData[32];
            #endregion
            #region Common Block: DVOD02 Initialization

            this._dvod02 = DVOD02;
            HU = DVOD02.doubleData[0];
            NCFN = DVOD02.intData[0];
            NETF = DVOD02.intData[1];
            NFE = DVOD02.intData[2];
            NJE = DVOD02.intData[3];
            NLU = DVOD02.intData[4];
            NNI = DVOD02.intData[5];
            NQU = DVOD02.intData[6];
            NST = DVOD02.intData[7];
            #endregion
            #endregion
        }
示例#16
0
文件: dvode.cs 项目: davidsiaw/neuron
        public DVSOL()
        {
            #region Initialization Common Blocks

            CommonBlock DVOD01 = new CommonBlock(48, 33, 0, 0);
            #endregion
            #region Dependencies (Initialization)

            DDOT ddot = new DDOT();
            DAXPY daxpy = new DAXPY();
            DGESL dgesl = new DGESL(ddot, daxpy);
            DGBSL dgbsl = new DGBSL(ddot, daxpy);
            #endregion
            #region Set Dependencies

            this._dgesl = dgesl; this._dgbsl = dgbsl;
            #endregion
            #region Data Initialization

            //ONE/1.0D0
            ONE = 1.0E0;
            //ZERO/0.0D0
            ZERO = 0.0E0;
            #endregion
            #region Common varaible Initialization

            #region Common Block: DVOD01 Initialization

            this._dvod01 = DVOD01;
            ACNRM = DVOD01.doubleData[0];
            CCMXJ = DVOD01.doubleData[1];
            CONP = DVOD01.doubleData[2];
            CRATE = DVOD01.doubleData[3];
            DRC = DVOD01.doubleData[4];
            //Start Array: EL  StartIndex: 5  LastIndex: 17
            EL = DVOD01.doubleData;
            o_el = 4;  //o_ = StartIndex-1
            //End Array: EL
            ETA = DVOD01.doubleData[18];
            ETAMAX = DVOD01.doubleData[19];
            H = DVOD01.doubleData[20];
            HMIN = DVOD01.doubleData[21];
            HMXI = DVOD01.doubleData[22];
            HNEW = DVOD01.doubleData[23];
            HSCAL = DVOD01.doubleData[24];
            PRL1 = DVOD01.doubleData[25];
            RC = DVOD01.doubleData[26];
            RL1 = DVOD01.doubleData[27];
            //Start Array: TAU  StartIndex: 28  LastIndex: 40
            TAU = DVOD01.doubleData;
            o_tau = 27;  //o_ = StartIndex-1
            //End Array: TAU
            //Start Array: TQ  StartIndex: 41  LastIndex: 45
            TQ = DVOD01.doubleData;
            o_tq = 40;  //o_ = StartIndex-1
            //End Array: TQ
            TN = DVOD01.doubleData[46];
            UROUND = DVOD01.doubleData[47];
            ICF = DVOD01.intData[0];
            INIT = DVOD01.intData[1];
            IPUP = DVOD01.intData[2];
            JCUR = DVOD01.intData[3];
            JSTART = DVOD01.intData[4];
            JSV = DVOD01.intData[5];
            KFLAG = DVOD01.intData[6];
            KUTH = DVOD01.intData[7];
            L = DVOD01.intData[8];
            LMAX = DVOD01.intData[9];
            LYH = DVOD01.intData[10];
            LEWT = DVOD01.intData[11];
            LACOR = DVOD01.intData[12];
            LSAVF = DVOD01.intData[13];
            LWM = DVOD01.intData[14];
            LIWM = DVOD01.intData[15];
            LOCJS = DVOD01.intData[16];
            MAXORD = DVOD01.intData[17];
            METH = DVOD01.intData[18];
            MITER = DVOD01.intData[19];
            MSBJ = DVOD01.intData[20];
            MXHNIL = DVOD01.intData[21];
            MXSTEP = DVOD01.intData[22];
            N = DVOD01.intData[23];
            NEWH = DVOD01.intData[24];
            NEWQ = DVOD01.intData[25];
            NHNIL = DVOD01.intData[26];
            NQ = DVOD01.intData[27];
            NQNYH = DVOD01.intData[28];
            NQWAIT = DVOD01.intData[29];
            NSLJ = DVOD01.intData[30];
            NSLP = DVOD01.intData[31];
            NYH = DVOD01.intData[32];
            #endregion
            #endregion
        }
示例#17
0
文件: dvode.cs 项目: davidsiaw/neuron
        public DVSET(CommonBlock DVOD01)
        {
            #region Data Initialization

            //CORTES/0.1D0
            CORTES = 0.1E0;
            //ONE/1.0D0
            ONE = 1.0E0;
            //SIX/6.0D0
            SIX = 6.0E0;
            //TWO/2.0D0
            TWO = 2.0E0;
            //ZERO/0.0D0
            ZERO = 0.0E0;
            #endregion
            #region Common varaible Initialization

            #region Common Block: DVOD01 Initialization

            this._dvod01 = DVOD01;
            ACNRM = DVOD01.doubleData[0];
            CCMXJ = DVOD01.doubleData[1];
            CONP = DVOD01.doubleData[2];
            CRATE = DVOD01.doubleData[3];
            DRC = DVOD01.doubleData[4];
            //Start Array: EL  StartIndex: 5  LastIndex: 17
            EL = DVOD01.doubleData;
            o_el = 4;  //o_ = StartIndex-1
            //End Array: EL
            ETA = DVOD01.doubleData[18];
            ETAMAX = DVOD01.doubleData[19];
            H = DVOD01.doubleData[20];
            HMIN = DVOD01.doubleData[21];
            HMXI = DVOD01.doubleData[22];
            HNEW = DVOD01.doubleData[23];
            HSCAL = DVOD01.doubleData[24];
            PRL1 = DVOD01.doubleData[25];
            RC = DVOD01.doubleData[26];
            RL1 = DVOD01.doubleData[27];
            //Start Array: TAU  StartIndex: 28  LastIndex: 40
            TAU = DVOD01.doubleData;
            o_tau = 27;  //o_ = StartIndex-1
            //End Array: TAU
            //Start Array: TQ  StartIndex: 41  LastIndex: 45
            TQ = DVOD01.doubleData;
            o_tq = 40;  //o_ = StartIndex-1
            //End Array: TQ
            TN = DVOD01.doubleData[46];
            UROUND = DVOD01.doubleData[47];
            ICF = DVOD01.intData[0];
            INIT = DVOD01.intData[1];
            IPUP = DVOD01.intData[2];
            JCUR = DVOD01.intData[3];
            JSTART = DVOD01.intData[4];
            JSV = DVOD01.intData[5];
            KFLAG = DVOD01.intData[6];
            KUTH = DVOD01.intData[7];
            L = DVOD01.intData[8];
            LMAX = DVOD01.intData[9];
            LYH = DVOD01.intData[10];
            LEWT = DVOD01.intData[11];
            LACOR = DVOD01.intData[12];
            LSAVF = DVOD01.intData[13];
            LWM = DVOD01.intData[14];
            LIWM = DVOD01.intData[15];
            LOCJS = DVOD01.intData[16];
            MAXORD = DVOD01.intData[17];
            METH = DVOD01.intData[18];
            MITER = DVOD01.intData[19];
            MSBJ = DVOD01.intData[20];
            MXHNIL = DVOD01.intData[21];
            MXSTEP = DVOD01.intData[22];
            N = DVOD01.intData[23];
            NEWH = DVOD01.intData[24];
            NEWQ = DVOD01.intData[25];
            NHNIL = DVOD01.intData[26];
            NQ = DVOD01.intData[27];
            NQNYH = DVOD01.intData[28];
            NQWAIT = DVOD01.intData[29];
            NSLJ = DVOD01.intData[30];
            NSLP = DVOD01.intData[31];
            NYH = DVOD01.intData[32];
            #endregion
            #endregion
        }
示例#18
0
文件: dvode.cs 项目: davidsiaw/neuron
        public DVODE()
        {
            #region Initialization Common Blocks

            CommonBlock DVOD01 = new CommonBlock(48, 33, 0, 0);
            CommonBlock DVOD02 = new CommonBlock(1, 8, 0, 0);
            #endregion
            #region Dependencies (Initialization)

            DVNORM dvnorm = new DVNORM();
            DCOPY dcopy = new DCOPY();
            DSCAL dscal = new DSCAL();
            IDAMAX idamax = new IDAMAX();
            DAXPY daxpy = new DAXPY();
            DDOT ddot = new DDOT();
            DUMSUM dumsum = new DUMSUM();
            DEWSET dewset = new DEWSET();
            IUMACH iumach = new IUMACH();
            DVSET dvset = new DVSET(DVOD01);
            DGEFA dgefa = new DGEFA(idamax, dscal, daxpy);
            DACOPY dacopy = new DACOPY(dcopy);
            DGBFA dgbfa = new DGBFA(idamax, dscal, daxpy);
            DVJAC dvjac = new DVJAC(dvnorm, dcopy, dscal, dgefa, dacopy, dgbfa, DVOD01, DVOD02);
            DGESL dgesl = new DGESL(ddot, daxpy);
            DGBSL dgbsl = new DGBSL(ddot, daxpy);
            DVSOL dvsol = new DVSOL(dgesl, dgbsl, DVOD01);
            DVNLSD dvnlsd = new DVNLSD(dvnorm, dcopy, dvjac, dvsol, dscal, daxpy, DVOD01, DVOD02);
            DUMACH dumach = new DUMACH(dumsum);
            DVHIN dvhin = new DVHIN(dvnorm);
            IXSAV ixsav = new IXSAV(iumach);
            XERRWD xerrwd = new XERRWD(ixsav);
            DVINDY dvindy = new DVINDY(dscal, xerrwd, DVOD01, DVOD02);
            DVJUST dvjust = new DVJUST(daxpy, DVOD01);
            DVSTEP dvstep = new DVSTEP(dvnorm, dvjust, dscal, dvset, daxpy, dcopy, DVOD01, DVOD02);
            #endregion
            #region Set Dependencies

            this._dvnlsd = dvnlsd; this._dumach = dumach; this._dvnorm = dvnorm; this._dcopy = dcopy; this._dewset = dewset;
            this._dvhin = dvhin;this._dscal = dscal; this._dvindy = dvindy; this._xerrwd = xerrwd; this._dvstep = dvstep;
            #endregion
            #region Data Initialization

            //MORD(1)/12
            MORD[1 + o_mord] = 12;
            //MORD(2)/5
            MORD[2 + o_mord] = 5;
            //MXSTP0/500
            MXSTP0 = 500;
            //MXHNL0/10
            MXHNL0 = 10;
            //ZERO/0.0D0
            ZERO = 0.0E0;
            //ONE/1.0D0
            ONE = 1.0E0;
            //TWO/2.0D0
            TWO = 2.0E0;
            //FOUR/4.0D0
            FOUR = 4.0E0;
            //PT2/0.2D0
            PT2 = 0.2E0;
            //HUN/100.0D0
            HUN = 100.0E0;
            #endregion
            #region Common varaible Initialization

            #region Common Block: DVOD01 Initialization

            this._dvod01 = DVOD01;
            ACNRM = DVOD01.doubleData[0];
            CCMXJ = DVOD01.doubleData[1];
            CONP = DVOD01.doubleData[2];
            CRATE = DVOD01.doubleData[3];
            DRC = DVOD01.doubleData[4];
            //Start Array: EL  StartIndex: 5  LastIndex: 17
            EL = DVOD01.doubleData;
            o_el = 4;  //o_ = StartIndex-1
            //End Array: EL
            ETA = DVOD01.doubleData[18];
            ETAMAX = DVOD01.doubleData[19];
            H = DVOD01.doubleData[20];
            HMIN = DVOD01.doubleData[21];
            HMXI = DVOD01.doubleData[22];
            HNEW = DVOD01.doubleData[23];
            HSCAL = DVOD01.doubleData[24];
            PRL1 = DVOD01.doubleData[25];
            RC = DVOD01.doubleData[26];
            RL1 = DVOD01.doubleData[27];
            //Start Array: TAU  StartIndex: 28  LastIndex: 40
            TAU = DVOD01.doubleData;
            o_tau = 27;  //o_ = StartIndex-1
            //End Array: TAU
            //Start Array: TQ  StartIndex: 41  LastIndex: 45
            TQ = DVOD01.doubleData;
            o_tq = 40;  //o_ = StartIndex-1
            //End Array: TQ
            TN = DVOD01.doubleData[46];
            UROUND = DVOD01.doubleData[47];
            ICF = DVOD01.intData[0];
            INIT = DVOD01.intData[1];
            IPUP = DVOD01.intData[2];
            JCUR = DVOD01.intData[3];
            JSTART = DVOD01.intData[4];
            JSV = DVOD01.intData[5];
            KFLAG = DVOD01.intData[6];
            KUTH = DVOD01.intData[7];
            L = DVOD01.intData[8];
            LMAX = DVOD01.intData[9];
            LYH = DVOD01.intData[10];
            LEWT = DVOD01.intData[11];
            LACOR = DVOD01.intData[12];
            LSAVF = DVOD01.intData[13];
            LWM = DVOD01.intData[14];
            LIWM = DVOD01.intData[15];
            LOCJS = DVOD01.intData[16];
            MAXORD = DVOD01.intData[17];
            METH = DVOD01.intData[18];
            MITER = DVOD01.intData[19];
            MSBJ = DVOD01.intData[20];
            MXHNIL = DVOD01.intData[21];
            MXSTEP = DVOD01.intData[22];
            N = DVOD01.intData[23];
            NEWH = DVOD01.intData[24];
            NEWQ = DVOD01.intData[25];
            NHNIL = DVOD01.intData[26];
            NQ = DVOD01.intData[27];
            NQNYH = DVOD01.intData[28];
            NQWAIT = DVOD01.intData[29];
            NSLJ = DVOD01.intData[30];
            NSLP = DVOD01.intData[31];
            NYH = DVOD01.intData[32];
            #endregion
            #region Common Block: DVOD02 Initialization

            this._dvod02 = DVOD02;
            HU = DVOD02.doubleData[0];
            NCFN = DVOD02.intData[0];
            NETF = DVOD02.intData[1];
            NFE = DVOD02.intData[2];
            NJE = DVOD02.intData[3];
            NLU = DVOD02.intData[4];
            NNI = DVOD02.intData[5];
            NQU = DVOD02.intData[6];
            NST = DVOD02.intData[7];
            #endregion
            #endregion
        }
示例#19
0
文件: TN.cs 项目: davidsiaw/neuron
        public GTIMS()
        {
            #region Initialization Common Blocks

            CommonBlock SUBSCR = new CommonBlock(0, 15, 0, 0);
            #endregion
            #region Common varaible Initialization

            #region Common Block: SUBSCR Initialization

            this._subscr = SUBSCR;
            LGV = SUBSCR.intData[0];
            LZ1 = SUBSCR.intData[1];
            LZK = SUBSCR.intData[2];
            LV = SUBSCR.intData[3];
            LSK = SUBSCR.intData[4];
            LYK = SUBSCR.intData[5];
            LDIAGB = SUBSCR.intData[6];
            LSR = SUBSCR.intData[7];
            LYR = SUBSCR.intData[8];
            LHYR = SUBSCR.intData[9];
            LHG = SUBSCR.intData[10];
            LHYK = SUBSCR.intData[11];
            LPK = SUBSCR.intData[12];
            LEMAT = SUBSCR.intData[13];
            LWTEST = SUBSCR.intData[14];
            #endregion
            #endregion
        }
示例#20
0
文件: TN.cs 项目: davidsiaw/neuron
        public INITPC()
        {
            #region Initialization Common Blocks

            CommonBlock SUBSCR = new CommonBlock(0, 15, 0, 0);
            #endregion
            #region Dependencies (Initialization)

            DDOT ddot = new DDOT();
            DCOPY dcopy = new DCOPY();
            INITP3 initp3 = new INITP3(ddot, dcopy);
            #endregion
            #region Set Dependencies

            this._initp3 = initp3;
            #endregion
            #region Common varaible Initialization

            #region Common Block: SUBSCR Initialization

            this._subscr = SUBSCR;
            LGV = SUBSCR.intData[0];
            LZ1 = SUBSCR.intData[1];
            LZK = SUBSCR.intData[2];
            LV = SUBSCR.intData[3];
            LSK = SUBSCR.intData[4];
            LYK = SUBSCR.intData[5];
            LDIAGB = SUBSCR.intData[6];
            LSR = SUBSCR.intData[7];
            LYR = SUBSCR.intData[8];
            LHYR = SUBSCR.intData[9];
            LHG = SUBSCR.intData[10];
            LHYK = SUBSCR.intData[11];
            LPK = SUBSCR.intData[12];
            LEMAT = SUBSCR.intData[13];
            LWTEST = SUBSCR.intData[14];
            #endregion
            #endregion
        }
示例#21
0
        public DECOMR(DEC dec, DECB decb, ELMHES elmhes, DECH dech, CommonBlock LINAL)
        {
            #region Set Dependencies

            this._dec = dec; this._decb = decb; this._elmhes = elmhes; this._dech = dech;
            #endregion
            #region Common varaible Initialization

            #region Common Block: LINAL Initialization

            this._linal = LINAL;
            MLE = LINAL.intData[0];
            MUE = LINAL.intData[1];
            MBJAC = LINAL.intData[2];
            MBB = LINAL.intData[3];
            MDIAG = LINAL.intData[4];
            MDIFF = LINAL.intData[5];
            MBDIAG = LINAL.intData[6];
            #endregion
            #endregion
        }
示例#22
0
文件: TN.cs 项目: davidsiaw/neuron
        public LMQNBC(DDOT ddot, DNRM2 dnrm2, STEP1 step1, CRASH crash, SETPAR setpar, CHKUCP chkucp, SETUCR setucr, ZTIME ztime, MONIT monit, MODLNP modlnp
            , DCOPY dcopy, STPMAX stpmax, LINDER linder, MODZ modz, CNVTST cnvtst, DXPY dxpy, CommonBlock SUBSCR)
        {
            #region Set Dependencies

            this._ddot = ddot; this._dnrm2 = dnrm2; this._step1 = step1; this._crash = crash; this._setpar = setpar;
            this._chkucp = chkucp;this._setucr = setucr; this._ztime = ztime; this._monit = monit; this._modlnp = modlnp;
            this._dcopy = dcopy;this._stpmax = stpmax; this._linder = linder; this._modz = modz; this._cnvtst = cnvtst;
            this._dxpy = dxpy;
            #endregion
            #region Common varaible Initialization

            #region Common Block: SUBSCR Initialization

            this._subscr = SUBSCR;
            LGV = SUBSCR.intData[0];
            LZ1 = SUBSCR.intData[1];
            LZK = SUBSCR.intData[2];
            LV = SUBSCR.intData[3];
            LSK = SUBSCR.intData[4];
            LYK = SUBSCR.intData[5];
            LDIAGB = SUBSCR.intData[6];
            LSR = SUBSCR.intData[7];
            LYR = SUBSCR.intData[8];
            LOLDG = SUBSCR.intData[9];
            LHG = SUBSCR.intData[10];
            LHYK = SUBSCR.intData[11];
            LPK = SUBSCR.intData[12];
            LEMAT = SUBSCR.intData[13];
            LWTEST = SUBSCR.intData[14];
            #endregion
            #endregion
        }
示例#23
0
        public DECOMR()
        {
            #region Initialization Common Blocks

            CommonBlock LINAL = new CommonBlock(0, 7, 0, 0);
            #endregion
            #region Dependencies (Initialization)

            DEC dec = new DEC();
            DECB decb = new DECB();
            ELMHES elmhes = new ELMHES();
            DECH dech = new DECH();
            #endregion
            #region Set Dependencies

            this._dec = dec; this._decb = decb; this._elmhes = elmhes; this._dech = dech;
            #endregion
            #region Common varaible Initialization

            #region Common Block: LINAL Initialization

            this._linal = LINAL;
            MLE = LINAL.intData[0];
            MUE = LINAL.intData[1];
            MBJAC = LINAL.intData[2];
            MBB = LINAL.intData[3];
            MDIAG = LINAL.intData[4];
            MDIFF = LINAL.intData[5];
            MBDIAG = LINAL.intData[6];
            #endregion
            #endregion
        }
示例#24
0
文件: TN.cs 项目: davidsiaw/neuron
        public MSOLVE(MSLV mslv, CommonBlock SUBSCR)
        {
            #region Set Dependencies

            this._mslv = mslv;
            #endregion
            #region Common varaible Initialization

            #region Common Block: SUBSCR Initialization

            this._subscr = SUBSCR;
            LGV = SUBSCR.intData[0];
            LZ1 = SUBSCR.intData[1];
            LZK = SUBSCR.intData[2];
            LV = SUBSCR.intData[3];
            LSK = SUBSCR.intData[4];
            LYK = SUBSCR.intData[5];
            LDIAGB = SUBSCR.intData[6];
            LSR = SUBSCR.intData[7];
            LYR = SUBSCR.intData[8];
            LHYR = SUBSCR.intData[9];
            LHG = SUBSCR.intData[10];
            LHYK = SUBSCR.intData[11];
            LPK = SUBSCR.intData[12];
            LEMAT = SUBSCR.intData[13];
            LWTEST = SUBSCR.intData[14];
            #endregion
            #endregion
        }
示例#25
0
文件: dvode.cs 项目: davidsiaw/neuron
        public DVINDY()
        {
            #region Initialization Common Blocks

            CommonBlock DVOD01 = new CommonBlock(48, 33, 0, 0);
            CommonBlock DVOD02 = new CommonBlock(1, 8, 0, 0);
            #endregion
            #region Dependencies (Initialization)

            DSCAL dscal = new DSCAL();
            IUMACH iumach = new IUMACH();
            IXSAV ixsav = new IXSAV(iumach);
            XERRWD xerrwd = new XERRWD(ixsav);
            #endregion
            #region Set Dependencies

            this._dscal = dscal; this._xerrwd = xerrwd;
            #endregion
            #region Data Initialization

            //HUN/100.0D0
            HUN = 100.0E0;
            //ZERO/0.0D0
            ZERO = 0.0E0;
            #endregion
            #region Common varaible Initialization

            #region Common Block: DVOD01 Initialization

            this._dvod01 = DVOD01;
            ACNRM = DVOD01.doubleData[0];
            CCMXJ = DVOD01.doubleData[1];
            CONP = DVOD01.doubleData[2];
            CRATE = DVOD01.doubleData[3];
            DRC = DVOD01.doubleData[4];
            //Start Array: EL  StartIndex: 5  LastIndex: 17
            EL = DVOD01.doubleData;
            o_el = 4;  //o_ = StartIndex-1
            //End Array: EL
            ETA = DVOD01.doubleData[18];
            ETAMAX = DVOD01.doubleData[19];
            H = DVOD01.doubleData[20];
            HMIN = DVOD01.doubleData[21];
            HMXI = DVOD01.doubleData[22];
            HNEW = DVOD01.doubleData[23];
            HSCAL = DVOD01.doubleData[24];
            PRL1 = DVOD01.doubleData[25];
            RC = DVOD01.doubleData[26];
            RL1 = DVOD01.doubleData[27];
            //Start Array: TAU  StartIndex: 28  LastIndex: 40
            TAU = DVOD01.doubleData;
            o_tau = 27;  //o_ = StartIndex-1
            //End Array: TAU
            //Start Array: TQ  StartIndex: 41  LastIndex: 45
            TQ = DVOD01.doubleData;
            o_tq = 40;  //o_ = StartIndex-1
            //End Array: TQ
            TN = DVOD01.doubleData[46];
            UROUND = DVOD01.doubleData[47];
            ICF = DVOD01.intData[0];
            INIT = DVOD01.intData[1];
            IPUP = DVOD01.intData[2];
            JCUR = DVOD01.intData[3];
            JSTART = DVOD01.intData[4];
            JSV = DVOD01.intData[5];
            KFLAG = DVOD01.intData[6];
            KUTH = DVOD01.intData[7];
            L = DVOD01.intData[8];
            LMAX = DVOD01.intData[9];
            LYH = DVOD01.intData[10];
            LEWT = DVOD01.intData[11];
            LACOR = DVOD01.intData[12];
            LSAVF = DVOD01.intData[13];
            LWM = DVOD01.intData[14];
            LIWM = DVOD01.intData[15];
            LOCJS = DVOD01.intData[16];
            MAXORD = DVOD01.intData[17];
            METH = DVOD01.intData[18];
            MITER = DVOD01.intData[19];
            MSBJ = DVOD01.intData[20];
            MXHNIL = DVOD01.intData[21];
            MXSTEP = DVOD01.intData[22];
            N = DVOD01.intData[23];
            NEWH = DVOD01.intData[24];
            NEWQ = DVOD01.intData[25];
            NHNIL = DVOD01.intData[26];
            NQ = DVOD01.intData[27];
            NQNYH = DVOD01.intData[28];
            NQWAIT = DVOD01.intData[29];
            NSLJ = DVOD01.intData[30];
            NSLP = DVOD01.intData[31];
            NYH = DVOD01.intData[32];
            #endregion
            #region Common Block: DVOD02 Initialization

            this._dvod02 = DVOD02;
            HU = DVOD02.doubleData[0];
            NCFN = DVOD02.intData[0];
            NETF = DVOD02.intData[1];
            NFE = DVOD02.intData[2];
            NJE = DVOD02.intData[3];
            NLU = DVOD02.intData[4];
            NNI = DVOD02.intData[5];
            NQU = DVOD02.intData[6];
            NST = DVOD02.intData[7];
            #endregion
            #endregion
        }
示例#26
0
文件: TN.cs 项目: davidsiaw/neuron
        public SETPAR()
        {
            #region Initialization Common Blocks

            CommonBlock SUBSCR = new CommonBlock(0, 15, 0, 0);
            #endregion
            #region Common varaible Initialization

            #region Common Block: SUBSCR Initialization

            this._subscr = SUBSCR;
            //Start Array: LSUB  StartIndex: 0  LastIndex: 13
            LSUB = SUBSCR.intData;
            o_lsub = -1;  //o_ = StartIndex-1
            //End Array: LSUB
            LWTEST = SUBSCR.intData[14];
            #endregion
            #endregion
        }
示例#27
0
        public ESTRAD(SOL sol, SOLB solb, SOLH solh, CommonBlock LINAL)
        {
            #region Set Dependencies

            this._sol = sol; this._solb = solb; this._solh = solh;
            #endregion
            #region Common varaible Initialization

            #region Common Block: LINAL Initialization

            this._linal = LINAL;
            MLE = LINAL.intData[0];
            MUE = LINAL.intData[1];
            MBJAC = LINAL.intData[2];
            MBB = LINAL.intData[3];
            MDIAG = LINAL.intData[4];
            MDIFF = LINAL.intData[5];
            MBDIAG = LINAL.intData[6];
            #endregion
            #endregion
        }
示例#28
0
        public COBYLATEST()
        {
            #region Initialization Common Blocks

            CommonBlock Default = new CommonBlock(0, 1, 0, 0);
            #endregion
            #region Dependencies (Initialization)

            TRSTLP trstlp = new TRSTLP();
            CALCFC calcfc = new CALCFC(Default);
            COBYLB cobylb = new COBYLB(calcfc, trstlp);
            COBYLA cobyla = new COBYLA(cobylb);
            #endregion
            #region Set Dependencies

            this._cobyla = cobyla;
            #endregion
            #region Common varaible Initialization

            #region Common Block: Default Initialization

            this._default = Default;
            NPROB = Default.intData[0];
            #endregion
            #endregion
        }
示例#29
0
        public ESTRAD()
        {
            #region Initialization Common Blocks

            CommonBlock LINAL = new CommonBlock(0, 7, 0, 0);
            #endregion
            #region Dependencies (Initialization)

            SOL sol = new SOL();
            SOLB solb = new SOLB();
            SOLH solh = new SOLH();
            #endregion
            #region Set Dependencies

            this._sol = sol; this._solb = solb; this._solh = solh;
            #endregion
            #region Common varaible Initialization

            #region Common Block: LINAL Initialization

            this._linal = LINAL;
            MLE = LINAL.intData[0];
            MUE = LINAL.intData[1];
            MBJAC = LINAL.intData[2];
            MBB = LINAL.intData[3];
            MDIAG = LINAL.intData[4];
            MDIFF = LINAL.intData[5];
            MBDIAG = LINAL.intData[6];
            #endregion
            #endregion
        }
示例#30
0
        public DECOMC()
        {
            #region Initialization Common Blocks

            CommonBlock LINAL = new CommonBlock(0, 7, 0, 0);
            #endregion
            #region Dependencies (Initialization)

            DECC decc = new DECC();
            DECBC decbc = new DECBC();
            DECHC dechc = new DECHC();
            #endregion
            #region Set Dependencies

            this._decc = decc; this._decbc = decbc; this._dechc = dechc;
            #endregion
            #region Common varaible Initialization

            #region Common Block: LINAL Initialization

            this._linal = LINAL;
            MLE = LINAL.intData[0];
            MUE = LINAL.intData[1];
            MBJAC = LINAL.intData[2];
            MBB = LINAL.intData[3];
            MDIAG = LINAL.intData[4];
            MDIFF = LINAL.intData[5];
            MBDIAG = LINAL.intData[6];
            #endregion
            #endregion
        }
示例#31
0
文件: dvode.cs 项目: davidsiaw/neuron
        public DVNLSD()
        {
            #region Initialization Common Blocks

            CommonBlock DVOD01 = new CommonBlock(48, 33, 0, 0);
            CommonBlock DVOD02 = new CommonBlock(1, 8, 0, 0);
            #endregion
            #region Dependencies (Initialization)

            DVNORM dvnorm = new DVNORM();
            DCOPY dcopy = new DCOPY();
            DSCAL dscal = new DSCAL();
            IDAMAX idamax = new IDAMAX();
            DAXPY daxpy = new DAXPY();
            DDOT ddot = new DDOT();
            DGEFA dgefa = new DGEFA(idamax, dscal, daxpy);
            DACOPY dacopy = new DACOPY(dcopy);
            DGBFA dgbfa = new DGBFA(idamax, dscal, daxpy);
            DVJAC dvjac = new DVJAC(dvnorm, dcopy, dscal, dgefa, dacopy, dgbfa, DVOD01, DVOD02);
            DGESL dgesl = new DGESL(ddot, daxpy);
            DGBSL dgbsl = new DGBSL(ddot, daxpy);
            DVSOL dvsol = new DVSOL(dgesl, dgbsl, DVOD01);
            #endregion
            #region Set Dependencies

            this._dvnorm = dvnorm; this._dcopy = dcopy; this._dvjac = dvjac; this._dvsol = dvsol; this._dscal = dscal;
            this._daxpy = daxpy;
            #endregion
            #region Data Initialization

            //CCMAX/0.3D0
            CCMAX = 0.3E0;
            //CRDOWN/0.3D0
            CRDOWN = 0.3E0;
            //MAXCOR/3
            MAXCOR = 3;
            //MSBP/20
            MSBP = 20;
            //RDIV/2.0D0
            RDIV = 2.0E0;
            //ONE/1.0D0
            ONE = 1.0E0;
            //TWO/2.0D0
            TWO = 2.0E0;
            //ZERO/0.0D0
            ZERO = 0.0E0;
            #endregion
            #region Common varaible Initialization

            #region Common Block: DVOD01 Initialization

            this._dvod01 = DVOD01;
            ACNRM = DVOD01.doubleData[0];
            CCMXJ = DVOD01.doubleData[1];
            CONP = DVOD01.doubleData[2];
            CRATE = DVOD01.doubleData[3];
            DRC = DVOD01.doubleData[4];
            //Start Array: EL  StartIndex: 5  LastIndex: 17
            EL = DVOD01.doubleData;
            o_el = 4;  //o_ = StartIndex-1
            //End Array: EL
            ETA = DVOD01.doubleData[18];
            ETAMAX = DVOD01.doubleData[19];
            H = DVOD01.doubleData[20];
            HMIN = DVOD01.doubleData[21];
            HMXI = DVOD01.doubleData[22];
            HNEW = DVOD01.doubleData[23];
            HSCAL = DVOD01.doubleData[24];
            PRL1 = DVOD01.doubleData[25];
            RC = DVOD01.doubleData[26];
            RL1 = DVOD01.doubleData[27];
            //Start Array: TAU  StartIndex: 28  LastIndex: 40
            TAU = DVOD01.doubleData;
            o_tau = 27;  //o_ = StartIndex-1
            //End Array: TAU
            //Start Array: TQ  StartIndex: 41  LastIndex: 45
            TQ = DVOD01.doubleData;
            o_tq = 40;  //o_ = StartIndex-1
            //End Array: TQ
            TN = DVOD01.doubleData[46];
            UROUND = DVOD01.doubleData[47];
            ICF = DVOD01.intData[0];
            INIT = DVOD01.intData[1];
            IPUP = DVOD01.intData[2];
            JCUR = DVOD01.intData[3];
            JSTART = DVOD01.intData[4];
            JSV = DVOD01.intData[5];
            KFLAG = DVOD01.intData[6];
            KUTH = DVOD01.intData[7];
            L = DVOD01.intData[8];
            LMAX = DVOD01.intData[9];
            LYH = DVOD01.intData[10];
            LEWT = DVOD01.intData[11];
            LACOR = DVOD01.intData[12];
            LSAVF = DVOD01.intData[13];
            LWM = DVOD01.intData[14];
            LIWM = DVOD01.intData[15];
            LOCJS = DVOD01.intData[16];
            MAXORD = DVOD01.intData[17];
            METH = DVOD01.intData[18];
            MITER = DVOD01.intData[19];
            MSBJ = DVOD01.intData[20];
            MXHNIL = DVOD01.intData[21];
            MXSTEP = DVOD01.intData[22];
            N = DVOD01.intData[23];
            NEWH = DVOD01.intData[24];
            NEWQ = DVOD01.intData[25];
            NHNIL = DVOD01.intData[26];
            NQ = DVOD01.intData[27];
            NQNYH = DVOD01.intData[28];
            NQWAIT = DVOD01.intData[29];
            NSLJ = DVOD01.intData[30];
            NSLP = DVOD01.intData[31];
            NYH = DVOD01.intData[32];
            #endregion
            #region Common Block: DVOD02 Initialization

            this._dvod02 = DVOD02;
            HU = DVOD02.doubleData[0];
            NCFN = DVOD02.intData[0];
            NETF = DVOD02.intData[1];
            NFE = DVOD02.intData[2];
            NJE = DVOD02.intData[3];
            NLU = DVOD02.intData[4];
            NNI = DVOD02.intData[5];
            NQU = DVOD02.intData[6];
            NST = DVOD02.intData[7];
            #endregion
            #endregion
        }