  We demonstrate that the singular value decomposition algorithm in conjunction
with the fast Fourier transform or finite difference procedures provides a
straightforward and accurate method for rapidly propagating electric fields in
the one-way Helmholtz formalism.
