fixed bug in least squares subroutine. Added mpigetrank==0 so that it only prints the extrapolation once