A Structure-Preserving Scheme for Computing the Gradient of the Eikonal Equation
Abstract
We propose a first-order, structure-preserving method for directly computing the normalized traveltime gradient of the eikonal equation in smooth first-arrival regions. The method combines analytical derivatives of spherical linear interpolation (SLERP) in a Godunov-type upwind discretization with spherical forward Euler (SFE) pseudo-time integration. This construction gives a tangent residual and keeps every iterate on 𝕊2 without projection or normalization. The SLERP discretization exactly preserves a family of curved orientation equilibria. Removing the background discretization bias improves absolute orientation accuracy near these equilibria, allowing a prescribed error tolerance to be reached on a coarser grid. We prove the equilibrium-preservation property and, separately, a smaller angular error for SFE than normalized Euler in a one-cell relaxation under the stated step-size condition. Two- and three-dimensional experiments examine convergence, geometric preservation, and computational cost. Comparisons with gradient recovery from scalar traveltimes show more accurate directions on the same grid in the curved media considered. The method provides direction fields for wave propagation and geometric optics, with potential use in adjoint-based inverse problems.
Related
- A Structure-Preserving LOBPCG Method for Computing Several of the Largest Generalized Singular Triplets
- Uniqueness and structure-preserving numerical simulations for an Ultrahyperbolic Schrödinger equation
- GradAgent: A Knowledge-Guided Multi-Agent System for Structure-Preserving Gradient-Flow Computation with an Application to Multicomponent Vesicle Dynamics
- Wave Breaking and Structure-Preserving Fully Discrete Method for Stochastic Camassa--Holm Equation with Additive Noise
- Rate of convergence of a fully discrete structure-preserving midpoint scheme for the stochastic Landau--Lifshitz--Gilbert equation
BibTeX
@misc{leung2026structurepreserving,
title = {{A Structure-Preserving Scheme for Computing the Gradient of the Eikonal Equation}},
author = {Ka Lun Leung and Shingyu Leung},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.28692},
eprint = {2610.28692},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.28692},
}
RIS
TY - UNPB
TI - A Structure-Preserving Scheme for Computing the Gradient of the Eikonal Equation
AU - Ka Lun Leung
AU - Shingyu Leung
PY - 2026
DA - 2026/10/20
PB - arXiv
DO - 10.48550/arxiv.2610.28692
UR - https://arxiv.org/abs/2610.28692
AB - We propose a first-order, structure-preserving method for directly computing the normalized traveltime gradient of the eikonal equation in smooth first-arrival regions. The method combines analytical derivatives of spherical linear interpolation (SLERP) in a Godunov-type upwind discretization with spherical forward Euler (SFE) pseudo-time integration. This construction gives a tangent residual and keeps every iterate on 𝕊^2 without projection or normalization. The SLERP discretization exactly preserves a family of curved orientation equilibria. Removing the background discretization bias improves absolute orientation accuracy near these equilibria, allowing a prescribed error tolerance to be reached on a coarser grid. We prove the equilibrium-preservation property and, separately, a smaller angular error for SFE than normalized Euler in a one-cell relaxation under the stated step-size condition. Two- and three-dimensional experiments examine convergence, geometric preservation, and computational cost. Comparisons with gradient recovery from scalar traveltimes show more accurate directions on the same grid in the curved media considered. The method provides direction fields for wave propagation and geometric optics, with potential use in adjoint-based inverse problems.
ER -