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