[논문 리뷰] FastGPL: a C++ library for fast evaluation of generalized polylogarithms

Yuxuan Wang, Li Lin Yang외 1명|arXiv (Cornell University)|2021. 12. 08.

Particle physics theoretical and experimental studies인용 수 4

한 줄 요약

FastGPL은 일반화된 다이로그함수(GPLs)의 수치적 평가를 가속화하는 C++ 라이브러리로, 두 단계의 접근 방식을 사용한다: 외부 프로그램을 통해 최적화된 표현식을 사전 계산하고, 이를 하드코딩하여 런타임 평가를 빠르게 한다. 이는 기존 도구인 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와 같은 기존 도구의 한계를 극복한다. 이는 특정 GPLs에 대해 반복 계산이 중복되어 발생하는 느린 평가를 방지하기 위함이다.
• 급격한 수열 전개에서 발생하는 부동소수점 오버플로우를 방지하기 위해 철저한 산술 처리를 통해 정확도 손실을 방지한다.
• 몬테카를로 이벤트 생성에 적합한 효율적이고 정밀한 GPL 평가를 가능하게 하여, 다음다음주요순서(이중고차) 정확도에 대응한다.
• 기호 라이브러리인 GiNaC이나 Fortran 기반 도구인 handyG와 같은 기존 도구에 비해 생산용으로 사용 가능한 고성능 대안을 제공한다.

제안 방법

• Vollinga와 Weinzierl(2005)의 반복 알고리즘을 두 단계의 파ipel라인으로 구현한다: 첫 번째로 외부 프로그램을 사용해 최적화된 하드코딩된 표현식을 생성하고, 두 번째로 이를 C++ 라이브러리에 통합하여 빠른 평가를 구현한다.
• 인덱스 크기 기반 분류 체계(작은 값 vs. 큰 값)를 사용해 평가 경로를 사전에 계산하고 캐시하여 중복 계산을 줄인다.
• 수열 합산 중 오버플로우와 상쇄 오차를 최소화하기 위해 부동소수점 산술을 철저히 다룬다.
• GPL의 셔플 대수를 활용하여 GPL 곱의 평가를 단순화하고 최적화한다.
• 기존의 진폭 도구(FeynArts, Mathematica, FIRE6, LiteRed)와 통합하여 마스터 인tegr럴을 GPL 표현식으로 매핑한다.
• 다양한 GPL 구성, 특히 고중량 및 수치적으로 불안정한 케이스를 포함해 handyG와 GiNaC와의 성능을 비교 분석한다.

실험 결과

연구 질문

• RQ1 두 단계의 사전 계산 및 하드코딩 방식이, 단순 반복 구현에 비해 일반화된 다이로그함수(GPLs)의 평가 시간을 크게 줄일 수 있는가?
• RQ2 철저한 산술 설계를 통해 부동소수점 오버플로우 및 상쇄 오차로 인한 수치적 불안정성을 어느 정도 완화할 수 있는가?
• RQ3 다양한 GPL 구성에서 FastGPL의 성능는 handyG와 GiNaC에 비해 속도와 정확도 측면에서 어떻게 비교되는가?
• RQ4 FastGPL은 향후 e+e− 충돌기에서 정밀한 힉스 물리학의 이중고차 정밀도 계산을 효율적이고 정확하게 가능하게 할 수 있는가?
• RQ5 인덱스 분류(예: 작은 인덱스 vs. 큰 인덱스)는 GPL 평가의 최적화 및 성능에 어떤 영향을 미치는가?

주요 결과

• 특정 무게 4 GPLs, 예를 들어 G(1.0025,0.989,0.45,0.89+0.24i;1)와 같이 반복 단계가 많은 경우, FastGPL는 handyG에 비해 최대 400배 빠른 성능을 보이며 평가 시간이 약 2.2초에서 약 0.006초로 감소한다.
• 다양한 무게 4 GPL 카테고리에 걸쳐 평균적으로 FastGPL는 handyG보다 약 20배 빠르며, 인덱스 구성에 따라 약 14배에서 약 23배의 성능 향상이 관찰된다.
• 표현식 절단을 피하고 부동소수점 연산을 철저히 관리함으로써 고정밀도를 유지하여, 다른 구현에서 흔히 발생하는 정확도 손실을 방지한다.
• 실제 응용 사례에서 FastGPL는 e+e− 충돌기에서 350 GeV 중심질량 에너지에서 두 루프 혼합 QCD-EW 보정을 힉스 생성에 대해 계산하는 데 성공했으며, 한 개의 진폭당 평가 시간이 약 0.5초였다.
• 복잡하거나 집중된 인덱스를 가진 GPL에서 성능 향상이 가장 두드러지며, 단순 반복 계산이 반복적으로 발생하는 기존 구현 방식에서의 심각한 저하를 방지한다.
• 라이브러리는 무게 4까지의 GPL에 국한되지만, 이는 정밀한 힉스 물리학 및 전자기상 상호작용 물리학에서 NNLO 계산에 관련된 대부분의 진폭을 커버한다.