799 {
801
802 const size_t nb_integration_pts = AssemblyDomainEleOp::getGaussPts().size2();
803 const size_t nb_base_functions = data.getN().size2();
804
806
807 auto next = [&]() { ++t_res_c; };
808
809 auto t_w = AssemblyDomainEleOp::getFTensor0IntegrationWeight();
810 auto &nf = AssemblyDomainEleOp::locF;
811 auto t_base = data.getFTensor0N();
812 for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
813 const double alpha = AssemblyDomainEleOp::getMeasure() * t_w;
814 ++t_w;
815 const auto res = alpha * t_res_c;
816 next();
817
818 size_t bb = 0;
819 for (; bb != AssemblyDomainEleOp::nbRows; ++bb) {
820 nf[bb] += t_base * res;
821 ++t_base;
822 }
823 for (; bb < nb_base_functions; ++bb)
824 ++t_base;
825 }
826
828}
#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()
