OpAdJointTestOp< SPACE_DIM, IntegrationType::GAUSS, DomainBaseOp > Struct Template Reference

Adjoint test operator for elastic problems. More...

Inheritance diagram for OpAdJointTestOp< SPACE_DIM, IntegrationType::GAUSS, DomainBaseOp >:

Collaboration diagram for OpAdJointTestOp< SPACE_DIM, IntegrationType::GAUSS, DomainBaseOp >:

Public Types
using	OP = DomainBaseOp

Public Member Functions
	OpAdJointTestOp (const std::string field_name, boost::shared_ptr< HookeOps::CommonData > comm_ptr)

Protected Member Functions
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &data)
	Integration implementation for adjoint test operator.

Protected Attributes
boost::shared_ptr< HookeOps::CommonData >	commPtr

Detailed Description

template<int SPACE_DIM>
struct OpAdJointTestOp< SPACE_DIM, IntegrationType::GAUSS, DomainBaseOp >

Adjoint test operator for elastic problems.

This operator implements the adjoint (transpose) of the elastic stiffness operator. It's used in the adjoint equation K^T * λ = ∂f/∂u to compute sensitivities.

The operator computes: ∫_Ω ∇v : σ dΩ where v are test functions and σ is the Cauchy stress tensor.

This is the transpose of the standard elastic operator ∫_Ω ∇u : C : ∇v dΩ

Template Parameters

SPACE_DIM

Spatial dimension (2 or 3)

Definition at line 1553 of file adjoint.cpp.

Member Typedef Documentation

OP

template<int SPACE_DIM>

using OpAdJointTestOp< SPACE_DIM, IntegrationType::GAUSS, DomainBaseOp >::OP = DomainBaseOp

Definition at line 1556 of file adjoint.cpp.

Constructor & Destructor Documentation

OpAdJointTestOp()

template<int SPACE_DIM>

OpAdJointTestOp< SPACE_DIM, IntegrationType::GAUSS, DomainBaseOp >::OpAdJointTestOp ( const std::string  field_name, boost::shared_ptr< HookeOps::CommonData >  comm_ptr  )

inline

Definition at line 1558 of file adjoint.cpp.

      : OP(field_name, field_name, DomainEleOp::OPROW), commPtr(comm_ptr) {}

Member Function Documentation

iNtegrate()

template<int SPACE_DIM>

MoFEMErrorCode OpAdJointTestOp< SPACE_DIM, IntegrationType::GAUSS, DomainBaseOp >::iNtegrate ( EntitiesFieldData::EntData &  row_data )

protected

Integration implementation for adjoint test operator.

Computes the integral: ∫_Ω ∇v : σ dΩ where v are test functions and σ is Cauchy stress tensor

This forms the right-hand side of the adjoint equation K^T * λ = ∂f/∂u

Parameters

row_data

Test function data for current element

Returns

MoFEMErrorCode Success or error code

Definition at line 1580 of file adjoint.cpp.

                                        {
  MoFEMFunctionBegin;
  // Define tensor indices for Einstein notation
  FTENSOR_INDEX(SPACE_DIM, i);
  FTENSOR_INDEX(SPACE_DIM, j);
  FTENSOR_INDEX(SPACE_DIM, k);
  FTENSOR_INDEX(SPACE_DIM, l);
  FTENSOR_INDEX(SPACE_DIM, I);
  FTENSOR_INDEX(SPACE_DIM, J);
 
  // Get element geometric information
  const double vol = OP::getMeasure();
  // Get integration weights for Gauss quadrature
  auto t_w = OP::getFTensor0IntegrationWeight();
  // Get test function gradients ∇v
  auto t_row_grad = row_data.getFTensor1DiffN<SPACE_DIM>();
  // Get Cauchy stress tensor σ from forward solution
  auto t_cauchy_stress =
      getFTensor2SymmetricFromMat<SPACE_DIM>(*(commPtr->getMatCauchyStress()));
  
  // Loop over integration points for numerical integration
  for (int gg = 0; gg != OP::nbIntegrationPts; gg++) {
    // take into account Jacobian
    const double alpha = t_w * vol;
 
    // get rhs vector
    auto t_nf = OP::template getNf<SPACE_DIM>();
    // loop over rows base functions
    int rr = 0;
    for (; rr != OP::nbRows / SPACE_DIM; rr++) {
      // Variation due to change in geometry (diff base)
      t_nf(j) += alpha * t_row_grad(i) * t_cauchy_stress(i, j);
 
      ++t_row_grad;
      ++t_nf;
    }
    for (; rr < OP::nbRowBaseFunctions; ++rr) {
      ++t_row_grad;
    }
 
    ++t_cauchy_stress;
    ++t_w; // move to another integration weight
  }
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commPtr

template<int SPACE_DIM>

boost::shared_ptr<HookeOps::CommonData> OpAdJointTestOp< SPACE_DIM, IntegrationType::GAUSS, DomainBaseOp >::commPtr

protected

Definition at line 1563 of file adjoint.cpp.

---

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

• tutorials/vec-7_shape_optimisation/adjoint.cpp