The derivation of a semiclassical surface hopping procedure from a formally exact solution of the Schrodinger equation is discussed. The fact that the derivation proceeds from an exact solution guarantees that all phase equation is discussed. The fact that the derivation proceeds from an exact solution g terms are completely and accurately included. Numerical evidence shows the method to be highly accurate. A Monte Carlo implementation of this method is considered, and recent work to significantly improve the statistical accuracy of the Monte Carlo approach is discussed.