1024 {
1026
1027#ifndef NDEBUG
1030 "Topological data pointer is null");
1031#endif
1032
1033 const int nb_dofs = data.
getIndices().size();
1034 if (!nb_dofs)
1036
1037 const int nb_integration_pts = OP::getGaussPts().size2();
1038
1039 auto t_normal = OP::getFTensor1NormalsAtGaussPts();
1040 auto t_w = OP::getFTensor0IntegrationWeight();
1041 const int nb_base_functions = data.
getN().size2() /
SPACE_DIM;
1043 auto t_obj_dtraction =
1044 topoData->getFTensorObjDDisplacement(nb_integration_pts);
1045
1048
1049 for (int gg = 0; gg != nb_integration_pts; ++gg) {
1050 auto t_nf = getFTensor1FromPtr<SPACE_DIM>(&*OP::locF.begin());
1051 int bb = 0;
1052 for (; bb != nb_dofs /
SPACE_DIM; ++bb) {
1054 t_w * (t_row_base_fun(
j) * t_normal(
j)) * t_obj_dtraction(
i) * 0.5;
1055 ++t_nf;
1056 ++t_row_base_fun;
1057 }
1058 for (; bb != nb_base_functions; ++bb)
1059 ++t_row_base_fun;
1060
1061 ++t_w;
1062 ++t_normal;
1063 ++t_obj_dtraction;
1064 }
1065
1067}
#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.