1176 {
1178
1179 const auto nb_integration_pts = AssemblyDomainEleOp::getGaussPts().size2();
1180 const auto nb_row_base_functions = row_data.getN().size2();
1181
1182 auto t_res_c_dtau = getFTensor0FromVec(
commonDataPtr->resCdTau);
1183 auto next = [&]() { ++t_res_c_dtau; };
1184
1185 auto t_w = AssemblyDomainEleOp::getFTensor0IntegrationWeight();
1186 auto t_row_base = row_data.getFTensor0N();
1187 for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
1188 const double alpha = AssemblyDomainEleOp::getMeasure() * t_w;
1189 ++t_w;
1190
1191 const auto res = alpha * (t_res_c_dtau);
1192 next();
1193
1194 auto mat_ptr = AssemblyDomainEleOp::locMat.data().begin();
1195 size_t rr = 0;
1196 for (; rr != AssemblyDomainEleOp::nbRows; ++rr) {
1197 auto t_col_base = col_data.getFTensor0N(gg, 0);
1198 for (size_t cc = 0; cc != AssemblyDomainEleOp::nbCols; ++cc) {
1199 *mat_ptr += t_row_base * t_col_base * res;
1200 ++t_col_base;
1201 ++mat_ptr;
1202 }
1203 ++t_row_base;
1204 }
1205 for (; rr < nb_row_base_functions; ++rr)
1206 ++t_row_base;
1207 }
1208
1210}
#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()
