[論文レビュー] FastGPL: a C++ library for fast evaluation of generalized polylogarithms

Yuxuan Wang, Li Lin Yang+1 more|arXiv (Cornell University)|Dec 8, 2021

Particle physics theoretical and experimental studies被引用数 4

ひとこと要約

FastGPL は、外部プログラムを用いて事前計算された最適化された式をハードコードすることで、重み 4 までの一般化多リログラム（GPLs）の数値評価を高速化する C++ ライブラリである。2段階のアプローチを採用し、浮動小数点のオーバーフローとキャンセルエラーを回避することで、高精度を確保している。手頃なツール（例：handyG）と比較して、困難なケースでは最大 400 倍の高速化を達成している。

ABSTRACT

We present FastGPL, a C++ library for the fast evaluation of generalized polylogarithms which appear in many multi-loop Feynman integrals. We implement the iterative algorithm proposed by Vollinga and Weinzierl in a two-step approach, i.e., we generate concise expressions using an external program and hard-code them into the numeric library. This allows efficient and accurate numeric evaluations of generalized polylogarithms suitable for Monte Carlo integration and event generation. Floating-point arithmetics are carefully taken care of to avoid loss of accuracy. As an application and demonstration, we calculate the two-loop corrections for Higgs boson production in the vector boson fusion channel at electron-positron colliders. FastGPL is expected to be useful for event generators at the next-to-next-to-leading order accuracy.

研究の動機と目的

• 高精度な量子場の理論計算における一般化多リログラム（GPLs）の評価における性能ボトルネックを解消すること。
• handyG などの既存ツールに見られる、繰り返しの反復計算による特定 GPL の評価の遅さといった制限を克服すること。
• 系列展開における浮動小数点のオーバーフローを回避し、正確性の損失を防ぐために、注意深く算術処理を行うこと。
• 次々に次に至るオーダー（NNLO）でのモンテカルロイベント生成に適した、効率的かつ正確な GPL 評価を可能にすること。
• シンボリックライブラリ（例：GiNaC）や Fortran ベースのツール（例：handyG）の代替として、生産環境で使用可能な高性能な代替手段を提供すること。

提案手法

• Vollinga と Weinzierl（2005）の反復アルゴリズムを、2段階のパイプラインとして実装する：まず外部プログラムを用いて最適化されたハードコード式を生成し、次にそれを C++ ライブラリに埋め込んで高速評価を実現する。
• インデックスの絶対値の大小（小・大）に基づく分類スキームを用い、事前に評価パスを計算・キャッシュすることで、重複計算を削減する。
• 系列和算の過程でオーバーフローとキャンセルエラーを最小限に抑えるために、浮動小数点算術処理を注意深く設計する。
• GPL のシャッフル代数を活用し、GPL の積の評価を簡素化・最適化する。
• 既存の振幅ツール（例：FeynArts、Mathematica、FIRE6、LiteRed）と統合し、マスターリンテグラルを GPL 式にマッピングする。
• 多様な GPL 構成、特に高重みや数値的に不安定なケースを含む、handyG および GiNaC とのパフォーマンスベンチマークを実施する。

実験結果

リサーチクエスチョン

• RQ1 2段階の事前計算とハードコード化アプローチは、単純な反復実装と比較して、一般化多リログラムの評価時間を顕著に短縮できるか？
• RQ2 注意深く設計された算術処理により、浮動小数点のオーバーフローとキャンセルによる数値的不安定性は、どの程度軽減可能か？
• RQ3 さまざまな GPL 構成において、FastGPL の速度と正確性は、handyG や GiNaC と比較してどの程度優れているか？
• RQ4 FastGPL は、将来の e+e− コライダーにおける精密ヒッグス物理学における次々に次に至るオーダー（NNLO）計算を、効率的かつ正確に可能にするか？
• RQ5 インデックス分類（例：小インデックス vs. 大インデックス）は、GPL 評価の最適化とパフォーマンスにどのような影響を及ぼすか？

主な発見

• 特定の重み 4 GPL において、例として G(1.0025,0.989,0.45,0.89+0.24i;1) のような反復ステップが多い場合、FastGPL は handyG と比較して最大 400 倍の高速化を達成し、評価時間が約 2.2 秒から約 0.006 秒に短縮された。
• 平均して、FastGPL はさまざまな重み 4 GPL カテゴリにおいて、handyG より約 20 倍高速であり、インデックス構成に応じて約 14 倍から約 23 倍の高速化が見られた。
• 系列の切り捨てを回避し、浮動小数点演算を注意深く管理することで、高精度な数値的正確性を維持しており、他の実装でよく見られる正確性の損失を防いでいる。
• 実際の応用では、FastGPL を用いて、350 GeV の中心系エネルギーにおける e+e− コライダーでのベクトロン・ボソン融合におけるヒッグス生成の二ループ混合 QCD-EW 紛糾訂正を計算可能にした。1 振幅あたりの評価時間は約 0.5 秒であった。
• 複雑なまたは集中的なインデックスを持つ GPL において、パフォーマンスの向上が顕著であり、単純な反復実装では再計算が繰り返され、深刻な遅延が生じる。
• ライブラリは重み 4 までの GPL に限定されているが、これは精密ヒッグス物理学および電弱物理学における NNLO 計算に関連する大部分の振幅をカバーしている。