MoFEM::GetFTensor4FromMatImpl< Tensor_Dim0, Tensor_Dim1, Tensor_Dim2, Tensor_Dim3, S, DL, M > Struct Template Reference

#include "include/Templates.hpp"

Static Public Member Functions
static auto	get (M &data, int rr=0, int cc=0)

Detailed Description

template<int Tensor_Dim0, int Tensor_Dim1, int Tensor_Dim2, int Tensor_Dim3, int S, class DL, class M>
struct MoFEM::GetFTensor4FromMatImpl< Tensor_Dim0, Tensor_Dim1, Tensor_Dim2, Tensor_Dim3, S, DL, M >

Definition at line 1263 of file Templates.hpp.

Member Function Documentation

get()

template<int Tensor_Dim0, int Tensor_Dim1, int Tensor_Dim2, int Tensor_Dim3, int S, class DL , class M >

static auto MoFEM::GetFTensor4FromMatImpl< Tensor_Dim0, Tensor_Dim1, Tensor_Dim2, Tensor_Dim3, S, DL, M >::get ( M &  data, int  rr = 0, int  cc = 0  )

inlinestatic

Definition at line 1264 of file Templates.hpp.

                                                          {
    using T = std::remove_cv_t<std::remove_reference_t<decltype(data(0, 0))>>;
    constexpr int tensor_size =
        Tensor_Dim0 * Tensor_Dim1 * Tensor_Dim2 * Tensor_Dim3;
    constexpr int stride =
        S == -1 ? (DL::isGaussByCoeffs ? tensor_size : 1) : S;
    if constexpr (DL::isGaussByCoeffs) {
      static_assert(
          stride % tensor_size == 0,
          "getFTensor4FromMat: stride should be a multiple of tensor_size");
 
#ifndef NDEBUG
      if (data.size2() != tensor_size) {
        THROW_MESSAGE("getFTensor4FromMat<" +
                      boost::lexical_cast<std::string>(Tensor_Dim0) + "," +
                      boost::lexical_cast<std::string>(Tensor_Dim1) + "," +
                      boost::lexical_cast<std::string>(Tensor_Dim2) + "," +
                      boost::lexical_cast<std::string>(Tensor_Dim3) +
                      ">: wrong size of columns in data matrix, should be " +
                      boost::lexical_cast<std::string>(tensor_size) +
                      " but is " +
                      boost::lexical_cast<std::string>(data.size2()));
      }
      auto *first = &data(rr + 0, cc + 0);
      auto *last = &data(rr + 0, cc + 0 + tensor_size - 1);
      if (last - first != tensor_size - 1) {
        THROW_MESSAGE("getFTensor4FromMat<" +
                      boost::lexical_cast<std::string>(Tensor_Dim0) + "," +
                      boost::lexical_cast<std::string>(Tensor_Dim1) + "," +
                      boost::lexical_cast<std::string>(Tensor_Dim2) + "," +
                      boost::lexical_cast<std::string>(Tensor_Dim3) +
                      ">: row slice is not contiguous");
      }
#endif
      return FTensor::Tensor4<FTensor::CursorPtr<T *, stride>, Tensor_Dim0,
                              Tensor_Dim1, Tensor_Dim2, Tensor_Dim3>{
          &data(rr + 0, cc + 0)};
    } else if constexpr (DL::isCoeffsByGauss) {
      std::array<T *, tensor_size> ptrs;
#ifndef NDEBUG
      if (data.size1() != tensor_size) {
        THROW_MESSAGE("getFTensor4FromMat<" +
                      boost::lexical_cast<std::string>(Tensor_Dim0) + "," +
                      boost::lexical_cast<std::string>(Tensor_Dim1) + "," +
                      boost::lexical_cast<std::string>(Tensor_Dim2) + "," +
                      boost::lexical_cast<std::string>(Tensor_Dim3) +
                      ">: wrong size of rows in data matrix, should be " +
                      boost::lexical_cast<std::string>(tensor_size) +
                      " but is " +
                      boost::lexical_cast<std::string>(data.size1()));
      }
#endif
      for (auto i = 0; i != tensor_size; ++i)
        ptrs[i] = &data(rr + i, cc + 0);
      return makeFTensor4FromPtrArray<T, Tensor_Dim0, Tensor_Dim1, Tensor_Dim2,
                                      Tensor_Dim3, stride>(
          ptrs, std::make_index_sequence<tensor_size>{});
    } else {
      static_assert(!std::is_same<M, M>::value,
                    "Unsupported data layout for getFTensor4FromMatImpl");
    }
  }

---

The documentation for this struct was generated from the following file:

• include/Templates.hpp