1210                                        {
 1212 
 1215  constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
 
 1217 
 1218  const auto nb_integration_pts = AssemblyDomainEleOp::getGaussPts().size2();
 1219  const auto nb_row_base_functions = row_data.getN().size2();
 1220 
 1221  auto t_c_dstrain =
 1222      getFTensor2SymmetricFromMat<SPACE_DIM>(
commonDataPtr->resCdStrain);
 
 1223  auto t_c_dplastic_strain =
 1224      getFTensor2SymmetricFromMat<SPACE_DIM>(
commonDataPtr->resCdPlasticStrain);
 
 1225 
 1226  auto next = [&]() {
 1227    ++t_c_dstrain;
 1228    ++t_c_dplastic_strain;
 1229  };
 1230 
 1232 
 1233  auto t_w = AssemblyDomainEleOp::getFTensor0IntegrationWeight();
 1234  auto t_row_base = row_data.getFTensor0N();
 1235  for (auto gg = 0; gg != nb_integration_pts; ++gg) {
 1236    const double alpha = AssemblyDomainEleOp::getMeasure() * t_w;
 1237    ++t_w;
 1238 
 1240    t_res_vec(L) =
 1241        t_L(
i, 
j, L) * (t_c_dplastic_strain(
i, 
j) - t_c_dstrain(
i, 
j));
 
 1242    next();
 1243 
 1246    size_t rr = 0;
 1247    for (; rr != AssemblyDomainEleOp::nbRows; ++rr) {
 1248      const auto row_base = alpha * t_row_base;
 1249      auto t_col_base = col_data.getFTensor0N(gg, 0);
 1250      for (
size_t cc = 0; cc != AssemblyDomainEleOp::nbCols / 
size_symm; cc++) {
 
 1251        t_mat(L) += (row_base * t_col_base) * t_res_vec(L);
 1252        ++t_mat;
 1253        ++t_col_base;
 1254      }
 1255      ++t_row_base;
 1256    }
 1257    for (; rr != nb_row_base_functions; ++rr)
 1258      ++t_row_base;
 1259  }
 1260 
 1262}
#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< 'j', 3 > j
auto get_mat_scalar_dtensor_sym(MatrixDouble &mat, FTensor::Number< 2 >)
auto symm_L_tensor(FTensor::Number< DIM >)