bugfix: the correct energy is *half* the integral.

[[Imported from SVN: r5807]]

src/harmonicenergystiffness.hh

@@ -89,7 +89,7 @@ energy(const Entity& element,
 
     }
 
-    return energy;
+    return 0.5 * energy;
 }
 
 #endif