Open AccessEfficient parallel inversion using the Neighbourhood Algorithm
Author(s) -
Rickwood P.,
Sambridge M.
Publication year - 2006
Publication title -
geochemistry, geophysics, geosystems
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.928
H-Index - 136
ISSN - 1525-2027
DOI - 10.1029/2006gc001246
Subject(s) - scalability , inversion (geology) , computer science , neighbourhood (mathematics) , algorithm , parallel algorithm , implementation , parallel computing , mathematical optimization , mathematics , geology , paleontology , mathematical analysis , structural basin , database , programming language
Issues controlling efficient parallel implementations of a popular direct search inversion algorithm are analyzed and discussed. A naive parallelization of a particular method, the Neighbourhood parameter search algorithm, leads to inefficient use of parallel architecture through lack of scalability and intolerance to hardware faults. These factors are quantified, and their origins are explained. A reformulation of the algorithm leads to dramatically improved performance when the cost of the forward problem is low and the number of unknowns is high. Numerical examples are used to illustrate the main results. Factors in the original Neighbourhood Algorithm which lead to poor parallel performance are likely to be present in other ensemble‐based inversion or global optimization algorithms. Hence the algorithmic solutions proposed may have widespread application.