[bugfix] Fix jacobian part of Piola transformations

Since jacobian(Inverse)Transposed may be a sparse matrix (e.g. DiagonalMatrix for YaspGrid) we must not use jacobian(Inverse)Transposed[i][j]. Instead iterate over the sparse row entries using sparseRange.

This is now covered by extending the test of GlobalValuedFiniteElement to YaspGrid and higher order RT elements.