fixed poisson solver, I think. adaptive newton needs improvement.