Abstract:We analyze the decomposition problem of multivariate polynomial-exponential functions from truncated series and present new algorithms to compute their decomposition. Using the duality between polynomials and formal power series, we first show how the elements in the dual of an Artinian algebra correspond to polynomial-exponential functions. They are also the solutions of systems of partial differential equations with constant coefficients. We relate their representation to the inverse system of the roots of the characteristic variety. Using the properties of Hankel operators, we establish a correspondence between polynomial exponential series and Artinian Gorenstein algebras. We generalize Kronecker theorem to the multivariate case, by showing that the symbol of a Hankel operator of finite rank is a polynomial-exponential series and by connecting the rank of the Hankel operator with the decomposition of the symbol. A generalization of Prony's approach to multivariate decomposition problems is presented , exploiting eigenvector methods for solving polynomial equations. We show how to compute the frequencies and weights of a minimal polynomial-exponential decomposition , using the first coefficients of the series. A key ingredient of the approach is the flat extension criteria, which leads to a multivariate generalization of a rank condition for a Carath{é}odory-Fej{é}r decomposition of multivariate Hankel matrices. A new algorithm is given to compute a basis of the Artinian Gorenstein algebra, based on a Gram-Schmidt orthogonalization process and to decompose polynomial-exponential series. A general framework for the applications of this approach is described and illustrated in different problems. We provide Kronecker-type theorems for convolution operators, showing that a convolution operator (or a cross-correlation operator) is of finite rank, if and only if, its symbol is a polynomial-exponential function, and we relate its rank to the decomposition of its symbol. We also present Kronecker-type theorems for the reconstruction of measures as weighted sums of Dirac measures from moments and for the decomposition of polynomial-exponential functions from values. Finally, we describe an application of this method for the sparse interpolation of polylog functions from values.

How to compute the constant term of a power of a Laurent polynomial efficiently

Dwork's congruences for the constant terms of powers of a Laurent polynomial

On a power series involving classical orthogonal polynomials

Polynomial-exponential decomposition from moments

Improved Laguerre Spectral Methods with Less Round-off Errors and Better Stability

GPU Accelerated Newton for Taylor Series Solutions of Polynomial Homotopies in Multiple Double Precision

On the Complexity and Parallel Implementation of Hensel's Lemma and Weierstrass Preparation

Symbolic computation of hypergeometric type and non-holonomic power series

Polynomial formulations of Multivariable Arithmetic Progressions of Alternating Powers

Effective homology and periods of complex projective hypersurfaces

A numeric study of power expansions around singular points of algebraic functions, their radii of convergence, and accuracy profiles

Approximation by a new sequence of operators involving Laguerre polynomials

Polynomial approximation on disjoint segments and amplification of approximation

Fast Approximate Polynomial Multipoint Evaluation and Applications

An Accurate Numerical Algorithm Based on the Generalized Narayana Polynomials to Solve a Class of Caputo-Fabrizio and Liouville-Caputo Fractional-Order Delay Differential Equations

Fast and backward stable computation of the eigenvalues and eigenvectors of matrix polynomials

Fast interpolation of sparse multivariate polynomials

Laurent expansions of meromorphic modular forms

On the Representation of Non-Holonomic Power Series

On the Hadamard and Vandermonde Determinants and the Bernoulli–Euler–Lagrange–Aitken Method for Calculating the Roots of Polynomials