PrecisionLauricella: Package for numerical computation of Lauricella functions depending on a parameter

Bezuglov, M. A. and Kniehl, B. A. and Onishchenko, A. I. and Veretin, O. L. (2025) PrecisionLauricella: Package for numerical computation of Lauricella functions depending on a parameter. COMPUTER PHYSICS COMMUNICATIONS, 316: 109812. ISSN 0010-4655, 1879-2944

Full text not available from this repository. (Request a copy)

Abstract

We introduce the PrecisionLauricella package, a computational tool developed in Wolfram Mathematica for high-precision numerical evaluations of the Laurent expansion coefficients of Lauricella functions whose parameters depend linearly on a small regulator, s. In practical multi-loop calculations, Lauricella functions are required only as series around s = 0, and PrecisionLauricella is designed specifically to deliver such coefficients with arbitrary precision. The package leverages a method based on analytic continuation via Frobenius generalized power series, providing an efficient and accurate alternative to conventional approaches relying on multi-dimensional series expansions or Mellin-Barnes representations. This one-dimensional approach is particularly advantageous for high-precision calculations and facilitates further optimization through sdependent reconstruction from evaluations at specific numerical values, enabling efficient parallelization. The underlying mathematical framework for this method has been detailed in our previous work, while the current paper focuses on the design, implementation, and practical applications of the PrecisionLauricella package. Program summary Program Title: PrecisionLauricella CPC Library link to program files: https://doi.org/10.17632/6f958yz2dr.1 Developer's repository link: https://bitbucket.org/BezuglovMaxim/precisionlauricella-package/src/main/ Licensing provisions: GPLv3 Programming language: Wolfram Mathematica Supplementary material: PrecisionLauricella_Examples.nb Nature of problem: Lauricella functions, generalizations of hypergeometric functions, appearing in physics and mathematics, including Feynman integrals and string theory. When their indices depend linearly on a small parameter s, their numerical evaluation becomes challenging due to the complexity of high-dimensional series and singularities. Traditional methods, like hypergeometric re-expansion or Mellin-Barnes integrals, often lack efficiency and precision. Managing multi-dimensional sums exacerbates computational costs, especially for high-precision requirements, making these approaches unsuitable for many practical applications. Thus, there is a pressing need for efficient, scalable methods capable of maintaining numerical accuracy and effectively handling parameter dependencies. Solution method: Our method uses the Frobenius approach to achieve analytic continuations of Lauricella functions through generalized power series. Representing the functions as one-dimensional series simplifies high-precision numerical evaluations compared to traditional methods relying on multi-dimensional expansions or Mellin-Barnes integrals. We further optimize calculations by reconstructing s dependencies from evaluations at specific values, enabling efficient parallelization and reducing computational costs. A comprehensive mathematical exposition of the method is provided in our previous work [1]. References [1] M. Bezuglov, B. Kniehl, A. Onishchenko, O. Veretin, High-precision numerical evaluation of Lauricella functions, arXiv:2502.03276, 2 2025.

Item Type: Article
Uncontrolled Keywords: MATHEMATICA-BASED PACKAGES; EXPANDING HYPERGEOMETRIC-FUNCTIONS; ANALYTIC CONTINUATION; DIFFERENTIAL REDUCTION; TRANSCENDENTAL FUNCTIONS; EPSILON EXPANSION; FEYNMAN DIAGRAMS; F-D; HYPERDIRE; EQUATIONS;
Subjects: 500 Science > 530 Physics
Divisions: Physics > Institute of Theroretical Physics
Depositing User: Dr. Gernot Deinzer
Date Deposited: 06 May 2026 08:15
Last Modified: 06 May 2026 08:15
URI: https://pred.uni-regensburg.de/id/eprint/66929

Actions (login required)

View Item View Item