Poisson2DHomogeneousOperators::OpDomainRhsVectorF Struct Reference

#include "tutorials/scl-1/src/poisson_2d_homogeneous.hpp"

Inheritance diagram for Poisson2DHomogeneousOperators::OpDomainRhsVectorF:

Collaboration diagram for Poisson2DHomogeneousOperators::OpDomainRhsVectorF:

Public Member Functions
	OpDomainRhsVectorF (std::string field_name, ScalarFunc source_term_function)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &data)

Private Attributes
ScalarFunc	sourceTermFunc

Detailed Description

Examples

poisson_2d_homogeneous.cpp.

Definition at line 94 of file poisson_2d_homogeneous.hpp.

Constructor & Destructor Documentation

OpDomainRhsVectorF()

Poisson2DHomogeneousOperators::OpDomainRhsVectorF::OpDomainRhsVectorF ( std::string field_name, ScalarFunc source_term_function )

inline

Examples

poisson_2d_homogeneous.hpp.

Definition at line 96 of file poisson_2d_homogeneous.hpp.

      : AssemblyDomainEleOp(field_name, field_name, DomainEleOp::OPROW),
        sourceTermFunc(source_term_function) {}

Member Function Documentation

iNtegrate()

MoFEMErrorCode Poisson2DHomogeneousOperators::OpDomainRhsVectorF::iNtegrate ( EntitiesFieldData::EntData & data )

inline

Examples

poisson_2d_homogeneous.hpp.

Definition at line 100 of file poisson_2d_homogeneous.hpp.

                                                           {
    MoFEMFunctionBegin;
 
    const int nb_dofs = data.getIndices().size();
 
    this->locF.resize(nb_dofs, false);
    this->locF.clear();
 
    // get element area
    const double area = getMeasure();
 
    // get number of integration points
    const int nb_integration_points = getGaussPts().size2();
    // get integration weights
    auto t_w = getFTensor0IntegrationWeight();
    auto t_coords = getFTensor1CoordsAtGaussPts();
 
    // get base function
    auto t_base = data.getFTensor0N();
 
    // START THE LOOP OVER INTEGRATION POINTS TO CALCULATE LOCAL VECTOR
    for (int gg = 0; gg != nb_integration_points; gg++) {
      const double a = t_w * area;
 
      // get source term at the integration point coordinates
      double body_source =
          sourceTermFunc(t_coords(0), t_coords(1), t_coords(2));
      for (int rr = 0; rr != nb_dofs; rr++) {
        this->locF[rr] += t_base * body_source * a;
 
        // move to the next base function
        ++t_base;
      }
 
      // move to the weight of the next integration point
      ++t_w;
      // move to the coordinates of the next integration point
      ++t_coords;
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

sourceTermFunc

ScalarFunc Poisson2DHomogeneousOperators::OpDomainRhsVectorF::sourceTermFunc

private

Examples

poisson_2d_homogeneous.hpp.

Definition at line 144 of file poisson_2d_homogeneous.hpp.

---

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

• tutorials/scl-1/src/poisson_2d_homogeneous.hpp