Related papers: GPU Implementation of Zippel Method for Feynman In…
This article presents a strongly polynomial-time algorithm for the general linear programming problem. This algorithm is an implicit reduction procedure that works as follows. Primal and dual problems are combined into a special system of…
We study partial fraction decompositions (PFDs) in several variables using tools from commutative algebra. We give criteria for when a rational function with poles on a hyperplane arrangement has a desirable PFD. Our criteria are obtained…
Reconstructing a complete object from its parts is a fundamental problem in many scientific domains. The purpose of this article is to provide a systematic survey on this topic. The reassembly problem requires understanding the attributes…
We propose a new lifting and recombination scheme for rational bivariate polynomial factorization that takes advantage of the Newton polytope geometry. We obtain a deterministic algorithm that can be seen as a sparse version of an algorithm…
This paper presents a new major release of the program FIESTA (Feynman Integral Evaluation by a Sector decomposiTion Approach). The new release is mainly aimed at optimal performance at large scales when one is increasing the number of…
A modular application of the integration by fractional expansion (IBFE) method for evaluating Feynman diagrams is extended to diagrams that contain loop triangle subdiagrams in their geometry. The technique is based in the replacement of…
A refinement of manifold data is a computational process, which produces a denser set of discrete data from a given one. Such refinements are closely related to multiresolution representations of manifold data by pyramid transforms, and…
In this talk, we discuss how ideas from geometry help to improve Feynman integral reduction and the construction of $\varepsilon$-factorised differential equations. In particular, we outline a systematic procedure to obtain an…
The generalized minimal residual (GMRES) algorithm is applied to image reconstruction using linear computed tomography (CT) models. The GMRES algorithm iteratively solves square, non-symmetric linear systems and it has practical application…
We describe an algorithm for splitting permutation representations of finite group over fields of characteristic zero into irreducible components. The algorithm is based on the fact that the components of the invariant inner product in…
Iterative phase retrieval methods based on the Gerchberg-Saxton (GS) or Fienup algorithm require a large number of iterations to converge to a meaningful solution. For complex-valued or phase objects, these approaches also suffer from…
Optimization on the Stiefel manifold or with orthogonality constraints is an important problem in many signal processing and data analysis applications such as Sparse Principal Component Analysis (SPCA). Algorithms such as the Riemannian…
In this updated vesion of ZMCintegral, we have added the functionality of integrations with parameter scan on distributed Graphics Processing Units(GPUs). Given a large parameter grid (up to 10^{10} parameter points to be scanned), the code…
Phase space cuts are implemented by inserting Heaviside theta functions in the integrands of momentum-space Feynman integrals. By directly parametrizing theta functions and constructing integration-by-parts (IBP) identities in the…
The Floyd-Warshall(FW) algorithm, is an ancient but a largely important algorithm used to solve the all-pairs simple-paths(APSP) problem. While the algorithm is available for use in open-source graph optimization libraries such as NetworkX,…
Selected inversion is essential for applications such as Bayesian inference, electronic structure calculations, and inverse covariance estimation, where computing only specific elements of large sparse matrix inverses significantly reduces…
We present an efficient algorithm for calculating multiloop Feynman integrals perturbatively.
We present an alternative technique for evaluating multiloop Feynman diagrams, using the integration by fractional expansion method. Here we consider generic diagrams that contain propagators with radiative corrections which topologically…
We present a new method for the decomposition of multi-loop Euclidean Feynman integrals into quasi-finite Feynman integrals. These are defined in shifted dimensions with higher powers of the propagators, make explicit both infrared and…
With the rapidly increasing rate of microlensing planet detections, microlensing modeling software faces significant challenges in computation efficiency. Here, we develop the Twinkle code, an efficient and robust binary-lens modeling…