EshelbianPlasticity::TopologicalTAOCtxImpl Struct Reference

Inheritance diagram for EshelbianPlasticity::TopologicalTAOCtxImpl:

Collaboration diagram for EshelbianPlasticity::TopologicalTAOCtxImpl:

Public Member Functions
	TopologicalTAOCtxImpl (EshelbianCore *ep, ForcesAndSourcesCore::GaussHookFun set_integration_at_interior, ForcesAndSourcesCore::GaussHookFun set_integration_at_face, SmartPetscObj< TS > time_solver)
Public Member Functions inherited from EshelbianPlasticity::TopologicalTAOCtx
	TopologicalTAOCtx ()=default
virtual	~TopologicalTAOCtx ()=default

Private Member Functions
auto	getMaterialH1PositionsStructureRawPtr ()
auto	getVertexMaterialH1PositionsTag ()
Tag	getTag (moab::Interface &moab, std::string tag_name, int size)
Tag	getJdXTag ()

Private Attributes
EshelbianCore *	ep_ptr
ForcesAndSourcesCore::GaussHookFun	integrationAtInterior
ForcesAndSourcesCore::GaussHookFun	integrationAtFace
SmartPetscObj< TS >	timeSolver
CommInterface::EntitiesPetscVector	geometryVec
CommInterface::EntitiesPetscVector	gradJVec

Friends
MoFEMErrorCode	cohesiveEvaluateObjectiveAndGradient (Tao tao, Vec x, PetscReal *f, Vec g, void *ctx)
MoFEMErrorCode	testTopologicalDerivative (TopologicalTAOCtx *ctx_ptr, const double alpha, const double rho, const double alpha_omega)

Detailed Description

Definition at line 31 of file EshelbianTopologicalDerivative.cpp.

Constructor & Destructor Documentation

TopologicalTAOCtxImpl()

EshelbianPlasticity::TopologicalTAOCtxImpl::TopologicalTAOCtxImpl ( EshelbianCore *  ep, ForcesAndSourcesCore::GaussHookFun  set_integration_at_interior, ForcesAndSourcesCore::GaussHookFun  set_integration_at_face, SmartPetscObj< TS >  time_solver  )

inline

Definition at line 32 of file EshelbianTopologicalDerivative.cpp.

      : ep_ptr(ep), integrationAtInterior(set_integration_at_interior),
        integrationAtFace(set_integration_at_face), timeSolver(time_solver),
 
        geometryVec(CommInterface::createEntitiesPetscVector(
            ep->mField.get_comm(), ep->mField.get_moab(), 0, 3, Sev::noisy)),
        gradJVec(CommInterface::createEntitiesPetscVector(
            ep->mField.get_comm(), ep->mField.get_moab(), 0, 3, Sev::noisy)) {}

Member Function Documentation

getJdXTag()

Tag EshelbianPlasticity::TopologicalTAOCtxImpl::getJdXTag ( )

inlineprivate

Definition at line 73 of file EshelbianTopologicalDerivative.cpp.

{ return getTag(ep_ptr->mField.get_moab(), "DJ_DX", 3); }

getMaterialH1PositionsStructureRawPtr()

auto EshelbianPlasticity::TopologicalTAOCtxImpl::getMaterialH1PositionsStructureRawPtr ( )

inlineprivate

Definition at line 54 of file EshelbianTopologicalDerivative.cpp.

                                               {
    return ep_ptr->mField.get_field_structure(ep_ptr->materialH1Positions);
  }

getTag()

Tag EshelbianPlasticity::TopologicalTAOCtxImpl::getTag ( moab::Interface &  moab, std::string  tag_name, int  size  )

inlineprivate

Definition at line 63 of file EshelbianTopologicalDerivative.cpp.

                                                                {
    Tag tag;
    rval = moab.tag_get_handle(tag_name.c_str(), size, MB_TYPE_DOUBLE, tag,
                               MB_TAG_CREAT | MB_TAG_DENSE, NULL);
    if (rval == MB_ALREADY_ALLOCATED)
      rval = MB_SUCCESS;
    CHK_MOAB_THROW(rval, "Failed to get tag " + tag_name);
    return tag;
  };

getVertexMaterialH1PositionsTag()

auto EshelbianPlasticity::TopologicalTAOCtxImpl::getVertexMaterialH1PositionsTag ( )

inlineprivate

Definition at line 58 of file EshelbianTopologicalDerivative.cpp.

                                         {
    auto field_ptr = getMaterialH1PositionsStructureRawPtr();
    return field_ptr->th_FieldDataVerts;
  }

Friends And Related Symbol Documentation

cohesiveEvaluateObjectiveAndGradient

MoFEMErrorCode cohesiveEvaluateObjectiveAndGradient ( Tao  tao, Vec  x, PetscReal *  f, Vec  g, void *  ctx  )

friend

