2955 {
2958 const int nb_integration_pts = getGaussPtsSlave().size2();
2960 xi_grad_mat.resize(6, nb_integration_pts, false);
2961 if (type == MBVERTEX)
2962 xi_grad_mat.clear();
2963
2966
2967 if (nb_dofs) {
2968
2969 auto t_grad_pos_xi = getFTensor2FromMat<3, 2>(xi_grad_mat);
2970
2971 for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
2974 for (size_t bb = 0; bb != nb_dofs / 3; ++bb) {
2975 t_grad_pos_xi(
i,
I) += t_diff_base(
I) * t_data(
i);
2976 ++t_data;
2977 ++t_diff_base;
2978 }
2979 ++t_grad_pos_xi;
2980 }
2981 }
2982
2984}
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
FTensor::Index< 'i', SPACE_DIM > i
constexpr IntegrationType I
FTensor::Tensor1< FTensor::PackPtr< double *, Tensor_Dim >, Tensor_Dim > getFTensor1DiffN(const FieldApproximationBase base)
Get derivatives of base functions.
const VectorDouble & getFieldData() const
get dofs values
FTensor::Tensor1< FTensor::PackPtr< double *, Tensor_Dim >, Tensor_Dim > getFTensor1FieldData()
Return FTensor of rank 1, i.e. vector from filed data coeffinects.