Poisson solver using relaxation method (MPI version)