EdCaunt

Can we resurrect this? It fixes an issue which still exists (I accidentally replicated the fix in #2357)

This was also found to affect substitutions at forward and backward time steps (although this may have been fixed since?). However, a new MFE is: ``` import numpy as np...