Implemented basic mixed dirichlet-neumann iterative solver
Spectral solver is wrong, or my iterative solver is wrong...