Related papers: Iteratively Reduce Auxiliary Scalar Product in Mul…

The computational cost associated with reducing tensor integrals to scalar integrals using the Passarino-Veltman method is dominated by the diagonalisation of large systems of equations. These systems of equations are sized according to the…

We show that, under certain circumstances, it is possible to automatically compute Jacobian-inverse-vector and Jacobian-inverse-transpose-vector products about as efficiently as Jacobian-vector and Jacobian-transpose-vector products. The…

We present here an efficient method which systematically reduces the rank of the augmented space and thereby helps to implement augmented space recursion for any real calculation. Our method is based on the symmetry of the Hamiltonian in…

In this work, we present the package {\sc NeatIBP}, which automatically generates small-size integration-by-parts (IBP) identities for Feynman integrals. Based on the syzygy and module intersection techniques, the generated IBP identities'…

We describe the unitarity approach for the numerical computation of two-loop integral coefficients of scattering amplitudes. It is well known that the leading propagator singularities of an amplitude's integrand are related to products of…

We present a new method to construct integration-by-part (IBP) identities from the viewpoint of differential geometry. Vectors for generating IBP identities are reformulated as differential forms, via Poincar\'{e} duality. Using the tools…

A general method for reducing tensor form factors, that appear in one-loop calculations in dimensional regularization, to scalar integrals is presented. The method is an extension of the reduction scheme introduced by Passarino and Veltman…

We present an efficient graphical approach to construct projectors for the tensor reduction of multi-loop Feynman integrals with both Lorentz and spinor indices in $D$ dimensions. An ansatz for the projectors is constructed making use of…

We recast Byerly's formula for integrals of products of Legendre polynomials. Then we adopt the idea to the case of Jacobi polynomials. After that, we use the formula to derive an asymptotic formula for integrals of products of Jacobi…

Convergence problems in coupled-cluster iterations are discussed, and a new iteration scheme is proposed. Whereas the Jacobi method inverts only the diagonal part of the large matrix of equation coefficients, we invert a matrix which also…

We present a Mathematica package which finds a basis of master integrals for the Feynman integral reduction. In this basis the dependence on the dimensional regularization in the denominators factorizes in kinematic independent polynomials.

Based on the method in Refs.~{\tt [D.~Kreimer, Z.\ Phys.\ C {\bf 54} (1992) 667} and {\tt Int.\ J.\ Mod.\ Phys.\ A {\bf 8} (1993) 1797]}, we present analytic results for scalar one-loop four-point Feynman integrals with complex internal…

This paper proposes several explicit and implicit multistep frequency response optimized integrators considering first or second order derivative. A prediction-based method aiming at accelerating a novel power system transient simulation…

We study the problem of solving integration-by-parts recurrence relations for a given class of Feynman integrals which is characterized by an arbitrary polynomial in the numerator and arbitrary integer powers of propagators, {\it i.e.}, the…

For solving the continuous Sylvester equation, a class of the multiplicative splitting iteration method is presented. We consider two symmetric positive definite splittings for each coefficient matrix of the continuous Sylvester equations…

We study the symmetric outer product decomposition which decomposes a fully (partially) symmetric tensor into a sum of rank-one fully (partially) symmetric tensors. We present iterative algorithms for the third-order partially symmetric…

This paper develops new combinatorial approaches to analyze and compute special set partitions, called complementary set partitions, which are fundamental in the study of generalized cumulants. Moving away from traditional graph-based and…

We construct product formulas for exponentials of commutators and explore their applications. First, we directly construct a third-order product formula with six exponentials by solving polynomial equations obtained using the operator…

We extend the notion of generalized unitarity cuts to accommodate loop integrals with higher powers of propagators. Such integrals frequently arise in for example integration-by-parts identities, Schwinger parametrizations and Mellin-Barnes…

We propose new linear combinations of compositions of a basic second-order scheme with appropriately chosen coefficients to construct higher order numerical integrators for differential equations. They can be considered as a generalization…