NonlinearPoissonOps::OpDomainLhs Struct Reference

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

Inheritance diagram for NonlinearPoissonOps::OpDomainLhs:

Collaboration diagram for NonlinearPoissonOps::OpDomainLhs:

Public Member Functions
	OpDomainLhs (std::string row_field_name, std::string col_field_name, boost::shared_ptr< VectorDouble > field_vec, boost::shared_ptr< MatrixDouble > field_grad_mat)
MoFEMErrorCode	iNtegrate (EntData &row_data, EntData &col_data)

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

Detailed Description

Definition at line 39 of file nonlinear_poisson_2d.hpp.

Constructor & Destructor Documentation

OpDomainLhs()

NonlinearPoissonOps::OpDomainLhs::OpDomainLhs ( std::string row_field_name, std::string col_field_name, boost::shared_ptr< VectorDouble > field_vec, boost::shared_ptr< MatrixDouble > field_grad_mat )

inline

Definition at line 41 of file nonlinear_poisson_2d.hpp.

      : AssemblyDomainEleOp(row_field_name, col_field_name, 
      DomainEleOp::OPROWCOL),
      fieldVec(field_vec), fieldGradMat(field_grad_mat) {}

Member Function Documentation

iNtegrate()

MoFEMErrorCode NonlinearPoissonOps::OpDomainLhs::iNtegrate ( EntData & row_data, EntData & col_data )

inline

Definition at line 49 of file nonlinear_poisson_2d.hpp.

                                           {
    MoFEMFunctionBegin;
 
    auto &locLhs = AssemblyDomainEleOp::locMat;
 
    const int nb_row_dofs = row_data.getIndices().size();
    const int nb_col_dofs = col_data.getIndices().size();
    // 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 solution (field value) at integration points
    auto t_field = getFTensor0FromVec(*fieldVec);
    // get gradient of field at integration points
    auto t_field_grad = getFTensor1FromMat<2>(*fieldGradMat);
 
    // get base functions on row
    auto t_row_base = row_data.getFTensor0N();
    // get derivatives of base functions on row
    auto t_row_diff_base = row_data.getFTensor1DiffN<2>();
 
    // 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 base functions on column
        auto t_col_base = col_data.getFTensor0N(gg, 0);
        // get derivatives of base functions on column
        auto t_col_diff_base = col_data.getFTensor1DiffN<2>(gg, 0);
 
        for (int cc = 0; cc != nb_col_dofs; cc++) {
          locLhs(rr, cc) += (((1 + t_field * t_field) * t_row_diff_base(i) *
                              t_col_diff_base(i)) +
                              (2.0 * t_field * t_field_grad(i) *
                              t_row_diff_base(i) * t_col_base)) *
                            a;
 
          // move to the next base functions on column
          ++t_col_base;
          // move to the derivatives of the next base function on column
          ++t_col_diff_base;
        }
 
        // move to the next base function on row
        ++t_row_base;
        // move to the derivatives of the next base function on row
        ++t_row_diff_base;
      }
 
      // move to the weight of the next integration point
      ++t_w;
      // 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::OpDomainLhs::fieldGradMat

private

Definition at line 116 of file nonlinear_poisson_2d.hpp.

fieldVec

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

private

Definition at line 115 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