Title Benchmarking and extension of fast linear algebra algorithms. Version française

Location

INRIA, Project CAFÉ
BP 93, 06902 France

Information

Manuel Bronstein

Description

The Café project has developped or implemented recent asymptotically fast algorithms for exact linear algebra [1,2,3]. For example, there are now 6 different algorithms for computing the determinants of matrices of polynomials or the exact solutions of linear systems. Our first experimental results indicate that choosing the appropriate algorithm to use on a given matrix depends on several parameters, including the size of the matrix but also the type of the coefficients of the polynomial entries (integers, rational numbers, elements of algebraic numbers fields or finite fields), and probably also the distribution of the degrees inside the matrix. We therefore need a careful and exhaustive benchmark in order to determine a strategy for automatically selecting the most appropriate algorithm in each case. This study should also help understand the practical behavior of those algorithms and eventually to refine and improve them. Finally, we are also interested in extending those algorithms to compute related quantities in linear algebra such as resultants or cyclic vectors and/or to implement some of them on a distributed platform.
[1] J.Abbott, M.Bronstein & T.Mulders,  Fast Deterministic Computation of Determinants of Dense Matrices , Proceedings of ISSAC'99, ACM Press, 197-204 (1999).
[2] T.Mulders & A.Storjohann,  Rational Solutions of Singular Linear Systems , Proceedings of ISSAC'2000, ACM Press, 242-249 (2000).
[3] T.Mulders & A.Storjohann,  Triangular Factorization of Polynomial Matrices,  ISSAC'2000 Poster.
 

Goals

A first objective is to develop a large and varied collection of test examples in order to gain enough data to develop a strategy for selecting the most appropriate algorithm for a given matrix. Further objectives, as time permits, will be to extend selected algorithms to compute resultants and cyclic vectors (linear dependencies) as well as to write a distributed implementation on top of the high--level Pi^it communication library.

Tools

Unix workstation, language Aldor, libraries  and Pi^it.

Duration

3 - 4 months, possible prolongation as a doctoral thesis if desired.


Manuel.Bronstein@sophia.inria.fr

November 16, 2000