to the new scalable MatMatMult() that is done in PETSc? Is it possible to even compare the flop rates sensibly? http://arxiv.org/abs/1203.1692 This is what Jeff is planning on using in his SciDAC. Matt -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener