adjust reinitialization params
move front velocity methods to heat.jl
compute front velocity within a narrow band

All other values have the level set evolution set to the average of the front
values. This is probably okay for simple symmetric problems, but would be easy
to extend to higher order interpolation for other problems.
devectorization
update hill 2d test