Dépôt numérique

Task-parallel implementation of 3D shortest path raytracing for geophysical applications.

Giroux, Bernard; Larouche, Benoît (2013). Task-parallel implementation of 3D shortest path raytracing for geophysical applications. Computers and Geosciences , vol. 54 . p. 130-141. DOI: 10.1016/j.cageo.2012.12.005.

Ce document n'est pas hébergé sur EspaceINRS.


This paper discusses two variants of the shortest path method and their parallel implementation on a shared-memory system. One variant is designed to perform raytracing in models with stepwise distributions of interval velocity while the other is better suited for continuous velocity models. Both rely on a discretization scheme where primary nodes are located at the corners of cuboid cells and where secondary nodes are found on the edges and sides of the cells. The parallel implementations allow raytracing concurrently for different sources, providing an attractive framework for ray-based tomography. The accuracy and performance of the implementations were measured by comparison with the analytic solution for a layered model and for a vertical gradient model. Mean relative error less than 0.2% was obtained with 5 secondary nodes for the layered model and 9 secondary nodes for the gradient model. Parallel performance depends on the level of discretization refinement, on the number of threads, and on the problem size, with the most determinant variable being the level of discretization refinement (number of secondary nodes). The results indicate that a good trade-off between speed and accuracy is achieved with the number of secondary nodes equal to 5. The programs are written in C++ and rely on the Standard Template Library and OpenMP.

Type de document:
Mots-clés libres: parallel; raytracing; shared memory; shortest path
Centre: Centre Eau Terre Environnement
Date de dépôt: 05 déc. 2016 20:50
Dernière modification: 05 déc. 2016 20:50
URI: http://espace.inrs.ca/id/eprint/3440

Actions (Identification requise)

Modifier la notice Modifier la notice