915                                    {
  917 
  918  const size_t nb_integration_pts = AssemblyDomainEleOp::getGaussPts().size2();
  919  const size_t nb_base_functions = data.getN().size2();
  920 
  922 
  923  auto next = [&]() { ++t_res_c; };
  924 
  925  auto t_w = AssemblyDomainEleOp::getFTensor0IntegrationWeight();
  926  auto &nf = AssemblyDomainEleOp::locF;
  927  auto t_base = data.getFTensor0N();
  928  for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
  929    const double alpha = AssemblyDomainEleOp::getMeasure() * t_w;
  930    ++t_w;
  931    const auto res = alpha * t_res_c;
  932    next();
  933 
  934    size_t bb = 0;
  935    for (; bb != AssemblyDomainEleOp::nbRows; ++bb) {
  936      nf[bb] += t_base * res;
  937      ++t_base;
  938    }
  939    for (; bb < nb_base_functions; ++bb)
  940      ++t_base;
  941  }
  942 
  944}
#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()