#include <src/finite_elements/LinearFormsIntegratorsImpl.hpp>

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

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

Public Member Functions
	OpNormalMixVecTimesVectorFieldImpl (const std::string field_name, boost::shared_ptr< MatrixDouble > u_ptr, ScalarFun beta_coeff=[](double, double, double) constexpr { return 1;}, boost::shared_ptr< Range > ents_ptr=nullptr)

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

Protected Attributes
boost::shared_ptr< MatrixDouble >	uPtr

ScalarFun	betaCoeff

Detailed Description

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

Definition at line 438 of file LinearFormsIntegratorsImpl.hpp.

Constructor & Destructor Documentation

◆ OpNormalMixVecTimesVectorFieldImpl()

template<int SPACE_DIM, typename OpBase >

MoFEM::OpNormalMixVecTimesVectorFieldImpl< SPACE_DIM, GAUSS, OpBase >::OpNormalMixVecTimesVectorFieldImpl	(	const std::string	field_name,
		boost::shared_ptr< MatrixDouble >	u_ptr,
		ScalarFun	beta_coeff = `[](double, double, double) constexpr { return 1; }`,
		boost::shared_ptr< Range >	ents_ptr = `nullptr`
	)

inline

Definition at line 440 of file LinearFormsIntegratorsImpl.hpp.

                                                                   { return 1; },
       boost::shared_ptr<Range> ents_ptr = nullptr)
       : OpBase(field_name, field_name, OpBase::OPROW, ents_ptr), uPtr(u_ptr),
         betaCoeff(beta_coeff) {}

Member Function Documentation

◆ iNtegrate()

template<int SPACE_DIM, typename OpBase >

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

protected

Definition at line 1120 of file LinearFormsIntegratorsImpl.hpp.

                                         {
   MoFEMFunctionBegin;
   FTensor::Index<'i', 3> i;
   FTensor::Index<'J', SPACE_DIM> J;
  
   const size_t nb_base_functions = row_data.getN().size2() / 3;
   // get element volume
   // get integration weights
   auto t_w = OpBase::getFTensor0IntegrationWeight();
   // get base function gradient on rows
   auto t_row_base = row_data.getFTensor1N<3>();
   // get coordinate at integration points
   auto t_coords = OpBase::getFTensor1CoordsAtGaussPts();
   // get normal
   auto t_normal = OpBase::getFTensor1NormalsAtGaussPts();
   // get field
   auto t_u = getFTensor1FromMat<SPACE_DIM>(*uPtr);
   // loop over integration points
   auto a = OpBase::getMeasure();
   for (int gg = 0; gg != OpBase::nbIntegrationPts; gg++) {
     // take into account Jacobian
     auto l2 = std::sqrt(t_normal(i) * t_normal(i));
     const double alpha =
         t_w * betaCoeff(t_coords(0), t_coords(1), t_coords(2)) * (a / l2);
     // get rhs vector
     auto t_nf = OpBase::template getNf<SPACE_DIM>();
     // loop over rows base functions
     int rr = 0;
     for (; rr != OpBase::nbRows / SPACE_DIM; ++rr) {
       t_nf(J) += alpha * (t_row_base(i) * t_normal(i)) * t_u(J);
       ++t_row_base;
       ++t_nf;
     }
     for (; rr < nb_base_functions; ++rr)
       ++t_row_base;
     ++t_coords;
     ++t_w; // move to another integration weight
     ++t_normal;
     ++t_u;
   }
   MoFEMFunctionReturn(0);
 }