dune-elasticity issueshttps://git.imp.fu-berlin.de/agnumpde/dune-elasticity/-/issues2020-05-05T08:57:49Zhttps://git.imp.fu-berlin.de/agnumpde/dune-elasticity/-/issues/1LocalDensity uses the same type for coordinates and derivatives2020-05-05T08:57:49Zoliver.sander_at_tu-dresden.deoliver.sander@tu-dresden.deLocalDensity uses the same type for coordinates and derivativesThe signature of the main interface method of `LocalDensity` is
virtual field_type operator() (const FieldVector<field_type,dim>& x, const FieldMatrix<field_type,dim,dim>& gradient) const = 0;
Note how it uses the same number typ...The signature of the main interface method of `LocalDensity` is
virtual field_type operator() (const FieldVector<field_type,dim>& x, const FieldMatrix<field_type,dim,dim>& gradient) const = 0;
Note how it uses the same number type for `x` and `gradient`. This is unfortunate, because when using algorithmic differentiation, `gradient` will typically be an active type but `x` will not.Patrick Jaappatrick.jaap@tu-dresden.dePatrick Jaappatrick.jaap@tu-dresden.de