Horizon / Plein textes La base de ressources documentaires de l'IRD

IRD

Publications des scientifiques de l'IRD

Martin R., Monteiller V., Komatitsch D., Perrouty S., Jessell Mark, Bonvalot Sylvain, Lindsay M. (2013). Gravity inversion using wavelet-based compression on parallel hybrid CPU/GPU systems : application to southwest Ghana. Geophysical Journal International, 195 (3), 1594-1619. ISSN 0956-540X

Lien direct chez l'éditeur doi:10.1093/gji/ggt334

Titre
Gravity inversion using wavelet-based compression on parallel hybrid CPU/GPU systems : application to southwest Ghana
Année de publication2013
Type de documentArticle référencé dans le Web of Science WOS:000326964400014
AuteursMartin R., Monteiller V., Komatitsch D., Perrouty S., Jessell Mark, Bonvalot Sylvain, Lindsay M.
SourceGeophysical Journal International, 2013, 195 (3), p. 1594-1619. ISSN 0956-540X
RésuméWe solve the 3-D gravity inverse problem using a massively parallel voxel (or finite element) implementation on a hybrid multi-CPU/multi-GPU (graphics processing units/GPUs) cluster. This allows us to obtain information on density distributions in heterogeneous media with an efficient computational time. In a new software package called TOMOFAST3D, the inversion is solved with an iterative least-square or a gradient technique, which minimizes a hybrid L-1-/L-2-norm-based misfit function. It is drastically accelerated using either Haar or fourth-order Daubechies wavelet compression operators, which are applied to the sensitivity matrix kernels involved in the misfit minimization. The compression process behaves like a pre-conditioning of the huge linear system to be solved and a reduction of two or three orders of magnitude of the computational time can be obtained for a given number of CPU processor cores. The memory storage required is also significantly reduced by a similar factor. Finally, we show how this CPU parallel inversion code can be accelerated further by a factor between 3.5 and 10 using GPU computing. Performance levels are given for an application to Ghana, and physical information obtained after 3-D inversion using a sensitivity matrix with around 5.37 trillion elements is discussed. Using compression the whole inversion process can last from a few minutes to less than an hour for a given number of processor cores instead of tens of hours for a similar number of processor cores when compression is not used.
Plan de classementGéophysique interne [066] ; Sciences fondamentales / Techniques d'analyse et de recherche [020]
Descr. géo.GHANA
LocalisationFonds IRD
Identifiant IRDPAR00011166
Lien permanenthttp://www.documentation.ird.fr/hor/PAR00011166

Export des données

Disponibilité des documents

Télechargment fichier PDF téléchargeable

Lien sur le Web lien chez l'éditeur

Accès réservé en accès réservé

HAL en libre accès sur HAL


Accès aux documents originaux :

Le FDI est labellisé CollEx

Accès direct

Bureau du chercheur

Site de la documentation

Espace intranet IST (accès réservé)

Suivi des publications IRD (accès réservé)

Mentions légales

Services Horizon

Poser une question

Consulter l'aide en ligne

Déposer une publication (accès réservé)

S'abonner au flux RSS

Voir les tableaux chronologiques et thématiques

Centres de documentation

Bondy

Montpellier (centre IRD)

Montpellier (MSE)

Nouméa

Papeete

Niamey

Ouagadougou

Tunis

La Paz

Quito