#include <src/finite_elements/FormsBrokenSpaceConstraintImpl.hpp>

Inheritance diagram for MoFEM::OpBrokenSpaceConstrainDHybridImpl< FIELD_DIM, GAUSS, OpBase >:

Collaboration diagram for MoFEM::OpBrokenSpaceConstrainDHybridImpl< FIELD_DIM, GAUSS, OpBase >:

Public Types
using	OP = OpBase

Public Member Functions
	OpBrokenSpaceConstrainDHybridImpl (const std::string row_field, boost::shared_ptr< std::vector< BrokenBaseSideData >> broken_side_data_ptr, const double beta)

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

Private Attributes
double	scalarBeta

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

Detailed Description

template<int FIELD_DIM, typename OpBase>
struct MoFEM::OpBrokenSpaceConstrainDHybridImpl< FIELD_DIM, GAUSS, OpBase >

Definition at line 484 of file FormsBrokenSpaceConstraintImpl.hpp.

Member Typedef Documentation

◆ OP

template<int FIELD_DIM, typename OpBase >

using MoFEM::OpBrokenSpaceConstrainDHybridImpl< FIELD_DIM, GAUSS, OpBase >::OP = OpBase

Definition at line 487 of file FormsBrokenSpaceConstraintImpl.hpp.

Constructor & Destructor Documentation

◆ OpBrokenSpaceConstrainDHybridImpl()

template<int FIELD_DIM, typename OpBase >

MoFEM::OpBrokenSpaceConstrainDHybridImpl< FIELD_DIM, GAUSS, OpBase >::OpBrokenSpaceConstrainDHybridImpl	(	const std::string	row_field,
		boost::shared_ptr< std::vector< BrokenBaseSideData >>	broken_side_data_ptr,
		const double	beta
	)

inline

Definition at line 489 of file FormsBrokenSpaceConstraintImpl.hpp.

493 : OpBase(row_field, row_field, OpBase::OPROW),

494 brokenSideDataPtr(broken_side_data_ptr), scalarBeta(beta) {}

Member Function Documentation

◆ iNtegrate()

template<int FIELD_DIM, typename OpBase >

MoFEMErrorCode MoFEM::OpBrokenSpaceConstrainDHybridImpl< FIELD_DIM, GAUSS, OpBase >::iNtegrate ( EntitiesFieldData::EntData & row_data )

private

Definition at line 504 of file FormsBrokenSpaceConstraintImpl.hpp.

                                         {
   MoFEMFunctionBegin;
  
   FTENSOR_INDEX(FIELD_DIM, i);
   FTENSOR_INDEX(3, J);
  
   OP::locF.resize(row_data.getIndices().size(), false);
   OP::locF.clear();
  
   for (auto &bd : *brokenSideDataPtr) {
     auto t_w = this->getFTensor0IntegrationWeight();
     auto t_normal = OpBase::getFTensor1NormalsAtGaussPts();
     auto t_row_base = row_data.getFTensor0N();
     auto t_flux = getFTensor2FromMat<FIELD_DIM, 3>(bd.getFlux());
     auto nb_base_functions = row_data.getN().size2() / 3;
     auto sense = bd.getSense();
     for (size_t gg = 0; gg != OpBase::nbIntegrationPts; ++gg) {
       auto t_vec = getFTensor1FromPtr<FIELD_DIM>(&*OP::locF.data().begin());
       size_t rr = 0;
       for (; rr != row_data.getIndices().size() / FIELD_DIM; ++rr) {
         t_vec(i) += (sense * t_w) * t_row_base * t_normal(J) * t_flux(i, J);
         ++t_row_base;
         ++t_vec;
       }
       for (; rr < nb_base_functions; ++rr)
         ++t_row_base;
       ++t_w;
       ++t_normal;
       ++t_flux;
     }
   }
  
   OP::locF *= scalarBeta;
  
   MoFEMFunctionReturn(0);
 }