From 9719c84bf1ea2f42c83cb51d8839d75de613d479 Mon Sep 17 00:00:00 2001
From: Oliver Sander <oliver.sander@tu-dresden.de>
Date: Mon, 1 Jul 2024 15:43:53 +0200
Subject: [PATCH] Minor simplification of CosseratRodEnergy
---
dune/gfe/assemblers/cosseratrodenergy.hh | 6 ++----
1 file changed, 2 insertions(+), 4 deletions(-)
diff --git a/dune/gfe/assemblers/cosseratrodenergy.hh b/dune/gfe/assemblers/cosseratrodenergy.hh
index 5545ab41..2cb9fb43 100644
--- a/dune/gfe/assemblers/cosseratrodenergy.hh
+++ b/dune/gfe/assemblers/cosseratrodenergy.hh
@@ -5,7 +5,6 @@
#include <dune/common/fmatrix.hh>
#include <dune/common/version.hh>
-#include <dune/common/transpose.hh>
#include <dune/istl/matrix.hh>
#include <dune/geometry/quadraturerules.hh>
@@ -248,12 +247,11 @@ namespace Dune::GFE {
const Entity& element,
const FieldVector<double,1>& pos) const
{
- const auto jit = element.geometry().jacobianInverseTransposed(pos);
+ const auto geometryJacobianInverse = element.geometry().jacobianInverse(pos);
auto value = localInterpolation.evaluate(pos);
- auto referenceDerivative = localInterpolation.evaluateDerivative(pos);
- auto derivative = referenceDerivative * transpose(jit);
+ auto derivative = localInterpolation.evaluateDerivative(pos) * geometryJacobianInverse;
using Number = std::decay_t<decltype(derivative[0][0])>;
FieldVector<Number,3> r_s = {derivative[0][0], derivative[1][0], derivative[2][0]};
--
GitLab