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

IF 3.4 2区 物理与天体物理 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS
M.A. Bezuglov , B.A. Kniehl , A.I. Onishchenko , O.L. Veretin
{"title":"PrecisionLauricella: Package for numerical computation of Lauricella functions depending on a parameter","authors":"M.A. Bezuglov ,&nbsp;B.A. Kniehl ,&nbsp;A.I. Onishchenko ,&nbsp;O.L. Veretin","doi":"10.1016/j.cpc.2025.109812","DOIUrl":null,"url":null,"abstract":"<div><div>We introduce the <span>PrecisionLauricella</span> 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, <em>ε</em>. In practical multi-loop calculations, Lauricella functions are required only as series around <span><math><mi>ε</mi><mo>=</mo><mn>0</mn></math></span>, and <span>PrecisionLauricella</span> 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 <em>ε</em>-dependent 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 <span>PrecisionLauricella</span> package.</div></div><div><h3>Program summary</h3><div><em>Program Title:</em> PrecisionLauricella</div><div><em>CPC Library link to program files:</em> <span><span>https://doi.org/10.17632/6f958yz2dr.1</span><svg><path></path></svg></span></div><div><em>Developer's repository link:</em> <span><span>https://bitbucket.org/BezuglovMaxim/precisionlauricella-package/src/main/</span><svg><path></path></svg></span></div><div><em>Licensing provisions:</em> GPLv3</div><div><em>Programming language:</em> Wolfram Mathematica</div><div><em>Supplementary material:</em> PrecisionLauricella_Examples.nb</div><div><em>Nature of problem:</em> 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 <em>ε</em>, 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.</div><div>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.</div><div><em>Solution method:</em> 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.</div><div>We further optimize calculations by reconstructing <em>ε</em> dependencies from evaluations at specific values, enabling efficient parallelization and reducing computational costs.</div><div>A comprehensive mathematical exposition of the method is provided in our previous work [1].</div></div><div><h3>References</h3><div><ul><li><span>[1]</span><span><div>M. Bezuglov, B. Kniehl, A. Onishchenko, O. Veretin, High-precision numerical evaluation of Lauricella functions, <span><span>arXiv:2502.03276</span><svg><path></path></svg></span>, 2 2025.</div></span></li></ul></div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"316 ","pages":"Article 109812"},"PeriodicalIF":3.4000,"publicationDate":"2025-08-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computer Physics Communications","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0010465525003145","RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0

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, ε. In practical multi-loop calculations, Lauricella functions are required only as series around ε=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 ε-dependent 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 ε, 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 ε 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.
Lauricella precision:根据参数对Lauricella函数进行数值计算的软件包
我们介绍了PrecisionLauricella软件包,这是一个由Wolfram Mathematica开发的计算工具,用于高精度地计算Lauricella函数的Laurent展开系数,其参数线性依赖于一个小调节器ε。在实际的多环计算中,Lauricella函数只需要ε=0附近的级数,而PrecisionLauricella专门设计用于以任意精度提供这些系数。该软件包利用基于Frobenius广义幂级数的解析延拓方法,为依赖于多维级数展开或Mellin-Barnes表示的传统方法提供了一种高效而准确的替代方法。这种一维方法对高精度计算特别有利,并通过特定数值计算的ε依赖重构促进进一步优化,从而实现高效的并行化。该方法的基本数学框架已在我们之前的工作中详细介绍,而本文的重点是PrecisionLauricella软件包的设计、实现和实际应用。程序摘要程序标题:PrecisionLauricellaCPC库链接到程序文件:https://doi.org/10.17632/6f958yz2dr.1Developer's存储库链接:https://bitbucket.org/BezuglovMaxim/precisionlauricella-package/src/main/Licensing条款:gplv3编程语言:Wolfram mathematicas补充材料:PrecisionLauricella_Examples。问题性质:Lauricella函数,超几何函数的推广,出现在物理和数学中,包括费曼积分和弦理论。当它们的指标线性依赖于一个小参数ε时,由于高维级数和奇异性的复杂性,它们的数值计算变得困难。传统的方法,如超几何再展开或梅林-巴恩斯积分,往往缺乏效率和精度。管理多维求和增加了计算成本,特别是对于高精度要求,使得这些方法不适合许多实际应用。因此,迫切需要一种高效、可扩展的方法,能够保持数值精度并有效地处理参数依赖性。求解方法:采用Frobenius方法,通过广义幂级数实现Lauricella函数的解析延拓。与依赖于多维展开或Mellin-Barnes积分的传统方法相比,将函数表示为一维级数简化了高精度数值计算。我们通过从特定值的评估中重建ε依赖关系来进一步优化计算,从而实现高效的并行化并降低计算成本。在我们以前的工作中提供了对该方法的全面的数学阐述。张晓东,张晓东,张晓东,等。Lauricella函数的高精度数值计算,应用化学学报,34(2):357 - 357。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
求助全文
约1分钟内获得全文 求助全文
来源期刊
Computer Physics Communications
Computer Physics Communications 物理-计算机:跨学科应用
CiteScore
12.10
自引率
3.20%
发文量
287
审稿时长
5.3 months
期刊介绍: The focus of CPC is on contemporary computational methods and techniques and their implementation, the effectiveness of which will normally be evidenced by the author(s) within the context of a substantive problem in physics. Within this setting CPC publishes two types of paper. Computer Programs in Physics (CPiP) These papers describe significant computer programs to be archived in the CPC Program Library which is held in the Mendeley Data repository. The submitted software must be covered by an approved open source licence. Papers and associated computer programs that address a problem of contemporary interest in physics that cannot be solved by current software are particularly encouraged. Computational Physics Papers (CP) These are research papers in, but are not limited to, the following themes across computational physics and related disciplines. mathematical and numerical methods and algorithms; computational models including those associated with the design, control and analysis of experiments; and algebraic computation. Each will normally include software implementation and performance details. The software implementation should, ideally, be available via GitHub, Zenodo or an institutional repository.In addition, research papers on the impact of advanced computer architecture and special purpose computers on computing in the physical sciences and software topics related to, and of importance in, the physical sciences may be considered.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
确定
请完成安全验证×
copy
已复制链接
快去分享给好友吧!
我知道了
右上角分享
点击右上角分享
0
联系我们:info@booksci.cn Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。 Copyright © 2023 布克学术 All rights reserved.
京ICP备2023020795号-1
ghs 京公网安备 11010802042870号
Book学术文献互助
Book学术文献互助群
群 号:604180095
Book学术官方微信