1028 {
1030
1031#ifndef NDEBUG
1034 "Topological data pointer is null");
1035#endif
1036
1037 const int nb_dofs = data.
getIndices().size();
1038 if (!nb_dofs)
1040
1041 const int nb_integration_pts = OP::getGaussPts().size2();
1042
1043 auto t_normal = OP::getFTensor1NormalsAtGaussPts();
1044 auto t_w = OP::getFTensor0IntegrationWeight();
1045 const int nb_base_functions = data.
getN().size2() /
SPACE_DIM;
1047 auto t_obj_dtraction =
1048 topoData->getFTensorObjDDisplacement(nb_integration_pts);
1049
1052
1053 for (int gg = 0; gg != nb_integration_pts; ++gg) {
1054 auto t_nf = getFTensor1FromPtr<SPACE_DIM>(&*OP::locF.begin());
1055 int bb = 0;
1056 for (; bb != nb_dofs /
SPACE_DIM; ++bb) {
1058 t_w * (t_row_base_fun(
j) * t_normal(
j)) * t_obj_dtraction(
i) * 0.5;
1059 ++t_nf;
1060 ++t_row_base_fun;
1061 }
1062 for (; bb != nb_base_functions; ++bb)
1063 ++t_row_base_fun;
1064
1065 ++t_w;
1066 ++t_normal;
1067 ++t_obj_dtraction;
1068 }
1069
1071}
#define FTENSOR_INDEX(DIM, I)
#define MoFEMFunctionReturnHot(a)
Last executable line of each PETSc function used for error handling. Replaces return()
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
@ MOFEM_DATA_INCONSISTENCY
#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
boost::shared_ptr< TopologicalData > topoData
MatrixDouble & getN(const FieldApproximationBase base)
get base functions this return matrix (nb. of rows is equal to nb. of Gauss pts, nb....
auto getFTensor1N(FieldApproximationBase base)
Get base functions for Hdiv/Hcurl spaces.
const VectorInt & getIndices() const
Get global indices of degrees of freedom on entity.