On ray tracing for sharp changing media.

Ray tracing is an integral part of modern imaging and inversion techniques. Several methods have been proposed depending on the requirements of the application. Some algorithms are best applied to fast changing material properties, like an interface between two differing media, while others are well suited to media with gradually changing properties, like composite materials. In this paper, an enhanced numerical algorithm for ray tracing is presented. Focus is given to solutions involving ordinary differential equations with initial-value conditions. The proposed algorithm is the result of a combination of two classical implementations and the authors show that it is well suited for media with both sharp and gradual changes in the index of refraction. Additionally, the authors present an application of ray path computation by using the technique known as the shooting method.