Changed PGAPACK_MPI to MPI_PGAPACK_USEMPI so it is grouped under MPI (very useful!)
removed a random extra loop

ToDo: Determine why the code ignores the upper limits