Poisson2DiscontGalerkinOperators::OpCalculateSideData Struct Reference

Operator tp collect data from elements on the side of Edge/Face. More...

#include "tutorials/scl-11/src/PoissonDiscontinousGalerkin.hpp"

Inheritance diagram for Poisson2DiscontGalerkinOperators::OpCalculateSideData:

Collaboration diagram for Poisson2DiscontGalerkinOperators::OpCalculateSideData:

Public Member Functions
	OpCalculateSideData (std::string field_name, std::string col_field_name)
MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntData &row_data, EntData &col_data)

Detailed Description

Operator tp collect data from elements on the side of Edge/Face.

Definition at line 38 of file PoissonDiscontinousGalerkin.hpp.

Constructor & Destructor Documentation

OpCalculateSideData()

Poisson2DiscontGalerkinOperators::OpCalculateSideData::OpCalculateSideData ( std::string  field_name, std::string  col_field_name  )

inline

Definition at line 40 of file PoissonDiscontinousGalerkin.hpp.

      : FaceSideOp(field_name, col_field_name, FaceSideOp::OPROWCOL) {
 
    std::fill(&doEntities[MBVERTEX], &doEntities[MBMAXTYPE], false);
 
    for (auto t = moab::CN::TypeDimensionMap[SPACE_DIM].first;
         t <= moab::CN::TypeDimensionMap[SPACE_DIM].second; ++t)
      doEntities[t] = true;
  }

Member Function Documentation

doWork()

MoFEMErrorCode Poisson2DiscontGalerkinOperators::OpCalculateSideData::doWork ( int  row_side, int  col_side, EntityType  row_type, EntityType  col_type, EntData &  row_data, EntData &  col_data  )

inline

Examples

PoissonDiscontinousGalerkin.hpp.

Definition at line 50 of file PoissonDiscontinousGalerkin.hpp.

                                           {
    MoFEMFunctionBeginHot;
 
    // Note: THat for L2 base data rows, and columns are the same, so operator
    // above can be simpler operator for the right hand side, and data can be
    // stored only for rows, since for columns data are the same. However for
    // complex multi-physics problems that not necessary would be a case. For
    // generality, we keep generic case.
 
    if ((CN::Dimension(row_type) == SPACE_DIM) &&
        (CN::Dimension(col_type) == SPACE_DIM)) {
      const auto nb_in_loop = getFEMethod()->nInTheLoop;
      indicesRowSideMap[nb_in_loop] = row_data.getIndices();
      indicesColSideMap[nb_in_loop] = col_data.getIndices();
      rowBaseSideMap[nb_in_loop] = row_data.getN();
      colBaseSideMap[nb_in_loop] = col_data.getN();
      rowDiffBaseSideMap[nb_in_loop] = row_data.getDiffN();
      colDiffBaseSideMap[nb_in_loop] = col_data.getDiffN();
      areaMap[nb_in_loop] = getMeasure();
      senseMap[nb_in_loop] = getSkeletonSense();
      if (!nb_in_loop) {
        indicesRowSideMap[1].clear();
        indicesColSideMap[1].clear();
        rowBaseSideMap[1].clear();
        colBaseSideMap[1].clear();
        rowDiffBaseSideMap[1].clear();
        colDiffBaseSideMap[1].clear();
        areaMap[1] = 0;
        senseMap[1] = 0;
      }
    } else {
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "Should not happen");
    }
 
    MoFEMFunctionReturnHot(0);
  }

---

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

• tutorials/scl-11/src/PoissonDiscontinousGalerkin.hpp