#include <src/finite_elements/LinearFormsIntegratorsImpl.hpp>

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

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

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

	OpMixTensorTimesGradUImpl (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

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

Detailed Description

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

Definition at line 347 of file LinearFormsIntegratorsImpl.hpp.

Constructor & Destructor Documentation

◆ OpMixTensorTimesGradUImpl() [1/2]

template<int SPACE_DIM, typename OpBase >

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

inline

Definition at line 348 of file LinearFormsIntegratorsImpl.hpp.

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

◆ OpMixTensorTimesGradUImpl() [2/2]

template<int SPACE_DIM, typename OpBase >

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

inline

Definition at line 352 of file LinearFormsIntegratorsImpl.hpp.

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

356 betaCoeff(beta_fun) {}

Member Function Documentation

◆ iNtegrate()

template<int SPACE_DIM, typename OpBase >

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

protected

Definition at line 976 of file LinearFormsIntegratorsImpl.hpp.

                                         {
   MoFEMFunctionBegin;
  
   const size_t nb_base_functions = row_data.getN().size2() / 3;
   auto t_w = this->getFTensor0IntegrationWeight();
   auto t_coords = this->getFTensor1CoordsAtGaussPts();
   auto t_base = row_data.getFTensor1N<3>();
   auto t_grad = getFTensor2FromMat<SPACE_DIM, 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 * betaCoeff(t_coords(0), t_coords(1), t_coords(2)) *
                  t_base(j) * t_grad(i, j);
       ++t_nf;
       ++t_base;
     }
     for (; bb < nb_base_functions; ++bb)
       ++t_base;
  
     ++t_grad;
     ++t_coords;
     ++t_w;
   }
  
   MoFEMFunctionReturn(0);
 }