Concepts.h

Go to the documentation of this file.

 
#ifndef PBAT_MATH_CONCEPTS_H
#define PBAT_MATH_CONCEPTS_H
 
#include <concepts>
#include <pbat/Aliases.h>
#include <pbat/common/Concepts.h>
 
namespace pbat {
namespace math {

template <class Q>
concept CQuadratureRule = requires(Q q)
{
    requires std::integral<decltype(Q::kDims)>;
    requires common::CContiguousArithmeticRange<decltype(q.points)>;
    requires common::CContiguousArithmeticRange<decltype(q.weights)>;
    {
        q.points.size()
    } -> std::convertible_to<int>;
    {
        q.weights.size()
    } -> std::convertible_to<int>;
};

template <class Q>
concept CFixedPointQuadratureRule = requires(Q q)
{
    requires CQuadratureRule<Q>;
    requires std::is_integral_v<decltype(Q::kPoints)>;
    {q.points.size() / q.weights.size() == Q::kDims};
    {q.weights.size() == Q::kPoints};
};

template <class Q>
concept CPolynomialQuadratureRule = requires(Q q)
{
    requires CQuadratureRule<Q>;
    requires std::is_integral_v<decltype(Q::kOrder)>;
};

template <class Q>
concept CFixedPointPolynomialQuadratureRule =
    CFixedPointQuadratureRule<Q> and CPolynomialQuadratureRule<Q>;
 
} // namespace math
} // namespace pbat
 
#endif // PBAT_MATH_CONCEPTS_H