Skip to content
Snippets Groups Projects

[bugfix] Fix jacobian part of Piola transformations

Merged 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

Loading
Loading

Activity

Filter activity
  • Approvals
  • Assignees & reviewers
  • Comments (from bots)
  • Comments (from users)
  • Commits & branches
  • Edits
  • Labels
  • Lock status
  • Mentions
  • Merge request status
  • Tracking
Please register or sign in to reply
Loading