Code modifications to use PETSc exclusively for the momentum matrix, as well as reducing the PETSc preallocation memory usage somewhat.