quadrature.hh

#ifndef __DUNE_ACFEM_QUADRATURE_HH__
#define __DUNE_ACFEM_QUADRATURE_HH__
 
#include <dune/fem/quadrature/cachingquadrature.hh>
#include <dune/fem/quadrature/intersectionquadrature.hh>
#include <dune/fem/quadrature/lumpingquadrature.hh>
 
namespace Dune {
 
  namespace ACFem {
 
    template<class GridPart>
    struct DefaultQuadratureTraits
    {
      typedef GridPart GridPartType;
      typedef Fem::CachingQuadrature<GridPartType, 0> BulkQuadratureType;
      typedef Fem::CachingQuadrature<GridPartType, 1> FaceQuadratureType;
 
      typedef BulkQuadratureType BulkMassQuadratureType;
      typedef FaceQuadratureType FaceMassQuadratureType;
    };
 
    template<class GridPart, bool faceLumping = false>
    struct LumpingQuadratureTraits
      : DefaultQuadratureTraits<GridPart>
    {
     protected:
      typedef DefaultQuadratureTraits<GridPart> BaseType;
     public:
      typedef typename BaseType::GridPartType GridPartType;
      typedef typename BaseType::BulkQuadratureType BulkQuadratureType;
      typedef typename BaseType::FaceQuadratureType FaceQuadratureType;
 
      typedef Fem::CachingLumpingQuadrature<GridPartType, 0> BulkMassQuadratureType;
      // Is it necessary to lump the face-quadrature to get a discrete
      // comparison principle? I guess not ...
      typedef typename std::conditional<faceLumping,
                                        Fem::CachingLumpingQuadrature<GridPartType, 1>,
                                        typename BaseType::FaceMassQuadratureType>::type
      FaceMassQuadratureType;
    };
 
 
  } // namespace ACFem
 
} // namespace Dune
 
#endif // __DUNE_ACFEM_QUADRATURE_HH__