[Accueil bibliotech]
Accueil > Les thèses en ligne de l'INP

A study on block flexible iterative solvers with applications to Earth imaging problem in geophysics

Ferreira Lago, Rafael (2013) A study on block flexible iterative solvers with applications to Earth imaging problem in geophysics. (Étude de méthodes itératives par bloc avec application à l’imagerie sismique en géophysique.)

Texte intégral disponible au format :

PDF - Nécessite un logiciel de visualisation PDF comme GSview, Xpdf ou Adobe Acrobat Reader
6.06 Mo

Prix Léopold Escande 2013 (en savoir plus)

Résumé

Les travaux de ce doctorat concernent le développement de méthodes itératives pour la résolution de systèmes linéaires creux de grande taille comportant de nombreux seconds membres. L’application visée est la résolution d’un problème inverse en géophysique visant à reconstruire la vitesse de propagation des ondes dans le sous-sol terrestre. Lorsque de nombreuses sources émettrices sont utilisées, ce problème inverse nécessite la résolution de systèmes linéaires complexes non symétriques non hermitiens comportant des milliers de seconds membres. Dans le cas tridimensionnel ces systèmes linéaires sont reconnus comme difficiles à résoudre plus particulièrement lorsque des fréquences élevées sont considérées. Le principal objectif de cette thèse est donc d’étendre les développements existants concernant les méthodes de Krylov par bloc. Nous étudions plus particulièrement les techniques de déflation dans le cas multiples seconds membres et recyclage de sous-espace dans le cas simple second membre. Des gains substantiels sont obtenus en terme de temps de calcul par rapport aux méthodes existantes sur des applications réalistes dans un environnement parallèle distribué. ABSTRACT : This PhD thesis concerns the development of flexible Krylov subspace iterative solvers for the solution of large sparse linear systems of equations with multiple right-hand sides. Our target application is the solution of the acoustic full waveform inversion problem in geophysics associated with the phenomena of wave propagation through an heterogeneous model simulating the subsurface of Earth. When multiple wave sources are being used, this problem gives raise to large sparse complex non-Hermitian and nonsymmetric linear systems with thousands of right-hand sides. Specially in the three-dimensional case and at high frequencies, this problem is known to be difficult. The purpose of this thesis is to develop a flexible block Krylov iterative method which extends and improves techniques already available in the current literature to the multiple right-hand sides scenario. We exploit the relations between each right-hand side to accelerate the convergence of the overall iterative method. We study both block deflation and single right-hand side subspace recycling techniques obtaining substantial gains in terms of computational time when compared to other strategies published in the literature, on realistic applications performed in a parallel environment.

Département ou laboratoire:Centre Européen de Recherche et Formation Avancées en Calcul Scientifique - CERFACS (Toulouse, France)
Directeur de thèse:Gratton, Serge et Vasseur, Xavier
Mots-clés:Sous-espaces de Krylov - Méthodes itératives - Calcul de haute performance - Equation de Helmholtz - Imagerie sismique. KEYWORDS : Krylov subspace methods - Iterative methods - High performance computing - Helmholtz equation - Earth imaging
Sujets:Informatique > Programmation et systèmes
Mathématiques appliquées
Déposé le:07 Novembre 2013

Administrateur seulement : modifier cet enregistrement


Contacts | Infos légales | Plan du site | Intranet

(c)INP de Toulouse 2012 - Tous droits réservés. -  INP Communication