1002                                        {
 1004 
 1009  constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
 
 1012 
 1013  auto &locMat = AssemblyDomainEleOp::locMat;
 1014 
 1015  const auto nb_integration_pts = AssemblyDomainEleOp::getGaussPts().size2();
 1016  const auto nb_row_base_functions = row_data.getN().size2();
 1017 
 1018  auto t_res_flow_dstrain =
 1019      getFTensor4DdgFromMat<DIM, DIM>(
commonDataPtr->resFlowDstrain);
 
 1020  auto t_res_flow_dplastic_strain =
 1021      getFTensor4DdgFromMat<DIM, DIM>(
commonDataPtr->resFlowDstrainDot);
 
 1023 
 1024  auto next = [&]() {
 1025    ++t_res_flow_dstrain;
 1026    ++t_res_flow_dplastic_strain;
 1027  };
 1028 
 1029  auto t_w = AssemblyDomainEleOp::getFTensor0IntegrationWeight();
 1030  auto t_row_base = row_data.getFTensor0N();
 1031  for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
 1032    double alpha = AssemblyDomainEleOp::getMeasure() * t_w;
 1033    ++t_w;
 1034 
 1036    t_res_mat(O, L) =
 1037        alpha * (t_L(
i, 
j, O) * ((t_res_flow_dplastic_strain(
i, 
j, 
k, 
l) -
 
 1038                                  t_res_flow_dstrain(
i, 
j, 
k, 
l)) *
 
 1040    next();
 1041 
 1042    size_t rr = 0;
 1043    for (; rr != AssemblyDomainEleOp::nbRows / 
size_symm; ++rr) {
 
 1046      auto t_col_base = col_data.getFTensor0N(gg, 0);
 1047      for (
size_t cc = 0; cc != AssemblyDomainEleOp::nbCols / 
size_symm; ++cc) {
 
 1048        t_mat(O, L) += ((t_row_base * t_col_base) * t_res_mat(O, L));
 1049        ++t_mat;
 1050        ++t_col_base;
 1051      }
 1052 
 1053      ++t_row_base;
 1054    }
 1055 
 1056    for (; rr < nb_row_base_functions; ++rr)
 1057      ++t_row_base;
 1058  }
 1059 
 1061}
#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
FTensor::Index< 'l', 3 > l
FTensor::Index< 'j', 3 > j
FTensor::Index< 'k', 3 > k
static auto get_mat_tensor_sym_dtensor_sym(size_t rr, MatrixDouble &mat, FTensor::Number< 2 >)
auto symm_L_tensor(FTensor::Number< DIM >)