still working on pde-solver routines