OpObjective Struct Reference

Inheritance diagram for OpObjective:

Collaboration diagram for OpObjective:

Public Types
using	OP = ForcesAndSourcesCore::UserDataOperator

Public Member Functions
	OpObjective (boost::shared_ptr< ObjectiveFunctionData > python_ptr, boost::shared_ptr< HookeOps::CommonData > comm_ptr, boost::shared_ptr< MatrixDouble > jac_ptr, boost::shared_ptr< MatrixDouble > u_ptr, boost::shared_ptr< double > glob_objective_ptr)
MoFEMErrorCode	doWork (int side, EntityType type, EntitiesFieldData::EntData &data)
	Compute objective function contributions at element level.

Private Attributes
boost::shared_ptr< ObjectiveFunctionData >	pythonPtr
boost::shared_ptr< HookeOps::CommonData >	commPtr
boost::shared_ptr< MatrixDouble >	jacPtr
boost::shared_ptr< MatrixDouble >	uPtr
boost::shared_ptr< double >	globObjectivePtr

Detailed Description

Definition at line 831 of file gradient.cpp.

Member Typedef Documentation

OP

using OpObjective::OP = ForcesAndSourcesCore::UserDataOperator

Definition at line 832 of file gradient.cpp.

Constructor & Destructor Documentation

OpObjective()

OpObjective::OpObjective ( boost::shared_ptr< ObjectiveFunctionData >  python_ptr, boost::shared_ptr< HookeOps::CommonData >  comm_ptr, boost::shared_ptr< MatrixDouble >  jac_ptr, boost::shared_ptr< MatrixDouble >  u_ptr, boost::shared_ptr< double >  glob_objective_ptr  )

inline

Definition at line 834 of file gradient.cpp.

      : OP(NOSPACE, OP::OPSPACE), pythonPtr(python_ptr), commPtr(comm_ptr),
        jacPtr(jac_ptr), uPtr(u_ptr), globObjectivePtr(glob_objective_ptr) {}

Member Function Documentation

doWork()

MoFEMErrorCode OpObjective::doWork ( int  side, EntityType  type, EntitiesFieldData::EntData &  data  )

inline

Compute objective function contributions at element level.

Evaluates Python objective function with current displacement and stress state, and accumulates global objective value and gradients.

Definition at line 848 of file gradient.cpp.

                                                        {
    MoFEMFunctionBegin;
 
    auto nb_gauss_pts =
        getGaussPts().size2(); // number of gauss points in the element
 
    auto objective_ptr = boost::make_shared<MatrixDouble>(
        1, nb_gauss_pts); // objective function values at gauss points
 
    auto evaluate_python = [&]() {
      MoFEMFunctionBegin;
      auto &coords = OP::getCoordsAtGaussPts();
      CHKERR pythonPtr->evalInteriorObjectiveFunction(
          coords, uPtr, commPtr->getMatCauchyStress(), commPtr->getMatStrain(),
          objective_ptr);
 
      auto t_obj = getFTensor0FromMat(*objective_ptr);
 
      auto vol = OP::getMeasure();
      auto t_w = getFTensor0IntegrationWeight();
      for (auto gg = 0; gg != nb_gauss_pts; ++gg) {
 
        auto alpha = t_w * vol;
        (*globObjectivePtr) += alpha * t_obj;
 
        ++t_w;
 
        ++t_obj;
      }
      MoFEMFunctionReturn(0);
    };
 
    CHKERR evaluate_python();
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

commPtr

boost::shared_ptr<HookeOps::CommonData> OpObjective::commPtr

private

Definition at line 888 of file gradient.cpp.

globObjectivePtr

boost::shared_ptr<double> OpObjective::globObjectivePtr

private

Definition at line 891 of file gradient.cpp.

jacPtr

boost::shared_ptr<MatrixDouble> OpObjective::jacPtr

private

Definition at line 889 of file gradient.cpp.

pythonPtr

boost::shared_ptr<ObjectiveFunctionData> OpObjective::pythonPtr

private

Definition at line 887 of file gradient.cpp.

uPtr

boost::shared_ptr<MatrixDouble> OpObjective::uPtr

private

Definition at line 890 of file gradient.cpp.

---

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

• tutorials/vec-7_shape_optimisation/gradient.cpp