Removed "const" from zero() function (since it is not constant).
Halved memory consumption and execution time of petscsparsematrix::buildinternalstructure (using the symmetric nonzero-structure).
FEInterpolation::evaldNdx now returns the determinant of the jacobian. They are almost always needed at the same time, and computing the determinant in evaldNdx is almost for free.