GEOWARE: An optimised Matlab software for determination of high-frequency geoid model using relational operators

Abstract

We present an optimised software (GEOWARE) for determination of high-frequency geoid height using terrestrial gravity measurements. The optimisation of Stokes integral is based on the extraction of a local area with a radius of a few hundreds kilometres around the computation point which complies with the specified spherical cap sizes. The extraction step is highly important because it detaches the dispensable compartments of the grid which are far from the computation domain. That makes it convenient to avoid passing through the compartments of the entire grid to test whether the spherical distances comply with the truncated cap size or not. Matlab relational operators and vectorisation are powerful optimisation tools because they can replace conditional statements and nested loops efficiently. GEOWARE has been compared with a non- optimised code over different sizes of cap size and it shows a significant improvement in the performance. The run time of GEOWARE in all cap sizes has up to 5 times smaller than that of the code before optimisation. GEOWARE is also compatible with modified Stokes, Newton and Poisson kernels.