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>
 
#include "types.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<bool conforming>
      using IntersectionQuadrature = Fem::IntersectionQuadrature<FaceQuadratureType, conforming>;
 
      template<bool conforming>
      using IntersectionMassQuadrature = Fem::IntersectionQuadrature<FaceMassQuadratureType, conforming>;
 
      static const bool hasMassQuadrature = false;
 
      static const bool hasFaceMassQuadrature = false;
    };
 
    template<class GridPart, bool bulkLumping, bool faceLumping>
    struct MaybeLumpingQuadratureTraits
      : public DefaultQuadratureTraits<GridPart>
    {
     protected:
      typedef DefaultQuadratureTraits<GridPart> BaseType;
     public:
      typedef typename BaseType::GridPartType GridPartType;
      typedef typename BaseType::BulkQuadratureType BulkQuadratureType;
      typedef typename BaseType::FaceQuadratureType FaceQuadratureType;
 
      using BaseType::IntersectionQuadrature;
 
      using BulkMassQuadratureType =
        ConditionalType<bulkLumping,
                        Fem::CachingLumpingQuadrature<GridPartType, 0>,
                        typename BaseType::BulkMassQuadratureType>;
 
      using FaceMassQuadratureType =
        ConditionalType<faceLumping,
                        Fem::CachingLumpingQuadrature<GridPartType, 1>,
                        typename BaseType::FaceMassQuadratureType>;
 
      template<bool conforming>
      using IntersectionMassQuadrature = Fem::IntersectionQuadrature<FaceMassQuadratureType, conforming>;
 
      static const bool hasMassQuadrature = bulkLumping;
 
      static const bool hasFaceMassQuadrature = faceLumping;
    };
 
    template<class GridPart>
    using LumpingQuadratureTraits = MaybeLumpingQuadratureTraits<GridPart, true, true>;
 
    template<class GridPart>
    using BulkLumpingQuadratureTraits = MaybeLumpingQuadratureTraits<GridPart, true, false>;
 
    template<class GridPart>
    using FaceLumpingQuadratureTraits = MaybeLumpingQuadratureTraits<GridPart, false, true>;
 
    template<class GridPart>
    using QuadratureTraits = DefaultQuadratureTraits<GridPart>;
 
 
  } // namespace ACFem
 
} // namespace Dune
 
#endif // __DUNE_ACFEM_QUADRATURE_HH__