16,139 research outputs found
New recurrence relationships between orthogonal polynomials which lead to new Lanczos-type algorithms
Lanczos methods for solving Ax = b consist in constructing a sequence of vectors (Xk),k = 1,... such that rk = b-AXk= Pk(A)r0, where Pk is the orthogonal polynomial of degree at most k with respect to the linear functional c defined as c(εi) = (y, Air0). Let P(1)k be the regular monic polynomial of degree k belonging to the family of formal orthogonal polynomials (FOP) with respect to c(1) defined as c(1)(εi) = c(εi+1). All Lanczos-type algorithms are characterized by the choice of one or two recurrence relationships, one for Pk and one for P(1)k. We shall study some new recurrence relations involving these two polynomials and their possible combinations to obtain new Lanczos-type algorithms. We will show that some recurrence relations exist, but cannot be used to derive Lanczos-type algorithms, while others do not exist at all
Transition probabilities for general birth-death processes with applications in ecology, genetics, and evolution
A birth-death process is a continuous-time Markov chain that counts the
number of particles in a system over time. In the general process with
current particles, a new particle is born with instantaneous rate
and a particle dies with instantaneous rate . Currently no robust and
efficient method exists to evaluate the finite-time transition probabilities in
a general birth-death process with arbitrary birth and death rates. In this
paper, we first revisit the theory of continued fractions to obtain expressions
for the Laplace transforms of these transition probabilities and make explicit
an important derivation connecting transition probabilities and continued
fractions. We then develop an efficient algorithm for computing these
probabilities that analyzes the error associated with approximations in the
method. We demonstrate that this error-controlled method agrees with known
solutions and outperforms previous approaches to computing these probabilities.
Finally, we apply our novel method to several important problems in ecology,
evolution, and genetics
Fast algorithm for border bases of Artinian Gorenstein algebras
Given a multi-index sequence , we present a new efficient algorithm
to compute generators of the linear recurrence relations between the terms of
. We transform this problem into an algebraic one, by identifying
multi-index sequences, multivariate formal power series and linear functionals
on the ring of multivariate polynomials. In this setting, the recurrence
relations are the elements of the kerne l\sigma of the Hankel operator
$H$\sigma associated to . We describe the correspondence between
multi-index sequences with a Hankel operator of finite rank and Artinian
Gorenstein Algebras. We show how the algebraic structure of the Artinian
Gorenstein algebra \sigma\sigma yields the
structure of the terms $\sigma\alpha N nAK[x 1 ,. .. , xnIHIA$ and the tables of multiplication by the variables in these
bases. It is an extension of Berlekamp-Massey-Sakata (BMS) algorithm, with
improved complexity bounds. We present applications of the method to different
problems such as the decomposition of functions into weighted sums of
exponential functions, sparse interpolation, fast decoding of algebraic codes,
computing the vanishing ideal of points, and tensor decomposition. Some
benchmarks illustrate the practical behavior of the algorithm
Polynomial-Time Algorithms for Quadratic Isomorphism of Polynomials: The Regular Case
Let and be
two sets of nonlinear polynomials over
( being a field). We consider the computational problem of finding
-- if any -- an invertible transformation on the variables mapping
to . The corresponding equivalence problem is known as {\tt
Isomorphism of Polynomials with one Secret} ({\tt IP1S}) and is a fundamental
problem in multivariate cryptography. The main result is a randomized
polynomial-time algorithm for solving {\tt IP1S} for quadratic instances, a
particular case of importance in cryptography and somewhat justifying {\it a
posteriori} the fact that {\it Graph Isomorphism} reduces to only cubic
instances of {\tt IP1S} (Agrawal and Saxena). To this end, we show that {\tt
IP1S} for quadratic polynomials can be reduced to a variant of the classical
module isomorphism problem in representation theory, which involves to test the
orthogonal simultaneous conjugacy of symmetric matrices. We show that we can
essentially {\it linearize} the problem by reducing quadratic-{\tt IP1S} to
test the orthogonal simultaneous similarity of symmetric matrices; this latter
problem was shown by Chistov, Ivanyos and Karpinski to be equivalent to finding
an invertible matrix in the linear space of matrices over and to compute the square root in a matrix
algebra. While computing square roots of matrices can be done efficiently using
numerical methods, it seems difficult to control the bit complexity of such
methods. However, we present exact and polynomial-time algorithms for computing
the square root in for various fields (including
finite fields). We then consider \\#{\tt IP1S}, the counting version of {\tt
IP1S} for quadratic instances. In particular, we provide a (complete)
characterization of the automorphism group of homogeneous quadratic
polynomials. Finally, we also consider the more general {\it Isomorphism of
Polynomials} ({\tt IP}) problem where we allow an invertible linear
transformation on the variables \emph{and} on the set of polynomials. A
randomized polynomial-time algorithm for solving {\tt IP} when
is presented. From an algorithmic point
of view, the problem boils down to factoring the determinant of a linear matrix
(\emph{i.e.}\ a matrix whose components are linear polynomials). This extends
to {\tt IP} a result of Kayal obtained for {\tt PolyProj}.Comment: Published in Journal of Complexity, Elsevier, 2015, pp.3
Schubert polynomials and Arakelov theory of symplectic flag varieties
Let X be the flag variety of the symplectic group. We propose a theory of
combinatorially explicit Schubert polynomials which represent the Schubert
classes in the Borel presentation of the cohomology ring of X. We use these
polynomials to describe the arithmetic Schubert calculus on X. Moreover, we
give a method to compute the natural arithmetic Chern numbers on X, and show
that they are all rational numbers.Comment: 22 pages; final versio
- …