1292                                        {
 1294 
 1295  const auto nb_integration_pts = AssemblyDomainEleOp::getGaussPts().size2();
 1296  const auto nb_row_base_functions = row_data.getN().size2();
 1297 
 1298  auto t_res_c_dtau = getFTensor0FromVec(
commonDataPtr->resCdTau);
 
 1299  auto next = [&]() { ++t_res_c_dtau; };
 1300 
 1301  auto t_w = AssemblyDomainEleOp::getFTensor0IntegrationWeight();
 1302  auto t_row_base = row_data.getFTensor0N();
 1303  for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
 1304    const double alpha = AssemblyDomainEleOp::getMeasure() * t_w;
 1305    ++t_w;
 1306 
 1307    const auto res = alpha * (t_res_c_dtau);
 1308    next();
 1309 
 1310    auto mat_ptr = AssemblyDomainEleOp::locMat.data().begin();
 1311    size_t rr = 0;
 1312    for (; rr != AssemblyDomainEleOp::nbRows; ++rr) {
 1313      auto t_col_base = col_data.getFTensor0N(gg, 0);
 1314      for (size_t cc = 0; cc != AssemblyDomainEleOp::nbCols; ++cc) {
 1315        *mat_ptr += t_row_base * t_col_base * res;
 1316        ++t_col_base;
 1317        ++mat_ptr;
 1318      }
 1319      ++t_row_base;
 1320    }
 1321    for (; rr < nb_row_base_functions; ++rr)
 1322      ++t_row_base;
 1323  }
 1324 
 1326}
#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()