pith. sign in

arxiv: 2605.31553 · v1 · pith:N5BBTC2Hnew · submitted 2026-05-29 · 🧮 math-ph · math.MP

Numerical analytical continuation of multivariate hypergeometric functions

Pith reviewed 2026-06-28 19:53 UTC · model grok-4.3

classification 🧮 math-ph math.MP
keywords multivariate hypergeometric functionsPfaffian systemsLaporta reductionFrobenius methodanalytic continuationholonomic systemsnumerical evaluation
0
0 comments X

The pith

A framework evaluates multivariate hypergeometric functions to high precision and continues them analytically by building Pfaffian systems with Laporta reduction then applying the Frobenius method.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a general framework for the high-precision numerical evaluation of multivariate hypergeometric functions, which satisfy holonomic systems of partial differential equations. It adapts techniques from multi-loop Feynman integral calculations by constructing Pfaffian systems for arbitrary cases through application of the Laporta reduction algorithm to suitable differential relations. A numerical scheme based on the Frobenius method then produces local power-series solutions with controlled precision and transports those solutions along chosen paths in the space of variables. The work includes a systematic treatment of multivaluedness, showing how the Frobenius method can reach different Riemann sheets and track changes under analytic continuation around singular loci.

Core claim

The central claim is that Pfaffian systems for arbitrary multivariate hypergeometric functions can be obtained by Laporta reduction on differential relations, after which a Frobenius-based numerical scheme computes local power series to controlled precision, transports solutions along paths, and handles branch structure by accessing different Riemann sheets during continuation around singular loci.

What carries the argument

Pfaffian systems constructed via the Laporta reduction algorithm applied to systems of differential relations, together with the Frobenius method for generating and transporting local power-series solutions.

If this is right

  • High-precision numerical values become available for solutions of holonomic PDE systems that define multivariate hypergeometric functions.
  • Solutions can be transported along arbitrary paths in the space of variables while preserving controlled accuracy.
  • Different Riemann sheets can be reached systematically by the same series construction around singular points.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same reduction-plus-Frobenius pipeline might be applied to other classes of holonomic functions beyond the hypergeometric case.
  • Efficiency gains could follow from combining the Laporta step with existing computer-algebra tools for holonomic systems.
  • Direct numerical checks on known low-variable cases would confirm whether the constructed Pfaffian systems reproduce established special-function identities.

Load-bearing premise

The Laporta reduction algorithm can be applied to suitable systems of differential relations to construct Pfaffian systems for arbitrary multivariate hypergeometric functions.

What would settle it

A concrete multivariate hypergeometric function for which the Laporta-derived Pfaffian system yields series solutions whose precision cannot be controlled or whose branch tracking fails to match known analytic properties under continuation around a singular locus.

Figures

Figures reproduced from arXiv: 2605.31553 by A.I. Onishchenko, B.A. Kniehl, M.A. Bezuglov, O.L. Veretin.

