An OpenCL implementation of ellipsoidal harmonics
The result's identifiers
Result code in IS VaVaI
<a href="https://www.isvavai.cz/riv?ss=detail&h=RIV%2F00025615%3A_____%2F15%3A%230002188" target="_blank" >RIV/00025615:_____/15:#0002188 - isvavai.cz</a>
Result on the web
<a href="http://dx.doi.org/10.1007/1345_2015_59" target="_blank" >http://dx.doi.org/10.1007/1345_2015_59</a>
DOI - Digital Object Identifier
<a href="http://dx.doi.org/10.1007/1345_2015_59" target="_blank" >10.1007/1345_2015_59</a>
Alternative languages
Result language
angličtina
Original language name
An OpenCL implementation of ellipsoidal harmonics
Original language description
The technology progress today makes it possible to treat most of the problems of physical geodesy by means of numerical arrangements hardly imaginable earlier. Nevertheless, considering an evaluation of spheroidal (spherical and ellipsoidal) harmonic functions in our typical tasks, we still observe a huge performance gap between our demands and capabilities of common CPUs. Methods used for calculating associated Legendre functions are mostly recursive and thus sequential. Therefore, it is challenging, but feasible, to arrange the processing of Legendre functions in a way that reduces memory utilisation and admits massive parallelism. Following this aim, we developed a streaming-parallel algorithm for computing oblate spheroidal harmonic functions and their derivatives. The algorithm is free of assumptions concerning the func-tion arguments, maximal degree/order or number of computation points and can be utilised on any data type, like a vector or scalar float, double or even integer numbers. Besides, it solves floating-point issues in the numerical treatment of Legendre functions. We demonstrate its Open Computing Language (OpenCL) implementation on a general-purpose graphics processing unit (GPGPU), which is ideal for its inexpensive computational power of some TFlops. Added performance benchmarks lead to the conclusion that our implementation on a single GPGPU device substantially outperforms recent multi-core CPUs, free of any precision penalty. Furthermore, thanks to the OpenCL standard, we can benefit from an excellent portability and scalability over heterogeneous parallel platforms. Let us note finally, that the topic presented is a matter of importance in many other application fields, not only in physical geodesy.
Czech name
—
Czech description
—
Classification
Type
D - Article in proceedings
CEP classification
DE - Earth magnetism, geodesy, geography
OECD FORD branch
—
Result continuities
Project
Result was created during the realization of more than one project. More information in the Projects tab.
Continuities
P - Projekt vyzkumu a vyvoje financovany z verejnych zdroju (s odkazem do CEP)
Others
Publication year
2015
Confidentiality
S - Úplné a pravdivé údaje o projektu nepodléhají ochraně podle zvláštních právních předpisů
Data specific for result type
Article name in the collection
VIII Hotine-Marussi Symposium on Mathematical Geodesy
ISBN
978-3-319-24548-5
ISSN
0939-9585
e-ISSN
—
Number of pages
9
Pages from-to
—
Publisher name
Springer-Verlag
Place of publication
Berlin
Event location
Rome
Event date
Jun 17, 2013
Type of event by nationality
WRD - Celosvětová akce
UT code for WoS article
—