/dports/math/cmlib/cmlib-3.0_8/src/quadpkd/ |
H A D | dqage.f | 152 3 RESABS,RESULT,RLIST,UFLOW 217 IF(KEYF.EQ.1) CALL DQK15(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 218 IF(KEYF.EQ.2) CALL DQK21(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 219 IF(KEYF.EQ.3) CALL DQK31(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 220 IF(KEYF.EQ.4) CALL DQK41(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 221 IF(KEYF.EQ.5) CALL DQK51(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 222 IF(KEYF.EQ.6) CALL DQK61(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 233 IF(IER.NE.0.OR.(ABSERR.LE.ERRBND.AND.ABSERR.NE.RESABS) 259 IF(KEYF.EQ.1) CALL DQK15(F,A1,B1,AREA1,ERROR1,RESABS,DEFAB1) 260 IF(KEYF.EQ.2) CALL DQK21(F,A1,B1,AREA1,ERROR1,RESABS,DEFAB1) [all …]
|
H A D | dqk15.f | 1 SUBROUTINE DQK15(F,A,B,RESULT,ABSERR,RESABS,RESASC) argument 57 1 D1MACH,EPMACH,F,FC,FSUM,FVAL1,FVAL2,FV1,FV2,HLGTH,RESABS,RESASC, 139 RESABS = ABS(RESK) 150 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 161 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 169 RESABS = RESABS*DHLGTH 174 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 175 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
H A D | dqk21.f | 1 SUBROUTINE DQK21(F,A,B,RESULT,ABSERR,RESABS,RESASC) argument 57 1 D1MACH,EPMACH,F,FC,FSUM,FVAL1,FVAL2,FV1,FV2,HLGTH,RESABS,RESASC, 147 RESABS = ABS(RESK) 158 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 169 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 177 RESABS = RESABS*DHLGTH 182 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 183 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
H A D | dqk15w.f | 1 SUBROUTINE DQK15W(F,W,P1,P2,P3,P4,KP,A,B,RESULT,ABSERR,RESABS, argument 70 2 P1,P2,P3,P4,RESABS,RESASC,RESG,RESK,RESKH,RESULT,UFLOW,W,WG,WGK, 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 179 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX((EPMACH* 180 1 0.5D+02)*RESABS,ABSERR)
|
H A D | dqk31.f | 1 SUBROUTINE DQK31(F,A,B,RESULT,ABSERR,RESABS,RESASC) argument 57 1 D1MACH,EPMACH,F,FC,FSUM,FVAL1,FVAL2,FV1,FV2,HLGTH,RESABS,RESASC, 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 191 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 192 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
H A D | dqk41.f | 1 SUBROUTINE DQK41(F,A,B,RESULT,ABSERR,RESABS,RESASC) argument 58 1 D1MACH,EPMACH,F,FC,FSUM,FVAL1,FVAL2,FV1,FV2,HLGTH,RESABS,RESASC, 172 RESABS = ABS(RESK) 183 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 194 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 202 RESABS = RESABS*DHLGTH 207 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 208 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
/dports/math/cmlib/cmlib-3.0_8/src/quadpks/ |
H A D | qage.f | 150 3 RESABS,RESULT,RLIST,UFLOW 215 IF(KEYF.EQ.1) CALL QK15(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 216 IF(KEYF.EQ.2) CALL QK21(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 217 IF(KEYF.EQ.3) CALL QK31(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 218 IF(KEYF.EQ.4) CALL QK41(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 219 IF(KEYF.EQ.5) CALL QK51(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 220 IF(KEYF.EQ.6) CALL QK61(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 232 IF(IER.NE.0.OR.(ABSERR.LE.ERRBND.AND.ABSERR.NE.RESABS) 258 IF(KEYF.EQ.1) CALL QK15(F,A1,B1,AREA1,ERROR1,RESABS,DEFAB1) 259 IF(KEYF.EQ.2) CALL QK21(F,A1,B1,AREA1,ERROR1,RESABS,DEFAB1) [all …]
|
H A D | qk15.f | 1 SUBROUTINE QK15(F,A,B,RESULT,ABSERR,RESABS,RESASC) argument 55 1 FV1,FV2,HLGTH,RESABS,RESASC,RESG,RESK,RESKH,RESULT,R1MACH,UFLOW, 123 RESABS = ABS(RESK) 134 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 145 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 153 RESABS = RESABS*DHLGTH 159 IF(RESABS.GT.UFLOW/(0.5E+02*EPMACH)) ABSERR = AMAX1 160 1 ((EPMACH*0.5E+02)*RESABS,ABSERR)
|
H A D | qk21.f | 1 SUBROUTINE QK21(F,A,B,RESULT,ABSERR,RESABS,RESASC) argument 55 1 FV1,FV2,HLGTH,RESABS,RESG,RESK,RESKH,RESULT,R1MACH,UFLOW,WG,WGK, 133 RESABS = ABS(RESK) 144 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 155 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 163 RESABS = RESABS*DHLGTH 169 IF(RESABS.GT.UFLOW/(0.5E+02*EPMACH)) ABSERR = AMAX1 170 1 ((EPMACH*0.5E+02)*RESABS,ABSERR)
|
H A D | qk31.f | 1 SUBROUTINE QK31(F,A,B,RESULT,ABSERR,RESABS,RESASC) argument 55 1 FV1,FV2,HLGTH,RESABS,RESASC,RESG,RESK,RESKH,RESULT,R1MACH,UFLOW, 135 RESABS = ABS(RESK) 146 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 157 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 165 RESABS = RESABS*DHLGTH 171 IF(RESABS.GT.UFLOW/(0.5E+02*EPMACH)) ABSERR = AMAX1 172 1 ((EPMACH*0.5E+02)*RESABS,ABSERR)
|
/dports/math/maxima/maxima-5.43.2/src/numerical/slatec/fortran/ |
H A D | dqage.f | 162 3 RESABS,RESULT,RLIST,UFLOW local 227 IF(KEYF.EQ.1) CALL DQK15(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 228 IF(KEYF.EQ.2) CALL DQK21(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 229 IF(KEYF.EQ.3) CALL DQK31(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 230 IF(KEYF.EQ.4) CALL DQK41(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 231 IF(KEYF.EQ.5) CALL DQK51(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 232 IF(KEYF.EQ.6) CALL DQK61(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 243 IF(IER.NE.0.OR.(ABSERR.LE.ERRBND.AND.ABSERR.NE.RESABS) 269 IF(KEYF.EQ.1) CALL DQK15(F,A1,B1,AREA1,ERROR1,RESABS,DEFAB1) 270 IF(KEYF.EQ.2) CALL DQK21(F,A1,B1,AREA1,ERROR1,RESABS,DEFAB1) [all …]
|
H A D | dqk15.f | 2 SUBROUTINE DQK15 (F, A, B, RESULT, ABSERR, RESABS, RESASC) argument 64 1 D1MACH,EPMACH,F,FC,FSUM,FVAL1,FVAL2,FV1,FV2,HLGTH,RESABS,RESASC, local 147 RESABS = ABS(RESK) 158 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 169 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 177 RESABS = RESABS*DHLGTH 182 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 183 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
H A D | dqk21.f | 2 SUBROUTINE DQK21 (F, A, B, RESULT, ABSERR, RESABS, RESASC) argument 64 1 D1MACH,EPMACH,F,FC,FSUM,FVAL1,FVAL2,FV1,FV2,HLGTH,RESABS,RESASC, local 155 RESABS = ABS(RESK) 166 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 177 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 185 RESABS = RESABS*DHLGTH 190 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 191 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
H A D | dqk15w.f | 3 + RESABS, RESASC) argument 77 2 P1,P2,P3,P4,RESABS,RESASC,RESG,RESK,RESKH,RESULT,UFLOW,W,WG,WGK, local 148 RESABS = ABS(RESK) 161 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 174 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 182 RESABS = RESABS*DHLGTH 187 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX((EPMACH* 188 1 0.5D+02)*RESABS,ABSERR)
|
H A D | dqk31.f | 2 SUBROUTINE DQK31 (F, A, B, RESULT, ABSERR, RESABS, RESASC) argument 64 1 D1MACH,EPMACH,F,FC,FSUM,FVAL1,FVAL2,FV1,FV2,HLGTH,RESABS,RESASC, local 164 RESABS = ABS(RESK) 175 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 186 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 194 RESABS = RESABS*DHLGTH 199 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 200 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
/dports/math/slatec/src/ |
H A D | qage.f | 161 3 RESABS,RESULT,RLIST,UFLOW 226 IF(KEYF.EQ.1) CALL QK15(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 227 IF(KEYF.EQ.2) CALL QK21(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 228 IF(KEYF.EQ.3) CALL QK31(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 229 IF(KEYF.EQ.4) CALL QK41(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 230 IF(KEYF.EQ.5) CALL QK51(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 231 IF(KEYF.EQ.6) CALL QK61(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 243 IF(IER.NE.0.OR.(ABSERR.LE.ERRBND.AND.ABSERR.NE.RESABS) 269 IF(KEYF.EQ.1) CALL QK15(F,A1,B1,AREA1,ERROR1,RESABS,DEFAB1) 270 IF(KEYF.EQ.2) CALL QK21(F,A1,B1,AREA1,ERROR1,RESABS,DEFAB1) [all …]
|
H A D | dqage.f | 162 3 RESABS,RESULT,RLIST,UFLOW 227 IF(KEYF.EQ.1) CALL DQK15(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 228 IF(KEYF.EQ.2) CALL DQK21(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 229 IF(KEYF.EQ.3) CALL DQK31(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 230 IF(KEYF.EQ.4) CALL DQK41(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 231 IF(KEYF.EQ.5) CALL DQK51(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 232 IF(KEYF.EQ.6) CALL DQK61(F,A,B,RESULT,ABSERR,DEFABS,RESABS) 243 IF(IER.NE.0.OR.(ABSERR.LE.ERRBND.AND.ABSERR.NE.RESABS) 269 IF(KEYF.EQ.1) CALL DQK15(F,A1,B1,AREA1,ERROR1,RESABS,DEFAB1) 270 IF(KEYF.EQ.2) CALL DQK21(F,A1,B1,AREA1,ERROR1,RESABS,DEFAB1) [all …]
|
H A D | dqk15.f | 2 SUBROUTINE DQK15 (F, A, B, RESULT, ABSERR, RESABS, RESASC) argument 64 1 D1MACH,EPMACH,F,FC,FSUM,FVAL1,FVAL2,FV1,FV2,HLGTH,RESABS,RESASC, 147 RESABS = ABS(RESK) 158 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 169 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 177 RESABS = RESABS*DHLGTH 182 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 183 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
H A D | qk15.f | 2 SUBROUTINE QK15 (F, A, B, RESULT, ABSERR, RESABS, RESASC) argument 64 1 FV1,FV2,HLGTH,RESABS,RESASC,RESG,RESK,RESKH,RESULT,R1MACH,UFLOW, 133 RESABS = ABS(RESK) 144 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 155 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 163 RESABS = RESABS*DHLGTH 169 IF(RESABS.GT.UFLOW/(0.5E+02*EPMACH)) ABSERR = MAX 170 1 ((EPMACH*0.5E+02)*RESABS,ABSERR)
|
H A D | dqk21.f | 2 SUBROUTINE DQK21 (F, A, B, RESULT, ABSERR, RESABS, RESASC) argument 64 1 D1MACH,EPMACH,F,FC,FSUM,FVAL1,FVAL2,FV1,FV2,HLGTH,RESABS,RESASC, 155 RESABS = ABS(RESK) 166 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 177 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 185 RESABS = RESABS*DHLGTH 190 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 191 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
H A D | qk21.f | 2 SUBROUTINE QK21 (F, A, B, RESULT, ABSERR, RESABS, RESASC) argument 64 1 FV1,FV2,HLGTH,RESABS,RESG,RESK,RESKH,RESULT,R1MACH,UFLOW,WG,WGK, 143 RESABS = ABS(RESK) 154 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 165 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 173 RESABS = RESABS*DHLGTH 179 IF(RESABS.GT.UFLOW/(0.5E+02*EPMACH)) ABSERR = MAX 180 1 ((EPMACH*0.5E+02)*RESABS,ABSERR)
|
H A D | qk31.f | 2 SUBROUTINE QK31 (F, A, B, RESULT, ABSERR, RESABS, RESASC) argument 64 1 FV1,FV2,HLGTH,RESABS,RESASC,RESG,RESK,RESKH,RESULT,R1MACH,UFLOW, 145 RESABS = ABS(RESK) 156 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 167 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 175 RESABS = RESABS*DHLGTH 181 IF(RESABS.GT.UFLOW/(0.5E+02*EPMACH)) ABSERR = MAX 182 1 ((EPMACH*0.5E+02)*RESABS,ABSERR)
|
H A D | dqk31.f | 2 SUBROUTINE DQK31 (F, A, B, RESULT, ABSERR, RESABS, RESASC) argument 64 1 D1MACH,EPMACH,F,FC,FSUM,FVAL1,FVAL2,FV1,FV2,HLGTH,RESABS,RESASC, 164 RESABS = ABS(RESK) 175 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 186 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 194 RESABS = RESABS*DHLGTH 199 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX 200 1 ((EPMACH*0.5D+02)*RESABS,ABSERR)
|
H A D | dqk15w.f | 3 + RESABS, RESASC) 77 2 P1,P2,P3,P4,RESABS,RESASC,RESG,RESK,RESKH,RESULT,UFLOW,W,WG,WGK, 148 RESABS = ABS(RESK) 161 RESABS = RESABS+WGK(JTW)*(ABS(FVAL1)+ABS(FVAL2)) 174 RESABS = RESABS+WGK(JTWM1)*(ABS(FVAL1)+ABS(FVAL2)) 182 RESABS = RESABS*DHLGTH 187 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = MAX((EPMACH* 188 1 0.5D+02)*RESABS,ABSERR)
|
/dports/math/octave/octave-6.4.0/liboctave/external/quadpack/ |
H A D | dqk21.f | 1 SUBROUTINE DQK21(F,A,B,RESULT,ABSERR,RESABS,RESASC,IERR) argument 53 * D1MACH,EPMACH,FC,FSUM,FVAL1,FVAL2,FV1,FV2,HLGTH,RESABS,RESASC, local 145 RESABS = DABS(RESK) 158 RESABS = RESABS+WGK(JTW)*(DABS(FVAL1)+DABS(FVAL2)) 171 RESABS = RESABS+WGK(JTWM1)*(DABS(FVAL1)+DABS(FVAL2)) 179 RESABS = RESABS*DHLGTH 184 IF(RESABS.GT.UFLOW/(0.5D+02*EPMACH)) ABSERR = DMAX1 185 * ((EPMACH*0.5D+02)*RESABS,ABSERR)
|