Poisson2DHomogeneousOperators::OpDomainLhsMatrixK Struct Reference

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

Inheritance diagram for Poisson2DHomogeneousOperators::OpDomainLhsMatrixK:

Collaboration diagram for Poisson2DHomogeneousOperators::OpDomainLhsMatrixK:

Public Member Functions
	OpDomainLhsMatrixK (std::string row_field_name, std::string col_field_name)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

Detailed Description

Definition at line 40 of file poisson_2d_homogeneous.hpp.

Constructor & Destructor Documentation

OpDomainLhsMatrixK()

Poisson2DHomogeneousOperators::OpDomainLhsMatrixK::OpDomainLhsMatrixK ( std::string row_field_name, std::string col_field_name )

inline

Definition at line 42 of file poisson_2d_homogeneous.hpp.

      : AssemblyDomainEleOp(row_field_name, col_field_name,
                            DomainEleOp::OPROWCOL) {}

Member Function Documentation

iNtegrate()

MoFEMErrorCode Poisson2DHomogeneousOperators::OpDomainLhsMatrixK::iNtegrate ( EntitiesFieldData::EntData & row_data, EntitiesFieldData::EntData & col_data )

inline

Examples

poisson_2d_homogeneous.hpp.

Definition at line 46 of file poisson_2d_homogeneous.hpp.

                                                               {
    MoFEMFunctionBegin;
 
    const int nb_row_dofs = row_data.getIndices().size();
    const int nb_col_dofs = col_data.getIndices().size();
 
    this->locMat.resize(nb_row_dofs, nb_col_dofs, false);
    this->locMat.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();
 
    // get derivatives of base functions on row
    auto t_row_diff_base = row_data.getFTensor1DiffN<SPACE_DIM>();
 
    // START THE LOOP OVER INTEGRATION POINTS TO CALCULATE LOCAL MATRIX
    for (int gg = 0; gg != nb_integration_points; gg++) {
      const double a = t_w * area;
 
      for (int rr = 0; rr != nb_row_dofs; ++rr) {
        // get derivatives of base functions on column
        auto t_col_diff_base = col_data.getFTensor1DiffN<SPACE_DIM>(gg, 0);
 
        for (int cc = 0; cc != nb_col_dofs; cc++) {
          this->locMat(rr, cc) += t_row_diff_base(i) * t_col_diff_base(i) * a;
 
          // move to the derivatives of the next base functions on column
          ++t_col_diff_base;
        }
 
        // move to the derivatives of the next base functions on row
        ++t_row_diff_base;
      }
 
      // move to the weight of the next integration point
      ++t_w;
    }
 
    MoFEMFunctionReturn(0);
  }

---

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

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