From 3a6a555f8daa524aab617938803fd938ad151600 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Tue, 3 Sep 2013 16:29:35 +0000
Subject: [PATCH] Use frobenius_norm2() instead of frobenius_norm()
Saves us one call to sqrt().
[[Imported from SVN: r9400]]
---
dune/gfe/cosseratenergystiffness.hh | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/dune/gfe/cosseratenergystiffness.hh b/dune/gfe/cosseratenergystiffness.hh
index 7868e45e..5505e0a3 100644
--- a/dune/gfe/cosseratenergystiffness.hh
+++ b/dune/gfe/cosseratenergystiffness.hh
@@ -365,9 +365,9 @@ public:
RT curvatureEnergy(const Tensor3<double,3,3,3>& DR) const
{
#ifdef DONT_USE_CURL
- return mu_ * std::pow(L_c_ * DR.frobenius_norm(),q_);
+ return mu_ * std::pow(L_c_ * DR.frobenius_norm2(),q_/2.0);
#else
- return mu_ * std::pow(L_c_ * curl(DR).frobenius_norm(),q_);
+ return mu_ * std::pow(L_c_ * curl(DR).frobenius_norm2(),q_/2.0);
#endif
}
--
GitLab