Validated and numerically efficient Chebyshev spectral methods for linear ordinary differential equations

Florent Bréhard 1, 2, 3 Nicolas Brisebarre 4, 2, 3 Mioara Joldes 1, 4
1 LAAS-MAC - Équipe Méthodes et Algorithmes en Commande
LAAS - Laboratoire d'analyse et d'architecture des systèmes [Toulouse]
3 ARIC - Arithmetic and Computing
Inria Grenoble - Rhône-Alpes, LIP - Laboratoire de l'Informatique du Parallélisme
Abstract : In this work we develop a validated numerics method for the solution of linear ordinary differential equations (LODEs). A wide range of algorithms (i.e., Runge-Kutta, collocation, spectral methods) exist for numerically computing approximations of the solutions. Most of these come with proofs of asymptotic convergence, but usually, provided error bounds are non-constructive. However, in some domains like critical systems and computer-aided mathematical proofs, one needs validated effective error bounds. We focus on both the theoretical and practical complexity analysis of a so-called \emph{a posteriori} quasi-Newton validation method, which mainly relies on a fixed-point argument of a contracting map. Specifically, given a polynomial approximation, obtained by some numerical algorithm and expressed in Chebyshev basis, our algorithm efficiently computes an accurate and rigorous error bound. For this, we study theoretical properties like compactness, convergence, invertibility of associated linear integral operators and their truncations in a suitable coefficient space of Chebyshev series. Then, we analyze the almost-banded matrix structure of these operators, which allows for very efficient numerical algorithms for both numerical solutions of LODEs and rigorous computation of the approximation error. Finally, several representative examples show the advantages of our algorithms as well as their theoretical and practical limits.
Type de document :
Pré-publication, Document de travail
Rapport LAAS n° 17177. 2017
Liste complète des métadonnées

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

https://hal.archives-ouvertes.fr/hal-01526272
Contributeur : Mioara Joldes <>
Soumis le : vendredi 28 juillet 2017 - 16:43:46
Dernière modification le : vendredi 20 avril 2018 - 15:44:26

Fichier

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

Identifiants

  • HAL Id : hal-01526272, version 2

Citation

Florent Bréhard, Nicolas Brisebarre, Mioara Joldes. Validated and numerically efficient Chebyshev spectral methods for linear ordinary differential equations. Rapport LAAS n° 17177. 2017. 〈hal-01526272v2〉

Partager

Métriques

Consultations de la notice

466

Téléchargements de fichiers

174