diff --git a/dune/gfe/averagedistanceassembler.hh b/dune/gfe/averagedistanceassembler.hh index f6676be291f4694dfaef50f75a865dbf4cbb1080..9d085d7c6a1efb357c687993ed0ff4be6e6465ea 100644 --- a/dune/gfe/averagedistanceassembler.hh +++ b/dune/gfe/averagedistanceassembler.hh @@ -3,8 +3,6 @@ #include <vector> -#include "rotation.hh" - template <class TargetSpace> class AverageDistanceAssembler { @@ -98,55 +96,4 @@ public: }; - -template <> -class AverageDistanceAssembler<Rotation<3,double> > -{ - typedef Rotation<3,double> TargetSpace; - - static const int size = TargetSpace::TangentVector::size; - -public: - - AverageDistanceAssembler(const std::vector<TargetSpace>& coefficients, - const std::vector<double>& weights) - : coefficients_(coefficients), - weights_(weights) - {} - - double value(const TargetSpace& x) const { - - double result = 0; - for (size_t i=0; i<coefficients_.size(); i++) { - double dist = TargetSpace::distance(coefficients_[i], x); - result += 0.5*weights_[i]*dist*dist; - } - - return result; - } - - void assembleGradient(const TargetSpace& x, - TargetSpace::TangentVector& gradient) const - { - DUNE_THROW(Dune::NotImplemented, "assembleGradient"); - } - - void assembleHessianApproximation(const TargetSpace& x, - Dune::FieldMatrix<double,size,size>& matrix) const - { - DUNE_THROW(Dune::NotImplemented, "assembleHessianApproximation"); - } - - void assembleHessian(const TargetSpace& x, - Dune::FieldMatrix<double,size,size>& matrix) const - { - DUNE_THROW(Dune::NotImplemented, "assembleHessian"); - } - - const std::vector<TargetSpace> coefficients_; - - const std::vector<double> weights_; - -}; - #endif