Removed a function that didn't do anything.
Imposing Dirichlet and Neumann works for the 2D Poisson case at least.
Incorporated the new way of imposing BCs.