748 {
749
750 using T = std::remove_reference_t<decltype(data(0, 0))>;
751 constexpr int symm_size = (Tensor_Dim * (Tensor_Dim + 1)) / 2;
752 constexpr int stride = S == -1 ? symm_size : S;
753 static_assert(stride % symm_size == 0,
754 "getFTensor2SymmetricFromMat: stride should be a multiple of "
755 "symm_size");
756
757#ifndef NDEBUG
758 if (data.size2() != symm_size)
760 "getFTensor2SymmetricFromMat<" +
761 boost::lexical_cast<std::string>(Tensor_Dim) +
762 ">: wrong size of data matrix, number of columns should be " +
763 boost::lexical_cast<std::string>(symm_size) + " but is " +
764 boost::lexical_cast<std::string>(data.size2()));
765 auto *first = &data(rr + 0, cc + 0);
766 auto *last = &data(rr + 0, cc + 0 + symm_size - 1);
767 if (last - first != symm_size - 1)
769 boost::lexical_cast<std::string>(Tensor_Dim) +
770 ">: row slice is not contiguous");
771#endif
773 Tensor_Dim>(&data(rr + 0, cc + 0));
774 }
#define THROW_MESSAGE(msg)
Throw MoFEM exception.