#include <src/finite_elements/BiLinearFormsIntegratorsImpl.hpp>

Inheritance diagram for MoFEM::OpMixDivTimesScalarImpl< SPACE_DIM, GAUSS, OpBase >:

Collaboration diagram for MoFEM::OpMixDivTimesScalarImpl< SPACE_DIM, GAUSS, OpBase >:

Public Member Functions
	OpMixDivTimesScalarImpl (const std::string row_field_name, const std::string col_field_name, ConstantFun alpha_fun, const bool assemble_transpose=false, const bool only_transpose=false)

Protected Member Functions
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

Protected Attributes
FTensor::Index< 'i', SPACE_DIM >	i
	summit Index More...

ConstantFun	alphaConstant

Detailed Description

template<int SPACE_DIM, typename OpBase>
struct MoFEM::OpMixDivTimesScalarImpl< SPACE_DIM, GAUSS, OpBase >

Definition at line 257 of file BiLinearFormsIntegratorsImpl.hpp.

Constructor & Destructor Documentation

◆ OpMixDivTimesScalarImpl()

template<int SPACE_DIM, typename OpBase >

MoFEM::OpMixDivTimesScalarImpl< SPACE_DIM, GAUSS, OpBase >::OpMixDivTimesScalarImpl	(	const std::string	row_field_name,
		const std::string	col_field_name,
		ConstantFun	alpha_fun,
		const bool	assemble_transpose = `false`,
		const bool	only_transpose = `false`
	)

inline

Definition at line 258 of file BiLinearFormsIntegratorsImpl.hpp.

       : OpBase(row_field_name, col_field_name, OpBase::OPROWCOL),
         alphaConstant(alpha_fun) {
     this->assembleTranspose = assemble_transpose;
     this->onlyTranspose = only_transpose;
   }

Member Function Documentation

◆ iNtegrate()

template<int SPACE_DIM, typename OpBase >

MoFEMErrorCode MoFEM::OpMixDivTimesScalarImpl< SPACE_DIM, GAUSS, OpBase >::iNtegrate	(	EntitiesFieldData::EntData &	row_data,
		EntitiesFieldData::EntData &	col_data
	)

protected

Definition at line 1129 of file BiLinearFormsIntegratorsImpl.hpp.

                                         {
   MoFEMFunctionBegin;
  
   auto t_w = this->getFTensor0IntegrationWeight();
  
   size_t nb_base_functions = row_data.getN().size2() / 3;
   auto t_row_diff_base = row_data.getFTensor2DiffN<3, SPACE_DIM>();
  
   const double alpha_constant = alphaConstant();
  
   for (size_t gg = 0; gg != OpBase::nbIntegrationPts; ++gg) {
     const double alpha = alpha_constant * this->getMeasure() * t_w;
  
     size_t rr = 0;
     for (; rr != OpBase::nbRows; ++rr) {
       const double t_row_div_base = t_row_diff_base(i, i);
       auto t_col_base = col_data.getFTensor0N(gg, 0);
       for (size_t cc = 0; cc != OpBase::nbCols; ++cc) {
         this->locMat(rr, cc) += alpha * t_row_div_base * t_col_base;
         ++t_col_base;
       }
       ++t_row_diff_base;
     }
     for (; rr < nb_base_functions; ++rr)
       ++t_row_diff_base;
  
     ++t_w;
   }
  
   MoFEMFunctionReturn(0);
 }