  We test current numerical implementations of laser-matter interactions by
comparison with exact analytical results. Focussing on photon emission
processes, it is found that the numerics accurately reproduce analytical
emission spectra in all considered regimes, except for the harmonic structures
often singled out as the most significant high intensity/multi-photon effects.
We find that this discrepancy originates in the use of the locally constant
field approximation.
