Somewhat working Poisson example, modulo boundary wrinkle
Marginally working Poisson solver