Assemble constrains vector. More...

#include <users_modules/tutorials/cor-2to5/src/PoissonOperators.hpp>

Inheritance diagram for PoissonExample::Op_g:

Collaboration diagram for PoissonExample::Op_g:

Public Types
typedef boost::function< double(const double, const double, const double)>	FVal

Public Member Functions
	Op_g (FVal f_value, const string field_name="L", const double beta=1)

Public Member Functions inherited from PoissonExample::OpBaseRhs< FaceElementForcesAndSourcesCore::UserDataOperator >
	OpBaseRhs (const std::string field_name)

MoFEMErrorCode	doWork (int row_side, EntityType row_type, EntitiesFieldData::EntData &row_data)
	This function is called by finite element. More...

Protected Member Functions
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &data)
	Integrate local constrains vector. More...

MoFEMErrorCode	aSsemble (EntitiesFieldData::EntData &data)
	assemble constrains vectors More...

Protected Attributes
FTensor::Number< 0 >	NX
	x-direction index More...

FTensor::Number< 1 >	NY
	y-direction index More...

FTensor::Number< 2 >	NZ
	z-direction index More...

FVal	fValue
	Function pointer evaluating values of "U" at the boundary. More...

VectorDouble	locVec

const double	bEta

Protected Attributes inherited from PoissonExample::OpBaseRhs< FaceElementForcesAndSourcesCore::UserDataOperator >
int	nbRows
	< error code More...

int	nbIntegrationPts
	number of integration points More...

Detailed Description

Assemble constrains vector.

\[ \mathbf{g} = \int_{\partial\Omega} \boldsymbol\psi \overline{u} \textrm{d}\partial\Omega \]

Examples: analytical_poisson_field_split.cpp.

Definition at line 416 of file PoissonOperators.hpp.

Member Typedef Documentation

◆ FVal

typedef boost::function<double(const double, const double, const double)> PoissonExample::Op_g::FVal

Examples: PoissonOperators.hpp.

Definition at line 420 of file PoissonOperators.hpp.

Constructor & Destructor Documentation

◆ Op_g()

PoissonExample::Op_g::Op_g	(	FVal	f_value,
		const string	field_name = `"L"`,
		const double	beta = `1`
	)

inline

Examples: PoissonOperators.hpp.

Definition at line 422 of file PoissonOperators.hpp.

       : OpBaseRhs<FaceElementForcesAndSourcesCore::UserDataOperator>(
             field_name),
         fValue(f_value), bEta(beta) {}

Member Function Documentation

◆ aSsemble()

MoFEMErrorCode PoissonExample::Op_g::aSsemble ( EntitiesFieldData::EntData & data )

inlineprotectedvirtual

assemble constrains vectors

Implements PoissonExample::OpBaseRhs< FaceElementForcesAndSourcesCore::UserDataOperator >.

Examples: PoissonOperators.hpp.

Definition at line 477 of file PoissonOperators.hpp.

                                                           {
     MoFEMFunctionBegin;
     const int *indices = &*data.getIndices().data().begin();
     const double *vals = &*locVec.data().begin();
     Vec f = getFEMethod()->ksp_f != PETSC_NULL ? getFEMethod()->ksp_f
                                                : getFEMethod()->snes_f;
     CHKERR VecSetValues(f, nbRows, indices, &*vals, ADD_VALUES);
     MoFEMFunctionReturn(0);
   }

◆ iNtegrate()

MoFEMErrorCode PoissonExample::Op_g::iNtegrate ( EntitiesFieldData::EntData & data )

inlineprotectedvirtual

Integrate local constrains vector.

Implements PoissonExample::OpBaseRhs< FaceElementForcesAndSourcesCore::UserDataOperator >.

Reimplemented in PoissonExample::OpResF_Boundary, and PoissonExample::OpRes_g.

Examples: PoissonOperators.hpp.

Definition at line 439 of file PoissonOperators.hpp.

                                                            {
     MoFEMFunctionBegin;
     // set size to local vector
     locVec.resize(nbRows, false);
     // clear local vector
     locVec.clear();
     // get face area
     const double area = getArea() * bEta;
     // get integration weight
     auto t_w = getFTensor0IntegrationWeight();
     // get base function
     auto t_l = data.getFTensor0N();
     // get coordinates at integration point
     auto t_coords = getFTensor1CoordsAtGaussPts();
     // make loop over integration points
     for (int gg = 0; gg != nbIntegrationPts; gg++) {
       // evaluate function on boundary and scale it by area and integration
       // weight
       double alpha =
           area * t_w * fValue(t_coords(NX), t_coords(NY), t_coords(NZ));
       // get element of vector
       FTensor::Tensor0<FTensor::PackPtr<double *, 1>> t_a(
           &*locVec.data().begin());
       //
       for (int rr = 0; rr != nbRows; rr++) {
         t_a += alpha * t_l;
         ++t_a;
         ++t_l;
       }
       ++t_w;
       ++t_coords;
     }
     MoFEMFunctionReturn(0);
   }