PoissonExample::OpC Struct Reference

Calculate constrains matrix. More...

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

Inheritance diagram for PoissonExample::OpC:

Collaboration diagram for PoissonExample::OpC:

Public Member Functions
	OpC (const bool assemble_transpose)
MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

Private Member Functions
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)
	Integrate local constrains matrix. More...
MoFEMErrorCode	aSsemble (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)
	integrate local constrains matrix More...

Private Attributes
int	nbRows
	< error code More...
int	nbCols
	number of dofs on column More...
int	nbIntegrationPts
	number of integration points More...
const bool	assembleTranspose
	assemble transpose, i.e. CT if set to true More...
MatrixDouble	locMat
	local constrains matrix More...

Detailed Description

Calculate constrains matrix.

\[ \mathbf{C} = \int_{\partial\Omega} \boldsymbol\psi \boldsymbol\phi \textrm{d}\partial\Omega \]

where \(\lambda \) is base function on boundary

Definition at line 302 of file PoissonOperators.hpp.

Constructor & Destructor Documentation

OpC()

PoissonExample::OpC::OpC ( const bool  assemble_transpose )

inline

Definition at line 304 of file PoissonOperators.hpp.

      : FaceElementForcesAndSourcesCore::UserDataOperator("L", "U", OPROWCOL,
                                                          false),
        assembleTranspose(assemble_transpose) {}

Member Function Documentation

aSsemble()

MoFEMErrorCode PoissonExample::OpC::aSsemble ( EntitiesFieldData::EntData &  row_data, EntitiesFieldData::EntData &  col_data  )

inlineprivate

integrate local constrains matrix

Examples

PoissonOperators.hpp.

Definition at line 385 of file PoissonOperators.hpp.

                                                                     {
    MoFEMFunctionBegin;
    // get indices on row
    const int *row_indices = &*row_data.getIndices().data().begin();
    // get indices on column
    const int *col_indices = &*col_data.getIndices().data().begin();
    Mat B = getFEMethod()->ksp_B != PETSC_NULL ? getFEMethod()->ksp_B
                                               : getFEMethod()->snes_B;
    // assemble local matrix
    CHKERR MatSetValues(B, nbRows, row_indices, nbCols, col_indices,
                        &*locMat.data().begin(), ADD_VALUES);
    // cerr << locMat << endl;
    if (assembleTranspose) {
      // assemble transpose of local matrix
      locMat = trans(locMat);
      CHKERR MatSetValues(B, nbCols, col_indices, nbRows, row_indices,
                          &*locMat.data().begin(), ADD_VALUES);
    }
    MoFEMFunctionReturn(0);
  }

doWork()

MoFEMErrorCode PoissonExample::OpC::doWork ( int  row_side, int  col_side, EntityType  row_type, EntityType  col_type, EntitiesFieldData::EntData &  row_data, EntitiesFieldData::EntData &  col_data  )

inline

Examples

PoissonOperators.hpp.

Definition at line 309 of file PoissonOperators.hpp.

                                                            {
    MoFEMFunctionBegin;
    // get number of dofs on row
    nbRows = row_data.getIndices().size();
    // exit here if no dofs on row, nothing to do
    if (!nbRows)
      MoFEMFunctionReturnHot(0);
    // get number of dofs on column,
    nbCols = col_data.getIndices().size();
    // exit here if no dofs on roe, nothing to do
    if (!nbCols)
      MoFEMFunctionReturnHot(0);
    // get number of integration points
    nbIntegrationPts = getGaussPts().size2();
    // integrate local constrains matrix
    CHKERR iNtegrate(row_data, col_data);
    // assemble local constrains matrix
    CHKERR aSsemble(row_data, col_data);
    MoFEMFunctionReturn(0);
  }

iNtegrate()

MoFEMErrorCode PoissonExample::OpC::iNtegrate ( EntitiesFieldData::EntData &  row_data, EntitiesFieldData::EntData &  col_data  )

inlineprivate

Integrate local constrains matrix.

Examples

PoissonOperators.hpp.

Definition at line 345 of file PoissonOperators.hpp.

                                                                      {
    MoFEMFunctionBegin;
    // set size of local constrains matrix
    locMat.resize(nbRows, nbCols, false);
    // clear matrix
    locMat.clear();
    // get area of element
    const double area = getArea();
    // get integration weights
    auto t_w = getFTensor0IntegrationWeight();
    // get base functions on entity
    auto t_row = row_data.getFTensor0N();
    // run over integration points
    for (int gg = 0; gg != nbIntegrationPts; gg++) {
      const double alpha = area * t_w;
      // get element of local matrix
      FTensor::Tensor0<FTensor::PackPtr<double *, 1>> c(
          &*locMat.data().begin());
      // run over base functions on rows
      for (int rr = 0; rr != nbRows; rr++) {
        // get first base functions on column for integration point gg
        auto t_col = col_data.getFTensor0N(gg, 0);
        // run over base function on column
        for (int cc = 0; cc != nbCols; cc++) {
          // integrate element of constrains matrix
          c += alpha * t_row * t_col;
          ++t_col; // move to next base function on column
          ++c;     // move to next element of local matrix
        }
        ++t_row; // move to next base function on row
      }
      ++t_w; // move to next integrate weight
    }
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

assembleTranspose

const bool PoissonExample::OpC::assembleTranspose

private

assemble transpose, i.e. CT if set to true

Examples

PoissonOperators.hpp.

Definition at line 339 of file PoissonOperators.hpp.

locMat

MatrixDouble PoissonExample::OpC::locMat

private

local constrains matrix

Examples

PoissonOperators.hpp.

Definition at line 341 of file PoissonOperators.hpp.

nbCols

int PoissonExample::OpC::nbCols

private

number of dofs on column

Examples

PoissonOperators.hpp.

Definition at line 337 of file PoissonOperators.hpp.

nbIntegrationPts

int PoissonExample::OpC::nbIntegrationPts

private

number of integration points

Examples

PoissonOperators.hpp.

Definition at line 338 of file PoissonOperators.hpp.

nbRows

int PoissonExample::OpC::nbRows

private

< error code

number of dofs on row

Examples

PoissonOperators.hpp.

Definition at line 336 of file PoissonOperators.hpp.

---

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

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