Quadrature points are not matched to the places Lutz' code puts them. Removed a bunch of redundant includes of the MPI header. Added new Esys_MPI.hpp header which has the extern "C" in it.
branches/refine/buckley/src/BuckleyDomain.cc