Skip to content

[bugfix] Fix jacobian part of Piola transformations

Carsten Gräser requested to merge bugfix/fix-piola-transform into master

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.

Merge request reports