#include <src/finite_elements/LinearFormsIntegratorsImpl.hpp>

Inheritance diagram for MoFEM::OpSourceImpl< 1, FIELD_DIM, GAUSS, SourceFunctionSpecialization::S< OpBase > >:

Collaboration diagram for MoFEM::OpSourceImpl< 1, FIELD_DIM, GAUSS, SourceFunctionSpecialization::S< OpBase > >:

Public Member Functions
	OpSourceImpl (const std::string field_name, TimeFun time_fun, VectorFun< FIELD_DIM > source_fun, boost::shared_ptr< Range > ents_ptr=nullptr)
	Construct a new Op Source Impl object. More...

	OpSourceImpl (const std::string field_name, VectorFun< FIELD_DIM > source_fun, boost::shared_ptr< Range > ents_ptr=nullptr)
	Construct a new Op Source Impl object. More...

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

Protected Attributes
VectorFun< FIELD_DIM >	sourceFun

Detailed Description

template<int FIELD_DIM, typename OpBase>
struct MoFEM::OpSourceImpl< 1, FIELD_DIM, GAUSS, SourceFunctionSpecialization::S< OpBase > >

Definition at line 66 of file LinearFormsIntegratorsImpl.hpp.

Constructor & Destructor Documentation

◆ OpSourceImpl() [1/2]

template<int FIELD_DIM, typename OpBase >

MoFEM::OpSourceImpl< 1, FIELD_DIM, GAUSS, SourceFunctionSpecialization::S< OpBase > >::OpSourceImpl	(	const std::string	field_name,
		TimeFun	time_fun,
		VectorFun< FIELD_DIM >	source_fun,
		boost::shared_ptr< Range >	ents_ptr = `nullptr`
	)

inline

Construct a new Op Source Impl object.

Parameters

field_name
time_fun
source_fun
ents_ptr

Definition at line 77 of file LinearFormsIntegratorsImpl.hpp.

80 : OpBase(field_name, field_name, OpBase::OPROW, time_fun, ents_ptr) {}

◆ OpSourceImpl() [2/2]

template<int FIELD_DIM, typename OpBase >

MoFEM::OpSourceImpl< 1, FIELD_DIM, GAUSS, SourceFunctionSpecialization::S< OpBase > >::OpSourceImpl	(	const std::string	field_name,
		VectorFun< FIELD_DIM >	source_fun,
		boost::shared_ptr< Range >	ents_ptr = `nullptr`
	)

inline

Construct a new Op Source Impl object.

Parameters

field_name
source_fun
ents_ptr

Definition at line 89 of file LinearFormsIntegratorsImpl.hpp.

91 : OpBase(field_name, field_name, OpBase::OPROW, ents_ptr),

92 sourceFun(source_fun) {}

Member Function Documentation

◆ iNtegrate()

template<int FIELD_DIM, typename OpBase >

MoFEMErrorCode MoFEM::OpSourceImpl< 1, FIELD_DIM, GAUSS, SourceFunctionSpecialization::S< OpBase > >::iNtegrate ( EntitiesFieldData::EntData & data )

protected

Definition at line 526 of file LinearFormsIntegratorsImpl.hpp.

                                                   {
   FTensor::Index<'i', FIELD_DIM> i;
   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();
   // get coordinate at integration points
   auto t_coords = OpBase::getFTensor1CoordsAtGaussPts();
   // loop over integration points
   for (int gg = 0; gg != OpBase::nbIntegrationPts; gg++) {
     // source file
     auto t_source = sourceFun(t_coords(0), t_coords(1), t_coords(2));
     // take into account Jacobian
     const double alpha = t_w * vol;
     // loop over rows base functions
     auto t_nf = getFTensor1FromArray<FIELD_DIM, FIELD_DIM>(OpBase::locF);
     int rr = 0;
     for (; rr != OpBase::nbRows / FIELD_DIM; ++rr) {
       t_nf(i) += alpha * t_row_base * t_source(i);
       ++t_row_base;
       ++t_nf;
     }
     for (; rr < OpBase::nbRowBaseFunctions; ++rr)
       ++t_row_base;
     ++t_coords;
     ++t_w; // move to another integration weight
   }
   MoFEMFunctionReturn(0);
 }