A variational formulation of the fast marching eikonal solver |

Unstructured (triangulated) grids have computational advantages over rectangular ones in three common situations:

- When the number of grid points can be substantially reduced by putting them on an irregular grid. This situation corresponds to irregular distribution of details in the propagation medium.
- When the computational domain has irregular boundaries. One possible kind of boundary corresponds to geological interfaces and seismic reflector surfaces (Wiggins et al., 1993). Another type of irregular boundary, in application to traveltime computations, is that of seismic rays. The method of bounding the numerical eikonal solution by ray envelopes has been introduced recently by Abgrall and Benamou (1996).
- When the grid itself needs to be dynamically updated to maintain a certain level of accuracy in the computation.

Computational aspects of triangular grid generation are outlined in Appendix A. A three-dimensional application would follow the same algorithmic patterns.

A variational formulation of the fast marching eikonal solver |

2013-03-03