#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, boost::shared_ptr< double > beta_ptr, boost::shared_ptr< Range > ents_ptr=nullptr)

	OpBrokenSpaceConstrainDHybridImpl (const std::string row_field, boost::shared_ptr< std::vector< BrokenBaseSideData >> broken_side_data_ptr, double beta, boost::shared_ptr< Range > ents_ptr=nullptr)

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

Private Attributes
boost::shared_ptr< double >	scalarBetaPtr

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

Constructor & Destructor Documentation

◆ OpBrokenSpaceConstrainDHybridImpl() [1/2]

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,
		boost::shared_ptr< double >	beta_ptr,
		boost::shared_ptr< Range >	ents_ptr = `nullptr`
	)

inline

Definition at line 512 of file FormsBrokenSpaceConstraintImpl.hpp.

517 : OpBase(row_field, row_field, OpBase::OPROW, ents_ptr),

518 brokenSideDataPtr(broken_side_data_ptr), scalarBetaPtr(beta_ptr) {}

◆ OpBrokenSpaceConstrainDHybridImpl() [2/2]

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,
		double	beta,
		boost::shared_ptr< Range >	ents_ptr = `nullptr`
	)

inline

Definition at line 520 of file FormsBrokenSpaceConstraintImpl.hpp.

       : OpBrokenSpaceConstrainDHybridImpl(row_field, broken_side_data_ptr,
                                           boost::make_shared<double>(beta),
                                           ents_ptr) {}

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 536 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;
     }
   }
  
   if(scalarBetaPtr)
     OP::locF *= *scalarBetaPtr;
  
   MoFEMFunctionReturn(0);
 }