[논문 리뷰] A C++ interface to QCDNUM

Valerio Bertone, M. Botje|arXiv (Cornell University)|2017. 12. 21.

Particle physics theoretical and experimental studies참고 문헌 5인용 수 3

한 줄 요약

이 논문은 FORTRAN 기반의 QCDNUM 프로그램을 위한 C++ 인터페이스를 제시하며, 현대적인 C++ 기반 고에너지물리 소프트웨어에 QCD 파arton 밀도 진화를 원활하게 통합할 수 있도록 한다. 기존 FORTRAN 루틴을 C++ 호환 함수로 래핑함으로써, QCDNUM의 높은 성능를 유지하면서도 객체지향 프로그래밍, 동적 메모리 관리, 타입 안정성을 지원한다. 검증된 예제를 통해 정확한 NNLO 진화 및 구조 함수 계산이 가능함을 보여준다.

ABSTRACT

In this document we report on the recent development of a C++ interface to the FORTRAN-based evolution program QCDNUM. A short description of the interface is given with a few basic examples of its usage.

연구 동기 및 목표

• 기존 FORTRAN 기반의 QCDNUM와 현대적인 C++ 기반 HEP 소프트웨어 프레임워크 간 격차를 메우기 위해.
• QCDNUM의 고성능 계산 성능를 유지하면서도 객체지향 C++ 응용 프로그램에서 사용 가능하게 하기 위해.
• autotools와 네임스페이스 캡슐화를 사용하여 이식 가능하고 표준 준수 C++ 인터페이스를 제공하기 위해.
• NNLO PDF 진화, 구조 함수 계산, 동적 메모리 접근과 같은 QCDNUM의 핵심 기능을 지원하기 위해.
• C++ 네이티브 호출을 가능하게 하여 xFitter 및 Alpos와 같은 현대 도구와의 통합을 촉진하기 위해.

제안 방법

• 이름이 동일하지만 소문자로 변환되고 QCDNUM 네임스페이스에 래핑된 C++ 함수로 기존 FORTRAN 서브루틴을 래핑하여 이름 충돌을 방지하기 위해.
• 암시적 FORTRAN 타입 선언을 명시적 C++ 변수 선언으로 전환하며, 정수형 및 이중 정밀도 형식을 포함한 처리를 수행하기 위해.
• FORTRAN 배열 인덱싱(1기반, 애니티쿼크를 위해 음수 인덱스 사용)을 0기반 C++ 배열로 변환하기 위해 인덱스 오프셋 조정을 수행하기 위해.
• FORTRAN A(n,m) 배열을 효율적인 메모리 접근을 위해 1차원 C++ 배열로 변환하기 위한 2차원 배열 매핑 함수 k(i,j,n)를 구현하기 위해.
• 함수 포인터를 처리하기 위해 함수 인수를 C++ 함수 포인터로 전달하여 진화 루틴에서 사용자 정의 PDF 입력 함수를 가능하게 하기 위해.
• C++와 FORTRAN 간의 상호 운용성을 위해 이름 매핑 및 링커 지원을 포함한 autotools를 사용하여 라이브러리 빌드를 수행하기 위해.

실험 결과

연구 질문

• RQ1 고성능 FORTRAN 기반 QCD 진화 프로그램을 어떻게 현대적인 C++ HEP 응용 프로그램에 효과적으로 노출시킬 수 있는가?
• RQ2 타입 안정성, 메모리 관리, 이름 규칙 측면에서 레거시 FORTRAN 코드와 C++ 간의 인터페이스화 과정에서 발생하는 주요 과제는 무엇인가?
• RQ3 C++ 인터페이스는 QCDNUM의 계산 속도를 유지하면서도 객체지향 프로그래밍 기법을 지원할 수 있는가?
• RQ4 FORTRAN과 C++ 간의 배열 인덱싱 차이를 안정적이고 유지보수 용이한 방식으로 올바르게 처리할 수 있는가?
• RQ5 autotools와 네임스페이스 캡슐화 사용이 인터페이스의 이식성과 사용성에 미치는 영향은 무엇인가?

주요 결과

• C++ 인터페이스는 QCDNUM의 NNLO PDF 진화 및 구조 함수 계산 루틴을 C++ 코드에서 직접 호출할 수 있도록 하며, 완전한 타입 안정성과 모듈성을 보장한다.
• 원본 FORTRAN 코드베이스를 재사용함으로써 QCDNUM의 고성능를 유지하며 재구현으로 인한 성능 저하를 방지한다.
• autotools의 사용으로 다중 플랫폼 호환성과 적절한 링킹이 보장되며, qcdnum 릴리스 17-01/14 버전부터 인터페이스를 이용할 수 있다.
• 예제 코드는 변수-플레이버러 수치 체계에서 비극성 PDF의 정확한 진화를 보여주며, x = 1e-3 및 μ² = 1000 GeV²에서 쿼크 해의 밀도가 1.8708임을 도출한다.
• 인터페이스를 통해 강한 상호작용 상수 αs(mZ²)는 0.11807로 정확히 계산되며, 이는 기존 FORTRAN 구현과의 일관성을 확인한다.
• 인터페이스는 그리드 설정, 가중치 계산, LO, NLO, NNLO 순서의 진화, 그리고 임의의 x 및 μ²에서의 PDF 보간 기능을 포함한 핵심 기능을 지원한다.