This actually returned the transposed of the inverse which is a severe bug.
added bugfix label
It's very surprising that no one ever noticed this. This code seems to be badly tested or rarely used. Unless you're only using grids with uniform grids, this bug has a severe impact on solutions.
The method was just added a year ago I think and probably not used or tested at all. I for example still use the old JIT in my assembly code.
I suspected this.
We're actually using jacobianInverse
in a dune-functions example. But this uses a uniform YaspGrid
, where you don't notice the bug.
enabled an automatic merge when the pipeline for 7b5bff97 succeeds
approved this merge request
merged
mentioned in commit 2fac40ac