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