made the comment about np.abs()**2 a bit less opaque, and referenced MJ's discussion in optics.py
(#563)