Knowledge

QUADPACK

Source 📝

670: 449:
is the only general-purpose routine for infinite intervals, and maps the infinite interval onto the semi-open interval (0,1] using a transformation then uses the same approach as QAGS, except with 15-point rather than 21-point Gauss–Kronrod quadrature. For an integral over the whole real line, the
499: 214:
and the routine attempts to perform the integration with an error no larger than that requested. There are nine such automatic routines in QUADPACK, in addition to a number of non-automatic routines. All but one of the automatic routines use
672:
This is not the best approach for all integrands: another transformation may be appropriate, or one might prefer to break up the original interval and use QAGI only on the infinite part.
494: 1310: 396:
The two general-purpose routines most suitable for use without further analysis of the integrand are QAGS for integration over a finite interval and QAGI for
665:{\displaystyle \int _{-\infty }^{+\infty }f(x)dx=\int _{0}^{1}{dt \over t^{2}}\left(f\left({\frac {1-t}{t}}\right)+f\left(-{\frac {1-t}{t}}\right)\right)\;.} 384:
Each of the adaptive routines also have versions suffixed by E that have an extended parameter list that provides more information and allows more control.
1300: 1222: 1119: 958: 1181: 1029: 177: 428: 1305: 894: 332: 60: 1295: 436: 320: 1050: 1192: 409: 169: 184: 211: 1186: 372: 165: 1136: 146: 440: 191:
library offers a 100% type-kind-generic multi-precision implementation of QUADPACK library in modern
697:
similar to QAGS but allows user to specify locations of internal singularities, discontinuities etc.
432: 424: 275: 216: 208: 142: 91: 86: 1274: 1241: 453: 188: 1155: 1115: 1109: 954: 397: 361: 343: 945:
Piessens, Robert; de Doncker-Kapenga, Elise; Überhuber, Christoph W.; Kahaner, David (1983).
1266: 1231: 385: 98: 1167: 950: 1289: 224: 154: 103: 1245: 1278: 974:
Fong, Kirby W.; Jefferson, Thomas H.; Suyehiro, Tokihiko; Walton, Lee (July 1993).
114: 1025: 401: 139: 79: 25: 1028:. Numerical Integration, Nonlinear Equations & Software (NINES) Group, 1270: 1236: 1217: 913: 192: 1087: 1079: 1001: 975: 917: 158: 150: 1135:
R Development Core Team and contributors worldwide (October 2010).
1218:"Algorithm 691; Improving QUADPACK automatic integration routines" 1078:
Piessens, Robert; De Doncker, Elise; Kahaner, David (1984-04-17).
173: 997: 207:
integration routines in which the user inputs the problem and an
353:
Oscillatory weight function (cos or sin) over a finite interval
1137:"integrate {stats}: Integration of One-Dimensional Functions" 1257:: a package for automatic cubature; framework description". 947:
QUADPACK: A subroutine package for automatic integration
1051:"scipy.integrate.quad -- SciPy v0.14.0 Reference Guide" 388:
versions of all routines were released with prefix D.
502: 456: 331:
Specified points of local difficulty (singularities,
149:
of one-dimensional functions. It was included in the
153:
Common Mathematical Library and is therefore in the
1182:"17.4 QAGS adaptive integration with singularities" 306: 284: 260: 246: 157:. The individual subprograms are also available on 109: 97: 85: 75: 59: 40: 24: 664: 488: 1253:Cools, R.; Haegemans, A. (2003). "Algorithm 824: 976:"Guide to the SLATEC Common Mathematical Library" 1141:Documentation for package ‘stats’ version 2.13.0 677:Brief overview of the other automatic routines 8: 19: 1216:Favati, P.; Lotti, G.; Romani, F. (1991). 854:Cauchy principal value of the integral of 658: 18: 1259:ACM Transactions on Mathematical Software 1235: 1223:ACM Transactions on Mathematical Software 630: 594: 574: 560: 554: 549: 515: 507: 501: 478: 455: 1311:Public-domain software with source code 1020: 1018: 905: 168:reimplemented the QUADPACK routines in 1163: 1153: 940: 938: 936: 934: 1073: 1071: 1069: 1067: 398:integration over an infinite interval 7: 1301:Numerical integration (quadrature) 519: 511: 299:Weight function of specified form 14: 400:. These two routines are used in 203:The main focus of QUADPACK is on 180:interface to part of QUADPACK. 1030:Katholieke Universiteit Leuven 685:simple non-adaptive integrator 533: 527: 475: 463: 431:within each subinterval, with 1: 1114:. A K Peters. p. 255. 1108:Zwillinger, Daniel (1992). 895:List of numerical libraries 1327: 691:simple adaptive integrator 209:absolute or relative error 489:{\displaystyle x=(1-t)/t} 71: 55: 1193:Free Software Foundation 429:Gauss–Kronrod quadrature 392:General-purpose routines 1111:Handbook of integration 450:transformation used is 227:for automatic routines 32:Elise deDoncker-Kapenga 1187:GNU Scientific Library 735:over a finite interval 666: 490: 373:Cauchy principal value 166:GNU Scientific Library 34:Christoph W. Überhuber 1271:10.1145/838250.838253 1237:10.1145/108556.108580 667: 491: 147:numerical integration 500: 454: 873:for user-specified 559: 523: 425:adaptive quadrature 228: 217:adaptive quadrature 21: 1306:Numerical software 1166:has generic name ( 1080:"Subroutine QPDOC" 914:"quadpack/changes" 662: 545: 503: 486: 427:based on 21-point 312:Simple integrator 291:General integrand 222: 26:Original author(s) 1296:Fortran libraries 1121:978-0-86720-293-9 960:978-3-540-12553-2 741:Fourier transform 646: 610: 580: 441:epsilon algorithm 382: 381: 378: 377: 362:Fourier transform 344:Infinite interval 303: 302: 281: 280: 257: 256: 133: 132: 1318: 1282: 1249: 1239: 1204: 1203: 1201: 1199: 1178: 1172: 1171: 1165: 1161: 1159: 1151: 1149: 1147: 1132: 1126: 1125: 1105: 1099: 1098: 1096: 1094: 1075: 1062: 1061: 1059: 1057: 1047: 1041: 1040: 1038: 1036: 1022: 1013: 1012: 1010: 1008: 994: 988: 987: 985: 983: 971: 965: 964: 942: 929: 928: 926: 924: 910: 884: 878: 872: 848: 837: 826: 783: 777: 771: 765: 734: 718: 671: 669: 668: 663: 657: 653: 652: 648: 647: 642: 631: 615: 611: 606: 595: 581: 579: 578: 569: 561: 558: 553: 522: 514: 495: 493: 492: 487: 482: 415: 407: 386:Double precision 307: 285: 261: 247: 229: 221: 129: 126: 124: 122: 120: 118: 116: 51: 49: 22: 16:Software library 1326: 1325: 1321: 1320: 1319: 1317: 1316: 1315: 1286: 1285: 1252: 1215: 1212: 1210:Further reading 1207: 1197: 1195: 1180: 1179: 1175: 1162: 1152: 1145: 1143: 1134: 1133: 1129: 1122: 1107: 1106: 1102: 1092: 1090: 1077: 1076: 1065: 1055: 1053: 1049: 1048: 1044: 1034: 1032: 1024: 1023: 1016: 1006: 1004: 996: 995: 991: 981: 979: 973: 972: 968: 961: 951:Springer-Verlag 944: 943: 932: 922: 920: 912: 911: 907: 903: 891: 880: 874: 855: 839: 828: 785: 779: 773: 767: 748: 720: 704: 679: 632: 626: 622: 596: 590: 586: 582: 570: 562: 498: 497: 452: 451: 413: 405: 394: 333:discontinuities 201: 113: 67: 47: 45: 41:Initial release 35: 33: 31: 30:Robert Piessens 17: 12: 11: 5: 1324: 1322: 1314: 1313: 1308: 1303: 1298: 1288: 1287: 1284: 1283: 1265:(3): 287–296. 1250: 1230:(2): 218–232. 1211: 1208: 1206: 1205: 1173: 1127: 1120: 1100: 1063: 1042: 1014: 989: 966: 959: 930: 904: 902: 899: 898: 897: 890: 887: 886: 885: 852: 849: 784:is smooth and 745: 742: 739: 736: 701: 698: 695: 692: 689: 686: 683: 678: 675: 674: 673: 661: 656: 651: 645: 641: 638: 635: 629: 625: 621: 618: 614: 609: 605: 602: 599: 593: 589: 585: 577: 573: 568: 565: 557: 552: 548: 544: 541: 538: 535: 532: 529: 526: 521: 518: 513: 510: 506: 485: 481: 477: 474: 471: 468: 465: 462: 459: 447: 444: 421: 393: 390: 380: 379: 376: 375: 370: 366: 365: 359: 355: 354: 351: 347: 346: 341: 337: 336: 329: 325: 324: 318: 314: 313: 310: 304: 301: 300: 297: 293: 292: 289: 282: 279: 278: 273: 269: 268: 265: 258: 255: 254: 251: 243: 242: 239: 236: 233: 200: 197: 187:module of the 131: 130: 111: 107: 106: 101: 95: 94: 89: 83: 82: 77: 73: 72: 69: 68: 65: 63: 57: 56: 53: 52: 42: 38: 37: 28: 15: 13: 10: 9: 6: 4: 3: 2: 1323: 1312: 1309: 1307: 1304: 1302: 1299: 1297: 1294: 1293: 1291: 1280: 1276: 1272: 1268: 1264: 1260: 1256: 1251: 1247: 1243: 1238: 1233: 1229: 1225: 1224: 1219: 1214: 1213: 1209: 1194: 1190: 1188: 1183: 1177: 1174: 1169: 1164:|author= 1157: 1142: 1138: 1131: 1128: 1123: 1117: 1113: 1112: 1104: 1101: 1089: 1085: 1081: 1074: 1072: 1070: 1068: 1064: 1052: 1046: 1043: 1031: 1027: 1021: 1019: 1015: 1003: 999: 993: 990: 977: 970: 967: 962: 956: 952: 948: 941: 939: 937: 935: 931: 919: 915: 909: 906: 900: 896: 893: 892: 888: 883: 877: 870: 866: 862: 858: 853: 850: 846: 842: 835: 831: 824: 820: 816: 812: 808: 804: 800: 796: 792: 788: 782: 776: 770: 763: 759: 755: 751: 746: 743: 740: 737: 732: 728: 724: 716: 712: 708: 702: 699: 696: 693: 690: 687: 684: 681: 680: 676: 659: 654: 649: 643: 639: 636: 633: 627: 623: 619: 616: 612: 607: 603: 600: 597: 591: 587: 583: 575: 571: 566: 563: 555: 550: 546: 542: 539: 536: 530: 524: 516: 508: 504: 483: 479: 472: 469: 466: 460: 457: 448: 445: 442: 438: 434: 430: 426: 422: 419: 418: 417: 411: 408:command) and 403: 399: 391: 389: 387: 374: 371: 368: 367: 364:(cos or sin) 363: 360: 357: 356: 352: 349: 348: 345: 342: 339: 338: 334: 330: 327: 326: 322: 321:Singularities 319: 316: 315: 311: 309: 308: 305: 298: 295: 294: 290: 287: 286: 283: 277: 274: 271: 270: 267:Non-adaptive 266: 263: 262: 259: 252: 249: 248: 245: 244: 240: 237: 234: 231: 230: 226: 225:naming scheme 220: 218: 213: 210: 206: 198: 196: 194: 190: 186: 181: 179: 175: 171: 167: 162: 160: 156: 155:public domain 152: 148: 144: 141: 137: 128: 112: 108: 105: 104:Public domain 102: 100: 96: 93: 90: 88: 84: 81: 78: 74: 70: 64: 62: 61:Final release 58: 54: 44:May 1981 43: 39: 36:David Kahaner 29: 27: 23: 1262: 1258: 1254: 1227: 1221: 1196:. Retrieved 1189:-- Reference 1185: 1176: 1144:. Retrieved 1140: 1130: 1110: 1103: 1091:. Retrieved 1083: 1054:. Retrieved 1045: 1035:November 13, 1033:. Retrieved 1007:November 13, 1005:. Retrieved 992: 982:November 13, 980:. Retrieved 978:. netlib.org 969: 946: 923:November 16, 921:. Retrieved 908: 881: 875: 868: 864: 860: 856: 844: 840: 833: 829: 822: 818: 814: 810: 806: 802: 798: 794: 790: 786: 780: 774: 768: 761: 757: 753: 749: 747:integral of 730: 726: 722: 714: 710: 706: 703:integral of 433:acceleration 423:uses global 416:function). 395: 383: 204: 202: 182: 163: 135: 134: 66:May 1987 1198:16 November 1146:16 November 1093:16 November 253:Quadrature 241:4th letter 223:Summary of 185:pm_quadpack 176:provides a 1290:Categories 1026:"QUADPACK" 998:"quadpack" 901:References 437:Peter Wynn 402:GNU Octave 238:3rd letter 235:2nd letter 232:1st letter 140:FORTRAN 77 80:FORTRAN 77 76:Written in 637:− 628:− 601:− 547:∫ 520:∞ 512:∞ 509:− 505:∫ 470:− 414:integrate 212:tolerance 205:automatic 189:ParaMonte 125:/QUADPACK 123:/software 119:.kuleuven 1246:19675880 1156:cite web 1084:QUADPACK 889:See also 836:= 0 or 1 778:, where 323:handled 276:Adaptive 199:Routines 136:QUADPACK 20:QUADPACK 1279:6855610 1255:CUBPACK 847:> –1 827:, with 193:Fortran 143:library 110:Website 99:License 92:Library 48:1981-05 46: ( 1277:  1244:  1118:  1088:netlib 1056:1 July 1002:Netlib 957:  918:Netlib 817:) log( 809:) log( 178:Python 159:netlib 151:SLATEC 1275:S2CID 1242:S2CID 793:) = ( 766:from 412:(the 404:(the 174:SciPy 138:is a 115:nines 1200:2010 1168:help 1148:2010 1116:ISBN 1095:2010 1058:2014 1037:2010 1009:2010 984:2010 955:ISBN 925:2010 879:and 851:QAWC 838:and 744:QAWS 738:QAWF 721:sin( 705:cos( 700:QAWO 694:QAGP 446:QAGI 420:QAGS 406:quad 183:The 164:The 145:for 87:Type 1267:doi 1232:doi 863:)/( 801:) ( 772:to 719:or 688:QAG 682:QNG 439:'s 435:by 335:…) 121:.be 117:.cs 1292:: 1273:. 1263:29 1261:. 1240:. 1228:17 1226:. 1220:. 1191:. 1184:. 1160:: 1158:}} 1154:{{ 1139:. 1086:. 1082:. 1066:^ 1017:^ 1000:. 953:. 949:. 933:^ 916:. 843:, 832:, 756:) 725:) 723:ωx 709:) 707:ωx 496:: 219:. 195:. 172:. 161:. 1281:. 1269:: 1248:. 1234:: 1202:. 1170:) 1150:. 1124:. 1097:. 1060:. 1039:. 1011:. 986:. 963:. 927:. 882:f 876:c 871:) 869:c 867:– 865:x 861:x 859:( 857:f 845:β 841:α 834:l 830:k 825:) 823:x 821:– 819:b 815:a 813:– 811:x 807:x 805:– 803:b 799:a 797:– 795:x 791:x 789:( 787:w 781:f 775:b 769:a 764:) 762:x 760:( 758:f 754:x 752:( 750:w 733:) 731:x 729:( 727:f 717:) 715:x 713:( 711:f 660:. 655:) 650:) 644:t 640:t 634:1 624:( 620:f 617:+ 613:) 608:t 604:t 598:1 592:( 588:f 584:( 576:2 572:t 567:t 564:d 556:1 551:0 543:= 540:x 537:d 534:) 531:x 528:( 525:f 517:+ 484:t 480:/ 476:) 473:t 467:1 464:( 461:= 458:x 443:. 410:R 369:C 358:F 350:O 340:I 328:P 317:S 296:W 288:G 272:A 264:N 250:Q 170:C 127:/ 50:)

Index

Original author(s)
Final release
FORTRAN 77
Type
Library
License
Public domain
nines.cs.kuleuven.be/software/QUADPACK/
FORTRAN 77
library
numerical integration
SLATEC
public domain
netlib
GNU Scientific Library
C
SciPy
Python
pm_quadpack
ParaMonte
Fortran
absolute or relative error
tolerance
adaptive quadrature
naming scheme
Adaptive
Singularities
discontinuities
Infinite interval
Fourier transform

Text is available under the Creative Commons Attribution-ShareAlike License. Additional terms may apply.