Data and Code for "Fast computation of terrain-induced gravitational and magnetic effects on arbitrary undulating surfaces"

We propose a new algorithm for the fast computation of terrain-induced gravitational and magnetic effects on arbitrary undulating surfaces. The algorithm is based on a combination of CANDECOMP/PARAFAC (CP) tensor decomposition of gravitational integral kernels and 2D Fast Fourier Transform (FFT) evaluation of discrete convolutions. We validate our algorithm by computing the gravitational potential (GP), the gravitational vector (GV), the gravity gradient tensor (GGT), and magnetic fields caused by densely-sampled topographic and bathymetric digital elevation models (DEMs) of selected mountainous areas around the globe. Both constant and variable density/magnetization models are considered. Comparisons between our new method and space-domain rigorous solutions show that with modeling errors well below existing instrumentation error levels, the calculation speed is accelerated thousands of times in all numerical tests. We release the source code of our algorithm in MATLAB language.