MatrixHelper

petsc output still needs work:

• using the builtin print functions, whose format is not consistent with the other operator types.
• manual implementation of petsc print is involved, because petsc interface lacks easy global access functions