Definition at line 1332 of file EshelbianCohesive.cpp.

                                                               {
 
  MoFEMFunctionBegin;
  auto cohesive_ctx = static_cast<CohesiveTAOCtxImpl *>(ctx);
  auto &ep = *(cohesive_ctx->ep_ptr);
 
#ifndef NDEBUG
  CHKERR VecView(delta_kappa, PETSC_VIEWER_STDOUT_WORLD);
#endif
  // Set delta_kappa values to the tag
  CHKERR VecCopy(delta_kappa, cohesive_ctx->kappaVec.second);
  CHKERR VecGhostUpdateBegin(cohesive_ctx->kappaVec.second, INSERT_VALUES,
                             SCATTER_FORWARD);
  CHKERR VecGhostUpdateEnd(cohesive_ctx->kappaVec.second, INSERT_VALUES,
                           SCATTER_FORWARD);
  CHKERR CommInterface::setTagFromVector(
      ep.mField.get_moab(), cohesive_ctx->kappaVec,
      get_delta_kappa_tag(ep.mField.get_moab()));
 
  // solve primal problem
  CHKERR evaluatePrimalProblemCohesiveImpl(ep, cohesive_ctx->timeSolver,
                                           cohesive_ctx->kspSolVec);
  // solve adjoint problem to get lambda
  auto ksp = snesGetKSP(tsGetSNES(cohesive_ctx->timeSolver));
  CHKERR evaluateCohesiveLambdaImpl(ep, cohesive_ctx->setIntegrationAtFrontFace,
                                    ksp, cohesive_ctx->lambdaVec);
  // evaluate dissipation and its gradient
  CHKERR evaluateDissipationAndGradImpl(
      ep, cohesive_ctx->setIntegrationAtFrontFace, cohesive_ctx->lambdaVec,
      cohesive_ctx->dissipationVec, cohesive_ctx->gradDissipationVec);
 
  CHKERR VecSum(cohesive_ctx->dissipationVec.second, f);
  CHKERR VecCopy(cohesive_ctx->gradDissipationVec.second, g);
 
  MOFEM_LOG("EP", Sev::inform)
      << "Cohesive objective function (negative total dissipation): " << *f;
 
  MoFEMFunctionReturn(0);
}

testTopologicalDerivative

MoFEMErrorCode testTopologicalDerivative ( TopologicalTAOCtx *  ctx_ptr, const double  alpha, const double  rho, const double  alpha_omega  )

friend

Definition at line 100 of file EshelbianTopologicalDerivative.cpp.

                                                                   {
 
  MoFEMFunctionBegin;
  auto *ctx_impl_ptr = dynamic_cast<TopologicalTAOCtxImpl *>(ctx_ptr);
  if (!ctx_impl_ptr)
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "Invalid context pointer type");
  auto &ep = *ctx_impl_ptr->ep_ptr;
  auto &m_field = ep.mField;
  auto interior_integration_hook = ctx_impl_ptr->integrationAtInterior;
  auto boundary_integration_hook = ctx_impl_ptr->integrationAtFace;
  auto th_top_tag = ctx_impl_ptr->getVertexMaterialH1PositionsTag();
  auto th_grad_tag = ctx_impl_ptr->getJdXTag(m_field.get_moab());
 
  boost::shared_ptr<double> J_ptr = boost::make_shared<double>(0.);
  auto dJ_dX_vec = ctx_impl_ptr->gradJVec.second;
 
  auto fe_spatial =
      boost::make_shared<VolumeElementForcesAndSourcesCore>(m_field);
  auto dJ_dx_vec = pushTopologicalSpatialOps(
      ep, fe_spatial, interior_integration_hook, boundary_integration_hook);
  auto fe_material =
      boost::make_shared<VolumeElementForcesAndSourcesCore>(m_field);
  CHKERR pushTopologicalMaterialOps(ep, fe_material, interior_integration_hook,
                                    boundary_integration_hook, J_ptr,
                                    th_grad_tag);
 
  auto lambda_vec = createDMVector(ep.dmElastic);
  auto ksp = snesGetKSP(tsGetSNES(ctx_impl_ptr->timeSolver));
  CHKERR KSPSolveTranspose(ksp, dJ_dx_vec, lambda_vec);
 
  auto fe_adjoint =
      boost::make_shared<VolumeElementForcesAndSourcesCore>(m_field);
  auto dJ_adjoint_vec = pushTopologicalOps_dJ_adjont_gradient(
      ep, fe_adjoint, interior_integration_hook, boundary_integration_hook,
      lambda_vec, dJ_dX_vec, alpha, rho, alpha_omega);
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

ep_ptr

EshelbianCore* EshelbianPlasticity::TopologicalTAOCtxImpl::ep_ptr

private

Definition at line 46 of file EshelbianTopologicalDerivative.cpp.

geometryVec

CommInterface::EntitiesPetscVector EshelbianPlasticity::TopologicalTAOCtxImpl::geometryVec

private

Definition at line 51 of file EshelbianTopologicalDerivative.cpp.

gradJVec

CommInterface::EntitiesPetscVector EshelbianPlasticity::TopologicalTAOCtxImpl::gradJVec

private

Definition at line 52 of file EshelbianTopologicalDerivative.cpp.

integrationAtFace

ForcesAndSourcesCore::GaussHookFun EshelbianPlasticity::TopologicalTAOCtxImpl::integrationAtFace

private

Definition at line 48 of file EshelbianTopologicalDerivative.cpp.

integrationAtInterior

ForcesAndSourcesCore::GaussHookFun EshelbianPlasticity::TopologicalTAOCtxImpl::integrationAtInterior

private

Definition at line 47 of file EshelbianTopologicalDerivative.cpp.

timeSolver

SmartPetscObj<TS> EshelbianPlasticity::TopologicalTAOCtxImpl::timeSolver

private

Definition at line 49 of file EshelbianTopologicalDerivative.cpp.

---

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

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