Accès direct au contenu

 

logo du site ENS

Logo du site Ens Lyon

Recherche

Accueil  >  ENS de Lyon > Recherche > Soutenances de thèses > Thèses en sc. exactes et exp.

Imprimer

Contribution à l’analyse d’algorithmes en arithmétique à virgule flottante / Contribution to error analysis of algorithms in floating-point arithmetic

Soutenance de thèse
Soutenance de thèse de M. Antoine PLET du LIP sous la direction de M. Jean-Michel MULLER
7 juillet 2017
14h00

Lieu(x) :

Site Monod - 46 allée d'Italie
Amphi A - Site Monod - ENS de Lyon

L’arithmétique virgule flottante est une approximation de l’arithmétique réelle dans laquelle chaque opération peut introduire une erreur. La norme IEEE 754 requiert que les opérations élémentaires soient aussi précises que possible, mais au cours d’un calcul, les erreurs d’arrondi s’accumulent et peuvent conduire à des résultats totalement faussés. Cela arrive avec une expression aussi simple que ab + cd, pour laquelle l’algorithme naïf retourne parfois un résultat aberrant, avec une erreur relative largement supérieure à 1. Il est donc important d’analyser les algorithmes utilisés pour contrôler l’erreur commise. Je m’intéresse à l’analyse de briques élémentaires du calcul en cherchant des bornes fines sur l’erreur relative.
Pour des algorithmes suffisamment précis, en arithmétique de base β et de précision p, on arrive en général à prouver une borne sur l'erreur de la forme α·u + o(u²) où α > 0 et u = 1/2·β1-p est l'unité d'arrondi. Comme indication de la finesse d'une telle borne, on peut fournir des exemples numériques pour les précisions standards qui approchent cette borne, ou bien un exemple paramétré par la précision qui génère une erreur de la forme α·u + o(u²), prouvant ainsi l'optimalité asymptotique de la borne. J’ai travaillé sur la formalisation d’une arithmétique à virgule flottante symbolique, sur des nombres paramétrés par la précision, et à son implantation dans le logiciel de calcul formel Maple.
J’ai aussi obtenu une borne d'erreur très fine pour un algorithme d’inversion complexe en arithmétique flottante. Ce résultat suggère le calcul d'une division décrit par la formule x/y = (1/y)·x, par opposition à x/y = (x·y)/|y|². Quel que soit l'algorithme utilisé pour effectuer la multiplication, nous avons une borne d'erreur plus petite pour les algorithmes décrits par la première formule.
Ces travaux sont réalisés avec mes directeurs de thèse, en collaboration avec Claude-Pierre
Jeannerod (CR Inria dans AriC, au LIP).

Retour au haut de la page

Pôle Recherche
 
 
Mise à jour le 23 juin 2017
Ens de Lyon
15 parvis René Descartes - BP 7000 69342 Lyon Cedex 07 - FRANCE
Tél. : Site René Descartes (siège) : +33 (0) 4 37 37 60 00 / Site Jacques Monod : +33 (0) 4 72 72 80 00