OpDomainResidualVector Struct Reference

Integrate the domain residual vector (RHS) More...

Inheritance diagram for OpDomainResidualVector:

Collaboration diagram for OpDomainResidualVector:

Public Member Functions
	OpDomainResidualVector (std::string field_name, boost::shared_ptr< MatrixDouble > field_grad_mat)
MoFEMErrorCode	iNtegrate (EntData &data)

Private Attributes
boost::shared_ptr< MatrixDouble >	fieldGradMat

Detailed Description

Integrate the domain residual vector (RHS)

\[ \sum\limits_j {\left[ {\int\limits_{{\Omega _e}} {\left( {{a_n}\nabla {\phi _i} \cdot \nabla {\phi _j} - a_n^3\nabla {\phi _i}\left( {\nabla u \cdot \nabla {\phi _j}} \right)\nabla u} \right)d{\Omega _e}} } \right]\delta {U_j}} = \int\limits_{{\Omega _e}} {{\phi _i}fd{\Omega _e}} - \int\limits_{{\Omega _e}} {\nabla {\phi _i}{a_n}\nabla ud{\Omega _e}} \\ {a_n} = \frac{1}{{{{\left( {1 + {{\left| {\nabla u} \right|}^2}} \right)}^{\frac{1}{2}}}}} \]

Definition at line 118 of file minimal_surface_equation.cpp.

Constructor & Destructor Documentation

OpDomainResidualVector()

OpDomainResidualVector::OpDomainResidualVector ( std::string field_name, boost::shared_ptr< MatrixDouble > field_grad_mat )

inline

Definition at line 120 of file minimal_surface_equation.cpp.

      : AssemblyDomainEleOp(field_name, field_name, AssemblyDomainEleOp::OPROW),
        fieldGradMat(field_grad_mat) {}

Member Function Documentation

iNtegrate()

MoFEMErrorCode OpDomainResidualVector::iNtegrate ( EntData & data )

inline

Definition at line 125 of file minimal_surface_equation.cpp.

                                          {
    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 gradient of the field at integration points
    auto t_field_grad = getFTensor1FromMat<2>(*fieldGradMat);
 
    // get base functions
    auto t_base = data.getFTensor0N();
    // get derivatives of base functions
    auto t_diff_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;
      const double an = 1. / std::sqrt(1 + t_field_grad(i) * t_field_grad(i));
 
      for (int rr = 0; rr != AssemblyDomainEleOp::nbRows; rr++) {
 
        // calculate components of the local vector
        // remember to use -= here due to PETSc consideration of Residual Vec
        nf[rr] += (t_diff_base(i) * t_field_grad(i)) * an * a;
 
        // move to the next base function
        ++t_base;
        // move to the derivatives of the next base function
        ++t_diff_base;
      }
 
      // move to the weight of the next integration point
      ++t_w;
      // move to the gradient of the field at the next integration point
      ++t_field_grad;
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

fieldGradMat

boost::shared_ptr<MatrixDouble> OpDomainResidualVector::fieldGradMat

private

Definition at line 173 of file minimal_surface_equation.cpp.

---

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

• tutorials/scl-5/minimal_surface_equation.cpp