The description of the passage from wave optics to geometrical optics claims that light rays are the integral curves of a certain vector field (the Poynting vector direction, normalized to 1). Here are the details, could you fill the blanks:
The wavelength $\lambda$ is much smaller than all other characteristic lengths.
The setting is a "nice" medium (with spatially varying refraction index $n(x)$) through which an "almost plane" wave is propagating. The wave (if linearly polarized) turns out to be representable by $\vec{E}{(x,t)} = \vec{E}_0\exp(i(\chi{(x)}-\omega t))$ and $\vec{B}{(x,t)} = \vec{B}_0\exp(i(\chi{(x)}-\omega t))$ with constant $\vec{E}_0$ and $\vec{B}_0$.
Maxwell's equations imply $(\vec{\nabla}\chi)^2=\frac{n^2\omega^2}{c^2}$ and a time-averaged Poynting vector $\vec{S}=\frac{c}{n}\vec{s}$, where $\vec{s}$ is the unit vector $\vec{s}=\frac{\vec{\nabla}\chi}{n\omega/c}$.
The integral curves of the field of unit vectors $\vec{s}$ are the light rays.
Working through the equations results for this for a ray trajectory $X(\tau)$ (where $\tau$ is just a parameter): $\frac{d}{d\tau}(n\frac{d\vec{X}}{d\tau})=\vec{\nabla} n$
How does one prove the jump from 2 to 3. Why follow the unit vectors and not the Poynting vectors themselves? Or even, why should the light rays be tangential to the Poynting vectors at all (besides intuition like "light rays should transport energy")?
Could someone give me the proof of 3. or point me to a reference?