#include "src/finite_elements/FormsBrokenSpaceConstraintImpl.hpp"

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

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

Public Types
using	OP = OpBrokenTopoBase< OpBase >

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

Public Member Functions
	OpTopoDerivativeBrokenSpaceConstrainDFluxImpl (const std::string row_field, boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_side_data_ptr, boost::shared_ptr< MatrixDouble > lamgrange_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 >	lagrangePtr

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::OpTopoDerivativeBrokenSpaceConstrainDFluxImpl< 3, GAUSS, OpBase >

Definition at line 723 of file FormsBrokenSpaceConstraintImpl.hpp.

Member Typedef Documentation

◆ OP

template<typename OpBase >

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

Definition at line 726 of file FormsBrokenSpaceConstraintImpl.hpp.

Constructor & Destructor Documentation

◆ OpTopoDerivativeBrokenSpaceConstrainDFluxImpl()

template<typename OpBase >

MoFEM::OpTopoDerivativeBrokenSpaceConstrainDFluxImpl< 3, GAUSS, OpBase >::OpTopoDerivativeBrokenSpaceConstrainDFluxImpl	(	const std::string	row_field,
		boost::shared_ptr< std::vector< BrokenBaseSideData > >	broken_side_data_ptr,
		boost::shared_ptr< MatrixDouble >	lamgrange_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 728 of file FormsBrokenSpaceConstraintImpl.hpp.

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

Member Function Documentation

◆ iNtegrate()

template<typename OpBase >

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

private

Definition at line 751 of file FormsBrokenSpaceConstraintImpl.hpp.

                                        {
  MoFEMFunctionBegin;
 
  // That 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.getVarFlux()
  
  using DL = DataLayoutTraits<DataLayout::GaussByCoeffs>;
  auto get_lagrange_at_pts =
      MatrixSizeHelper<GetFTensor1FromMatType<3, -1, DL>, DL>::get(
          *lagrangePtr, OpBase::nbIntegrationPts);
 
  FTENSOR_INDEX(3, J);
  FTENSOR_INDEX(3, i);
  FTENSOR_INDEX(3, j);
  FTENSOR_INDEX(3, k);
 
  auto get_ftensor1 = [](MatrixDouble &m) {
    return FTensor::Tensor1<FTensor::PackPtr<double *, 3>, 3>(
        &m(0, 0), &m(0, 1), &m(0, 2));
  };
 
  double triangle_scale = 1.0;
  if (OP::getFEType() == MBTRI)
    triangle_scale = 0.5;
 
  double d_loc_J = 0.0;
 
  for (auto &bd : *brokenSideDataPtr) {
    auto sense = bd.getSense();
 
    if (dJPtr) {
      auto t_w = this->getFTensor0IntegrationWeight();
      auto t_t1 = get_ftensor1(*this->tangent1DiffPtr);
      auto t_t2 = get_ftensor1(*this->tangent2DiffPtr);
      auto t_lagrange_at_pts = get_lagrange_at_pts();
      auto t_adjoint_lambda = getFTensor2FromMat<3, 3>(bd.getVarFlux());
 
      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_lagrange_at_pts(i) *
                    (t_adjoint_lambda(i, J) * t_tmp_normal(J)));
        ++t_w;
        ++t_t1;
        ++t_t2;
        ++t_lagrange_at_pts;
        ++t_adjoint_lambda;
      }
    }
  }
 
  // The H(div) contravariant Piola trace P.n dA is material-map invariant.
  // The normal variation is cancelled by the variation of bd.getVarFlux().
  if (dJPtr) {
    if (scalarBetaPtr)
      d_loc_J *= *scalarBetaPtr;
    *dJPtr += d_loc_J;
  }
 
  MoFEMFunctionReturn(0);
}