Extend jacobian interface of geometry classes
- Jun 01, 2022
-
-
Carsten Gräser authored
-
Carsten Gräser authored
-
Carsten Gräser authored
-
Carsten Gräser authored
-
Carsten Gräser authored
-
This adds methods jacobian(local)
and jacobianInverse(local)
and corresponding typedefs Jacobian
and JacobianInverse
for the return types to the Geometry
interface. More specifically:
Geometry
implementations provided by dune-geometry. All implementations return either a FieldMatrix
or a DiagonalMatrix
by value, even if jacobianTransposed()
or jacobianInverseTransposed()
provide cached references. Since this only involves transposing the cached matrix, it is still cheap, while not requiring to cache another to matrices.checkgeometry.hh
by the new interface. This requires that all geometry implementations adopt the new interface as well.checkgeometry.hh
this implicitly extends the requirements for the grid geometries. But a default implementation is provided there by dune-grid!577 (merged). Hence the latter should be merged before the present MR.