693 {
695
696 const auto nb_integration_pts = AssemblyDomainEleOp::getGaussPts().size2();
697 const auto nb_row_base_functions = row_data.getN().size2();
698
699 auto t_res_c_dtemperature =
701 auto next = [&]() { ++t_res_c_dtemperature; };
702
703 auto t_w = AssemblyDomainEleOp::getFTensor0IntegrationWeight();
704 auto t_row_base = row_data.getFTensor0N();
705 for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
706 const double alpha = AssemblyDomainEleOp::getMeasure() * t_w;
707 ++t_w;
708
709 const auto res = alpha * (t_res_c_dtemperature);
710 next();
711
712 auto mat_ptr = AssemblyDomainEleOp::locMat.data().begin();
713 size_t rr = 0;
714 for (; rr != AssemblyDomainEleOp::nbRows; ++rr) {
715 auto t_col_base = col_data.getFTensor0N(gg, 0);
716 for (size_t cc = 0; cc != AssemblyDomainEleOp::nbCols; ++cc) {
717 *mat_ptr += t_row_base * t_col_base * res;
718 ++t_col_base;
719 ++mat_ptr;
720 }
721 ++t_row_base;
722 }
723 for (; rr < nb_row_base_functions; ++rr)
724 ++t_row_base;
725 }
726
728}
#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()