diff --git a/fem/l2projection.hh b/fem/l2projection.hh deleted file mode 100644 index 1bc6e20d41f1b2c9d2a07200e7a3dd328580b8d1..0000000000000000000000000000000000000000 --- a/fem/l2projection.hh +++ /dev/null @@ -1,74 +0,0 @@ -// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*- -// vi: set et ts=4 sw=2 sts=2: -#ifndef DUNE_L2_PROJECTION_HH -#define DUNE_L2_PROJECTION_HH - -#include <dune/quadrature/quadraturerules.hh> - - -namespace Dune -{ - - /** \brief Projects an analytical function orthogonally (in the L2-sense) - onto a FE-space - */ - template <class DiscreteFunctionType> - class L2Projection - { - typedef typename DiscreteFunctionType::FunctionSpaceType FunctionSpaceType; - - public: - /** \brief Do the projection - * - * \tparam polOrd The order of the quadrature scheme used - * \tparam FunctionType The type of the class implementing the analytical function - */ - template <int polOrd, class FunctionType> - void project (FunctionType &f, DiscreteFunctionType &discFunc) - { - const typename DiscreteFunctionType::FunctionSpaceType - & functionSpace_= discFunc.getFunctionSpace(); - - discFunc.clear(); - - typedef typename FunctionSpaceType::GridType GridType; - typedef typename FunctionSpaceType::IteratorType IteratorType; - typedef typename DiscreteFunctionType::LocalFunctionType LocalFuncType; - - const int dim = GridType::dimension; - - typename FunctionSpaceType::RangeType ret (0.0); - typename FunctionSpaceType::RangeType phi (0.0); - - IteratorType it = functionSpace_.begin(); - IteratorType endit = functionSpace_.end(); - - assert( it != endit ); - - // Get quadrature rule - const QuadratureRule<double, dim>& quad = QuadratureRules<double, dim>::rule(it->geometry().type(), polOrd); - - for( ; it != endit ; ++it) - { - LocalFuncType lf = discFunc.localFunction( *it ); - - const typename FunctionSpaceType::BaseFunctionSetType & set = - functionSpace_.getBaseFunctionSet(*it); - - for(int i=0; i<lf.numberOfDofs(); i++) - { - for(unsigned int qP = 0; qP < quad.size(); qP++) - { - double det = (*it).geometry().integrationElement(quad[qP].position()); - f.evaluate((*it).geometry().global( quad[qP].position() ), ret); - set.eval(i,quad[qP].position(),phi); - lf[i] += det * quad[qP].weight() * (ret * phi); - } - } - } - } - }; - -} // end namespace - -#endif diff --git a/fem/norms/Makefile.am b/fem/norms/Makefile.am deleted file mode 100644 index a4b7fb9a3068ee09a6cd97f5a5d5af70d06665e7..0000000000000000000000000000000000000000 --- a/fem/norms/Makefile.am +++ /dev/null @@ -1,6 +0,0 @@ -# $Id$ - -femnormsdir = $(includedir)/dune/fem/norms -femnorms_HEADERS = l2norm.hh norm.hh - -include $(top_srcdir)/am/global-rules diff --git a/solver/common/.gitignore b/solver/common/.gitignore deleted file mode 100644 index 55880268f2ad9f36919c5800c511bfd8ee056caa..0000000000000000000000000000000000000000 --- a/solver/common/.gitignore +++ /dev/null @@ -1,3 +0,0 @@ -Makefile -Makefile.in -semantic.cache \ No newline at end of file