From 353e29b5dda51bcbaa6f7acc06cd524fe193246a Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Wed, 13 Dec 2006 13:48:39 +0000
Subject: [PATCH] bugfix in computeEnergy()
[[Imported from SVN: r1086]]
---
src/rodassembler.cc | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/rodassembler.cc b/src/rodassembler.cc
index fe3e3bef..237ba991 100644
--- a/src/rodassembler.cc
+++ b/src/rodassembler.cc
@@ -663,7 +663,7 @@ computeEnergy(const std::vector<Configuration>& sol) const
FieldVector<double,blocksize> referenceStrain = getStrain(referenceConfiguration_, it, quadPos);
for (int i=0; i<3; i++)
- energy += weight * 0.5 * A_[i] * (strain[i] - referenceStrain[i]) * (strain[i] * referenceStrain[i]);
+ energy += weight * 0.5 * A_[i] * (strain[i] - referenceStrain[i]) * (strain[i] - referenceStrain[i]);
}
@@ -685,7 +685,7 @@ computeEnergy(const std::vector<Configuration>& sol) const
// Part II: the bending and twisting energy
for (int i=0; i<3; i++)
- energy += weight * 0.5 * K_[i] * (strain[i+3] - referenceStrain[i+3]) * (strain[i+3] * referenceStrain[i+3]);
+ energy += weight * 0.5 * K_[i] * (strain[i+3] - referenceStrain[i+3]) * (strain[i+3] - referenceStrain[i+3]);
}
--
GitLab