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 1024 of file EshelbianTopologicalDerivativeOperators.cpp.

1024 {
1026
1027#ifndef NDEBUG
1028 if (!topoData)
1029 SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
1030 "Topological data pointer is null");
1031#endif // NDEBUG
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;
1042 auto t_row_base_fun = data.getFTensor1N<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) {
1053 t_nf(i) +=
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)
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: