/dports/math/cmlib/cmlib-3.0_8/src/ddassl/ |
H A D | ddajac.f | 1 SUBROUTINE DDAJAC(NEQ,X,Y,YPRIME,DELTA,CJ,H, argument 78 DEL=DSIGN(DEL,H*YPRIME(I)) 81 YPSAVE=YPRIME(I) 83 YPRIME(I)=YPRIME(I)+CJ*DEL 91 YPRIME(I)=YPSAVE 127 WM(IPSAVE+K)=YPRIME(N) 130 DEL=DSIGN(DEL,H*YPRIME(N)) 133 510 YPRIME(N)=YPRIME(N)+CJ*DEL 134 CALL RES(X,Y,YPRIME,E,IRES,RPAR,IPAR) 139 YPRIME(N)=WM(IPSAVE+K) [all …]
|
H A D | ddaini.f | 1 SUBROUTINE DDAINI(X,Y,YPRIME,NEQ, argument 56 DIMENSION Y(*),YPRIME(*),WT(*) 84 100 PHI(I,2)=YPRIME(I) 100 250 Y(I)=Y(I)+H*YPRIME(I) 111 CALL RES(XNEW,Y,YPRIME,DELTA,IRES,RPAR,IPAR) 119 CALL DDAJAC(NEQ,XNEW,Y,YPRIME,DELTA,CJ,H, 147 330 YPRIME(I)=YPRIME(I)-CJ*DELTA(I) 193 420 YPRIME(I)=YPRIME(I)-CJ*DELTA(I) 232 610 YPRIME(I)=PHI(I,2)
|
/dports/math/cmlib/cmlib-3.0_8/src/sdassl/ |
H A D | sdajac.f | 1 SUBROUTINE SDAJAC(NEQ,X,Y,YPRIME,DELTA,CJ,H, argument 77 DEL=SIGN(DEL,H*YPRIME(I)) 80 YPSAVE=YPRIME(I) 82 YPRIME(I)=YPRIME(I)+CJ*DEL 90 YPRIME(I)=YPSAVE 126 WM(IPSAVE+K)=YPRIME(N) 129 DEL=SIGN(DEL,H*YPRIME(N)) 132 510 YPRIME(N)=YPRIME(N)+CJ*DEL 133 CALL RES(X,Y,YPRIME,E,IRES,RPAR,IPAR) 138 YPRIME(N)=WM(IPSAVE+K) [all …]
|
H A D | sdaini.f | 1 SUBROUTINE SDAINI(X,Y,YPRIME,NEQ, argument 55 DIMENSION Y(*),YPRIME(*),WT(*) 83 100 PHI(I,2)=YPRIME(I) 99 250 Y(I)=Y(I)+H*YPRIME(I) 110 CALL RES(XNEW,Y,YPRIME,DELTA,IRES,RPAR,IPAR) 118 CALL SDAJAC(NEQ,XNEW,Y,YPRIME,DELTA,CJ,H, 146 330 YPRIME(I)=YPRIME(I)-CJ*DELTA(I) 192 420 YPRIME(I)=YPRIME(I)-CJ*DELTA(I) 231 610 YPRIME(I)=PHI(I,2)
|
/dports/math/octave/octave-6.4.0/liboctave/external/daspk/ |
H A D | dmatd.f | 82 CALL JACD(X,Y,YPRIME,WM,CJ,RPAR,IPAR) 94 DEL=SIGN(DEL,H*YPRIME(I)) 97 YPSAVE=YPRIME(I) 99 YPRIME(I)=YPRIME(I)+CJ*DEL 108 YPRIME(I)=YPSAVE 128 CALL JACD(X,Y,YPRIME,WM,CJ,RPAR,IPAR) 148 WM(IPSAVE+K)=YPRIME(N) 151 DEL=SIGN(DEL,H*YPRIME(N)) 154 510 YPRIME(N)=YPRIME(N)+CJ*DEL 161 YPRIME(N)=WM(IPSAVE+K) [all …]
|
H A D | dyypnw.f | 5 SUBROUTINE DYYPNW (NEQ, Y, YPRIME, CJ, RL, P, ICOPT, ID, argument 38 DIMENSION Y(*), YPRIME(*), YNEW(*), YPNEW(*), ID(*), P(*) 44 YPNEW(I) = YPRIME(I) 47 YPNEW(I) = YPRIME(I) - RL*CJ*P(I) 53 YPNEW(I) = YPRIME(I)
|
H A D | dnedd.f | 5 SUBROUTINE DNEDD(X,Y,YPRIME,NEQ,RES,JACD,PDUM,H,WT, argument 133 DIMENSION Y(*),YPRIME(*),WT(*) 187 310 YPRIME(I)=0.0D0 191 320 YPRIME(I)=YPRIME(I)+GAMMA(J)*PHI(I,J) 199 CALL RES(X,Y,YPRIME,CJ,DELTA,IRES,RPAR,IPAR) 209 CALL DMATD(NEQ,X,Y,YPRIME,DELTA,CJ,H,IERJ,WT,E,WM,IWM, 220 CALL DNSD(X,Y,YPRIME,NEQ,RES,PDUM,WT,RPAR,IPAR,DUMSVR,
|
H A D | dnedk.f | 5 SUBROUTINE DNEDK(X,Y,YPRIME,NEQ,RES,JACK,PSOL, argument 131 DIMENSION Y(*),YPRIME(*),WT(*) 191 310 YPRIME(I)=0.0D0 195 320 YPRIME(I)=YPRIME(I)+GAMMA(J)*PHI(I,J) 203 CALL RES(X,Y,YPRIME,CJ,DELTA,IRES,RPAR,IPAR) 213 CALL JACK (RES, IRES, NEQ, X, Y, YPRIME, WT, DELTA, E, H, CJ, 223 CALL DNSK(X,Y,YPRIME,NEQ,RES,PSOL,WT,RPAR,IPAR,SAVR,
|
H A D | dnsk.f | 5 SUBROUTINE DNSK(X,Y,YPRIME,NEQ,RES,PSOL,WT,RPAR,IPAR, argument 94 DIMENSION Y(*),YPRIME(*),WT(*),DELTA(*),E(*),SAVR(*) 125 CALL DSLVK (NEQ, Y, X, YPRIME, SAVR, DELTA, WT, WM, IWM, 135 360 YPRIME(I) = YPRIME(I) - CJ*DELTA(I) 159 CALL RES(X,Y,YPRIME,CJ,DELTA,IRES,RPAR,IPAR)
|
H A D | ddasid.f | 5 SUBROUTINE DDASID(X,Y,YPRIME,NEQ,ICOPT,ID,RES,JACD,PDUM,H,WT, argument 93 DIMENSION Y(*),YPRIME(*),ID(*),WT(*),ICNSTR(*) 112 CALL RES(X,Y,YPRIME,CJ,DELTA,IRES,RPAR,IPAR) 130 CALL DMATD(NEQ,X,Y,YPRIME,DELTA,CJ,H,IERJ,WT,R, 136 CALL DNSID(X,Y,YPRIME,NEQ,ICOPT,ID,RES,WT,RPAR,IPAR,DELTA,R, 147 CALL RES(X,Y,YPRIME,CJ,DELTA,IRES,RPAR,IPAR)
|
H A D | datv.f | 5 SUBROUTINE DATV (NEQ, Y, TN, YPRIME, SAVR, V, WGHT, YPTEM, RES, argument 79 DOUBLE PRECISION Y, TN, YPRIME, SAVR, V, WGHT, YPTEM, Z, VTEM, local 81 DIMENSION Y(*), YPRIME(*), SAVR(*), V(*), WGHT(*), YPTEM(*), 99 YPTEM(I) = YPRIME(I) + VTEM(I)*CJ 118 CALL PSOL (NEQ, TN, Y, YPRIME, SAVR, YPTEM, CJ, WGHT, WP, IWP,
|
/dports/math/slatec/src/ |
H A D | ddajac.f | 98 DEL=SIGN(DEL,H*YPRIME(I)) 101 YPSAVE=YPRIME(I) 103 YPRIME(I)=YPRIME(I)+CJ*DEL 104 CALL RES(X,Y,YPRIME,E,IRES,RPAR,IPAR) 111 YPRIME(I)=YPSAVE 147 WM(IPSAVE+K)=YPRIME(N) 149 DEL=SIGN(DEL,H*YPRIME(N)) 152 510 YPRIME(N)=YPRIME(N)+CJ*DEL 153 CALL RES(X,Y,YPRIME,E,IRES,RPAR,IPAR) 158 YPRIME(N)=WM(IPSAVE+K) [all …]
|
H A D | sdajac.f | 97 DEL=SIGN(DEL,H*YPRIME(I)) 100 YPSAVE=YPRIME(I) 102 YPRIME(I)=YPRIME(I)+CJ*DEL 103 CALL RES(X,Y,YPRIME,E,IRES,RPAR,IPAR) 110 YPRIME(I)=YPSAVE 146 WM(IPSAVE+K)=YPRIME(N) 148 DEL=SIGN(DEL,H*YPRIME(N)) 151 510 YPRIME(N)=YPRIME(N)+CJ*DEL 152 CALL RES(X,Y,YPRIME,E,IRES,RPAR,IPAR) 157 YPRIME(N)=WM(IPSAVE+K) [all …]
|
H A D | ddaini.f | 2 SUBROUTINE DDAINI (X, Y, YPRIME, NEQ, RES, JAC, H, WT, IDID, RPAR, argument 58 * X, Y(*), YPRIME(*), H, WT(*), RPAR(*), PHI(NEQ,*), DELTA(*), 94 100 PHI(I,2)=YPRIME(I) 108 250 Y(I)=Y(I)+H*YPRIME(I) 119 CALL RES(X,Y,YPRIME,DELTA,IRES,RPAR,IPAR) 127 CALL DDAJAC(NEQ,X,Y,YPRIME,DELTA,CJ,H, 151 330 YPRIME(I)=YPRIME(I)-CJ*DELTA(I) 194 420 YPRIME(I)=YPRIME(I)-CJ*DELTA(I) 230 610 YPRIME(I)=PHI(I,2)
|
H A D | sdaini.f | 2 SUBROUTINE SDAINI (X, Y, YPRIME, NEQ, RES, JAC, H, WT, IDID, RPAR, argument 57 REAL X, Y(*), YPRIME(*), H, WT(*), RPAR(*), PHI(NEQ,*), DELTA(*), 92 100 PHI(I,2)=YPRIME(I) 106 250 Y(I)=Y(I)+H*YPRIME(I) 117 CALL RES(X,Y,YPRIME,DELTA,IRES,RPAR,IPAR) 125 CALL SDAJAC(NEQ,X,Y,YPRIME,DELTA,CJ,H, 149 330 YPRIME(I)=YPRIME(I)-CJ*DELTA(I) 192 420 YPRIME(I)=YPRIME(I)-CJ*DELTA(I) 228 610 YPRIME(I)=PHI(I,2)
|
H A D | dhstrt.f | 2 SUBROUTINE DHSTRT (DF, NEQ, A, B, Y, YPRIME, ETOL, MORDER, SMALL, argument 149 3 SRYDPB, TOLEXP, TOLMIN, TOLP, TOLSUM, Y, YDPB, YP, YPRIME 150 DIMENSION Y(*),YPRIME(*),ETOL(*),SPY(*),PV(*),YP(*), 175 YP(J) = SF(J) - YPRIME(J) 210 DELF = DHVNRM(YPRIME,NEQ) 215 SPY(J) = YPRIME(J) 216 YP(J) = YPRIME(J) 240 PV(J) = YP(J) - YPRIME(J)
|
H A D | ddassl.f | 1149 CALL DDAINI(TN,Y,YPRIME,NEQ, 1164 370 RWORK(ITEMP + I - 1) = H*YPRIME(I) 1197 CALL DDATRP(TN,TN,Y,YPRIME,NEQ,IWORK(LKOLD), 1226 CALL DDATRP(TN,TN,Y,YPRIME,NEQ,IWORK(LKOLD), 1311 CALL DDASTP(TN,Y,YPRIME,NEQ, 1331 CALL DDATRP(TN,TOUT,Y,YPRIME,NEQ, 1340 535 CALL DDATRP(TN,TOUT,Y,YPRIME,NEQ, 1347 CALL DDATRP(TN,TOUT,Y,YPRIME,NEQ, 1358 545 CALL DDATRP(TN,TSTOP,Y,YPRIME,NEQ, 1368 552 CALL DDATRP(TN,TSTOP,Y,YPRIME,NEQ, [all …]
|
/dports/math/octave/octave-6.4.0/liboctave/external/dassl/ |
H A D | ddajac.f | 1 SUBROUTINE DDAJAC (NEQ, X, Y, YPRIME, DELTA, CJ, H, argument 99 DEL=SIGN(DEL,H*YPRIME(I)) 102 YPSAVE=YPRIME(I) 104 YPRIME(I)=YPRIME(I)+CJ*DEL 112 YPRIME(I)=YPSAVE 148 WM(IPSAVE+K)=YPRIME(N) 150 DEL=SIGN(DEL,H*YPRIME(N)) 153 510 YPRIME(N)=YPRIME(N)+CJ*DEL 154 CALL RES(X,Y,YPRIME,E,IRES,RPAR,IPAR) 159 YPRIME(N)=WM(IPSAVE+K) [all …]
|
H A D | ddaini.f | 1 SUBROUTINE DDAINI (X, Y, YPRIME, NEQ, RES, JAC, H, WT, IDID, RPAR, argument 57 * X, Y(*), YPRIME(*), H, WT(*), RPAR(*), PHI(NEQ,*), DELTA(*), local 93 100 PHI(I,2)=YPRIME(I) 107 250 Y(I)=Y(I)+H*YPRIME(I) 118 CALL RES(X,Y,YPRIME,DELTA,IRES,RPAR,IPAR) 126 CALL DDAJAC(NEQ,X,Y,YPRIME,DELTA,CJ,H, 150 330 YPRIME(I)=YPRIME(I)-CJ*DELTA(I) 193 420 YPRIME(I)=YPRIME(I)-CJ*DELTA(I) 229 610 YPRIME(I)=PHI(I,2)
|
/dports/math/math77/MATH77/demo/ |
H A D | drsjacg2.f | 211 YPRIME(1)=ZERO 212 YPRIME(2)=ZERO 213 YPRIME(3)=ZERO 214 YPRIME(4)=-GRAVITY 215 YPRIME(5)=ZERO 350 YPRIME(1)=ZERO 351 CALL SCOPY(5,YPRIME,0,YPRIME,1) 480 YPRIME(1)=ZERO 481 YPRIME(2)=ZERO 482 YPRIME(3)=ZERO [all …]
|
H A D | drdjacg2.f | 211 YPRIME(1)=ZERO 212 YPRIME(2)=ZERO 213 YPRIME(3)=ZERO 214 YPRIME(4)=-GRAVITY 215 YPRIME(5)=ZERO 350 YPRIME(1)=ZERO 351 CALL DCOPY(5,YPRIME,0,YPRIME,1) 480 YPRIME(1)=ZERO 481 YPRIME(2)=ZERO 482 YPRIME(3)=ZERO [all …]
|
/dports/math/cmlib/cmlib-3.0_8/test/ |
H A D | sdassl.f | 119 YPRIME(1) = -10.0E0 120 YPRIME(2) = 10.0E0 188 CALL RES2(T,Y,DELTA,YPRIME,IRES,RPAR,IPAR) 227 DIMENSION Y(*),YPRIME(*),DELTA(*) 228 DELTA(1) = YPRIME(1) + 10.0E0*Y(1) 232 SUBROUTINE JAC1(T,Y,YPRIME,PD,CJ,RPAR,IPAR) argument 233 DIMENSION Y(*),YPRIME(*),PD(2,*) 241 DIMENSION Y(*),YPRIME(*),DELTA(*) 249 10 DELTA(K) = D - YPRIME(K) 252 SUBROUTINE JAC2(T,Y,YPRIME,PD,CJ,RPAR,IPAR) argument [all …]
|
H A D | ddassl.f | 120 YPRIME(1) = -10.0D0 121 YPRIME(2) = 10.0D0 189 CALL RES2(T,Y,DELTA,YPRIME,IRES,RPAR,IPAR) 229 DIMENSION Y(*),YPRIME(*),DELTA(*) 230 DELTA(1) = YPRIME(1) + 10.0D0*Y(1) 234 SUBROUTINE JAC1(T,Y,YPRIME,PD,CJ,RPAR,IPAR) argument 236 DIMENSION Y(*),YPRIME(*),PD(2,*) 245 DIMENSION Y(*),YPRIME(*),DELTA(*) 253 10 DELTA(K) = D - YPRIME(K) 256 SUBROUTINE JAC2(T,Y,YPRIME,PD,CJ,RPAR,IPAR) argument [all …]
|
/dports/math/scilab/scilab-6.1.1/scilab/modules/differential_equations/src/fortran/ |
H A D | ddassl.f | 94 100 PHI(I,2)=YPRIME(I) 159 330 YPRIME(I)=YPRIME(I)-CJ*DELTA(I) 202 420 YPRIME(I)=YPRIME(I)-CJ*DELTA(I) 238 610 YPRIME(I)=PHI(I,2) 369 YPSAVE=YPRIME(I) 371 YPRIME(I)=YPRIME(I)+CJ*DEL 380 YPRIME(I)=YPSAVE 422 510 YPRIME(N)=YPRIME(N)+CJ*DEL 2407 310 YPRIME(I)=0.0D0 2411 320 YPRIME(I)=YPRIME(I)+GAMMA(J)*PHI(I,J) [all …]
|
H A D | ddaskr.f | 1407 DIMENSION Y(*),YPRIME(*) 4799 310 YPRIME(I)=0.0D0 4803 320 YPRIME(I)=YPRIME(I)+GAMMA(J)*PHI(I,J) 5005 340 YPRIME(I)=YPRIME(I)-CJ*DELTA(I) 5145 YPSAVE=YPRIME(I) 5147 YPRIME(I)=YPRIME(I)+CJ*DEL 5156 YPRIME(I)=YPSAVE 5201 510 YPRIME(N)=YPRIME(N)+CJ*DEL 6090 310 YPRIME(I)=0.0D0 6094 320 YPRIME(I)=YPRIME(I)+GAMMA(J)*PHI(I,J) [all …]
|