From 11db41d05fc8b753d96837b8d5c344913d486fbc Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Mon, 22 Mar 2010 10:45:38 +0000
Subject: [PATCH] add a few missing methods. I can compute harmonic maps for
RealTuples now
[[Imported from SVN: r5783]]
---
src/realtuple.hh | 22 +++++++++++++++++++++-
1 file changed, 21 insertions(+), 1 deletion(-)
diff --git a/src/realtuple.hh b/src/realtuple.hh
index 8681b3c2..a18c7ed9 100644
--- a/src/realtuple.hh
+++ b/src/realtuple.hh
@@ -15,6 +15,8 @@ class RealTuple
{
public:
+ typedef double ctype;
+
typedef Dune::FieldVector<double,N> EmbeddedTangentVector;
typedef Dune::FieldVector<double,N> TangentVector;
@@ -29,6 +31,21 @@ public:
data_ = v;
}
+ /** \brief Copy constructor */
+ RealTuple(const RealTuple<N>& other)
+ : data_(other.data_)
+ {}
+
+ /** \brief Constructor from FieldVector*/
+ RealTuple(const Dune::FieldVector<double,N>& other)
+ : data_(other)
+ {}
+
+ RealTuple& operator=(const Dune::FieldVector<double,N>& other) {
+ data_ = other;
+ return *this;
+ }
+
/** \brief The exponention map */
static RealTuple exp(const RealTuple& p, const TangentVector& v) {
return RealTuple(p.data_+v);
@@ -53,7 +70,10 @@ public:
Unlike the distance itself the squared distance is differentiable at zero
*/
static EmbeddedTangentVector derivativeOfDistanceSquaredWRTSecondArgument(const RealTuple& a, const RealTuple& b) {
- return -2*(a.data_ - b.data_);
+ EmbeddedTangentVector result = a.data_;
+ result -= b.data_;
+ result *= -2;
+ return result;
}
/** \brief Compute the Hessian of the squared distance function keeping the first argument fixed
--
GitLab