Lhs for H dU. More...

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

Inheritance diagram for FreeSurfaceOps::OpLhsH_dU:

Collaboration diagram for FreeSurfaceOps::OpLhsH_dU:

Public Member Functions
	OpLhsH_dU (const std::string h_field_name, const std::string u_field_name, boost::shared_ptr< MatrixDouble > grad_h_ptr)

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

Private Attributes
boost::shared_ptr< MatrixDouble >	gradHPtr

Detailed Description

Lhs for H dU.

Parameters

h_field_name	Name of the field associated with the row operator (H)
u_field_name	Name of the field associated with the column operator (U)
grad_h_ptr	Pointer to the free surface height gradient matrix (∇H)

Returns: MoFEMErrorCode

(Jacobian block ∂R_H/∂U) Linearization of the conserved phase equation (Lovric 3.1c) w.r.t. U

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

Definition at line 1091 of file FreeSurfaceOps.hpp.

Constructor & Destructor Documentation

◆ OpLhsH_dU()

FreeSurfaceOps::OpLhsH_dU::OpLhsH_dU	(	const std::string	h_field_name,
		const std::string	u_field_name,
		boost::shared_ptr< MatrixDouble >	grad_h_ptr
	)

inline

Definition at line 1092 of file FreeSurfaceOps.hpp.

      : AssemblyDomainEleOp(h_field_name, u_field_name,
                            AssemblyDomainEleOp::OPROWCOL),
        gradHPtr(grad_h_ptr) {
    sYmm = false;
    assembleTranspose = false;
  }

Member Function Documentation

◆ iNtegrate()

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

inline

Definition at line 1100 of file FreeSurfaceOps.hpp.

                                                               {
    MoFEMFunctionBegin;
 
    const double vol = getMeasure();
    auto t_grad_h = getFTensor1FromMat<SPACE_DIM>(*gradHPtr);
    auto t_coords = getFTensor1CoordsAtGaussPts();
 
    auto t_row_base = row_data.getFTensor0N();
    auto t_w = getFTensor0IntegrationWeight();
 
    for (int gg = 0; gg != nbIntegrationPts; gg++) {
 
      const auto r = t_coords(0);
      const auto alpha = t_w * vol * cylindrical(r);
      auto t_mat = getFTensor1FromPtr<U_FIELD_DIM>(&locMat(0, 0));
      FTensor::Tensor1<double, SPACE_DIM> t_row;
 
      int rr = 0;
      for (; rr != nbRows; ++rr) {
        t_row(i) = (alpha * t_row_base) * t_grad_h(i);
        auto t_col_base = col_data.getFTensor0N(gg, 0);
        for (int cc = 0; cc != nbCols / U_FIELD_DIM; ++cc) {
          t_mat(i) += t_row(i) * t_col_base;
          ++t_mat;
          ++t_col_base;
        }
        ++t_row_base;
      }
 
      for (; rr < nbRowBaseFunctions; ++rr)
        ++t_row_base;
 
      ++t_grad_h;
      ++t_w;
      ++t_coords;
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

◆ gradHPtr

boost::shared_ptr<MatrixDouble> FreeSurfaceOps::OpLhsH_dU::gradHPtr

private

Definition at line 1142 of file FreeSurfaceOps.hpp.

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

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

Public Member Functions