768 {
770
773 constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
775
776 const auto nb_integration_pts = AssemblyDomainEleOp::getGaussPts().size2();
777 const size_t nb_row_base_functions = row_data.getN().size2();
778 auto &locMat = AssemblyDomainEleOp::locMat;
779
780 auto t_res_flow_dtemp =
782
784
785 auto next = [&]() { ++t_res_flow_dtemp; };
786
787 auto t_w = AssemblyDomainEleOp::getFTensor0IntegrationWeight();
788 auto t_row_base = row_data.getFTensor0N();
789 for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
790 double alpha = AssemblyDomainEleOp::getMeasure() * t_w;
791 ++t_w;
793 t_res_vec(L) = alpha * (t_res_flow_dtemp(
i,
j) * t_L(
i,
j, L));
794 next();
795
796 size_t rr = 0;
797 for (; rr != AssemblyDomainEleOp::nbRows /
size_symm; ++rr) {
800 auto t_col_base = col_data.getFTensor0N(gg, 0);
801 for (size_t cc = 0; cc != AssemblyDomainEleOp::nbCols; cc++) {
802 t_mat(L) += t_row_base * t_col_base * t_res_vec(L);
803 ++t_mat;
804 ++t_col_base;
805 }
806 ++t_row_base;
807 }
808 for (; rr != nb_row_base_functions; ++rr)
809 ++t_row_base;
810 }
811
813}
#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< 'j', 3 > j
static auto get_mat_tensor_sym_dscalar(size_t rr, MatrixDouble &mat, FTensor::Number< 2 >)
auto symm_L_tensor(FTensor::Number< DIM >)