We present a scalar beam propagation method (BPM) that is tailored for step
-index waveguides. For many structures, especially periodic ones such as st
ar couplers, it is many times faster than and avoids the spatial quantizati
on of existing BPMs. It can handle large propagation angles but not very la
rge effective index steps and does not perfectly preserve energy. We provid
e the two-dimensional version in C-language code.