PoissonExample::Op_g Struct Reference

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. More...
virtual MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &data)=0
	Class dedicated to integrate operator. More...
virtual MoFEMErrorCode	aSsemble (EntitiesFieldData::EntData &data)=0
	Class dedicated to assemble operator to global system vector. 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 419 of file PoissonOperators.hpp.

Member Typedef Documentation

FVal

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

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

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

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_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 )

inlineprotected

Integrate local constrains vector.

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);
  }

Member Data Documentation

bEta

const double PoissonExample::Op_g::bEta

protected

Definition at line 437 of file PoissonOperators.hpp.

fValue

FVal PoissonExample::Op_g::fValue

protected

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

Definition at line 434 of file PoissonOperators.hpp.

locVec

VectorDouble PoissonExample::Op_g::locVec

protected

Definition at line 436 of file PoissonOperators.hpp.

NX

FTensor::Number<0> PoissonExample::Op_g::NX

protected

x-direction index

Definition at line 431 of file PoissonOperators.hpp.

NY

FTensor::Number<1> PoissonExample::Op_g::NY

protected

y-direction index

Definition at line 432 of file PoissonOperators.hpp.

NZ

FTensor::Number<2> PoissonExample::Op_g::NZ

protected

z-direction index

Definition at line 433 of file PoissonOperators.hpp.

---

The documentation for this struct was generated from the following file:

• tutorials/cor-2to5/src/PoissonOperators.hpp