add exact solution (as initial) and work on FD formulas; note FIXMEs remain