Related papers: Reducing differential equations for multiloop mast…
The evolution of a continuous time Markov process with a finite number of states is usually calculated by the Master equation - a linear differential equations with a singular generator matrix. We derive a general method for reducing the…
We describe the application of differential reduction algorithms for Feynman Diagram calculation. We illustrate the procedure in the context of generalized hypergeometric functions, and give an example for a type of q-loop bubble diagram.
An algorithmic method to exploit a general class of infinitesimal symmetries for reducing stochastic differential equations is presented and a natural definition of reconstruction, inspired by the classical reconstruction by quadratures, is…
In this paper we derive and analyze an algorithm for inverting quaternion matrices. The algorithm is an analogue of the Frobenius algorithm for the complex matrix inversion. On the theory side, we prove that our algorithm is more efficient…
We consider deformations of $2\times2$ and $3\times3$ matrix linear ODEs with rational coefficients with respect to singular points of Fuchsian type which don't satisfy the well-known system of Schlesinger equations (or its natural…
We present a new method for the numerical evaluation of arbitrary loop integrals in dimensional regularization. We first derive Mellin-Barnes integral representations and apply an algorithmic technique, based on the Cauchy theorem, to…
We evaluate a four-loop conformal integral, i.e. an integral over four four-dimensional coordinates, by turning to its dimensionally regularized version and applying differential equations for the set of the corresponding 213 master…
We present in detail two resummation methods emerging from the application of the Simplified Differential Equations approach to a canonical basis of master integrals. The first one is a method which allows for an easy determination of the…
We consider the numerical solution of large-scale symmetric differential matrix Riccati equations. Under certain hypotheses on the data, reduced order methods have recently arisen as a promising class of solution strategies, by forming…
In the context of linear inverse problems, we propose and study a general iterative regularization method allowing to consider large classes of regularizers and data-fit terms. The algorithm we propose is based on a primal-dual diagonal…
A method for calculating phase-space master integrals for the decay process $1 \to n$ massless partons in QCD using integration-by-parts and differential equations techniques is discussed. The method is based on the appropriate choice of…
This paper describes a method of calculating the transforms, currently obtained via Fourier and reverse Fourier transforms. The method allows calculating efficiently the transforms of a signal having an arbitrary dimension of the digital…
Linear algebraic expressions are the essence of many computationally intensive problems, including scientific simulations and machine learning applications. However, translating high-level formulations of these expressions to efficient…
For a system of partial differential equations that has an extended Kovalevskaya form, a reduction procedure is presented that allows one to use a local (point, contact, or higher) symmetry of a system and a symmetry-invariant conservation…
Joint diagonalization, the process of finding a shared set of approximate eigenvectors for a collection of matrices, arises in diverse applications such as multidimensional harmonic analysis or quantum information theory. This task is…
In this paper, we present numerical procedures to compute solutions of partial differential equations posed on fractals. In particular, we consider the strong form of the equation using standard graph Laplacian matrices and also weak forms…
We extent the standard approach of dimensional regularization of Feynman diagrams: we replace the transition to lower dimensions by a 'natural' cut-off regulator. Introducing an external regulator of mass Lambda^(2e), we regain in the limit…
We propose a new algorithm for efficiently solving the damped Fisher matrix in large-scale scenarios where the number of parameters significantly exceeds the number of available samples. This problem is fundamental for natural gradient…
In the paper we deal with linear fractional control problems with constant delays in the state. Single-order systems with fractional derivative in Caputo sense of orders between 0 and 1 are considered. The aim is to introduce a new…
We define linearly reducible elliptic Feynman integrals, and we show that they can be algorithmically solved up to arbitrary order of the dimensional regulator in terms of a 1-dimensional integral over a polylogarithmic integrand, which we…