Concepts.h

Go to the documentation of this file.

 
#ifndef PBAT_FEM_CONCEPTS_H
#define PBAT_FEM_CONCEPTS_H
 
#include "pbat/Aliases.h"
#include "pbat/common/Concepts.h"
 
#include <concepts>
 
namespace pbat::fem {

template <class T>
concept CElement = requires(T t)
{
    typename T::AffineBaseType;
    {T::bHasConstantJacobian}->std::convertible_to<bool>;
    // Should be valid for argument > 1 as well, but we don't check that.
    typename T::template QuadratureType<1, Scalar>;
    {T::kOrder}->std::convertible_to<int>;
    {T::kDims}->std::convertible_to<int>;
    {T::kNodes}->std::convertible_to<int>;
    requires common::CContiguousIndexRange<decltype(T::Coordinates)>;
    requires common::CContiguousIndexRange<decltype(T::Vertices)>;
    {t.N(Vector<T::kDims>{})}->std::convertible_to<Vector<T::kNodes>>;
    {t.GradN(Vector<T::kDims>{})}->std::convertible_to<Matrix<T::kNodes, T::kDims>>;
};

template <class M>
concept CMesh = requires(M m)
{
    requires CElement<typename M::ElementType>;
    requires common::CFloatingPoint<typename M::ScalarType>;
    requires common::CIndex<typename M::IndexType>;
    {M::kDims}->std::convertible_to<int>;
    {M::kOrder}->std::convertible_to<int>;
    {m.X} /*->std::
         convertible_to<Eigen::Matrix<typename M::ScalarType, Eigen::Dynamic, Eigen::Dynamic>>*/
    ;
    {m.E} /*->std::convertible_to<Eigen::Matrix<typename M::IndexType, Eigen::Dynamic,
             Eigen::Dynamic>>*/
    ;
};
 
} // namespace pbat::fem
 
#endif // PBAT_FEM_CONCEPTS_H