1d eikonal solver uses implicit numerical diffusion for better stability
remove accidental factor of 0.5 from numerical diffusion function
semi-implicit solver for heat equation

this allows larger timesteps, plus CN has better asymptotic error properties
than Euler

fixes #4