#include <src/finite_elements/LinearFormsIntegratorsImpl.hpp>

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

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

Public Member Functions
	OpMixVecTimesDivLambdaImpl (const std::string field_name, boost::shared_ptr< MatrixDouble > mat_vals)

	OpMixVecTimesDivLambdaImpl (const std::string field_name, boost::shared_ptr< MatrixDouble > mat_vals, ScalarFun beta_fun)

Protected Member Functions
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &data)

Protected Attributes
boost::shared_ptr< MatrixDouble >	matVals

ScalarFun	betaCoeff = [](double, double, double) constexpr { return 1; }

FTensor::Index< 'i', SPACE_DIM >	i

Detailed Description

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

Definition at line 313 of file LinearFormsIntegratorsImpl.hpp.

Constructor & Destructor Documentation

◆ OpMixVecTimesDivLambdaImpl() [1/2]

template<int SPACE_DIM, typename OpBase >

MoFEM::OpMixVecTimesDivLambdaImpl< SPACE_DIM, GAUSS, OpBase >::OpMixVecTimesDivLambdaImpl	(	const std::string	field_name,
		boost::shared_ptr< MatrixDouble >	mat_vals
	)

inline

Definition at line 314 of file LinearFormsIntegratorsImpl.hpp.

316 : OpBase(field_name, field_name, OpBase::OPROW), matVals(mat_vals) {}

◆ OpMixVecTimesDivLambdaImpl() [2/2]

template<int SPACE_DIM, typename OpBase >

MoFEM::OpMixVecTimesDivLambdaImpl< SPACE_DIM, GAUSS, OpBase >::OpMixVecTimesDivLambdaImpl	(	const std::string	field_name,
		boost::shared_ptr< MatrixDouble >	mat_vals,
		ScalarFun	beta_fun
	)

inline

Definition at line 318 of file LinearFormsIntegratorsImpl.hpp.

321 : OpBase(field_name, field_name, OpBase::OPROW), matVals(mat_vals),

322 betaCoeff(beta_fun) {}

Member Function Documentation

◆ iNtegrate()

template<int SPACE_DIM, typename OpBase >

MoFEMErrorCode MoFEM::OpMixVecTimesDivLambdaImpl< SPACE_DIM, GAUSS, OpBase >::iNtegrate ( EntitiesFieldData::EntData & data )

protected

Definition at line 1010 of file LinearFormsIntegratorsImpl.hpp.

                                         {
   MoFEMFunctionBegin;
  
   const size_t nb_base_functions = row_data.getN().size2();
   auto t_w = this->getFTensor0IntegrationWeight();
   auto t_coords = this->getFTensor1CoordsAtGaussPts();
   auto t_base = row_data.getFTensor0N();
   auto t_div = getFTensor1FromMat<SPACE_DIM>(*(matVals));
  
   for (size_t gg = 0; gg != OpBase::nbIntegrationPts; ++gg) {
     const double alpha = this->getMeasure() * t_w;
     auto t_nf = OpBase::template getNf<SPACE_DIM>();
  
     size_t bb = 0;
     for (; bb != this->nbRows / SPACE_DIM; ++bb) {
       t_nf(i) += alpha * t_base *
                  betaCoeff(t_coords(0), t_coords(1), t_coords(2)) * t_div(i);
       ++t_nf;
       ++t_base;
     }
     for (; bb < nb_base_functions; ++bb)
       ++t_base;
  
     ++t_coords;
     ++t_div;
     ++t_w;
   }
  
   MoFEMFunctionReturn(0);
 }