lumpingquadraturerule.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 LUMPINGQUADRATURERULE_HH
#define LUMPINGQUADRATURERULE_HH
 
#include <dune/geometry/referenceelements.hh>
#include <dune/geometry/quadraturerules.hh>
 
template<class ct, int dim>
class LumpingQuadratureRule :
    public Dune::QuadratureRule<ct, dim>
{
        typedef typename Dune::QuadratureRule<ct, dim> Base;
        typedef typename Dune::QuadraturePoint<ct, dim> QuadPoint;
 
    public:
 
        LumpingQuadratureRule(const Dune::GeometryType& gt) :
            Base(gt, 1)
        {
            auto refElement = Dune::ReferenceElements<ct, dim>::general(this->type());
            int size = refElement.size(dim);
            ct weight = refElement.volume()/size;
            this->reserve(size);
            for(int i=0; i<size; ++i)
                this->push_back(QuadPoint(refElement.position(i,dim), weight));
        }
 
};
 
 
#endif