TY - GEN
T1 - Rapid inversion of large airborne AEM data datasets utilizing massively parallel coprocessors
AU - Kirkegaard, C.
AU - Andersen, K.
AU - Christiansen, A. V.
AU - Auken, E.
AU - Boesen, T.
PY - 2015
Y1 - 2015
N2 - While full 2D and 3D inversion schemes are emerging for airborne TEM datasets, the workhorse of large scale surveying still remains the well established 1D forward formulation. Modern airborne TEM surveys typically span several thousand line kilometers of data, which can be very time consuming to invert even within the framework of 1D forward modeling. Here, we present how we have modified our existing inversion code to offload its 1D forward computations to massively parallel Intel Xeon Phi co-processors. A prerequisite for good performance on this type of next generation technology is for a code to provide not only good parallel scaling, but also make efficient use of vector instructions. The latter is not possible without modification to the established framework of 1D forward modeling and we demonstrate how this problem can be overcome in a straight forward manner. We show how the modified algorithm provides virtually ideal parallel scaling and almost ideal use of vector instructions on both multi-core processors and massively parallel Intel Xeon Phi co-processors. The use of vector instructions alone provide a speedup of almost 8x on the co-processor, allowing for full inversion of several thousand line kilometers of data per hour.
AB - While full 2D and 3D inversion schemes are emerging for airborne TEM datasets, the workhorse of large scale surveying still remains the well established 1D forward formulation. Modern airborne TEM surveys typically span several thousand line kilometers of data, which can be very time consuming to invert even within the framework of 1D forward modeling. Here, we present how we have modified our existing inversion code to offload its 1D forward computations to massively parallel Intel Xeon Phi co-processors. A prerequisite for good performance on this type of next generation technology is for a code to provide not only good parallel scaling, but also make efficient use of vector instructions. The latter is not possible without modification to the established framework of 1D forward modeling and we demonstrate how this problem can be overcome in a straight forward manner. We show how the modified algorithm provides virtually ideal parallel scaling and almost ideal use of vector instructions on both multi-core processors and massively parallel Intel Xeon Phi co-processors. The use of vector instructions alone provide a speedup of almost 8x on the co-processor, allowing for full inversion of several thousand line kilometers of data per hour.
UR - http://www.scopus.com/inward/record.url?scp=84956667482&partnerID=8YFLogxK
U2 - 10.3997/2214-4609.201413867
DO - 10.3997/2214-4609.201413867
M3 - Conference article in proceedings
AN - SCOPUS:84956667482
T3 - 1st European Airborne Electromagnetics Conference - Held at Near Surface Geoscience 2015
SP - 36
EP - 40
BT - 1st European Airborne Electromagnetics Conference - Held at Near Surface Geoscience 2015
PB - European Association of Geoscientists and Engineers
T2 - 1st European Airborne Electromagnetics Conference
Y2 - 6 September 2015 through 10 September 2015
ER -