22 research outputs found
Evaluating parametric holonomic sequences using rectangular splitting
We adapt the rectangular splitting technique of Paterson and Stockmeyer to
the problem of evaluating terms in holonomic sequences that depend on a
parameter. This approach allows computing the -th term in a recurrent
sequence of suitable type using "expensive" operations at the cost
of an increased number of "cheap" operations.
Rectangular splitting has little overhead and can perform better than either
naive evaluation or asymptotically faster algorithms for ranges of
encountered in applications. As an example, fast numerical evaluation of the
gamma function is investigated. Our work generalizes two previous algorithms of
Smith.Comment: 8 pages, 2 figure
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
A Comparative Study of Two Real Root Isolation Methods
Recent progress in polynomial elimination has rendered the computation of the real roots of ill-conditioned polynomials of high degree (over 1000) with huge coefficients (several thousand digits) a critical operation in computer algebra.
To rise to the occasion, the only method-candidate that has been considered by various authors for modification and improvement has been the Collins-Akritas bisection method [1], which is a based on a variation of Vincentâs theorem [2].
The most recent example is the paper by Rouillier and Zimmermann [3], where the authors present â... a new algorithm, which is optimal in terms of memory usage and as fast as both Collins and Akritasâ algorithm and Krandick variant ...â [3]
In this paper we compare our own continued fractions method CF [4] (which is directly based on Vincentâs theorem) with the best bisection method REL described in [3]. Experimentation with the data presented in [3] showed that, with respect to time, our continued fractions method CF is by far superior to REL, whereas the two are about equal with respect to space
Parallel Integer Polynomial Multiplication
We propose a new algorithm for multiplying dense polynomials with integer
coefficients in a parallel fashion, targeting multi-core processor
architectures. Complexity estimates and experimental comparisons demonstrate
the advantages of this new approach
A Fast Algorithm for Computing the p-Curvature
We design an algorithm for computing the -curvature of a differential
system in positive characteristic . For a system of dimension with
coefficients of degree at most , its complexity is \softO (p d r^\omega)
operations in the ground field (where denotes the exponent of matrix
multiplication), whereas the size of the output is about . Our
algorithm is then quasi-optimal assuming that matrix multiplication is
(\emph{i.e.} ). The main theoretical input we are using is the
existence of a well-suited ring of series with divided powers for which an
analogue of the Cauchy--Lipschitz Theorem holds.Comment: ISSAC 2015, Jul 2015, Bath, United Kingdo
On the Complexity of Real Root Isolation
We introduce a new approach to isolate the real roots of a square-free
polynomial with real coefficients. It is assumed that
each coefficient of can be approximated to any specified error bound. The
presented method is exact, complete and deterministic. Due to its similarities
to the Descartes method, we also consider it practical and easy to implement.
Compared to previous approaches, our new method achieves a significantly better
bit complexity. It is further shown that the hardness of isolating the real
roots of is exclusively determined by the geometry of the roots and not by
the complexity or the size of the coefficients. For the special case where
has integer coefficients of maximal bitsize , our bound on the bit
complexity writes as which improves the best bounds
known for existing practical algorithms by a factor of . The crucial
idea underlying the new approach is to run an approximate version of the
Descartes method, where, in each subdivision step, we only consider
approximations of the intermediate results to a certain precision. We give an
upper bound on the maximal precision that is needed for isolating the roots of
. For integer polynomials, this bound is by a factor lower than that of
the precision needed when using exact arithmetic explaining the improved bound
on the bit complexity
Products of Ordinary Differential Operators by Evaluation and Interpolation
It is known that multiplication of linear differential operators over ground
fields of characteristic zero can be reduced to a constant number of matrix
products. We give a new algorithm by evaluation and interpolation which is
faster than the previously-known one by a constant factor, and prove that in
characteristic zero, multiplication of differential operators and of matrices
are computationally equivalent problems. In positive characteristic, we show
that differential operators can be multiplied in nearly optimal time.
Theoretical results are validated by intensive experiments
Continued Fraction Expansion of Real Roots of Polynomial Systems
We present a new algorithm for isolating the real roots of a system of
multivariate polynomials, given in the monomial basis. It is inspired by
existing subdivision methods in the Bernstein basis; it can be seen as
generalization of the univariate continued fraction algorithm or alternatively
as a fully analog of Bernstein subdivision in the monomial basis. The
representation of the subdivided domains is done through homographies, which
allows us to use only integer arithmetic and to treat efficiently unbounded
regions. We use univariate bounding functions, projection and preconditionning
techniques to reduce the domain of search. The resulting boxes have optimized
rational coordinates, corresponding to the first terms of the continued
fraction expansion of the real roots. An extension of Vincent's theorem to
multivariate polynomials is proved and used for the termination of the
algorithm. New complexity bounds are provided for a simplified version of the
algorithm. Examples computed with a preliminary C++ implementation illustrate
the approach.Comment: 10 page
Fast algorithms for differential equations in positive characteristic
We address complexity issues for linear differential equations in
characteristic : resolution and computation of the -curvature. For
these tasks, our main focus is on algorithms whose complexity behaves well with
respect to . We prove bounds linear in on the degree of polynomial
solutions and propose algorithms for testing the existence of polynomial
solutions in sublinear time , and for determining a whole
basis of the solution space in quasi-linear time ; the
notation indicates that we hide logarithmic factors. We show that
for equations of arbitrary order, the -curvature can be computed in
subquadratic time , and that this can be improved to
for first order equations and to for classes of
second order equations