High accuracy computational methods for the semiclassical Schrödinger equation