#include <src/finite_elements/LinearFormsIntegratorsImpl.hpp>

Inheritance diagram for MoFEM::OpMixDivTimesUImpl< 3, 1, SPACE_DIM, GAUSS, OpBase, CoordSys >:

Collaboration diagram for MoFEM::OpMixDivTimesUImpl< 3, 1, SPACE_DIM, GAUSS, OpBase, CoordSys >:

Public Member Functions
	OpMixDivTimesUImpl (const std::string field_name, boost::shared_ptr< VectorDouble > vec_vals, ScalarFun beta=[](double, double, double) constexpr { return 1;}, boost::shared_ptr< Range > ents_ptr=nullptr)

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

Protected Attributes
ScalarFun	betaConst

boost::shared_ptr< VectorDouble >	vecVals

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

Detailed Description

template<int SPACE_DIM, typename OpBase, CoordinateTypes CoordSys>
struct MoFEM::OpMixDivTimesUImpl< 3, 1, SPACE_DIM, GAUSS, OpBase, CoordSys >

Definition at line 265 of file LinearFormsIntegratorsImpl.hpp.

Constructor & Destructor Documentation

◆ OpMixDivTimesUImpl()

template<int SPACE_DIM, typename OpBase , CoordinateTypes CoordSys>

MoFEM::OpMixDivTimesUImpl< 3, 1, SPACE_DIM, GAUSS, OpBase, CoordSys >::OpMixDivTimesUImpl	(	const std::string	field_name,
		boost::shared_ptr< VectorDouble >	vec_vals,
		ScalarFun	beta = `[](double, double, double) constexpr { return 1; }`,
		boost::shared_ptr< Range >	ents_ptr = `nullptr`
	)

inline

Definition at line 267 of file LinearFormsIntegratorsImpl.hpp.

                                                             { return 1; },
       boost::shared_ptr<Range> ents_ptr = nullptr)
       : OpBase(field_name, field_name, OpBase::OPROW, ents_ptr),
         vecVals(vec_vals), betaConst(beta) {}

Member Function Documentation

◆ iNtegrate()

template<int SPACE_DIM, typename OpBase , CoordinateTypes CoordSys>

MoFEMErrorCode MoFEM::OpMixDivTimesUImpl< 3, 1, SPACE_DIM, GAUSS, OpBase, CoordSys >::iNtegrate ( EntitiesFieldData::EntData & data )

protected

Definition at line 882 of file LinearFormsIntegratorsImpl.hpp.

                                         {
   MoFEMFunctionBegin;
  
   const size_t nb_base_functions = row_data.getN().size2() / 3;
   auto t_w = this->getFTensor0IntegrationWeight();
   // get coordinate at integration points
   auto t_coords = OpBase::getFTensor1CoordsAtGaussPts();
   auto t_diff_base = row_data.getFTensor2DiffN<3, SPACE_DIM>();
   auto t_u = getFTensor0FromVec(*(vecVals));
  
   for (size_t gg = 0; gg != OpBase::nbIntegrationPts; ++gg) {
  
     const double alpha = this->getMeasure() * t_w *
                          betaConst(t_coords(0), t_coords(1), t_coords(2));
     ;
  
     size_t bb = 0;
     for (; bb != this->nbRows; ++bb) {
       const double t_div_base = t_diff_base(j, j);
       OpBase::locF[bb] += alpha * t_div_base * t_u;
       ++t_diff_base;
     }
     for (; bb < nb_base_functions; ++bb)
       ++t_diff_base;
  
     ++t_u;
     ++t_w;
     ++t_coords;
   }
  
   MoFEMFunctionReturn(0);
 }