[논문 리뷰] CHAPLIN - Complex Harmonic Polylogarithms in Fortran

Stephan Buehler, Claude Duhr|arXiv (Cornell University)|2011. 06. 28.

Advanced Mathematical Identities참고 문헌 24인용 수 14

한 줄 요약

이 논문은 임의의 복소수 변수에 대해 무게 4까지의 모든 조화 다이로그함수(HPLs)를 수치적으로 평가하기 위한 고성능 Fortran 라이브러리인 Chaplin을 제시한다. HPLs를 수렴이 빠른 최소 32개의 기저 함수로 감소시켜 원 안에서는 최적화된 급수 전개, 원 밖에서는 역행 관계를 적용함으로써 빠르고 정확한 결과를 도출하며, 이는 고에너지 물리학의 복잡한 운동학 조건에서의 고차원 루프 계산에 효과적으로 활용될 수 있다.

ABSTRACT

Abstract We present a new Fortran library to evaluate all harmonic polylogarithms up to weight four numerically for any complex argument. The algorithm is based on a reduction of harmonic polylogarithms up to weight four to a minimal set of basis functions that are computed numerically using series expansions allowing for fast and reliable numerical results. Title of program: Chaplin Catalogue Id: AETC_v1_0 Nature of problem Numerical evaluation of harmonic polylogarithms. Versions of this program held in the CPC repository in Mendeley Data AETC_v1_0; Chaplin; 10.1016/j.cpc.2014.05.022 This program has been imported from the CPC Program Library held at Queen's University Belfast (1969-2018)

연구 동기 및 목표

• 임의의 복소수 변수에 대해 무게 4까지의 조화 다이로그함수(HPLs)를 신속하고 신뢰성 있게 수치적으로 평가할 수 있는 라이브러리를 개발하는 것.
• 양자장이론에서 다음 다음 주요순서(Next-to-Next-to-Leading Order, NNLO) 계산에서 HPL 평가의 효율성 문제를 해결하는 것.
• 기존 도구보다 빠른 속도를 확보하면서도 복소평면 전역에서 수치 정확도를 유지하는 Fortran 기반 솔루션을 제공하는 것.
• 질량이 있는 입자를 포함한 일반적인 운동학 구성 조건과 복소 질량 체계에서 HPL의 정확한 계산을 가능하게 하는 것.

제안 방법

• 라이브러리는 무게 4까지의 모든 HPLs를 고전적 다이로그함수와 두 개의 새로운 무게 4 함수를 포함한 최소 32개의 기저 함수 집합으로 감소시킨다.
• 원 내부에서는 수렴이 매우 빠른 적절히 선택된 급수 전개를 사용하여 기저 함수를 계산한다.
• 원 외부에서는 변수를 원 내부로 매핑하기 위해 역행 관계를 적용하여 평가한다.
• HPL의 셔플 대수 구조를 활용하여 고중량 함수를 저중량 성분으로 체계적으로 감소시킨다.
• 수치적 안정성과 정밀도는 최적화된 급수 전개와 로그함수 및 다이로그함수 특이점의 철저한 처리를 통해 확보된다.
• 구현은 고에너지 물리학의 성능이 중요한 응용 분야에 적합하도록 Fortran 77을 선택하였다.

실험 결과

연구 질문

• RQ1 무게 4까지의 조화 다이로그함수는 임의의 복소수 변수에 대해 어떻게 효율적이고 정확하게 평가할 수 있는가?
• RQ2 수렴이 매우 빠른 수치 평가를 가능하게 하는 최소 기저 집합으로의 HPL 최적 감소는 무엇인가?
• RQ3 원 내외에서 수렴이 매우 빠르도록 구성된 급수 전개를 어떻게 설계할 수 있는가?
• RQ4 기존의 기호 계산 또는 C++ 기반 도구와 비교해 복소수 기반 구현의 성능는 어떻게 되는가?
• RQ5 다양한 복소 운동학 구성 조건에서 이 방법의 수치 정확도와 신뢰성은 어떠한가?

주요 결과

• Chaplin 라이브러리는 호출당 0.2에서 400 마이크로초의 시간이 소요되어, 무게 벡터와 변수에 따라 달라진다.
• 메서드는 매우 높은 수치 정밀도를 확보하여, z = 2.0 + 2.0i와 같은 테스트 포인트에서 HPL(Mathematica) 및 GiNaC와 비교해 최소 14자리 소수까지 일치한다.
• 기존의 두 개의 새로운 무게 4 함수를 포함한 32개의 기저 함수로의 감소는 체계적이고 효율적인 평가 전략을 가능하게 한다.
• 특히 대부분의 평가가 수행되는 원 내부에서 수렴이 매우 빠른 급수 전개가 구성되어 있다.
• 역행 관계의 사용은 원 외부의 변수에 대해서도 신뢰할 수 있는 평가를 가능하게 하여, 복소 질량 체계와 일반적인 루프 적분으로의 적용 범위를 확장한다.
• hplog(Fortran)에 비해 뛰어난 성능를 보이며, GiNaC와 같은 기호 계산 도구와도 경쟁 가능한 속도를 확보하고 있으며, 네이티브 컴파일링과 이식성도 확보하고 있다.