1122                                        {
 1124 
 1127  constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
 
 1129 
 1130  const auto nb_integration_pts = AssemblyDomainEleOp::getGaussPts().size2();
 1131  const size_t nb_row_base_functions = row_data.getN().size2();
 1132  auto &locMat = AssemblyDomainEleOp::locMat;
 1133 
 1134  auto t_res_flow_dtau =
 1135      getFTensor2SymmetricFromMat<DIM>(
commonDataPtr->resFlowDtau);
 
 1136 
 1138 
 1139  auto next = [&]() { ++t_res_flow_dtau; };
 1140 
 1141  auto t_w = AssemblyDomainEleOp::getFTensor0IntegrationWeight();
 1142  auto t_row_base = row_data.getFTensor0N();
 1143  for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
 1144    double alpha = AssemblyDomainEleOp::getMeasure() * t_w;
 1145    ++t_w;
 1147    t_res_vec(L) = alpha * (t_res_flow_dtau(
i, 
j) * t_L(
i, 
j, L));
 
 1148    next();
 1149 
 1150    size_t rr = 0;
 1151    for (; rr != AssemblyDomainEleOp::nbRows / 
size_symm; ++rr) {
 
 1152      auto t_mat =
 1154      auto t_col_base = col_data.getFTensor0N(gg, 0);
 1155      for (size_t cc = 0; cc != AssemblyDomainEleOp::nbCols; cc++) {
 1156        t_mat(L) += t_row_base * t_col_base * t_res_vec(L);
 1157        ++t_mat;
 1158        ++t_col_base;
 1159      }
 1160      ++t_row_base;
 1161    }
 1162    for (; rr != nb_row_base_functions; ++rr)
 1163      ++t_row_base;
 1164  }
 1165 
 1167}
#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
static auto get_mat_tensor_sym_dscalar(size_t rr, MatrixDouble &mat, FTensor::Number< 2 >)
auto symm_L_tensor(FTensor::Number< DIM >)