Efficient and accurate computation of upper bounds of approximation errors - ENS de Lyon - École normale supérieure de Lyon Accéder directement au contenu
Pré-Publication, Document De Travail Année : 2010

Efficient and accurate computation of upper bounds of approximation errors

John Harrison
  • Fonction : Auteur
  • PersonId : 865917
Mioara Maria Joldes
Christoph Lauter

Résumé

For purposes of actual evaluation, mathematical functions f are commonly replaced by approximation polynomials p. Examples include floating-point implementations of elementary functions, quadrature or more theoretical proof work involving transcendental functions. Replacing f by p induces a relative error epsilon =p/f - 1. In order to ensure the validity of the use of p instead of f, the maximum error, i.e. the supremum norm of epsilon must be safely bounded above. Numerical algorithms for supremum norms are efficient but cannot offer the required safety. Previous validated approaches often require tedious manual intervention. If they are automated, they have several drawbacks, such as the lack of quality guarantees. In this article a novel, automated supremum norm algorithm with a priori quality is proposed. It focuses on the validation step and paves the way for formally certified supremum norms. Key elements are the use of intermediate approximation polynomials with bounded truncation remainder and a non-negativity test based on a Sum-of-squares expression of polynomials. The new algorithm was implemented in the tool Sollya. The article includes experimental results on real-life examples.
Fichier principal
Vignette du fichier
RRLIP2010-2.pdf (1.03 Mo) Télécharger le fichier
Origine : Fichiers produits par l'(les) auteur(s)

Dates et versions

ensl-00445343 , version 1 (08-01-2010)
ensl-00445343 , version 2 (14-07-2010)

Identifiants

  • HAL Id : ensl-00445343 , version 1

Citer

Sylvain Chevillard, John Harrison, Mioara Maria Joldes, Christoph Lauter. Efficient and accurate computation of upper bounds of approximation errors. 2010. ⟨ensl-00445343v1⟩
614 Consultations
1353 Téléchargements

Partager

Gmail Facebook X LinkedIn More