12#include <range/v3/range/conversion.hpp>
13#include <range/v3/view/transform.hpp>
20template <
typename ShapeFunction>
22 boost::mp11::mp_identity<ShapeFunction>)
26 using MeshElement =
typename ShapeFunction::MeshElement;
28 auto const& integration_method =
29 NumLib::IntegrationMethodRegistry::getIntegrationMethod<MeshElement>(
32 using ShapeMatrixPolicy =
40 constexpr int dim = 3;
50 ranges::views::transform([](
auto& sm) {
return std::move(sm.N); }) |
51 ranges::to<std::vector>;
60 using Indices = mp_iota<mp_size<ETLs>>;
63 [
this, integration_order](
auto index)
65 using ETL = mp_at_c<ETLs, index>;
70 integration_order, mp_identity<ShapeFunctionHigherOrder>{});
76 integration_order, mp_identity<ShapeFunctionLowerOrder>{});
MeshElementType const element
boost::mp11::mp_at< ShapeFunctionsLowerOrder, boost::mp11::mp_find< MeshElements, MeshElement > > ShapeFunctionLowerOrder
ShapeMatrixCache(unsigned const integration_order)
boost::mp11::mp_at< ShapeFunctionsHigherOrder, boost::mp11::mp_find< MeshElements, MeshElement > > ShapeFunctionHigherOrder
ShapeMatrixVectorsLowerOrder_N Nss_lower_order_
ShapeMatrixVectorsHigherOrder_N Nss_higher_order_
typename ElementTraitsLagrange::ShapeFunction GetShapeFunctionHigherOrder
typename ElementTraitsLagrange::LowerOrderShapeFunction GetShapeFunctionLowerOrder
std::vector< typename ShapeMatricesType::ShapeMatrices, Eigen::aligned_allocator< typename ShapeMatricesType::ShapeMatrices > > initShapeMatrices(MeshLib::Element const &e, bool const is_axially_symmetric, IntegrationMethod const &integration_method)
static auto initShapeMatrices(unsigned const integration_order, boost::mp11::mp_identity< ShapeFunction >)