Parallel Implementation of Interval Matrix Multiplication - Inria - Institut national de recherche en sciences et technologies du numérique Accéder directement au contenu
Article Dans Une Revue Reliable Computing Année : 2013

Parallel Implementation of Interval Matrix Multiplication

Résumé

Two main and not necessarily compatible objectives when implementing the product of two dense matrices with interval coefficients are accuracy and efficiency. In this work, we focus on an implementation on multicore architectures. One direction successfully explored to gain performance in execution time is the representation of intervals by their midpoints and radii rather than the classical representation by endpoints. Computing with the midpoint-radius representation enables the use of optimized floating-point BLAS and consequently the performances benefit from the performances of the BLAS routines. Several variants of interval matrix multiplication have been proposed, that correspond to various trade-offs between accuracy and efficiency, including some efficient ones proposed by Rump in 2012. However, in order to guarantee that the computed result encloses the exact one, these efficient algorithms rely on an assumption on the order of execution of floating-point operations which is not verified by most implementations of BLAS. In this paper, an algorithm for interval matrix product is proposed that verifies this assumption. Furthermore, several optimizations are proposed and the implementation on a multicore architecture compares reasonably well with a non-guaranteed implementation based on MKL, the optimized BLAS of Intel: the overhead is most of the time less than 2 and never exceeds 3. This implementation also exhibits a good scalability.
Fichier principal
Vignette du fichier
parallel-igemm.pdf (440.38 Ko) Télécharger le fichier
Origine : Fichiers produits par l'(les) auteur(s)
Loading...

Dates et versions

hal-00801890 , version 1 (18-03-2013)
hal-00801890 , version 2 (11-12-2013)

Identifiants

  • HAL Id : hal-00801890 , version 2

Citer

Nathalie Revol, Philippe Théveny. Parallel Implementation of Interval Matrix Multiplication. Reliable Computing, 2013, 19 (1), pp.91-106. ⟨hal-00801890v2⟩
293 Consultations
458 Téléchargements

Partager

Gmail Facebook X LinkedIn More