EshelbianPlasticity::OpTopologicalObjectivePythonImpl Struct Reference

Inheritance diagram for EshelbianPlasticity::OpTopologicalObjectivePythonImpl:

Collaboration diagram for EshelbianPlasticity::OpTopologicalObjectivePythonImpl:

Public Types
using	OP = ForcesAndSourcesCore::UserDataOperator

Public Member Functions
	OpTopologicalObjectivePythonImpl (boost::shared_ptr< DataAtIntegrationPts > data_at_pts_ptr, boost::shared_ptr< TopologicalData > topo_p, boost::shared_ptr< ObjectiveFunctionData > python_ptr)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)

Private Attributes
boost::shared_ptr< DataAtIntegrationPts >	dataAtPts
boost::shared_ptr< TopologicalData >	topoData
boost::shared_ptr< ObjectiveFunctionData >	pythonPtr

Detailed Description

Definition at line 202 of file EshelbianTopologicalDerivativeOperators.cpp.

Member Typedef Documentation

OP

using EshelbianPlasticity::OpTopologicalObjectivePythonImpl::OP = ForcesAndSourcesCore::UserDataOperator

Definition at line 204 of file EshelbianTopologicalDerivativeOperators.cpp.

Constructor & Destructor Documentation

OpTopologicalObjectivePythonImpl()

EshelbianPlasticity::OpTopologicalObjectivePythonImpl::OpTopologicalObjectivePythonImpl ( boost::shared_ptr< DataAtIntegrationPts >  data_at_pts_ptr, boost::shared_ptr< TopologicalData >  topo_p, boost::shared_ptr< ObjectiveFunctionData >  python_ptr  )

inline

Definition at line 206 of file EshelbianTopologicalDerivativeOperators.cpp.

      : OP(NOSPACE, OP::OPSPACE), dataAtPts(data_at_pts_ptr), topoData(topo_p),
        pythonPtr(python_ptr) {}

Member Function Documentation

doWork()

MoFEMErrorCode EshelbianPlasticity::OpTopologicalObjectivePythonImpl::doWork	(	int	side,
		EntityType	type,
		EntData &	data
	)

Definition at line 221 of file EshelbianTopologicalDerivativeOperators.cpp.

                                                                       {
  MoFEMFunctionBegin;
 
#ifndef NDEBUG
  if (!dataAtPts)
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "DataAtIntegrationPts pointer is null");
  if (!topoData)
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "Topological data pointer is null");
  if (!pythonPtr)
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "ObjectiveFunctionData pointer is null");
#endif // NDEBUG
 
  constexpr int size_full = SPACE_DIM * SPACE_DIM;
  const int nb_gauss_pts = getGaussPts().size2();
  if (!nb_gauss_pts)
    MoFEMFunctionReturnHot(0);
 
  auto stress_full_ptr = boost::make_shared<MatrixDouble>();
  stress_full_ptr->resize(size_full, nb_gauss_pts, false);
  stress_full_ptr->clear();
  auto strain_full_ptr = boost::make_shared<MatrixDouble>();
  strain_full_ptr->resize(size_full, nb_gauss_pts, false);
  strain_full_ptr->clear();
 
  auto t_stress = getFTensor2FromMat<3, 3>(*stress_full_ptr);
  auto t_strain = getFTensor2FromMat<3, 3>(*strain_full_ptr);
 
  auto t_P = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(dataAtPts->approxPAtPts);
  auto t_log_u =
      getFTensor2SymmetricFromMat<3>(dataAtPts->logStretchTensorAtPts);
 
  auto next = [&]() {
    ++t_stress;
    ++t_strain;
    ++t_P;
    ++t_log_u;
  };
 
  FTENSOR_INDEX(SPACE_DIM, i);
  FTENSOR_INDEX(SPACE_DIM, j);
 
  for (auto gg = 0; gg != nb_gauss_pts; ++gg) {
    // we have to handle all variants, that will render how the Jacobian
    // gradient is evaluated.
    t_stress(i, j) = t_P(i, j);
    t_strain(i, j) = t_log_u(i, j);
    next();
  }
 
  auto &coords = OP::getCoordsAtGaussPts();
  CHKERR pythonPtr->evalInteriorObjectiveFunction(
      coords, dataAtPts->getSmallWL2AtPts(), stress_full_ptr, strain_full_ptr,
      topoData->getObjAtPts(), false);
 
  CHKERR pythonPtr->evalInteriorObjectiveGradientStrain(
      coords, dataAtPts->getSmallWL2AtPts(), stress_full_ptr, strain_full_ptr,
      topoData->getObjDStrainAtPts(), false);
 
  CHKERR pythonPtr->evalInteriorObjectiveGradientU(
      coords, dataAtPts->getSmallWL2AtPts(), stress_full_ptr, strain_full_ptr,
      topoData->getObjDDisplacementAtPts(), false);
 
  CHKERR pythonPtr->evalInteriorObjectiveGradientStress(
      coords, dataAtPts->getSmallWL2AtPts(), stress_full_ptr, strain_full_ptr,
      topoData->getObjDStressAtPts(), false);
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

dataAtPts

boost::shared_ptr<DataAtIntegrationPts> EshelbianPlasticity::OpTopologicalObjectivePythonImpl::dataAtPts

private

Definition at line 216 of file EshelbianTopologicalDerivativeOperators.cpp.

pythonPtr

boost::shared_ptr<ObjectiveFunctionData> EshelbianPlasticity::OpTopologicalObjectivePythonImpl::pythonPtr

private

Definition at line 218 of file EshelbianTopologicalDerivativeOperators.cpp.

topoData

boost::shared_ptr<TopologicalData> EshelbianPlasticity::OpTopologicalObjectivePythonImpl::topoData

private

Definition at line 217 of file EshelbianTopologicalDerivativeOperators.cpp.

---

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

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