Semi-Automatic Floating-Point Implementation of Special Functions

Christoph Lauter 1, * Marc Mezzarobba 1, *
* Auteur correspondant
1 PEQUAN - Performance et Qualité des Algorithmes Numériques
LIP6 - Laboratoire d'Informatique de Paris 6
Abstract : This work introduces an approach to the computer-assisted implementation of mathematical functions geared toward special functions such as those occurring in mathematical physics. The general idea is to start with an exact symbolic representation of a function and automate as much as possible of the process of implementing it. In order to deal with a large class of special functions, our symbolic representation is an implicit one: the input is a linear differential equation with polynomial coefficients along with initial values. The output is a C program to evaluate the solution of the equation using domain splitting, argument reduction and polynomial approximations in double-precision arithmetic, in the usual style of mathematical libraries. Our generation method combines symbolic-numeric manipulations of linear ODEs with interval-based tools for the floating-point implementation of "black-box" functions. We describe a prototype code generator that can automatically produce implementations on moderately large intervals. Implementations on the whole real line are possible in some cases but require manual tool setup and code integration. Due to this limitation and as some heuristics remain, we refer to our method as "semi-automatic" at this stage. Along with other examples, we present an implementation of the Voigt profile with fixed parameters that may be of independent interest.
Type de document :
Communication dans un congrès
IEEE 22nd Symposium on Computer Arithmetic, ARITH 22, Jun 2015, Lyon, France. IEEE, pp.58-65, 2015, 〈10.1109/ARITH.2015.12〉
Liste complète des métadonnées

Littérature citée [28 références]  Voir  Masquer  Télécharger

https://hal.archives-ouvertes.fr/hal-01137953
Contributeur : Marc Mezzarobba <>
Soumis le : jeudi 2 avril 2015 - 09:21:12
Dernière modification le : jeudi 20 juillet 2017 - 09:28:21
Document(s) archivé(s) le : vendredi 3 juillet 2015 - 10:07:14

Fichiers

frankenstein.pdf
Fichiers produits par l'(les) auteur(s)

Licence


Distributed under a Creative Commons Paternité 4.0 International License

Identifiants

Collections

Citation

Christoph Lauter, Marc Mezzarobba. Semi-Automatic Floating-Point Implementation of Special Functions. IEEE 22nd Symposium on Computer Arithmetic, ARITH 22, Jun 2015, Lyon, France. IEEE, pp.58-65, 2015, 〈10.1109/ARITH.2015.12〉. 〈hal-01137953〉

Partager

Métriques

Consultations de la notice

216

Téléchargements de fichiers

191