/dports/math/cmlib/cmlib-3.0_8/src/quadpks/ |
H A D | qk15w.f | 1 SUBROUTINE QK15W(F,W,P1,P2,P3,P4,KP,A,B,RESULT,ABSERR,RESABS, argument 68 2 HLGTH,P1,P2,P3,P4,RESABS,RESASC,RESG,RESK,RESKH,RESULT,UFLOW, 140 RESABS = ABS(RESK) 153 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 166 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 174 RESABS = RESABS*DHLGTH 180 IF(RESABS.GT.UFLOW/(0.5E+02*EPMACH)) ABSERR = AMAX1((EPMACH* 181 1 0.5E+02)*RESABS,ABSERR)
|
H A D | qk41.f | 1 SUBROUTINE QK41(F,A,B,RESULT,ABSERR,RESABS,RESASC) argument 56 1 FV1,FV2,HLGTH,RESABS, 146 RESABS = ABS(RESK) 157 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 168 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 176 RESABS = RESABS*DHLGTH 182 IF(RESABS.GT.UFLOW/(0.5E+02*EPMACH)) ABSERR = AMAX1 183 1 ((EPMACH*0.5E+02)*RESABS,ABSERR)
|
H A D | qk51.f | 1 SUBROUTINE QK51(F,A,B,RESULT,ABSERR,RESABS,RESASC) argument 55 1 FV1,FV2,HLGTH,RESABS,RESASC,RESG,RESK,RESKH,RESULT,R1MACH,UFLOW, 153 RESABS = ABS(RESK) 164 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 175 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 183 RESABS = RESABS*DHLGTH 189 IF(RESABS.GT.UFLOW/(0.5E+02*EPMACH)) ABSERR = AMAX1 190 1 ((EPMACH*0.5E+02)*RESABS,ABSERR)
|
H A D | qk61.f | 1 SUBROUTINE QK61(F,A,B,RESULT,ABSERR,RESABS,RESASC) argument 55 1 FV1,FV2,HLGTH,RESABS,RESASC,RESG,RESK,RESKH,RESULT,R1MACH,UFLOW, 163 RESABS = ABS(RESK) 174 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 185 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 193 RESABS = RESABS*DHLGTH 199 IF(RESABS.GT.UFLOW/(0.5E+02*EPMACH)) ABSERR = AMAX1 200 1 ((EPMACH*0.5E+02)*RESABS,ABSERR)
|
H A D | qk15i.f | 1 SUBROUTINE QK15I(F,BOUN,INF,A,B,RESULT,ABSERR,RESABS,RESASC) argument 75 2 FV2,HLGTH,RESABS,RESASC,RESG,RESK,RESKH,RESULT,TABSC1,TABSC2, 155 RESABS = ABS(RESK) 173 RESABS = RESABS+WGK(J)*(ABS(FVAL1)+ABS(FVAL2)) 182 RESABS = RESABS*HLGTH 186 IF(RESABS.GT.UFLOW/(0.5E+02*EPMACH)) ABSERR = AMAX1 187 1 ((EPMACH*0.5E+02)*RESABS,ABSERR)
|
H A D | qng.f | 81 2 RES87,RESABS,RESASC,RESKH,R1MACH,SAVFUN,UFLOW,W10,W21A,W43A, 251 RESABS = W21B(6)*ABS(FCENTR) 259 RESABS = RESABS+W21A(K)*(ABS(FVAL1)+ABS(FVAL2)) 272 RESABS = RESABS+W21B(K)*(ABS(FVAL1)+ABS(FVAL2)) 281 RESABS = RESABS*DHLGTH 329 IF (RESABS.GT.UFLOW/(0.5E+02*EPMACH)) ABSERR = AMAX1 330 1 ((EPMACH*0.5E+02)*RESABS,ABSERR)
|
/dports/math/maxima/maxima-5.43.2/src/numerical/slatec/fortran/ |
H A D | dqk41.f | 2 SUBROUTINE DQK41 (F, A, B, RESULT, ABSERR, RESABS, RESASC) argument 65 1 D1MACH,EPMACH,F,FC,FSUM,FVAL1,FVAL2,FV1,FV2,HLGTH,RESABS,RESASC, local 180 RESABS = ABS(RESK) 191 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 202 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 210 RESABS = RESABS*DHLGTH 215 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 216 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
H A D | dqk51.f | 2 SUBROUTINE DQK51 (F, A, B, RESULT, ABSERR, RESABS, RESASC) argument 65 1 D1MACH,EPMACH,F,FC,FSUM,FVAL1,FVAL2,FV1,FV2,HLGTH,RESABS,RESASC, local 193 RESABS = ABS(RESK) 204 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 215 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 223 RESABS = RESABS*DHLGTH 228 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 229 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
H A D | dqk61.f | 2 SUBROUTINE DQK61 (F, A, B, RESULT, ABSERR, RESABS, RESASC) argument 64 1 D1MACH,EPMACH,F,FC,FSUM,FVAL1,FVAL2,FV1,FV2,HLGTH,RESABS,RESASC, local 203 RESABS = ABS(RESK) 214 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 225 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 233 RESABS = RESABS*DHLGTH 238 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 239 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
H A D | dqk15i.f | 2 SUBROUTINE DQK15I (F, BOUN, INF, A, B, RESULT, ABSERR, RESABS, argument 85 2 RESABS,RESASC,RESG,RESK,RESKH,RESULT,TABSC1,TABSC2,UFLOW,WG,WGK, local 173 RESABS = ABS(RESK) 191 RESABS = RESABS+WGK(J)*(ABS(FVAL1)+ABS(FVAL2)) 200 RESABS = RESABS*HLGTH 204 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 205 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
H A D | dqng.f | 91 2 FV3,FV4,HLGTH,RESULT,RES10,RES21,RES43,RES87,RESABS,RESASC, local 301 RESABS = W21B(6)*ABS(FCENTR) 309 RESABS = RESABS+W21A(K)*(ABS(FVAL1)+ABS(FVAL2)) 322 RESABS = RESABS+W21B(K)*(ABS(FVAL1)+ABS(FVAL2)) 331 RESABS = RESABS*DHLGTH 378 IF (RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 379 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
/dports/math/slatec/src/ |
H A D | qk15w.f | 3 + RESABS, RESASC) 77 2 HLGTH,P1,P2,P3,P4,RESABS,RESASC,RESG,RESK,RESKH,RESULT,UFLOW, 150 RESABS = ABS(RESK) 163 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 176 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 184 RESABS = RESABS*DHLGTH 190 IF(RESABS.GT.UFLOW/(0.5E+02*EPMACH)) ABSERR = MAX((EPMACH* 191 1 0.5E+02)*RESABS,ABSERR)
|
H A D | qk41.f | 2 SUBROUTINE QK41 (F, A, B, RESULT, ABSERR, RESABS, RESASC) argument 65 1 FV1,FV2,HLGTH,RESABS, 156 RESABS = ABS(RESK) 167 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 178 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 186 RESABS = RESABS*DHLGTH 192 IF(RESABS.GT.UFLOW/(0.5E+02*EPMACH)) ABSERR = MAX 193 1 ((EPMACH*0.5E+02)*RESABS,ABSERR)
|
H A D | dqk41.f | 2 SUBROUTINE DQK41 (F, A, B, RESULT, ABSERR, RESABS, RESASC) argument 65 1 D1MACH,EPMACH,F,FC,FSUM,FVAL1,FVAL2,FV1,FV2,HLGTH,RESABS,RESASC, 180 RESABS = ABS(RESK) 191 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 202 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 210 RESABS = RESABS*DHLGTH 215 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 216 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
H A D | qk51.f | 2 SUBROUTINE QK51 (F, A, B, RESULT, ABSERR, RESABS, RESASC) argument 64 1 FV1,FV2,HLGTH,RESABS,RESASC,RESG,RESK,RESKH,RESULT,R1MACH,UFLOW, 163 RESABS = ABS(RESK) 174 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 185 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 193 RESABS = RESABS*DHLGTH 199 IF(RESABS.GT.UFLOW/(0.5E+02*EPMACH)) ABSERR = MAX 200 1 ((EPMACH*0.5E+02)*RESABS,ABSERR)
|
H A D | dqk61.f | 2 SUBROUTINE DQK61 (F, A, B, RESULT, ABSERR, RESABS, RESASC) argument 64 1 D1MACH,EPMACH,F,FC,FSUM,FVAL1,FVAL2,FV1,FV2,HLGTH,RESABS,RESASC, 203 RESABS = ABS(RESK) 214 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 225 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 233 RESABS = RESABS*DHLGTH 238 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 239 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
H A D | dqk51.f | 2 SUBROUTINE DQK51 (F, A, B, RESULT, ABSERR, RESABS, RESASC) argument 65 1 D1MACH,EPMACH,F,FC,FSUM,FVAL1,FVAL2,FV1,FV2,HLGTH,RESABS,RESASC, 193 RESABS = ABS(RESK) 204 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 215 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 223 RESABS = RESABS*DHLGTH 228 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 229 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
H A D | dqk15i.f | 2 SUBROUTINE DQK15I (F, BOUN, INF, A, B, RESULT, ABSERR, RESABS, argument 85 2 RESABS,RESASC,RESG,RESK,RESKH,RESULT,TABSC1,TABSC2,UFLOW,WG,WGK, 164 RESABS = ABS(RESK) 182 RESABS = RESABS+WGK(J)*(ABS(FVAL1)+ABS(FVAL2)) 191 RESABS = RESABS*HLGTH 195 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 196 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
H A D | qk61.f | 2 SUBROUTINE QK61 (F, A, B, RESULT, ABSERR, RESABS, RESASC) argument 64 1 FV1,FV2,HLGTH,RESABS,RESASC,RESG,RESK,RESKH,RESULT,R1MACH,UFLOW, 173 RESABS = ABS(RESK) 184 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 195 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 203 RESABS = RESABS*DHLGTH 209 IF(RESABS.GT.UFLOW/(0.5E+02*EPMACH)) ABSERR = MAX 210 1 ((EPMACH*0.5E+02)*RESABS,ABSERR)
|
H A D | qk15i.f | 2 SUBROUTINE QK15I (F, BOUN, INF, A, B, RESULT, ABSERR, RESABS, argument 85 2 FV2,HLGTH,RESABS,RESASC,RESG,RESK,RESKH,RESULT,TABSC1,TABSC2, 166 RESABS = ABS(RESK) 184 RESABS = RESABS+WGK(J)*(ABS(FVAL1)+ABS(FVAL2)) 193 RESABS = RESABS*HLGTH 197 IF(RESABS.GT.UFLOW/(0.5E+02*EPMACH)) ABSERR = MAX 198 1 ((EPMACH*0.5E+02)*RESABS,ABSERR)
|
/dports/math/cmlib/cmlib-3.0_8/src/quadpkd/ |
H A D | dqk51.f | 1 SUBROUTINE DQK51(F,A,B,RESULT,ABSERR,RESABS,RESASC) argument 57 1 D1MACH,EPMACH,F,FC,FSUM,FVAL1,FVAL2,FV1,FV2,HLGTH,RESABS,RESASC, 184 RESABS = ABS(RESK) 195 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 206 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 214 RESABS = RESABS*DHLGTH 219 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 220 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
H A D | dqk61.f | 1 SUBROUTINE DQK61(F,A,B,RESULT,ABSERR,RESABS,RESASC) argument 57 1 D1MACH,EPMACH,F,FC,FSUM,FVAL1,FVAL2,FV1,FV2,HLGTH,RESABS,RESASC, 195 RESABS = ABS(RESK) 206 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 217 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 225 RESABS = RESABS*DHLGTH 230 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 231 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
H A D | dqk15i.f | 1 SUBROUTINE DQK15I(F,BOUN,INF,A,B,RESULT,ABSERR,RESABS,RESASC) argument 77 2 RESABS,RESASC,RESG,RESK,RESKH,RESULT,TABSC1,TABSC2,UFLOW,WG,WGK, 155 RESABS = ABS(RESK) 173 RESABS = RESABS+WGK(J)*(ABS(FVAL1)+ABS(FVAL2)) 182 RESABS = RESABS*HLGTH 186 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 187 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
H A D | dqng.f | 83 2 FV3,FV4,HLGTH,RESULT,RES10,RES21,RES43,RES87,RESABS,RESASC, 292 RESABS = W21B(6)*ABS(FCENTR) 300 RESABS = RESABS+W21A(K)*(ABS(FVAL1)+ABS(FVAL2)) 313 RESABS = RESABS+W21B(K)*(ABS(FVAL1)+ABS(FVAL2)) 322 RESABS = RESABS*DHLGTH 369 IF (RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 370 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
/dports/math/octave/octave-6.4.0/liboctave/external/quadpack/ |
H A D | dqk15i.f | 1 SUBROUTINE DQK15I(F,BOUN,INF,A,B,RESULT,ABSERR,RESABS,RESASC, argument 74 * RESABS,RESASC,RESG,RESK,RESKH,RESULT,TABSC1,TABSC2,UFLOW,WG,WGK, local 167 RESABS = DABS(RESK) 195 RESABS = RESABS+WGK(J)*(DABS(FVAL1)+DABS(FVAL2)) 204 RESABS = RESABS*HLGTH 208 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = DMAX1 209 * ((EPMACH*0.5D+02)*RESABS,ABSERR)
|