IntegrationGaussLegendreRegular.h

Go to the documentation of this file.

 
#pragma once
 
#include <array>
#include <cmath>
 
#include "BaseLib/Error.h"
#include "MathLib/Integration/GaussLegendre.h"
#include "MathLib/WeightedPoint.h"
 
namespace NumLib
{
template <unsigned N_DIM>
class IntegrationGaussLegendreRegular
{
public:
    explicit IntegrationGaussLegendreRegular(unsigned order = 2) : _order(order)
    {
        this->setIntegrationOrder(order);
    }
 
    void setIntegrationOrder(unsigned order)
    {
        this->_n_sampl_pt = static_cast<unsigned>(std::pow(order, N_DIM));
        this->_order = order;
    }
 
    unsigned getIntegrationOrder() const { return _order; }
 
    unsigned getNumberOfPoints() const { return _n_sampl_pt; }
 
    MathLib::WeightedPoint getWeightedPoint(unsigned const igp) const
    {
        return getWeightedPoint(getIntegrationOrder(), igp);
    }
 
    static std::array<unsigned, N_DIM> getPositionIndices(unsigned order,
                                                          unsigned igp);
 
    static MathLib::WeightedPoint getWeightedPoint(unsigned const order,
                                                   unsigned const igp);
 
private:
    template <typename Method>
    static MathLib::WeightedPoint getWeightedPoint(
        std::array<unsigned, N_DIM> const& pos);
 
private:
    unsigned _order;
    unsigned _n_sampl_pt{0};
};
 
}  // namespace NumLib
 
#include "IntegrationGaussLegendreRegular-impl.h"