731 {
733
738 constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
740
741 const auto nb_integration_pts = AssemblyDomainEleOp::getGaussPts().size2();
742 const auto nb_base_functions = data.getN().size2();
743
744 auto t_res_flow = getFTensor2SymmetricFromMat<DIM>(
commonDataPtr->resFlow);
745
746 auto t_L = FTensor::symm_l_tensor<double, DIM>();
747
748 auto next = [&]() { ++t_res_flow; };
749
750 auto t_w = AssemblyDomainEleOp::getFTensor0IntegrationWeight();
751 auto t_base = data.getFTensor0N();
752 auto &nf = AssemblyDomainEleOp::locF;
753 for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
754 double alpha = AssemblyDomainEleOp::getMeasure() * t_w;
755 ++t_w;
756
758 t_rhs(L) = alpha * (t_res_flow(
i,
j) * t_L(
i,
j, L));
759 next();
760
761 auto t_nf = getFTensor1FromArray<size_symm, size_symm>(nf);
762 size_t bb = 0;
763 for (; bb != AssemblyDomainEleOp::nbRows /
size_symm; ++bb) {
764 t_nf(L) += t_base * t_rhs(L);
765 ++t_base;
766 ++t_nf;
767 }
768 for (; bb < nb_base_functions; ++bb)
769 ++t_base;
770 }
771
773}
#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()
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Index< 'l', 3 > l
FTensor::Index< 'j', 3 > j
FTensor::Index< 'k', 3 > k
