EshelbianPlasticity::dJ_duGammaImpl Struct Reference

Inheritance diagram for EshelbianPlasticity::dJ_duGammaImpl:

Collaboration diagram for EshelbianPlasticity::dJ_duGammaImpl:

Public Member Functions
MoFEMErrorCode	integrate (EntData &data)
Public Member Functions inherited from EshelbianPlasticity::OpAssembleTopologicalObjectiveDerivativeImplBase< OpAssembleFace >
	OpAssembleTopologicalObjectiveDerivativeImplBase (const std::string &field_name, boost::shared_ptr< DataAtIntegrationPts > data_ptr, boost::shared_ptr< TopologicalData > topo_ptr, boost::shared_ptr< double > J_ptr, SmartPetscObj< Vec > assemble_vec, Tag topo_tag)
MoFEMErrorCode	assemble (EntData &data) override
Public Member Functions inherited from OpAssembleBasic< FaceUserDataOperator >
	OpAssembleBasic (const std::string &field_name, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const char type)
	OpAssembleBasic (std::string row_field, std::string col_field, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const char type, const bool assemble_symmetry, ScaleOff scale_off=[]() { return 1;})
	OpAssembleBasic (const FieldSpace space)
virtual MoFEMErrorCode	integrate (int row_side, EntityType row_type, EntData &data)
virtual MoFEMErrorCode	integrate (EntData &row_data, EntData &col_data)
virtual MoFEMErrorCode	assemble (int row_side, EntityType row_type, EntData &data)
virtual MoFEMErrorCode	assemble (int row_side, int col_side, EntityType row_type, EntityType col_type, EntData &row_data, EntData &col_data)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)
MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntData &row_data, EntData &col_data)

Additional Inherited Members
Public Types inherited from EshelbianPlasticity::OpAssembleTopologicalObjectiveDerivativeImplBase< OpAssembleFace >
using	OP = OpAssembleFace
Public Types inherited from OpAssembleBasic< FaceUserDataOperator >
using	ScaleOff = boost::function< double()>
Public Attributes inherited from OpAssembleBasic< FaceUserDataOperator >
const bool	assembleSymmetry
boost::shared_ptr< DataAtIntegrationPts >	dataAtPts
	data at integration pts
VectorDouble	nF
	local right hand side vector
MatrixDouble	K
	local tangent matrix
MatrixDouble	transposeK
ScaleOff	scaleOff
Protected Attributes inherited from EshelbianPlasticity::OpAssembleTopologicalObjectiveDerivativeImplBase< OpAssembleFace >
double	locJ
boost::shared_ptr< double >	JPtr
SmartPetscObj< Vec >	assembleVec
Tag	topoTag
boost::shared_ptr< TopologicalData >	topoData

Detailed Description

Definition at line 170 of file EshelbianTopologicalDerivativeOperators.cpp.

Member Function Documentation

integrate()

MoFEMErrorCode EshelbianPlasticity::dJ_duGammaImpl::integrate ( EntData &  data )

virtual

Reimplemented from OpAssembleBasic< FaceUserDataOperator >.

Definition at line 589 of file EshelbianTopologicalDerivativeOperators.cpp.

                                                      {
  MoFEMFunctionBegin;
 
#ifndef NDEBUG
  if (!topoData)
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "Topological data pointer is null");
#endif // NDEBUG
 
  const int nb_dofs = data.getIndices().size();
  if (!nb_dofs)
    MoFEMFunctionReturnHot(0);
 
  const int nb_integration_pts = data.getN().size1();
 
#ifndef NDEBUG
  const int disp_rows = topoData->objDDisplacementAtPts.size1();
  const int disp_cols = topoData->objDDisplacementAtPts.size2();
  if (disp_cols != nb_integration_pts)
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "objDDisplacementAtPts columns (%d) != nb integration points (%d)",
            disp_cols, nb_integration_pts);
 
  if (disp_rows != SPACE_DIM) {
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "objDDisplacementAtPts row size is %d, expected %d", disp_rows,
            SPACE_DIM);
  }
#endif // NDEBUG
 
  auto t_w = getFTensor0IntegrationWeight();
  const int nb_base_functions = data.getN().size2();
  auto t_row_base_fun = data.getFTensor0N();
  auto t_obj_du_gamma =
      getFTensor1FromMat<SPACE_DIM>(topoData->objDDisplacementAtPts);
 
  FTENSOR_INDEX(SPACE_DIM, i);
 
  auto get_ftensor1 = [](auto &v) {
    return FTensor::Tensor1<FTensor::PackPtr<double *, SPACE_DIM>, SPACE_DIM>(
        &v[0], &v[1], &v[2]);
  };
 
  for (int gg = 0; gg != nb_integration_pts; ++gg) {
    const double a = t_w * getMeasure();
    auto t_nf = get_ftensor1(nF);
 
    int bb = 0;
    for (; bb != nb_dofs / SPACE_DIM; ++bb) {
      t_nf(i) += a * t_row_base_fun * t_obj_du_gamma(i);
      ++t_nf;
      ++t_row_base_fun;
    }
    for (; bb != nb_base_functions; ++bb)
      ++t_row_base_fun;
 
    ++t_w;
    ++t_obj_du_gamma;
  }
 
  MoFEMFunctionReturn(0);
}

---

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

• users_modules/eshelbian_plasticity/src/impl/EshelbianTopologicalDerivativeOperators.cpp