Loading article…
| 原作者 | ロベルト・ ピーセンス エリーズ・デドンカー=カペンガ クリストフ・W・ユーバーフーバー デヴィッド・カハーナー |
|---|---|
| 初回リリース | 1981年5月 |
| 最終リリース | 1987年5月[1]
|
| 書かれた | フォートラン77 |
| タイプ | 図書館 |
| ライセンス | パブリックドメイン |
| Webサイト | nines.cs.kuleuven.be/software/QUADPACK/ |
QUADPACKは、1次元関数の数値積分のためのFORTRAN 77 ライブラリです。 [2]これはSLATEC Common Mathematical Libraryに含まれていたため、パブリックドメインとなっています。[3]個々のサブプログラムはnetlibでも入手できます。[4]
GNU Scientific LibraryはQUADPACKルーチンをCで再実装しました。SciPyはQUADPACKの一部にPythonインターフェースを提供します。 [5] [6]
ParaMonte ライブラリの pm_quadpack モジュールは、最新のFortranで QUADPACK ライブラリの 100% 型種別汎用多精度実装を提供します。
ルーチン
QUADPACKの主な焦点は自動積分ルーチンにあり、ユーザーは問題と絶対または相対誤差 許容値を入力し、ルーチンは要求された誤差を超えない範囲で積分を実行しようとします。QUADPACKには、このような自動ルーチンが9つあり、さらに非自動ルーチンもいくつかあります。自動ルーチンのうち1つを除いてすべて適応求積法を使用します。[7]
各適応ルーチンには、より多くの情報を提供し、より詳細な制御を可能にする拡張パラメータ リストを持つ E という接尾辞が付いたバージョンもあります。すべてのルーチンの倍精度バージョンは、接頭辞 D でリリースされました。
汎用ルーチン
積分関数をさらに解析せずに使用するのに最適な2つの汎用ルーチンは、有限区間での積分を行うQAGSと無限区間での積分を行うQAGIです。[7]これら2つのルーチンは、 GNU Octave(quadコマンド)[5]とR(関数) [9]で使用されていますintegrate。
- 品質保証
- 各サブ区間内で21点ガウス・クロンロッド積分法に基づくグローバル適応積分法を使用し、ピーター・ウィンのイプシロンアルゴリズムによって加速される。 [7] [10]
- カギ
- は無限区間に対する唯一の汎用ルーチンであり、変換を使用して無限区間を半開区間 (0,1] にマッピングし、21 点ではなく 15 点のガウス・クロンロッド積分法を使用することを除いて QAGS と同じアプローチを使用します。[2]実数直線全体の積分の場合、使用される変換は次のとおりです。[2]これはすべての積分関数に最適なアプローチではありません。別の変換が適切である可能性があります。または、元の区間を分割して無限部分に対してのみ QAGI を使用することをお勧めします。[7]
その他の自動ルーチンの概要
- QNG
- 単純な非適応積分器
- グアグ
- 単純な適応積分器
- グアム
- QAGS に似ていますが、ユーザーが内部特異点、不連続点などの位置を指定できます。
- カウォ
- 有限区間におけるcos( ωx ) f ( x )またはsin( ωx ) f ( x )の積分
- クアラルンプール
- フーリエ変換
- QAWS
- w ( x ) f ( x )をaからbまで積分したもの。ここでfは滑らかで、w ( x ) = ( x – a ) α ( b – x ) β log k ( x – a ) log l ( b – x )であり、k、l = 0 または 1、α、β > –1である。
- クアラルンプール
- ユーザ指定のcとfに対するf ( x )/( x - c )の積分のコーシー主値[2]
参照
参考文献
- ^ "quadpack/changes". Netlib . 2010年11月16日閲覧。
- ^ abcd ピーセンス、ロバート;デ・ドンカー=カペンガ、エリーゼ。ユーバーフーバー、クリストフ W.デビッド・カハナー (1983)。QUADPACK: 自動統合のためのサブルーチン パッケージ。スプリンガー・フェルラーグ。ISBN 978-3-540-12553-2。
- ^ Fong, Kirby W.; Jefferson, Thomas H.; Suyehiro, Tokihiko; Walton, Lee (1993 年 7 月). 「SLATEC 共通数学ライブラリガイド」. netlib.org . 2010 年11 月 13 日閲覧。
- ^ "quadpack". Netlib . 2010年11月13日閲覧。
- ^ ab 「クアッドパック」。数値積分、非線形方程式およびソフトウェア (NINES) グループ、ルーヴェン カトリック大学。2010 年11 月 13 日に取得。
- ^ 「scipy.integrate.quad -- SciPy v0.14.0 リファレンス ガイド」 。2014年7 月 1 日閲覧。
- ^ abcd ピーセンス、ロバート;デ・ドンカー、エリーゼ。デヴィッド・カハナー (1984-04-17)。 「サブルーチンQPDOC」。クアッドパック。ネットライブラリ。2010 年11 月 16 日に取得。
- ^ ツヴィリンガー、ダニエル (1992).統合のハンドブック。 AKピーターズ。 p. 255.ISBN 978-0-86720-293-9。
- ^ R 開発コアチームおよび世界中の貢献者 (2010 年 10 月)。「integrate {stats}: 1 次元関数の積分」。パッケージ 'stats' バージョン 2.13.0 のドキュメント。2010年11 月 16 日に取得。
{{cite web}}:|author=一般的な名前があります (ヘルプ) - ^ 「17.4 特異点を持つ QAGS 適応積分」。GNU科学ライブラリ- リファレンス。フリーソフトウェア財団。2010年11 月 16 日閲覧。
さらに読む
- Favati, P.; Lotti, G.; Romani, F. (1991). 「アルゴリズム 691; QUADPACK 自動積分ルーチンの改善」. ACM Transactions on Mathematical Software . 17 (2): 218–232. doi : 10.1145/108556.108580 . S2CID 19675880.
- Cools, R.; Haegemans, A. (2003). 「アルゴリズム 824: CUBPACK : 自動立体化パッケージ; フレームワークの説明」. ACM Transactions on Mathematical Software . 29 (3): 287–296. doi :10.1145/838250.838253. S2CID 6855610.
