added limits, so each part of the matrix is only computed once. Not tested, as I don't have mpi on this machine. Will test on Kratos later tonight