compositequadraturerule.hh

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright © DUNE-FUFEM Project contributors, see file AUTHORS.md
// SPDX-License-Identifier: LicenseRef-GPL-2.0-only-with-DUNE-exception OR LGPL-3.0-or-later
 
#ifndef COMPOSITE_QUADRATURE_RULE_HH
#define COMPOSITE_QUADRATURE_RULE_HH
 
#include <dune/geometry/quadraturerules.hh>
#include <dune/geometry/virtualrefinement.hh>
 
 
template <class ct, int dim>
class CompositeQuadratureRule:
    public Dune::QuadratureRule<ct,dim>
{
    public:
        CompositeQuadratureRule(const Dune::QuadratureRule<ct,dim>& quad, int refinement)
        {
            
            typedef Dune::StaticRefinement<Dune::GeometryTypes::simplex(dim).id(),
                                           ct,
                                           Dune::GeometryTypes::simplex(dim).id(),
                                           dim> Refinement;
 
            int numberOfSubelements = (1<<(dim*refinement));
            //ct volume = Dune::ReferenceElements<ct,dim>::general(quad.type()).volume()/numberOfSubelements;
            ct volumeFraction = ct(1.0)/numberOfSubelements;
            
            auto eSubEnd = Refinement::eEnd(Dune::refinementLevels(refinement));
            auto eSubIt  = Refinement::eBegin(Dune::refinementLevels(refinement));
            
            for (; eSubIt != eSubEnd; ++eSubIt) {
                
                for (size_t i=0; i<quad.size(); i++) {
                    
                    const Dune::FieldVector<ct,dim>& quadPos = quad[i].position();
             
                    this->push_back(Dune::QuadraturePoint<ct,dim>(eSubIt.geometry().global(quadPos), volumeFraction*quad[i].weight()));
                    
                }
                
            }
                
        }
 
};
 
#endif