#include <src/finite_elements/LinearFormsIntegratorsImpl.hpp>

Inheritance diagram for MoFEM::OpConvectiveTermRhsImpl< 1, FIELD_DIM, SPACE_DIM, GAUSS, OpBase >:

Collaboration diagram for MoFEM::OpConvectiveTermRhsImpl< 1, FIELD_DIM, SPACE_DIM, GAUSS, OpBase >:

Public Member Functions
	OpConvectiveTermRhsImpl (const std::string field_name, boost::shared_ptr< MatrixDouble > u_ptr, boost::shared_ptr< MatrixDouble > y_grad_ptr, ConstantFun source_fun=[]() constexpr { return 1;})

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

Protected Attributes
boost::shared_ptr< MatrixDouble >	uPtr

boost::shared_ptr< MatrixDouble >	yGradPtr

ConstantFun	alphaConstant

Detailed Description

template<int FIELD_DIM, int SPACE_DIM, typename OpBase>
struct MoFEM::OpConvectiveTermRhsImpl< 1, FIELD_DIM, SPACE_DIM, GAUSS, OpBase >

Definition at line 474 of file LinearFormsIntegratorsImpl.hpp.

Constructor & Destructor Documentation

◆ OpConvectiveTermRhsImpl()

template<int FIELD_DIM, int SPACE_DIM, typename OpBase >

MoFEM::OpConvectiveTermRhsImpl< 1, FIELD_DIM, SPACE_DIM, GAUSS, OpBase >::OpConvectiveTermRhsImpl	(	const std::string	field_name,
		boost::shared_ptr< MatrixDouble >	u_ptr,
		boost::shared_ptr< MatrixDouble >	y_grad_ptr,
		ConstantFun	source_fun = `[]() constexpr { return 1; }`
	)

inline

Definition at line 476 of file LinearFormsIntegratorsImpl.hpp.

                                               { return 1; })
       : OpBase(field_name, field_name, OpBase::OPROW), uPtr(u_ptr),
         yGradPtr(y_grad_ptr), alphaConstant(source_fun) {}

Member Function Documentation

◆ iNtegrate()

template<int FIELD_DIM, int SPACE_DIM, typename OpBase >

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

protected

Definition at line 1203 of file LinearFormsIntegratorsImpl.hpp.

                                         {
   MoFEMFunctionBegin;
  
   auto t_w = this->getFTensor0IntegrationWeight();
   auto t_base = row_data.getFTensor0N();
  
   auto t_u = getFTensor1FromMat<SPACE_DIM>(*uPtr);
   auto t_grad_y = getFTensor2FromMat<FIELD_DIM, SPACE_DIM>(*yGradPtr);
  
   FTensor::Index<'i', SPACE_DIM> i;
   FTensor::Index<'J', FIELD_DIM> J;
   const double alpha_constant = alphaConstant();
   for (size_t gg = 0; gg != OpBase::nbIntegrationPts; ++gg) {
  
     // get element volume
     const double vol = OpBase::getMeasure();
  
     FTensor::Tensor1<double, FIELD_DIM> t_c;
     t_c(J) = (t_grad_y(J, i) * t_u(i)) * (t_w * vol * alpha_constant);
  
     auto t_nf = getFTensor1FromArray<FIELD_DIM, FIELD_DIM>(OpBase::locF);
     int rr = 0;
     for (; rr != OpBase::nbRows / FIELD_DIM; ++rr) {
       t_nf(J) += t_base * t_c(J);
       ++t_base;
       ++t_nf;
     }
     for (; rr < OpBase::nbRowBaseFunctions; ++rr)
       ++t_base;
     ++t_w; // move to another integration weight
  
     ++t_u;
     ++t_grad_y;
   }
  
   MoFEMFunctionReturn(0);
 }