Abstract: We
propose an intermediate 'Semi-Lagrangian'
formulation of Geometric Optics type problems which
stands between the classical 'Lagrangian'
Hamiltonian systems and its 'Eulerian'
Hamilton-Jacobi Partial Differential Equations
counterpart. The goal is to design a numerical
method which takes the best of both world in order
to compute numerically the so called 'multi-valued'
solution. We present a numerical algorithm in the
'paraxial' 2-D case, discuss its limitation and its
possible extension to higher dimensions.
|