// This file contains a simple prototypical Poisson solver for the unit square domain. // The PDE to be solved reads: // -Laplace(u) = f in the domain [0,1]x[0,1 ...