#include "src/finite_elements/FormsBrokenSpaceConstraintImpl.hpp"

Inheritance diagram for MoFEM::OpTopoDerivativeBrokenSpaceConstrainDHybridImpl< 3, GAUSS, OpBase >:

Collaboration diagram for MoFEM::OpTopoDerivativeBrokenSpaceConstrainDHybridImpl< 3, GAUSS, OpBase >:

Public Types
using	OP = OpBrokenTopoBase< OpBase >

Public Types inherited from MoFEM::OpBrokenTopoBase< OpBase >
using	OP = OpBase

Public Member Functions
	OpTopoDerivativeBrokenSpaceConstrainDHybridImpl (const std::string row_field, boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_side_data_ptr, boost::shared_ptr< MatrixDouble > adjoint_hybrid_ptr, boost::shared_ptr< MatrixDouble > tangent1_ptr, boost::shared_ptr< MatrixDouble > tangent2_ptr, boost::shared_ptr< double > beta_ptr, SmartPetscObj< Vec > assemble_vec, Tag th, boost::shared_ptr< double > dJ_ptr=nullptr, boost::shared_ptr< Range > ents_ptr=nullptr)

Private Member Functions
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data)

Private Attributes
boost::shared_ptr< std::vector< BrokenBaseSideData > >	brokenSideDataPtr

boost::shared_ptr< MatrixDouble >	adjointHybridPtr

boost::shared_ptr< double >	scalarBetaPtr

boost::shared_ptr< double >	dJPtr

Additional Inherited Members
Protected Member Functions inherited from MoFEM::OpBrokenTopoBase< OpBase >
	OpBrokenTopoBase (const std::string row_field, boost::shared_ptr< MatrixDouble > tangent1_diff_ptr, boost::shared_ptr< MatrixDouble > tangent2_diff_ptr, SmartPetscObj< Vec > assemble_vec, Tag th, boost::shared_ptr< Range > ents_ptr=nullptr)

MoFEMErrorCode	aSsemble (EntitiesFieldData::EntData &row_data)

Protected Attributes inherited from MoFEM::OpBrokenTopoBase< OpBase >
boost::shared_ptr< MatrixDouble >	tangent1DiffPtr

boost::shared_ptr< MatrixDouble >	tangent2DiffPtr

SmartPetscObj< Vec >	assembleVec

Tag	thGradTag

Detailed Description

template<typename OpBase>
struct MoFEM::OpTopoDerivativeBrokenSpaceConstrainDHybridImpl< 3, GAUSS, OpBase >

Definition at line 817 of file FormsBrokenSpaceConstraintImpl.hpp.

Member Typedef Documentation

◆ OP

template<typename OpBase >

using MoFEM::OpTopoDerivativeBrokenSpaceConstrainDHybridImpl< 3, GAUSS, OpBase >::OP = OpBrokenTopoBase<OpBase>

Definition at line 820 of file FormsBrokenSpaceConstraintImpl.hpp.

Constructor & Destructor Documentation

◆ OpTopoDerivativeBrokenSpaceConstrainDHybridImpl()

template<typename OpBase >

MoFEM::OpTopoDerivativeBrokenSpaceConstrainDHybridImpl< 3, GAUSS, OpBase >::OpTopoDerivativeBrokenSpaceConstrainDHybridImpl	(	const std::string	row_field,
		boost::shared_ptr< std::vector< BrokenBaseSideData > >	broken_side_data_ptr,
		boost::shared_ptr< MatrixDouble >	adjoint_hybrid_ptr,
		boost::shared_ptr< MatrixDouble >	tangent1_ptr,
		boost::shared_ptr< MatrixDouble >	tangent2_ptr,
		boost::shared_ptr< double >	beta_ptr,
		SmartPetscObj< Vec >	assemble_vec,
		Tag	th,
		boost::shared_ptr< double >	dJ_ptr = `nullptr`,
		boost::shared_ptr< Range >	ents_ptr = `nullptr`
	)

inline

Definition at line 822 of file FormsBrokenSpaceConstraintImpl.hpp.

      : OP(row_field, tangent1_ptr, tangent2_ptr, assemble_vec, th, ents_ptr),
        brokenSideDataPtr(broken_side_data_ptr),
        adjointHybridPtr(adjoint_hybrid_ptr), scalarBetaPtr(beta_ptr),
        dJPtr(dJ_ptr) {}

Member Function Documentation

◆ iNtegrate()

template<typename OpBase >

MoFEMErrorCode MoFEM::OpTopoDerivativeBrokenSpaceConstrainDHybridImpl< 3, GAUSS, OpBase >::iNtegrate ( EntitiesFieldData::EntData & row_data )

private

Definition at line 846 of file FormsBrokenSpaceConstraintImpl.hpp.

                                                                   {
  MoFEMFunctionBegin;
 
  using DL = DataLayoutTraits<DataLayout::GaussByCoeffs>;
  auto get_adjoint_lambda_at_pts =
      MatrixSizeHelper<GetFTensor1FromMatType<3, -1, DL>, DL>::get(
          *adjointHybridPtr, OpBase::nbIntegrationPts);
 
  FTENSOR_INDEX(3, i);
  FTENSOR_INDEX(3, j);
  FTENSOR_INDEX(3, k);
  FTENSOR_INDEX(3, J);
 
  OP::locF.resize(row_data.getIndices().size(), false);
  OP::locF.clear();
 
  auto get_ftensor1 = [](MatrixDouble &m) {
    return FTensor::Tensor1<FTensor::PackPtr<double *, 3>, 3>(
        &m(0, 0), &m(0, 1), &m(0, 2));
  };
 
  double d_loc_J = 0;
 
  double triangle_scale = 1.0;
  if (OP::getFEType() == MBTRI)
    triangle_scale = 0.5;
 
  for (auto &bd : *brokenSideDataPtr) {
    auto sense = bd.getSense();
 
    // This is used for testing the H(div) contravariant Piola trace P.n dA invariance
    // with respect to material map variations. The variation of the normal is
    // cancelled by the variation of bd.getFlux() in the adjoint lambda at pts
 
    if (dJPtr) {
      auto t_w = this->getFTensor0IntegrationWeight();
      auto t_t1 = get_ftensor1(*this->tangent1DiffPtr);
      auto t_t2 = get_ftensor1(*this->tangent2DiffPtr);
      auto t_flux = getFTensor2FromMat<3, 3>(bd.getFlux());
      auto t_adjoint_lambda_at_pts = get_adjoint_lambda_at_pts();
      for (size_t gg = 0; gg != OpBase::nbIntegrationPts; ++gg) {
        FTensor::Tensor1<double, 3> t_tmp_normal;
        t_tmp_normal(j) = FTensor::levi_civita(i, j, k) * t_t1(k) * t_t2(i);
        d_loc_J +=
            (sense * triangle_scale * t_w) *
            (t_adjoint_lambda_at_pts(i) * (t_flux(i, J) * t_tmp_normal(J)));
        ++t_w;
        ++t_t1;
        ++t_t2;
        ++t_flux;
        ++t_adjoint_lambda_at_pts;
      }
    }
  }
 
  // The L2 hybrid multiplier is not material-map transformed, while the H(div)
  // contravariant Piola trace P.n dA is invariant.
  if (dJPtr) {
    if (scalarBetaPtr)
      d_loc_J *= *scalarBetaPtr;
    *dJPtr += d_loc_J;
  }
 
  MoFEMFunctionReturn(0);
}