From a1a00afde0714373c7b327ce57e2d95a7209eea5 Mon Sep 17 00:00:00 2001
From: Lisa Julia Nebel <lisa_julia.nebel@tu-dresden.de>
Date: Thu, 8 Oct 2020 10:30:09 +0200
Subject: [PATCH] Correct typo in the calculation of the alternator tensor
---
dune/gfe/nonplanarcosseratshellenergy.hh | 2 +-
dune/gfe/surfacecosseratenergy.hh | 2 +-
2 files changed, 2 insertions(+), 2 deletions(-)
diff --git a/dune/gfe/nonplanarcosseratshellenergy.hh b/dune/gfe/nonplanarcosseratshellenergy.hh
index 85e97321..48506937 100644
--- a/dune/gfe/nonplanarcosseratshellenergy.hh
+++ b/dune/gfe/nonplanarcosseratshellenergy.hh
@@ -230,7 +230,7 @@ energy(const typename Basis::LocalView& localView,
for (int alpha=0; alpha<2; alpha++)
for (int beta=0; beta<2; beta++)
- c += sqrt(aScalar) * eps[alpha][beta] * Dune::GFE::dyadicProduct(aContravariant[alpha], aContravariant[beta]);
+ c += aScalar * eps[alpha][beta] * Dune::GFE::dyadicProduct(aContravariant[alpha], aContravariant[beta]);
// Second fundamental form
// The derivative of the normal field
diff --git a/dune/gfe/surfacecosseratenergy.hh b/dune/gfe/surfacecosseratenergy.hh
index 1aad21af..b7b5d122 100644
--- a/dune/gfe/surfacecosseratenergy.hh
+++ b/dune/gfe/surfacecosseratenergy.hh
@@ -245,7 +245,7 @@ RT energy(const typename Basis::LocalView& localView,
for (int alpha=0; alpha<2; alpha++)
for (int beta=0; beta<2; beta++)
- c += sqrt(aScalar) * eps[alpha][beta] * Dune::GFE::dyadicProduct(aContravariant[alpha], aContravariant[beta]);
+ c += aScalar * eps[alpha][beta] * Dune::GFE::dyadicProduct(aContravariant[alpha], aContravariant[beta]);
// Second fundamental form
// The derivative of the normal field
--
GitLab