Abstract:Efficient computation of the geopotential gradient is essential for numerical propagators, particularly in scenarios involving low Earth orbits. Conventional geopotential calculations are based on spherical harmonics series, which become computationally demanding as the degree/order increases. This computational burden can be mitigated by means of parallelized algorithms. Additionally, certain situations lend themselves to high parallelization, such as the propagation of space debris catalogs, satellite mega-constellations, or the dispersion of particles resulting from a space collision event. This paper introduces an optimized Graphics Processing Unit (GPU) implementation designed to facilitate extensive parallelization in the geopotential gradient calculation. The formulation developed in this study is not specific to any GPU. However, to illustrate the low-level optimizations necessary for an efficient implementation, we selected the Compute Unified Device Architecture (CUDA) as the dominant and de facto standard in parallel computing. Nevertheless, most of the concepts and optimizations presented in this paper are also valid for other GPU architectures. Built upon the spherical harmonic expansion using the Cunningham formulation, which is well-suited for GPU computations, our implementation offers several variants with different tradeoffs between speed and accuracy. Besides GPU double precision, we introduced a mixed precision arithmetic –a hybrid between single and double precision– that exploits GPU capabilities with a low penalty in accuracy. The proposed algorithm was implemented as a software reusable module, and its performance was evaluated against GMAT, GODOT, and Orekit astrodynamic codes. The algorithm's accuracy in double precision is comparable to such codes. The mixed precision version showed enough accuracy for LEO satellite propagation, with around 1 m difference in four days. Testing across different CUDA architectures revealed very high speed-up factors compared to a single CPU, reaching a speed-up of 645 for the mixed precision variant and 450 for the double precision one in the propagation of about 3200 objects with a geopotential of degree/order 126 × 126 using an A100 GPU device.

Efficient method for parallel computation of geodesic transformation on CPU

HyGrid: A CPU-GPU Hybrid Convolution-Based Gridding Algorithm in Radio Astronomy.

Large Image Fast Compositing Based on Uniform Sample

Fast algorithm for dilation and erosion using arbitrary flat structuring element: Improvement of urbach and Wilkinson's algorithm to GPU computing

Parallel stream computing implementation of GIST algorithm on heterogeneous platform

GPU Accelerated Convex Hull Computation.

Technical Section: Real-time B-spline Free-Form Deformation via GPU acceleration

Parallel Implementation of Dense Optical Flow Computation on Many-Core Processor.

Parallel Implementation and Optimization of High Definition Video Real-Time Dehazing

Balancing Cpu And Gpu: Real-Time Visualization Of Large Scale 3d Scanning Models

A Systematic Parallel Strategy for Generating Contours from Large-Scale DEM Data Using Collaborative CPUs and GPUs

A Graphics Processing Unit Implementation and Optimization for Parallel Double-Difference Seismic Tomography

Efficient deformable geometry image-maps.

The Accelerated Fusion Algorithm of Remote Sensing Image Data Based on GPU

Accelerating Spatial Cross-Matching on CPU-GPU Hybrid Platform With CUDA and OpenACC

Parallel Spatial Hashing for Collision Detection of Deformable Surfaces

A Strategy for Parallelising Polygon Rasterisation Algorithms Using Multi-Core Cpus

Parallel and Scalable Heat Methods for Geodesic Distance Computation

Efficient computation of the geopotential gradient in Graphic Processing Units

Efficient Same-Dimensional Implicit Time Advancement Parallel Scheme and Optimization Methods for the Iteration Parameters Using a Graphics-Processing Unit

Parallel Smoothers for Matrix-Based Geometric Multigrid Methods on Locally Refined Meshes Using Multicore CPUs and GPUs