EshelbianTopologicalDerivative.hpp

Go to the documentation of this file.

/**
 * @file EshelbianTopologicalDerivative.hpp
 * @brief 
 * @version 0.1
 * @date 2026-02-11
 * 
 * @copyright Copyright (c) 2026
 * 
 */
 
#ifndef ESHELBIAN_TOPOLOGICAL_DERIVATIVE_HPP
#define ESHELBIAN_TOPOLOGICAL_DERIVATIVE_HPP
 
namespace EshelbianPlasticity {
 
struct TopologicalData : boost::enable_shared_from_this<TopologicalData> {
  MatrixDouble jacobianAtPts;
  MatrixDouble invJacobianAtPtr;
  VectorDouble detJacobianAtPts;
  MatrixDouble tangent1AtPts;
  MatrixDouble tangent2AtPts;
  MatrixDouble normalAtPts;
  VectorDouble objAtPts;
  MatrixDouble objDStrainAtPts;
  MatrixDouble objDStressAtPts;
  MatrixDouble objDDisplacementAtPts;
  boost::shared_ptr<MatrixDouble> getJacobianAtPts() {
    return boost::shared_ptr<MatrixDouble>(shared_from_this(), &jacobianAtPts);
  }
  boost::shared_ptr<MatrixDouble> getInvJacobianAtPts() {
    return boost::shared_ptr<MatrixDouble>(shared_from_this(),
                                           &invJacobianAtPtr);
  }
  boost::shared_ptr<VectorDouble> getDetJacobianAtPts() {
    return boost::shared_ptr<VectorDouble>(shared_from_this(),
                                           &detJacobianAtPts);
  }
  boost::shared_ptr<MatrixDouble> getTangent1AtPts() {
    return boost::shared_ptr<MatrixDouble>(shared_from_this(), &tangent1AtPts);
  }
  boost::shared_ptr<MatrixDouble> getTangent2AtPts() {
    return boost::shared_ptr<MatrixDouble>(shared_from_this(), &tangent2AtPts);
  }
  boost::shared_ptr<MatrixDouble> getNormalAtPts() {
    return boost::shared_ptr<MatrixDouble>(shared_from_this(), &normalAtPts);
  }
  boost::shared_ptr<VectorDouble> getObjAtPts() {
    return boost::shared_ptr<VectorDouble>(shared_from_this(), &objAtPts);
  }
  boost::shared_ptr<MatrixDouble> getObjDStrainAtPts() {
    return boost::shared_ptr<MatrixDouble>(shared_from_this(),
                                           &objDStrainAtPts);
  }
  boost::shared_ptr<MatrixDouble> getObjDStressAtPts() {
    return boost::shared_ptr<MatrixDouble>(shared_from_this(),
                                           &objDStressAtPts);
  }
  boost::shared_ptr<MatrixDouble> getObjDDisplacementAtPts() {
    return boost::shared_ptr<MatrixDouble>(shared_from_this(),
                                           &objDDisplacementAtPts);
  }
};
 
struct TopologicalTAOCtx {
  TopologicalTAOCtx() = default;
  virtual ~TopologicalTAOCtx() = default;
};
 
boost::shared_ptr<TopologicalTAOCtx> createTopologicalTAOCtx(
    EshelbianCore *ep,
    ForcesAndSourcesCore::GaussHookFun set_integration_at_interior,
    ForcesAndSourcesCore::GaussHookFun set_integration_at_face,
    SmartPetscObj<TS> time_solver);
 
MoFEMErrorCode topologicalEvaluateObjectiveAndGradient(Tao tao, Vec sol,
                                                       PetscReal *f, Vec g,
                                                       void *ctx);
 
MoFEMErrorCode testTopologicalDerivative(TopologicalTAOCtx *ctx_ptr, Vec sol,
                                         PetscReal *f, Vec g);
 
}
 
#endif // ESHELBIAN_TOPOLOGICAL_DERIVATIVE_HPP