adding forgotten compareDDS and compareminusDS functions
mpi.m and pic.cc and pi.cc all agree now on DDS, minusDS and phi on iteration 1