[petsc-users] Element by element comparison of matrices

PETSc Users, I want to compare two Jacobians matrices (one with derivatives calculated analytically; the other numerically). I want to apply relatives and/or absolute tolerances. Does anyone know if such capability is built into PETSc? I cannot find anything other the MatEqual() with compares down to the bit. Otherwise, I plan to use MatGetValues() and compare the elements individually. Just hoping there is something more convenient hidden somewhere. Thanks, Glenn