#include <src/finite_elements/TriLinearFormsIntegrators.hpp>

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

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

Public Member Functions
	OpConvectiveTermLhsDuImpl (const std::string field_name_row, const std::string field_name_col, boost::shared_ptr< MatrixDouble > y_grad_ptr, ConstantFun alpha_fun=[]() { return 1;})

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

Protected Attributes
ConstantFun	alphaConstant

boost::shared_ptr< MatrixDouble >	yGradPtr

Detailed Description

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

Definition at line 64 of file TriLinearFormsIntegrators.hpp.

Constructor & Destructor Documentation

◆ OpConvectiveTermLhsDuImpl()

template<int FIELD_DIM, int SPACE_DIM, typename OpBase >

MoFEM::OpConvectiveTermLhsDuImpl< 1, FIELD_DIM, SPACE_DIM, GAUSS, OpBase >::OpConvectiveTermLhsDuImpl	(	const std::string	field_name_row,
		const std::string	field_name_col,
		boost::shared_ptr< MatrixDouble >	y_grad_ptr,
		ConstantFun	alpha_fun = `[]() { return 1; }`
	)

inline

Definition at line 66 of file TriLinearFormsIntegrators.hpp.

                                    { return 1; })
       : OpBase(field_name_row, field_name_col, OpBase::OPROWCOL),
         yGradPtr(y_grad_ptr), alphaConstant(alpha_fun) {
  
     this->assembleTranspose = false;
     this->onlyTranspose = false;
     this->sYmm = false;
   }

Member Function Documentation

◆ iNtegrate()

template<int FIELD_DIM, int SPACE_DIM, typename OpBase >

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

protected

Definition at line 239 of file TriLinearFormsIntegrators.hpp.

                                         {
   MoFEMFunctionBegin;
  
   // get element volume
   const double vol = OpBase::getMeasure();
   // get integration weights
   auto t_w = OpBase::getFTensor0IntegrationWeight();
   // get base function gradient on rows
   auto t_row_base = row_data.getFTensor0N();
  
   auto t_grad_y = getFTensor2FromMat<FIELD_DIM, SPACE_DIM>(*yGradPtr);
   FTensor::Index<'I', FIELD_DIM> I;
   FTensor::Index<'k', SPACE_DIM> k;
  
   auto get_t_mat = [&](const int rr) {
     std::array<double *, FIELD_DIM * SPACE_DIM> ptrs;
     int s = 0;
     for (int j = 0; j != FIELD_DIM; ++j)
       for (auto i = 0; i != SPACE_DIM; ++i, ++s)
         ptrs[s] = &OpBase::locMat(rr + j, i);
     return FTensor::Tensor2<FTensor::PackPtr<double *, SPACE_DIM>, FIELD_DIM,
                             SPACE_DIM>(ptrs);
   };
  
   const double alpha_constant = alphaConstant();
   // loop over integration points
   for (int gg = 0; gg != OpBase::nbIntegrationPts; gg++) {
     // take into account Jacobian
     const double alpha = t_w * vol * alpha_constant;
  
     // loop over rows base functions
     int rr = 0;
     for (; rr != OpBase::nbRows / FIELD_DIM; rr++) {
       // get column base functions gradient at gauss point gg
       auto t_col_base = col_data.getFTensor0N(gg, 0);
       // get mat
       auto t_mat = get_t_mat(FIELD_DIM * rr);
       // loop over columns
       for (int cc = 0; cc != OpBase::nbCols / SPACE_DIM; cc++) {
         // calculate element of local matrix
         t_mat(I, k) += (alpha * t_row_base * t_col_base) * t_grad_y(I, k);
         ++t_col_base;
         ++t_mat;
       }
       ++t_row_base;
     }
     for (; rr < OpBase::nbRowBaseFunctions; ++rr)
       ++t_row_base;
  
     ++t_grad_y;
     ++t_w; // move to another integration weight
   }
   MoFEMFunctionReturn(0);
 };