AGL_Quadrature.h

Go to the documentation of this file.

#pragma once

#include "AdaptiveQuadrature.h"
#include "GaussLegendreRules.h"

#include <cassert>

namespace UA_CoMP {
namespace Num {

template <int N>
class GaussLegendreSelect
{
};

template <>
class GaussLegendreSelect<2>
{
public:
        using Rule = GaussLegendre2Rule;
};

template <>
class GaussLegendreSelect<4>
{
public:
        using Rule = GaussLegendre2Rule;
};

template <>
class GaussLegendreSelect<8>
{
public:
        using Rule = GaussLegendre8Rule;
};

template <>
class GaussLegendreSelect<16>
{
public:
        using Rule = GaussLegendre16Rule;
};

template <>
class GaussLegendreSelect<32>
{
public:
        using Rule = GaussLegendre32Rule;
};

template <int N1, int N2, typename ConvergencePolicy = ComboDifference, typename CellCountPolicy = OnViolationThrow>
class AGL_Quadrature
    : public AdaptiveQuadrature<typename GaussLegendreSelect<N1>::Rule, typename GaussLegendreSelect<N2>::Rule,
                                ConvergencePolicy, CellCountPolicy>
{
public:
        AGL_Quadrature(double convergenceTolerance, unsigned int initialCellCount, unsigned int cellCountLimit)
            : AdaptiveQuadrature<typename GaussLegendreSelect<N1>::Rule, typename GaussLegendreSelect<N2>::Rule,
                                 ConvergencePolicy, CellCountPolicy>(convergenceTolerance, initialCellCount,
                                                                     cellCountLimit)
        {
                assert(N1 != N2);
        }

private:
};

} // namespace Num
} // namespace UA_CoMP