Works on one proc. Currently uses identity matrices rather than the RNG.