v0.16.0
Loading...
Searching...
No Matches
Public Member Functions | List of all members
EshelbianPlasticity::dJ_dTractionImpl Struct Reference
Inheritance diagram for EshelbianPlasticity::dJ_dTractionImpl:
[legend]
Collaboration diagram for EshelbianPlasticity::dJ_dTractionImpl:
[legend]

Public Member Functions

MoFEMErrorCode iNtegrate (EntData &data)
 
- Public Member Functions inherited from EshelbianPlasticity::OpAssembleBrokenFaceTopologicalDerivativeImplBase
 OpAssembleBrokenFaceTopologicalDerivativeImplBase (boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_base_side_data, boost::shared_ptr< TopologicalData > topo_ptr, boost::shared_ptr< double > J_ptr, SmartPetscObj< Vec > assemble_vec, Tag topo_tag, boost::shared_ptr< Range > ents_ptr=nullptr)
 
MoFEMErrorCode aSsemble (EntData &data) override
 

Additional Inherited Members

- Public Types inherited from EshelbianPlasticity::OpAssembleBrokenFaceTopologicalDerivativeImplBase
using OP = typename FormsIntegrators< FaceUserDataOperator >::Assembly< A >::OpBrokenBase
 
- Protected Attributes inherited from EshelbianPlasticity::OpAssembleBrokenFaceTopologicalDerivativeImplBase
boost::shared_ptr< doubleJPtr
 
SmartPetscObj< Vec > assembleVec
 
Tag topoTag
 
boost::shared_ptr< TopologicalDatatopoData
 

Detailed Description

Definition at line 174 of file EshelbianTopologicalDerivativeOperators.cpp.

Member Function Documentation

◆ iNtegrate()

MoFEMErrorCode EshelbianPlasticity::dJ_dTractionImpl::iNtegrate ( EntData data)

Definition at line 1028 of file EshelbianTopologicalDerivativeOperators.cpp.

1028 {
1030
1031#ifndef NDEBUG
1032 if (!topoData)
1033 SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
1034 "Topological data pointer is null");
1035#endif // NDEBUG
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;
1046 auto t_row_base_fun = data.getFTensor1N<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) {
1057 t_nf(i) +=
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)
constexpr int SPACE_DIM
#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
Definition definitions.h:31
#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
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.

The documentation for this struct was generated from the following file: