Previously, the matrix was always assembled using a P1-basis, even if we used
it to measure errors of vectors pertaining to P2 or other spaces.  Nobody noticed
this, because the h1SemiNorm routine would not fail if the vector was larger
than the matrix.   Sigh...

[[Imported from SVN: r9770]]