diff --git a/src/localgeodesicfestiffness.hh b/src/localgeodesicfestiffness.hh index 7b9f5b5394451f428fc024024016b93ac63c5198..d22b91818f616d8a00eacaec1898b0dbb2979a61 100644 --- a/src/localgeodesicfestiffness.hh +++ b/src/localgeodesicfestiffness.hh @@ -52,6 +52,12 @@ class LocalGeodesicFEStiffness (i==2)*eps)); } + static void infinitesimalVariation(Rotation<2,double>& c, double eps, int i) + { + Dune::FieldVector<double,1> v(eps); + c = Rotation<2,double>::exp(c,v); + } + static void infinitesimalVariation(RealTuple<2>& c, double eps, int i) { Dune::FieldVector<double,2> v(0);