Moved to dune-solvers. We don't need it in ag-common.

[[Imported from SVN: r2809]]

dune-solvers/computeenergy.hh

+#ifndef CONTACT_COMPUTE_ENERGY_HH
+#define CONTACT_COMPUTE_ENERGY_HH
+
+template <class M, class V>
+double computeEnergy(const M& mat, const V& x, const V& rhs)
+{
+    V tmp(x.size());
+    tmp = 0;
+    mat.umv(x, tmp);
+
+//     std::cout << "tmp*x: " << (tmp*x) << "    rhs*x: " << (rhs*x) << std::endl;
+//     exit(0);
+    return 0.5*(tmp*x) - rhs*x;
+}
+
+template <class M, class V>
+double computeEnergy(const M& mat, const V& x)
+{
+    V tmp(x.size());
+    tmp = 0;
+    mat.umv(x, tmp);
+    
+    return 0.5*(tmp*x);
+}
+
+#endif