MoFEM::GetFTensor2FromMatImpl< Tensor_Dim1, Tensor_Dim2, 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_Dim1, int Tensor_Dim2, int S, class DL, class M>
struct MoFEM::GetFTensor2FromMatImpl< Tensor_Dim1, Tensor_Dim2, S, DL, M >

Definition at line 482 of file Templates.hpp.

Member Function Documentation

get()

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

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

inlinestatic

Definition at line 483 of file Templates.hpp.

                                                          {
    // Preserve const matrix elements in the non-owning tensor view.
    using T = std::remove_reference_t<decltype(data(0, 0))>;
    constexpr int stride =
        S == -1 ? (DL::isGaussByCoeffs ? Tensor_Dim1 * Tensor_Dim2 : 1) : S;
    if constexpr (DL::isGaussByCoeffs) {
      static_assert(stride % (Tensor_Dim1 * Tensor_Dim2) == 0,
                    "getFTensor2FromMat: stride should be a multiple of "
                    "Tensor_Dim1 * Tensor_Dim2");
 
#ifndef NDEBUG
      if (data.size2() != Tensor_Dim1 * Tensor_Dim2) {
        THROW_MESSAGE(
            "getFTensor2FromMat<" +
            boost::lexical_cast<std::string>(Tensor_Dim1) + "," +
            boost::lexical_cast<std::string>(Tensor_Dim2) +
            ">: wrong size of columns in data matrix, should be " +
            boost::lexical_cast<std::string>(Tensor_Dim1 * Tensor_Dim2) +
            " but is " + boost::lexical_cast<std::string>(data.size2()));
      }
      auto *first = &data(rr + 0, cc + 0);
      auto *last = &data(rr + 0, cc + 0 + Tensor_Dim1 * Tensor_Dim2 - 1);
      if (last - first != Tensor_Dim1 * Tensor_Dim2 - 1) {
        THROW_MESSAGE("getFTensor2FromMat<" +
                      boost::lexical_cast<std::string>(Tensor_Dim1) + "," +
                      boost::lexical_cast<std::string>(Tensor_Dim2) +
                      ">: row slice is not contiguous");
      }
#endif
      return FTensor::Tensor2<FTensor::CursorPtr<T *, stride>, Tensor_Dim1,
                              Tensor_Dim2>(&data(rr + 0, cc + 0));
    } else if constexpr (DL::isCoeffsByGauss) {
      std::array<T *, Tensor_Dim1 * Tensor_Dim2> ptrs;
#ifndef NDEBUG
      if (data.size1() != Tensor_Dim1 * Tensor_Dim2) {
        THROW_MESSAGE(
            "getFTensor2FromMat<" +
            boost::lexical_cast<std::string>(Tensor_Dim1) + "," +
            boost::lexical_cast<std::string>(Tensor_Dim2) +
            ">: wrong size of rows in data matrix, should be " +
            boost::lexical_cast<std::string>(Tensor_Dim1 * Tensor_Dim2) +
            " but is " + boost::lexical_cast<std::string>(data.size1()));
      }
#endif
      for (auto i = 0; i != Tensor_Dim1 * Tensor_Dim2; ++i)
        ptrs[i] = &data(rr + i, cc + 0);
      return FTensor::Tensor2<FTensor::PackPtr<T *, stride>, Tensor_Dim1,
                              Tensor_Dim2>(ptrs);
    } else {
      static_assert(!std::is_same<M, M>::value,
                    "Unsupported data layout for getFTensor2FromMatImpl");
    }
  }

---

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

• include/Templates.hpp