6,575 research outputs found
A simple and fast algorithm for computing exponentials of power series
International audienceAs was initially shown by Brent, exponentials of truncated power series can be computed using a constant number of polynomial multiplications. This note gives a relatively simple algorithm with a low constant factor
Fast computation of power series solutions of systems of differential equations
We propose new algorithms for the computation of the first N terms of a
vector (resp. a basis) of power series solutions of a linear system of
differential equations at an ordinary point, using a number of arithmetic
operations which is quasi-linear with respect to N. Similar results are also
given in the non-linear case. This extends previous results obtained by Brent
and Kung for scalar differential equations of order one and two
Software for Exact Integration of Polynomials over Polyhedra
We are interested in the fast computation of the exact value of integrals of
polynomial functions over convex polyhedra. We present speed ups and extensions
of the algorithms presented in previous work. We present the new software
implementation and provide benchmark computations. The computation of integrals
of polynomials over polyhedral regions has many applications; here we
demonstrate our algorithmic tools solving a challenge from combinatorial voting
theory.Comment: Major updat
Efficient sum-of-exponentials approximations for the heat kernel and their applications
In this paper, we show that efficient separated sum-of-exponentials
approximations can be constructed for the heat kernel in any dimension. In one
space dimension, the heat kernel admits an approximation involving a number of
terms that is of the order for any x\in\bbR and
, where is the desired precision. In all
higher dimensions, the corresponding heat kernel admits an approximation
involving only terms for fixed accuracy
. These approximations can be used to accelerate integral
equation-based methods for boundary value problems governed by the heat
equation in complex geometry. The resulting algorithms are nearly optimal. For
points in the spatial discretization and time steps, the cost is
in terms of both memory and CPU time for
fixed accuracy . The algorithms can be parallelized in a
straightforward manner. Several numerical examples are presented to illustrate
the accuracy and stability of these approximations.Comment: 23 pages, 5 figures, 3 table
Computing the Exponential of Large Block-Triangular Block-Toeplitz Matrices Encountered in Fluid Queues
The Erlangian approximation of Markovian fluid queues leads to the problem of
computing the matrix exponential of a subgenerator having a block-triangular,
block-Toeplitz structure. To this end, we propose some algorithms which exploit
the Toeplitz structure and the properties of generators. Such algorithms allow
to compute the exponential of very large matrices, which would otherwise be
untreatable with standard methods. We also prove interesting decay properties
of the exponential of a generator having a block-triangular, block-Toeplitz
structure
Efficient implementation of the Hardy-Ramanujan-Rademacher formula
We describe how the Hardy-Ramanujan-Rademacher formula can be implemented to
allow the partition function to be computed with softly optimal
complexity and very little overhead. A new implementation
based on these techniques achieves speedups in excess of a factor 500 over
previously published software and has been used by the author to calculate
, an exponent twice as large as in previously reported
computations.
We also investigate performance for multi-evaluation of , where our
implementation of the Hardy-Ramanujan-Rademacher formula becomes superior to
power series methods on far denser sets of indices than previous
implementations. As an application, we determine over 22 billion new
congruences for the partition function, extending Weaver's tabulation of 76,065
congruences.Comment: updated version containing an unconditional complexity proof;
accepted for publication in LMS Journal of Computation and Mathematic
- …