#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< 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

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 687 of file FormsBrokenSpaceConstraintImpl.hpp.

Member Typedef Documentation

◆ OP

template<typename OpBase >

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

Definition at line 690 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< Range >	ents_ptr = `nullptr`
	)

inline

Definition at line 692 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) {}

Member Function Documentation

◆ iNtegrate()

template<typename OpBase >

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

private

Definition at line 713 of file FormsBrokenSpaceConstraintImpl.hpp.

                                        {
  MoFEMFunctionBegin;
 
  FTENSOR_INDEX(3, I);
  FTENSOR_INDEX(3, J);
  FTENSOR_INDEX(3, i);
  FTENSOR_INDEX(3, j);
  FTENSOR_INDEX(3, k);
  FTENSOR_INDEX(3, m);
  FTensor::Number<0> N0;
  FTensor::Number<1> N1;
 
  auto get_ftensor1 = [](MatrixDouble &m) {
    return FTensor::Tensor1<FTensor::PackPtr<double *, 3>, 3>(
        &m(0, 0), &m(0, 1), &m(0, 2));
  };
 
  for (auto &bd : *brokenSideDataPtr) {
 
    auto t_w = this->getFTensor0IntegrationWeight();
    auto t_lagrange = getFTensor1FromMat<3>(*lagrangePtr);
    auto t_adjoint_lambda = getFTensor2FromMat<3, 3>(bd.getVarFlux());
    auto t_t1 = get_ftensor1(*this->tangent1DiffPtr);
    auto t_t2 = get_ftensor1(*this->tangent2DiffPtr);
 
    auto t_diff_base = row_data.getFTensor1DiffN<2>();
    auto nb_base_functions = row_data.getN().size2();
 
    constexpr auto t_kd = FTensor::Kronecker_Delta<double>();
 
    for (size_t gg = 0; gg != OpBase::nbIntegrationPts; ++gg) {
 
      FTensor::Tensor2<double, 3, 3> t_normal_diff;
      t_normal_diff(j, m) = FTensor::levi_civita(i, j, k) * t_t1(k) *
                            (t_kd(i, m) * t_diff_base(N1));
      t_normal_diff(j, m) += FTensor::levi_civita(i, j, k) *
                             (t_kd(k, m) * t_diff_base(N0)) * t_t2(i);
      auto t_vec = getFTensor1FromPtr<3>(&*OP::locF.data().begin());
      size_t rr = 0;
      for (; rr != row_data.getIndices().size() / 3; ++rr) {
        t_vec(m) += (0.5 * t_w) * (t_adjoint_lambda(i, J) * t_lagrange(i)) *
                    t_normal_diff(J, m);
        ++t_diff_base;
        ++t_vec;
      }
      for (; rr < nb_base_functions; ++rr)
        ++t_diff_base;
      ++t_w;
      ++t_lagrange;
      ++t_adjoint_lambda;
      ++t_t1;
      ++t_t2;
    }
  }
 
  if (scalarBetaPtr)
    OP::locF *= *scalarBetaPtr;
 
  MoFEMFunctionReturn(0);
}