FreeSurfaceOps::OpLhsG_dG Struct Reference

Lhs for G dG. More...

#include "tutorials/vec-5_free_surface/src/FreeSurfaceOps.hpp"

Inheritance diagram for FreeSurfaceOps::OpLhsG_dG:

Collaboration diagram for FreeSurfaceOps::OpLhsG_dG:

Public Member Functions
	OpLhsG_dG (const std::string field_name)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

Detailed Description

Lhs for G dG.

Parameters

field_name

Name of the field associated with the operator

Returns

MoFEMErrorCode

(Jacobian block ∂R_G/∂G) Linearization of the chemical potential equation (Lovric 3.1d) w.r.t. G

Examples

mofem/tutorials/vec-5_free_surface/free_surface.cpp.

Definition at line 1537 of file FreeSurfaceOps.hpp.

Constructor & Destructor Documentation

OpLhsG_dG()

FreeSurfaceOps::OpLhsG_dG::OpLhsG_dG ( const std::string  field_name )

inline

Definition at line 1539 of file FreeSurfaceOps.hpp.

      : AssemblyDomainEleOp(field_name, field_name,
                            AssemblyDomainEleOp::OPROWCOL) {
    sYmm = true;
  }

Member Function Documentation

iNtegrate()

MoFEMErrorCode FreeSurfaceOps::OpLhsG_dG::iNtegrate ( EntitiesFieldData::EntData &  row_data, EntitiesFieldData::EntData &  col_data  )

inline

Definition at line 1545 of file FreeSurfaceOps.hpp.

                                                               {
    MoFEMFunctionBegin;
 
    const double vol = getMeasure();
 
    auto t_row_base = row_data.getFTensor0N();
    auto t_w = getFTensor0IntegrationWeight();
    auto t_coords = getFTensor1CoordsAtGaussPts();
 
    for (int gg = 0; gg != nbIntegrationPts; ++gg) {
 
      const double r = t_coords(0);
      const double alpha = t_w * vol * cylindrical(r);
 
      int rr = 0;
      for (; rr != nbRows; ++rr) {
        auto t_col_base = col_data.getFTensor0N(gg, 0);
        const double beta = alpha * t_row_base;
        for (int cc = 0; cc != nbCols; ++cc) {
          locMat(rr, cc) += (t_col_base * beta);
          ++t_col_base;
        }
 
        ++t_row_base;
      }
 
      for (; rr < nbRowBaseFunctions; ++rr) {
        ++t_row_base;
      }
 
      ++t_w;
      ++t_coords;
    }
 
    MoFEMFunctionReturn(0);
  }

---

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

• tutorials/vec-5_free_surface/src/FreeSurfaceOps.hpp