diff --git a/dune/gfe/hyperbolichalfspacepoint.hh b/dune/gfe/hyperbolichalfspacepoint.hh index ebab7b2703bd24bb7928970b97070e248af68305..0c3a5bdc4df206408b4c1a450fae039d73ec8833 100644 --- a/dune/gfe/hyperbolichalfspacepoint.hh +++ b/dune/gfe/hyperbolichalfspacepoint.hh @@ -527,6 +527,12 @@ public: return Dune::FieldMatrix<T,N,N>(result); } + + /** \brief Scalar product of two tangent vectors */ + T metric(const TangentVector& v, const TangentVector& w) const + { + return v*w/(data_[N-1]*data_[N-1]); + } /** \brief Write unit vector object to output stream */ friend std::ostream& operator<< (std::ostream& s, const HyperbolicHalfspacePoint& p)