Efficient computation with structured matrices and arithmetic expressions

Christophe Mouilleron 1
1 ARENAIRE - Computer arithmetic
Inria Grenoble - Rhône-Alpes, LIP - Laboratoire de l'Informatique du Parallélisme
Abstract : Designing efficient code in practice for a given computation is a hard task. In this thesis, we tackle this issue in two different situations. The first part of the thesis introduces some algorithmic improvements in structured linear algebra. We first show how to extend an algorithm by Cardinal for inverting Cauchy-like matrices to the other common structures. This approach, which mainly relies on products of the type "structured matrix × matrix", leads to a theoretical speed-up of a factor up to 7 that we also observe in practice. Then, we extend some works on Toeplitz-like matrices and prove that, for any of the common structures, the product of an n×n structured matrix of displacement rank α by an n×α matrix can be computed in Õ(α^(ω-1)n). This leads to direct inversion algorithms in Õ(α^(ω-1)n) , that do not rely on a reduction to the Toeplitz-like case. The second part of the thesis deals with the implementation of arithmetic expressions. This topic raises several issues like finding the minimum number of operations, and maximizing the speed or the accuracy when using some finite-precision arithmetic. Making use of the inductive nature of arithmetic expressions enables the design of algorithms that help to answer such questions. We thus present a set of algorithms for generating evaluation schemes, counting them, and optimizing them according to one or several criteria. These algorithms are part of a library that we have developed and used, among other things, in order to decrease the running time of a code generator for a mathematical library, and to study optimality issues about the evaluation of a small degree polynomial with scalar coefficients at a matrix point.
Document type :
Theses
Other [cs.OH]. Ecole normale supérieure de lyon - ENS LYON, 2011. English. < NNT : 2011ENSL0652 >
Liste complète des métadonnées


https://tel.archives-ouvertes.fr/tel-00688388
Contributor : Abes Star <>
Submitted on : Tuesday, April 17, 2012 - 3:17:35 PM
Last modification on : Friday, June 16, 2017 - 11:28:28 AM
Document(s) archivé(s) le : Wednesday, July 18, 2012 - 2:31:54 AM

File

MOUILLERON_Christophe_2011_The...
Version validated by the jury (STAR)

Identifiers

  • HAL Id : tel-00688388, version 1

Collections

Citation

Christophe Mouilleron. Efficient computation with structured matrices and arithmetic expressions. Other [cs.OH]. Ecole normale supérieure de lyon - ENS LYON, 2011. English. < NNT : 2011ENSL0652 >. <tel-00688388>

Share

Metrics

Record views

642

Document downloads

936