NonlinearPoissonOps::OpDomainRhs Struct Reference

#include "tutorials/scl-4/src/nonlinear_poisson_2d.hpp"

Inheritance diagram for NonlinearPoissonOps::OpDomainRhs:

Collaboration diagram for NonlinearPoissonOps::OpDomainRhs:

Public Member Functions
	OpDomainRhs (std::string field_name, ScalarFunc source_term_function, boost::shared_ptr< VectorDouble > field_vec, boost::shared_ptr< MatrixDouble > field_grad_mat)
MoFEMErrorCode	iNtegrate (EntData &data)

Private Attributes
ScalarFunc	sourceTermFunc
boost::shared_ptr< VectorDouble >	fieldVec
boost::shared_ptr< MatrixDouble >	fieldGradMat

Detailed Description

Definition at line 119 of file nonlinear_poisson_2d.hpp.

Constructor & Destructor Documentation

OpDomainRhs()

NonlinearPoissonOps::OpDomainRhs::OpDomainRhs ( std::string field_name, ScalarFunc source_term_function, boost::shared_ptr< VectorDouble > field_vec, boost::shared_ptr< MatrixDouble > field_grad_mat )

inline

Definition at line 121 of file nonlinear_poisson_2d.hpp.

      : AssemblyDomainEleOp(field_name, field_name, AssemblyDomainEleOp::OPROW),
        sourceTermFunc(source_term_function), fieldVec(field_vec),
        fieldGradMat(field_grad_mat) {}

Member Function Documentation

iNtegrate()

MoFEMErrorCode NonlinearPoissonOps::OpDomainRhs::iNtegrate ( EntData & data )

inline

Definition at line 129 of file nonlinear_poisson_2d.hpp.

                                          {
    MoFEMFunctionBegin;
 
    auto &nf = AssemblyDomainEleOp::locF;
 
 
    // 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 coordinates of the integration point
    auto t_coords = getFTensor1CoordsAtGaussPts();
    // get solution (field value) at integration point
    auto t_field = getFTensor0FromVec(*fieldVec);
    // get gradient of field value of integration point
    auto t_field_grad = getFTensor1FromMat<2>(*fieldGradMat);
 
    // get base function
    auto t_base = data.getFTensor0N();
    // get derivatives of base function
    auto t_grad_base = data.getFTensor1DiffN<2>();
 
    // 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;
      double body_source =
          sourceTermFunc(t_coords(0), t_coords(1), t_coords(2));
 
      // calculate the local vector
      for (int rr = 0; rr != AssemblyDomainEleOp::nbRows; rr++) {
        nf[rr] +=
            (-t_base * body_source +
              t_grad_base(i) * t_field_grad(i) * (1 + t_field * t_field)) *
            a;
 
        // move to the next base function
        ++t_base;
        // move to the derivatives of the next base function
        ++t_grad_base;
      }
 
      // move to the weight of the next integration point
      ++t_w;
      // move to the coordinates of the next integration point
      ++t_coords;
      // move to the solution (field value) at the next integration point
      ++t_field;
      // move to the gradient of field value at the next integration point
      ++t_field_grad;
 
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

fieldGradMat

boost::shared_ptr<MatrixDouble> NonlinearPoissonOps::OpDomainRhs::fieldGradMat

private

Definition at line 190 of file nonlinear_poisson_2d.hpp.

fieldVec

boost::shared_ptr<VectorDouble> NonlinearPoissonOps::OpDomainRhs::fieldVec

private

Definition at line 189 of file nonlinear_poisson_2d.hpp.

sourceTermFunc

ScalarFunc NonlinearPoissonOps::OpDomainRhs::sourceTermFunc

private

Definition at line 188 of file nonlinear_poisson_2d.hpp.

---

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

• tutorials/scl-4/src/nonlinear_poisson_2d.hpp