OpStateSensitivity Struct Reference

Inheritance diagram for OpStateSensitivity:

Collaboration diagram for OpStateSensitivity:

Public Types
using	OP = DomainBaseOp

Public Member Functions
	OpStateSensitivity (const std::string field_name, boost::shared_ptr< ObjectiveFunctionData > python_ptr, boost::shared_ptr< HookeOps::CommonData > comm_ptr, boost::shared_ptr< MatrixDouble >u_ptr)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &data)

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

Detailed Description

Examples

mofem/tutorials/vec-7_shape_optimisation/adjoint.cpp.

Definition at line 1476 of file adjoint.cpp.

Member Typedef Documentation

OP

using OpStateSensitivity::OP = DomainBaseOp

Definition at line 1478 of file adjoint.cpp.

Constructor & Destructor Documentation

OpStateSensitivity()

OpStateSensitivity::OpStateSensitivity ( const std::string  field_name, boost::shared_ptr< ObjectiveFunctionData >  python_ptr, boost::shared_ptr< HookeOps::CommonData >  comm_ptr, boost::shared_ptr< MatrixDouble >  u_ptr  )

inline

Definition at line 1479 of file adjoint.cpp.

      : OP(field_name, field_name, DomainEleOp::OPROW), pythonPtr(python_ptr), commPtr(comm_ptr),
        uPtr(u_ptr){}

Member Function Documentation

iNtegrate()

MoFEMErrorCode OpStateSensitivity::iNtegrate ( EntitiesFieldData::EntData &  data )

inline

Examples

mofem/tutorials/vec-7_shape_optimisation/adjoint.cpp.

Definition at line 1486 of file adjoint.cpp.

  {
    MoFEMFunctionBegin;
    FTENSOR_INDEX(SPACE_DIM, i);
    FTENSOR_INDEX(SPACE_DIM, j);
    FTENSOR_INDEX(SPACE_DIM, k);
    FTENSOR_INDEX(SPACE_DIM, l);
 
    constexpr auto symm_size = (SPACE_DIM * (SPACE_DIM + 1)) / 2; 
    auto nb_gauss_pts = getGaussPts().size2(); 
 
    auto objective_dstress = boost::make_shared<MatrixDouble>(symm_size, nb_gauss_pts);
    auto objective_dstrain = boost::make_shared<MatrixDouble>(symm_size, nb_gauss_pts);
    auto objective_du = boost::make_shared<MatrixDouble>(SPACE_DIM, nb_gauss_pts);
 
    
 
    auto evaluate_python = [&]() 
    {
      MoFEMFunctionBegin;
      auto &coords = OP::getCoordsAtGaussPts();
      CHKERR pythonPtr->evalInteriorObjectiveGradientStress(
          coords, uPtr, commPtr->getMatCauchyStress(), commPtr->getMatStrain(),
          objective_dstress);
      CHKERR pythonPtr->evalInteriorObjectiveGradientStrain(
          coords, uPtr, commPtr->getMatCauchyStress(), commPtr->getMatStrain(),
          objective_dstrain);
      CHKERR pythonPtr->evalInteriorObjectiveGradientU(
          coords, uPtr, commPtr->getMatCauchyStress(), commPtr->getMatStrain(),
          objective_du);
 
      auto vol = OP::getMeasure();
      auto t_w = OP::getFTensor0IntegrationWeight();
 
      auto t_D = getFTensor4DdgFromMat<SPACE_DIM, SPACE_DIM, 0>(*(commPtr->matDPtr));
      auto t_row_grad = data.getFTensor1DiffN<SPACE_DIM>();
      auto t_row_base = data.getFTensor0N();
 
      auto t_obj_dstress = getFTensor2SymmetricFromMat<SPACE_DIM>(*objective_dstress);
      auto t_obj_dstrain = getFTensor2SymmetricFromMat<SPACE_DIM>(*objective_dstrain);
      auto t_obj_du = getFTensor1FromMat<SPACE_DIM>(*objective_du);
      
 
      for(int gg = 0; gg != nb_gauss_pts; ++gg)
      {
        const double alpha = t_w * vol;
        FTensor::Tensor2<double, SPACE_DIM, SPACE_DIM> t_adjoint_stress;
        t_adjoint_stress(i,j) = t_D(i, j, k, l) * t_obj_dstress(k, l) + t_obj_dstrain(i, j);
 
        auto t_nf = OP::template getNf<SPACE_DIM>();
        int rr = 0;
        for (; rr != OP::nbRows / SPACE_DIM; rr++)
         {
          t_nf(j) += alpha * t_row_grad(i) * t_adjoint_stress(i, j);
          t_nf(j) += alpha * t_row_base * t_obj_du(j);
 
          ++t_row_grad;
          ++t_row_base;
          ++t_nf;
         }
 
        for (; rr < OP::nbRowBaseFunctions; ++rr)
         {
          ++t_row_grad;
          ++t_row_base;
         }
         ++t_obj_dstrain;
         ++t_obj_dstress;
         ++t_obj_du;
         ++t_w;
 
      }
      MoFEMFunctionReturn(0);
 
      
    };
    CHKERR evaluate_python();
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

commPtr

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

private

Examples

mofem/tutorials/vec-7_shape_optimisation/adjoint.cpp.

Definition at line 1569 of file adjoint.cpp.

pythonPtr

boost::shared_ptr<ObjectiveFunctionData> OpStateSensitivity::pythonPtr

private

Examples

mofem/tutorials/vec-7_shape_optimisation/adjoint.cpp.

Definition at line 1568 of file adjoint.cpp.

uPtr

boost::shared_ptr<MatrixDouble> OpStateSensitivity::uPtr

private

Examples

mofem/tutorials/vec-7_shape_optimisation/adjoint.cpp.

Definition at line 1570 of file adjoint.cpp.

---

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

• tutorials/vec-7_shape_optimisation/adjoint.cpp