FMeikonal - Fast marching eikonal solver
FMeikonal < VelocityCube.H >TTCube.H pars
file with multiple shots (n2 - number of shots, n1=3)
[0., o2 + 0.5*(n2-1)*d2, o3 + 0.5*(n3-1)*d3 shot location
[1,1,1] constant-velocity box around the source
[d1,d2,d3] constant-velocity box around the source
[1] velocity (1) or slowness (0)
velocity (or slowness) at the source [vel[0]]
[2] first (1), second (2) or third (3) order
Fast marching eikonal solver
the Hwt3d manpage, the fastmarch manpage, the Gfgradz manpage
Be careful with sources outside the domain. You should smooth the velocity field if you want high-order accuracy. Stability problems *may* occur if traveltime accuracy is beyond machine precision. The program uses constant velocity ray-tracing inside an initial box, and until it detects a velocity variation, whichever is greater.
seis/travel