Assemble constrains vector. More...

#include "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.

virtual MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &data)=0
	Class dedicated to integrate operator.

virtual MoFEMErrorCode	aSsemble (EntitiesFieldData::EntData &data)=0
	Class dedicated to assemble operator to global system vector.

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

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

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

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

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

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

VectorDouble	locVec

const double	bEta

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

int	nbIntegrationPts
	number of integration points

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 419 of file PoissonOperators.hpp.

Member Typedef Documentation

◆ FVal

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

Examples: PoissonOperators.hpp.

Definition at line 423 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 425 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 )

inlineprotected

assemble constrains vectors

Examples: PoissonOperators.hpp.

Definition at line 480 of file PoissonOperators.hpp.

                                                          {
    MoFEMFunctionBegin;
    const int *indices = &*data.getIndices().data().begin();
    const double *vals = &*locVec.data().begin();
    Vec f = getFEMethod()->ksp_f != PETSC_NULLPTR ? 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 )

inlineprotected

Integrate local constrains vector.

Examples: PoissonOperators.hpp.

Definition at line 442 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);
  }