Two point seismic ray tracing in layered media using bending
We develop an efficient algorithm for two point seismic ray tracing in layered media by means of the bending method based on the Fermat's principle of stationary traveltime. We use the fact that the matrix of the second derivatives of traveltime function with respect to the points of intersection at the interfaces has a block tridiagonal structure. This feature enables efficient matrix inversion and leads to a fast implementation of the Newton Raphson method for finding a ray path with stationary traveltime.