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()