Précédent : Inversion sismique Remonter : Résultats nouveaux Suivant : Inversion électromagnétique
Précédent : Inversion sismique
Remonter : Résultats nouveaux Suivant :
Inversion
électromagnétique
Participants : Michel Kern, Guillaume Dumont.
Nous avons repris le code de migration sismique 3D par équations paraxiales écrit par M. Kern pour des machines vectorielles afin de le porter sur les ordinateurs à mémoire distribuée (Cray T3E, SGI Origin 2000). Nous avons voulu tirer profit de cette opération pour évaluer la bibliothèque PETSc du Argonne National Laboratory. PETSc fournit des routines et des structures de données permettant la résolution d'équations aux dérivées partielles sur calculateurs parallèles. En pratique, l'utilisateur prend en charge tout ce qui concerne la physique particulière de son problème, et PETSc fournit les solveurs algébriques.
Le code écrit par Guillaume Dumonttourne actuellement sur l'Origin 2000 du Centre Charles Hermite à Nancy. PETSc a permis une mise au point rapide du code, avec une performance raisonnable. Toutefois, certains systèmes (correspondant aux basses fréquences) se sont avérés très difficiles à résoudre, et l'absence de préconditionneurs robustes en parallèle s'est rapidement fait sentir. Cela montre à la fois les possibilités et les limites de ce type d'outils.