Figure 1
Figure 1. Figure 1: Graphical representation of the linear system for the Horn function [PITH_FULL_IMAGE:figures/full_fig_p011_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Identification of master derivatives for the function [PITH_FULL_IMAGE:figures/full_fig_p012_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: One of the possible paths with complex singular points is shown. The singular points [PITH_FULL_IMAGE:figures/full_fig_p015_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: This figure illustrates the path of the analytic continuation of the system in Eq. ( [PITH_FULL_IMAGE:figures/full_fig_p017_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Examples of Feynman integrals expressible in terms of multivariate hypergeometric [PITH_FULL_IMAGE:figures/full_fig_p018_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Schematic picture of analytic continuation in the [PITH_FULL_IMAGE:figures/full_fig_p023_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Geometric picture of monodromy for the Appell function [PITH_FULL_IMAGE:figures/full_fig_p026_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Riemann sphere representation of the monodromy at infinity. The north pole corre [PITH_FULL_IMAGE:figures/full_fig_p028_8.png] view at source ↗
read the original abstract

We present a general framework for the high-precision numerical evaluation of multivariate hypergeometric functions defined as solutions of holonomic systems of partial differential equations. Our approach adapts and extends methods originally developed for multi-loop Feynman integrals to the setting of hypergeometric functions of many variables. In particular, we construct Pfaffian systems for arbitrary multivariate hypergeometric functions by applying the Laporta reduction algorithm to suitable systems of differential relations. Next, we construct a numerical scheme based on the Frobenius method, which allows us to compute local power-series solutions with controlled precision and to transport them along prescribed paths in the space of variables. A central part of the paper is devoted to a systematic analysis of multivaluedness and branch structure: we show how the Frobenius method can be used to access different Riemann sheets in a controlled way and to track changes of the solution under analytic continuation around singular loci.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 2 minor

Summary. The paper presents a general framework for high-precision numerical evaluation of multivariate hypergeometric functions defined as solutions of holonomic systems of PDEs. It adapts Laporta reduction (originally for Feynman integrals) to construct Pfaffian systems from suitable differential relations, then applies the Frobenius method to obtain local power-series solutions that are transported along paths in variable space while systematically tracking multivaluedness, branch structure, and monodromy around singular loci.

Significance. If the described procedure holds, the work supplies an explicit, parameter-free algorithmic scheme for arbitrary multivariate hypergeometric functions, with controllable precision and explicit monodromy tracking. This extends existing single-variable or special-case methods and directly addresses a practical need in mathematical physics for reproducible high-precision continuation of multivalued functions.

minor comments (2)
  1. The abstract refers to 'suitable systems of differential relations' without naming the generating set; the main text should include an explicit example of the input relations for a non-trivial multivariate case (e.g., Appell F1 or Lauricella FD).
  2. Notation for the Pfaffian matrix and the connection form should be introduced once with a clear definition before the first algorithmic step.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive summary, significance assessment, and recommendation to accept the manuscript. No major comments were raised.

Circularity Check

0 steps flagged

No significant circularity

full rationale

The manuscript describes an explicit algorithmic pipeline: Laporta reduction applied to a generating set of differential relations from the holonomic ideal yields a Pfaffian system, after which the Frobenius method produces local series solutions that are continued along paths with monodromy tracking. All steps are presented as standard, parameter-free operations on the input holonomic system; no equation is shown to be defined in terms of its own output, no fitted parameter is relabeled as a prediction, and no load-bearing uniqueness claim rests on a self-citation chain. The derivation is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review; no explicit free parameters, axioms, or invented entities are stated or derivable from the given text.

pith-pipeline@v0.9.1-grok · 5689 in / 1058 out tokens · 20464 ms · 2026-06-28T19:53:26.045682+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

109 extracted references · 67 canonical work pages · 37 internal anchors

  1. [1]

    Weinzierl, Feynman Integrals

    S. Weinzierl, Feynman Integrals. 1 2022, 2201.03593

  2. [2]

    Mellin-Barnes Integrals: A Primer on Particle Physics Applications,

    I. Dubovyk, J. Gluza, and G. Somogyi, “Mellin-Barnes Integrals: A Primer on Particle Physics Applications,” Lect. Notes Phys. , vol. 1008, p. pp., 2022, 2211.13733

  3. [3]

    Smirnov, Feynman Integral Calculus

    V. Smirnov, Feynman Integral Calculus . Springer Berlin Heidelberg, 2006

  4. [4]

    Hypergeometric representation of the two-loop equal mass sunrise diagram

    O. V. Tarasov, “Hypergeometric representation of the two-loop equal mass sunrise dia- gram,” Phys. Lett. B , vol. 638, pp. 195–201, 2006, hep-ph/0603227

  5. [5]

    Space-time dimensionality D as complex variable: calculating loop integrals using dimensional recurrence relation and analytical properties with respect to D

    R. N. Lee, “Space-time dimensionality D as complex variable: Calculating loop integrals using dimensional recurrence relation and analytical properties with respect to D,” Nucl. Phys. B , vol. 830, pp. 474–492, 2010, 0911.0252

  6. [6]

    DRA method: Powerful tool for the calculation of the loop integrals

    R. N. Lee, “DRA method: Powerful tool for the calculation of the loop integrals,” J. Phys. Conf. Ser. , vol. 368, p. 012050, 2012, 1203.4868

  7. [7]

    Functional reduction of one-loop Feynman integrals with arbitrary masses,

    O. V. Tarasov, “Functional reduction of one-loop Feynman integrals with arbitrary masses,” JHEP, vol. 06, p. 155, 2022, 2203.00143

  8. [8]

    On Series and Integral Represen- tations of Some NRQCD Master Integrals,

    M. A. Bezuglov, A. V. Kotikov, and A. I. Onishchenko, “On Series and Integral Represen- tations of Some NRQCD Master Integrals,” JETP Lett. , vol. 116, no. 1, pp. 61–69, 2022, 2205.14115

  9. [9]

    Non-planar elliptic vertex,

    M. A. Bezuglov and A. I. Onishchenko, “Non-planar elliptic vertex,” JHEP, vol. 04, p. 045, 2022, 2112.05096

  10. [10]

    HypergeometricstructuresinFeyn- man integrals

    J. Blümlein, M. Saragnese, and C. Schneider, “Hypergeometric Structures in Feynman Integrals,” Ann. Math. Artif. Intell. , vol. 91, pp. 591–649, 2023, 2111.15501. 34

  11. [11]

    Feynman amplitudes as power series,

    D. S. Kershaw, “Feynman amplitudes as power series,” Phys. Rev. D, vol. 8, pp. 2708–2713, 1973

  12. [12]

    One loop integrals revisited. 1. The Two point functions,

    D. Kreimer, “One loop integrals revisited. 1. The Two point functions,” Z. Phys. C , vol. 54, pp. 667–672, 1992

  13. [13]

    One loop integrals revisited. 2. The Three point functions,

    D. Kreimer, “One loop integrals revisited. 2. The Three point functions,” Int. J. Mod. Phys. A , vol. 8, pp. 1797–1814, 1993

  14. [14]

    Loop Integrals, R Functions and their Analytic Continuation

    L. Brucher, J. Franzkowski, and D. Kreimer, “Loop integrals, R functions and their ana- lytic continuation,” Mod. Phys. Lett. A , vol. 9, pp. 2335–2346, 1994, hep-th/9307055

  15. [15]

    A Method of evaluating massive Feynman integrals,

    E. E. Boos and A. I. Davydychev, “A Method of evaluating massive Feynman integrals,” Theor. Math. Phys. , vol. 89, pp. 1052–1063, 1991

  16. [16]

    General results for massive N point Feynman diagrams with different masses,

    A. I. Davydychev, “General results for massive N point Feynman diagrams with different masses,” J. Math. Phys. , vol. 33, pp. 358–369, 1992

  17. [17]

    Two-loop two-point functions with masses: asymptotic expansions and Taylor series, in any dimension

    D. J. Broadhurst, J. Fleischer, and O. V. Tarasov, “Two loop two point functions with masses: Asymptotic expansions and Taylor series, in any dimension,” Z. Phys. C , vol. 60, pp. 287–302, 1993, hep-ph/9304303

  18. [18]

    Hypergeometric Functions and Feynman Diagrams

    M. Kalmykov, V. Bytev, B. A. Kniehl, S.-O. Moch, B. F. L. Ward, and S. A. Yost, “Hyper- geometric Functions and Feynman Diagrams,” in Antidifferentiation and the Calculation of Feynman Amplitudes , 12 2020, 2012.14492

  19. [19]

    HypExp, a Mathematica package for expanding hypergeometric functions around integer-valued parameters

    T. Huber and D. Maitre, “HypExp: A Mathematica package for expanding hypergeometric functions around integer-valued parameters,” Comput. Phys. Commun. , vol. 175, pp. 122– 144, 2006, hep-ph/0507094

  20. [20]

    HypExp 2, Expanding Hypergeometric Functions about Half-Integer Parameters

    T. Huber and D. Maitre, “HypExp 2, Expanding Hypergeometric Functions about Half- Integer Parameters,” Comput. Phys. Commun. , vol. 178, pp. 755–776, 2008, 0708.2443

  21. [21]

    - XSummer - Transcendental Functions and Symbolic Summation in Form

    S. Moch and P. Uwer, “XSummer: Transcendental functions and symbolic summation in form,” Comput. Phys. Commun. , vol. 174, pp. 759–770, 2006, math-ph/0508008

  22. [22]

    Symbolic Expansion of Transcendental Functions

    S. Weinzierl, “Symbolic expansion of transcendental functions,” Comput. Phys. Commun. , vol. 145, pp. 357–370, 2002, math-ph/0201011

  23. [23]

    Analytic and Algorithmic Aspects of Generalized Harmonic Sums and Polylogarithms

    J. Ablinger, J. Blümlein, and C. Schneider, “Analytic and Algorithmic Aspects of Gen- eralized Harmonic Sums and Polylogarithms,” J. Math. Phys. , vol. 54, p. 082301, 2013, 1302.0378

  24. [24]

    NumExp: Numerical epsilon expansion of hypergeometric functions

    Z.-W. Huang and J. Liu, “NumExp: Numerical epsilon expansion of hypergeometric func- tions,” Comput. Phys. Commun. , vol. 184, pp. 1973–1980, 2013, 1209.3971

  25. [25]

    MultiHypExp: A Mathematica package for expanding multivariate hy- pergeometric functions in terms of multiple polylogarithms

    S. Bera, “MultiHypExp: A Mathematica package for expanding multivariate hypergeo- metric functions in terms of multiple polylogarithms,” Comput. Phys. Commun. , vol. 297, p. 109060, 2024, 2306.11718. 35

  26. [26]

    Expansion of hypergeometric functions in terms of polylogarithms with a nontrivial change of variables

    M. A. Bezuglov and A. I. Onishchenko, “Expansion of hypergeometric functions in terms of polylogarithms with a nontrivial change of variables,” Theor. Math. Phys. , vol. 219, no. 3, pp. 871–896, 2024, 2312.06242

  27. [27]

    A new approach to the epsilon expansion of generalized hypergeometric functions

    D. Greynat and J. Sesma, “A new approach to the epsilon expansion of generalized hyper- geometric functions,” Comput. Phys. Commun. , vol. 185, pp. 472–478, 2014, 1302.2423

  28. [28]

    Derivatives of the Pochhammer and reciprocal Pochhammer symbols and their use in epsilon-expansions of Appell and Kampé de Fériet functions,

    D. Greynat, J. Sesma, and G. Vulvert, “Derivatives of the Pochhammer and reciprocal Pochhammer symbols and their use in epsilon-expansions of Appell and Kampé de Fériet functions,” J. Math. Phys. , vol. 55, p. 043501, 2014

  29. [29]

    pySecDec: a toolbox for the numerical evaluation of multi-scale integrals

    S. Borowka, G. Heinrich, S. Jahn, S. P. Jones, M. Kerner, J. Schlenk, and T. Zirke, “pySecDec: a toolbox for the numerical evaluation of multi-scale integrals,” Comput. Phys. Commun., vol. 222, pp. 313–326, 2018, 1703.09692

  30. [30]

    Fonctions hypergéométriques et hypersphériques. Polynômes d’Hermite

    P. Appell and J. Kampé de Fériet, “Fonctions hypergéométriques et hypersphériques. Polynômes d’Hermite.. ” Paris: Gauthier-Villars. vii, 434 p. (1926)., 1926

  31. [31]

    Hypergeometric functions of two variables,

    A. Erdélyi, “Hypergeometric functions of two variables,” Acta Math., vol. 83, pp. 131–164, 1950

  32. [32]

    Integration of the partial differential equations for the hypergeometric functions F1 and FD of two and more variables,

    P. O. M. Olsson, “Integration of the partial differential equations for the hypergeometric functions F1 and FD of two and more variables,” J. Math. Phys. , vol. 5, pp. 420–430, 1964

  33. [33]

    Multiple hypergeometric functions and applications

    H. Exton, “Multiple hypergeometric functions and applications. ” Mathematics and Its Applications. Chichester: Ellis Horwood Limited, Publisher; New York etc.: Halsted Press, a division of John Wiley & Sons, Inc. 312 p., 1976

  34. [34]

    Olsson.wl: A Mathematica Pack- age for the Computation of Linear Transformations of Multivariable Hypergeometric Func- tions,

    B. Ananthanarayan, S. Bera, S. Friot, and T. Pathak, “Olsson.wl: A Mathematica Pack- age for the Computation of Linear Transformations of Multivariable Hypergeometric Func- tions,” Springer Proc. Phys. , vol. 304, pp. 738–740, 2024

  35. [35]

    Analytic continuations and numerical evaluation of the Appell F1, F3, Lauricella FD(3) and Lauricella-Saran FS(3) and their application to Feynman integrals

    S. Bera and T. Pathak, “Analytic continuations and numerical evaluation of the Appell F1, F3, Lauricella FD(3) and Lauricella-Saran FS(3) and their application to Feynman integrals,” Comput. Phys. Commun. , vol. 306, p. 109386, 2025, 2403.02237

  36. [36]

    Analytic continuation of appell’s hypergeometric series f 2 to the vicinity of the singular point x= 1, y= 1,

    G. Hahne, “Analytic continuation of appell’s hypergeometric series f 2 to the vicinity of the singular point x= 1, y= 1,” Journal of Mathematical Physics , vol. 10, no. 3, pp. 524–531, 1969

  37. [37]

    Analytic continuation formulas and Jacobi-type relations for Lauricella function,

    S. I. Bezrodnykh, “Analytic continuation formulas and Jacobi-type relations for Lauricella function,” Dokl. Math. , vol. 93, no. 2, pp. 129–134, 2016

  38. [38]

    On the analytic continuation of the Lauricella function F (N ) D ,

    S. I. Bezrodnykh, “On the analytic continuation of the Lauricella function F (N ) D ,” Math. Notes, vol. 100, no. 2, pp. 318–324, 2016

  39. [39]

    Analytic continuation of the appell function f 1 and integration of the associated system of equations in the logarithmic case,

    S. I. Bezrodnykh, “Analytic continuation of the appell function f 1 and integration of the associated system of equations in the logarithmic case,” Computational Mathematics and Mathematical Physics , vol. 57, pp. 559–589, 2017. 36

  40. [40]

    Analytic continuation of the Lauricella function F (N ) D with arbitrary number of variables,

    S. I. Bezrodnykh, “Analytic continuation of the Lauricella function F (N ) D with arbitrary number of variables,” Integral Transforms Spec. Funct., vol. 29, no. 1, pp. 21–42, 2018

  41. [41]

    Analytic continuation of Lauricella’s functions F (N ) A , F (N ) B and F (N ) D ,

    S. I. Bezrodnykh, “Analytic continuation of Lauricella’s functions F (N ) A , F (N ) B and F (N ) D ,” Integral Transforms Spec. Funct., vol. 31, no. 11, pp. 921–940, 2020

  42. [42]

    The Lauricella hypergeometric function F (N ) D , the Riemann-Hilbert problem, and some applications,

    S. I. Bezrodnykh, “The Lauricella hypergeometric function F (N ) D , the Riemann-Hilbert problem, and some applications,” Russ. Math. Surv. , vol. 73, no. 6, pp. 941–1031, 2018

  43. [43]

    HYPERDIRE: HYPERgeometric functions DIfferential REduction: MATHEMATICA based packages for differential reduction of generalized hypergeometric functions pFq, F1,F2,F3,F4

    V. V. Bytev, M. Y. Kalmykov, and B. A. Kniehl, “HYPERDIRE, HYPERgeometric func- tions DIfferential REduction: MATHEMATICA-based packages for differential reduction of generalized hypergeometric functions pFp−1, F1,F2,F3,F4,” Comput. Phys. Commun. , vol. 184, pp. 2332–2342, 2013, 1105.3565

  44. [44]

    V. V. Bytev and B. A. Kniehl, “HYPERDIRE HYPERgeometric functions DIfferential REduction: Mathematica-based packages for the differential reduction of generalized hy- pergeometric functions: Horn-type hypergeometric functions of two variables,” Comput. Phys. Commun. , vol. 189, pp. 128–154, 2015, 1309.2806

  45. [45]

    V. V. Bytev, M. Y. Kalmykov, and S.-O. Moch, “HYPERgeometric functions DIfferential REduction (HYPERDIRE): MATHEMATICA based packages for differential reduction of generalized hypergeometric functions: FD and FS Horn-type hypergeometric functions of three variables,” Comput. Phys. Commun. , vol. 185, pp. 3041–3058, 2014, 1312.5777

  46. [46]

    V. V. Bytev and B. A. Kniehl, “HYPERDIRE—HYPERgeometric functions DIfferen- tial REduction: Mathematica-based packages for the differential reduction of generalized hypergeometric functions: Lauricella function Fc of three variables,” Comput. Phys. Com- mun., vol. 206, pp. 78–83, 2016, 1602.00917

  47. [47]

    Saito, B

    M. Saito, B. Sturmfels, and N. Takayama, Gröbner deformations of hypergeometric dif- ferential equations, vol. 6 of Algorithms Comput. Math. Berlin: Springer, 2000

  48. [48]

    Differential equations method: New technique for massive Feynman dia- grams calculation,

    A. V. Kotikov, “Differential equations method: New technique for massive Feynman dia- grams calculation,” Phys. Lett. , vol. B254, pp. 158–164, 1991

  49. [49]

    New method of massive Feynman diagrams calculation,

    A. V. Kotikov, “New method of massive Feynman diagrams calculation,” Mod. Phys. Lett., vol. A6, pp. 677–692, 1991

  50. [50]

    Differential equations method: The Calculation of vertex type Feynman diagrams,

    A. V. Kotikov, “Differential equations method: The Calculation of vertex type Feynman diagrams,” Phys. Lett. , vol. B259, pp. 314–322, 1991

  51. [51]

    Differential equation method: The Calculation of N point Feynman dia- grams,

    A. V. Kotikov, “Differential equation method: The Calculation of N point Feynman dia- grams,” Phys. Lett., vol. B267, pp. 123–127, 1991. [Erratum: Phys. Lett.B295,409(1992)]

  52. [52]

    Differential Equations for Feynman Graph Amplitudes,

    E. Remiddi, “Differential equations for Feynman graph amplitudes,” Nuovo Cim. , vol. A110, pp. 1435–1452, 1997, hep-th/9711188

  53. [53]

    Multiloop integrals in dimensional regularization made simple

    J. M. Henn, “Multiloop integrals in dimensional regularization made simple,” Phys. Rev. Lett., vol. 110, p. 251601, 2013, 1304.1806. 37

  54. [54]

    Reducing differential equations for multiloop master integrals

    R. N. Lee, “Reducing differential equations for multiloop master integrals,” JHEP, vol. 04, p. 108, 2015, 1411.0911

  55. [55]

    A Theorem on Analytical Calculability of Four Loop Renormalization Group Functions,

    F. V. Tkachov, “A Theorem on Analytical Calculability of Four Loop Renormalization Group Functions,” Phys. Lett. , vol. 100B, pp. 65–68, 1981

  56. [56]

    Integration by Parts: The Algorithm to Calculate beta Functions in 4 Loops,

    K. G. Chetyrkin and F. V. Tkachov, “Integration by Parts: The Algorithm to Calculate beta Functions in 4 Loops,” Nucl. Phys. , vol. B192, pp. 159–204, 1981

  57. [57]

    High- precision numerical evaluation of Lauricella functions

    M. A. Bezuglov, B. A. Kniehl, A. I. Onishchenko, and O. L. Veretin, “High-precision numerical evaluation of Lauricella functions,” Nucl. Phys. B , vol. 1018, p. 116994, 2025, 2502.03276

  58. [58]

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

    M. A. Bezuglov, B. A. Kniehl, A. I. Onishchenko, and O. L. Veretin, “PrecisionLauricella: Package for numerical computation of Lauricella functions depending on a parameter,” Comput. Phys. Commun. , vol. 316, p. 109812, 2025, 2502.07935

  59. [59]

    High-precision calculation of multi-loop Feynman integrals by difference equations

    S. Laporta, “High-precision calculation of multiloop Feynman integrals by difference equa- tions,” Int. J. Mod. Phys. A , vol. 15, pp. 5087–5159, 2000, hep-ph/0102033

  60. [60]

    Aomoto and M

    K. Aomoto and M. Kita, Theory of Hypergeometric Functions . Springer Monographs in Mathematics, Springer, 2011

  61. [61]

    Kira - A Feynman Integral Reduction Program

    P. Maierhöfer, J. Usovitsch, and P. Uwer, “Kira—A Feynman integral reduction program,” Comput. Phys. Commun. , vol. 230, pp. 99–112, 2018, 1705.05610

  62. [62]

    Kira 1.2 Release Notes

    P. Maierhöfer and J. Usovitsch, “Kira 1.2 Release Notes,” 12 2018, 1812.01491

  63. [63]

    Lange, J

    F. Lange, J. Usovitsch, and Z. Wu, “Kira 3: integral reduction with efficient seeding and optimized equation selection,” 5 2025, 2505.20197

  64. [64]

    Automatic Integral Reduction for Higher Order Perturbative Calculations

    C. Anastasiou and A. Lazopoulos, “Automatic integral reduction for higher order pertur- bative calculations,” JHEP, vol. 07, p. 046, 2004, hep-ph/0404258

  65. [65]

    Blade: A package for block-triangular form improved Feynman integrals decomposition,

    X. Guan, X. Liu, Y.-Q. Ma, and W.-H. Wu, “Blade: A package for block-triangular form improved Feynman integrals decomposition,” Comput. Phys. Commun. , vol. 310, p. 109538, 2025, 2405.14621

  66. [66]

    Algorithm FIRE -- Feynman Integral REduction

    A. V. Smirnov, “Algorithm FIRE – Feynman Integral REduction,” JHEP, vol. 10, p. 107, 2008, 0807.3243

  67. [67]

    FIRE4, LiteRed and accompanying tools to solve integration by parts relations

    A. V. Smirnov and V. A. Smirnov, “FIRE4, LiteRed and accompanying tools to solve integration by parts relations,” Comput. Phys. Commun. , vol. 184, pp. 2820–2827, 2013, 1302.5885

  68. [68]

    FIRE5: a C++ implementation of Feynman Integral REduction

    A. V. Smirnov, “FIRE5: A C++ implementation of Feynman Integral REduction,” Com- put. Phys. Commun. , vol. 189, pp. 182–191, 2015, 1408.2372

  69. [69]

    Reduze - Feynman Integral Reduction in C++

    C. Studerus, “Reduze – Feynman integral reduction in C++,” Comput. Phys. Commun. , vol. 181, pp. 1293–1300, 2010, 0912.2546. 38

  70. [70]

    Reduze 2 - Distributed Feynman Integral Reduction

    A. von Manteuffel and C. Studerus, “Reduze 2 - Distributed Feynman Integral Reduction,” 1 2012, 1201.4330

  71. [71]

    On the computation of finite bottom-quark mass effects in Higgs boson production

    R. Mueller and D. G. Öztürk, “On the computation of finite bottom-quark mass effects in Higgs boson production,” JHEP, vol. 08, p. 055, 2016, arXiv:1512.08570

  72. [72]

    Two-loop $gg \to Hg$ amplitude mediated by a nearly massless quark

    K. Melnikov, L. Tancredi, and C. Wever, “Two-loop gg ! Hg amplitude mediated by a nearly massless quark,” JHEP, vol. 11, p. 104, 2016, arXiv:1610.03747

  73. [73]

    Three-loop massive tadpoles and polylogarithms through weight six

    B. A. Kniehl, A. F. Pikelner, and O. L. Veretin, “Three-loop massive tadpoles and poly- logarithms through weight six,” JHEP, vol. 08, p. 024, 2017, arXiv:1705.05136

  74. [74]

    Solving differential equations for Feynman integrals by expansions near singular points

    R. N. Lee, A. V. Smirnov, and V. A. Smirnov, “Solving differential equations for Feynman integrals by expansions near singular points,” JHEP, vol. 03, p. 008, 2018, arXiv:1709.07525

  75. [75]

    Evaluating `elliptic' master integrals at special kinematic values: using differential equations and their solutions via expansions near singular points

    R. N. Lee, A. V. Smirnov, and V. A. Smirnov, “Evaluating ‘elliptic’ master integrals at special kinematic values: using differential equations and their solutions via expansions near singular points,” JHEP, vol. 07, p. 102, 2018, arXiv:1805.00227

  76. [76]

    Two-loop dia- grams in non-relativistic QCD with elliptics,

    B. A. Kniehl, A. V. Kotikov, A. I. Onishchenko, and O. L. Veretin, “Two-loop dia- grams in non-relativistic QCD with elliptics,” Nucl. Phys. , vol. B948, p. 114780, 2019, arXiv:1907.04638

  77. [77]

    Generalised power series expansions for the elliptic planar families of Higgs + jet production at two loops

    F. Moriello, “Generalised power series expansions for the elliptic planar families of Higgs + jet production at two loops,” JHEP, vol. 01, p. 150, 2020, 1907.13234

  78. [78]

    Evaluating a family of two-loop non-planar master integrals for Higgs + jet production with full heavy-quark mass dependence,

    R. Bonciani, V. Del Duca, H. Frellesvig, J. M. Henn, M. Hidding, L. Maestri, F. Moriello, G. Salvatori, and V. A. Smirnov, “Evaluating a family of two-loop non-planar master integrals for Higgs + jet production with full heavy-quark mass dependence,” JHEP, vol. 01, p. 132, 2020, 1907.13156

  79. [79]

    The complete set of two-loop master integrals for Higgs + jet production in QCD,

    H. Frellesvig, M. Hidding, L. Maestri, F. Moriello, and G. Salvatori, “The complete set of two-loop master integrals for Higgs + jet production in QCD,” JHEP, vol. 06, p. 093, 2020, 1911.06308

  80. [80]

    DiffExp, a Mathematica package for computing Feynman integrals in terms of one-dimensional series expansions

    M. Hidding, “DiffExp, a Mathematica package for computing Feynman integrals in terms of one-dimensional series expansions,” Comput. Phys. Commun. , vol. 269, p. 108125, 2021, 2006.05510

Showing first 80 references.