From b068e2839e0801390e9876caaded46e5a0453d7e Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Thu, 29 Apr 2010 13:23:27 +0000
Subject: [PATCH] add methods projectOntoTangentSpace() and orthonormalFrame()
[[Imported from SVN: r5993]]
---
dune/gfe/realtuple.hh | 19 +++++++++++++++++++
1 file changed, 19 insertions(+)
diff --git a/dune/gfe/realtuple.hh b/dune/gfe/realtuple.hh
index a45e027e..c4abeed2 100644
--- a/dune/gfe/realtuple.hh
+++ b/dune/gfe/realtuple.hh
@@ -86,10 +86,29 @@ public:
return result;
}
+ /** \brief Project tangent vector of R^n onto the tangent space */
+ EmbeddedTangentVector projectOntoTangentSpace(const EmbeddedTangentVector& v) const {
+ return v;
+ }
+
/** \brief The global coordinates, if you really want them */
const Dune::FieldVector<double,N>& globalCoordinates() const {
return data_;
}
+
+ /** \brief Compute an orthonormal basis of the tangent space of S^n.
+
+ This basis is of course not globally continuous.
+ */
+ Dune::FieldMatrix<double,N,N> orthonormalFrame() const {
+
+ Dune::FieldMatrix<double,N,N> result;
+
+ for (int i=0; i<N; i++)
+ for (int j=0; j<N; j++)
+ result[i][j] = (i==j);
+ return result;
+ }
/** \brief Write LocalKey object to output stream */
friend std::ostream& operator<< (std::ostream& s, const RealTuple& realTuple)
--
GitLab