check if a photoionization edge is crossed (from Doppler shifts) along path element ds

See #73