OpTauStabilizationDispLhsBc Struct Reference

#include "users_modules/eshelbian_plasticity/src/EshelbianOperators.hpp"

Inheritance diagram for OpTauStabilizationDispLhsBc:

Collaboration diagram for OpTauStabilizationDispLhsBc:

Public Types
using	OP = OpBrokenBaseBrokenBase
Public Types inherited from OpBrokenBaseBrokenBase
using	OP = OpStabBrokenBaseImpl< OpMassVectorFace >
Public Types inherited from OpStabBrokenBaseImpl< OpMassVectorFace >
using	BASE = OpBrokenBaseImpl< OpMassVectorFace >

Public Member Functions
	OpTauStabilizationDispLhsBc (boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_base_side_data, boost::shared_ptr< BcDispVec > &bc_disp_ptr, ScalarFun tau_coeff, boost::shared_ptr< Range > ents_ptr=nullptr)
Public Member Functions inherited from OpBrokenBaseBrokenBase
	OpBrokenBaseBrokenBase (boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_base_side_data, ScalarFun beta, boost::shared_ptr< Range > ents_ptr=nullptr)
MoFEMErrorCode	doWork (int row_side, EntityType row_type, EntitiesFieldData::EntData &row_data)
Public Member Functions inherited from OpStabBrokenBaseImpl< OpMassVectorFace >
MoFEMErrorCode	doWork (int row_side, EntityType row_type, EntitiesFieldData::EntData &row_data)

Protected Member Functions
MoFEMErrorCode	iNtegrate (EntData &row_data, EntData &col_data)

Protected Attributes
boost::shared_ptr< BcDispVec >	bcDispPtr
Protected Attributes inherited from OpStabBrokenBaseImpl< OpMassVectorFace >
boost::weak_ptr< MatrixDouble >	fluxMatPtr

Detailed Description

Examples

/home/lk58p/mofem_install/vanilla_dev_release/mofem-cephas/mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianPlasticity.cpp.

Definition at line 1320 of file EshelbianOperators.hpp.

Member Typedef Documentation

OP

using OpTauStabilizationDispLhsBc::OP = OpBrokenBaseBrokenBase

Definition at line 1322 of file EshelbianOperators.hpp.

Constructor & Destructor Documentation

OpTauStabilizationDispLhsBc()

OpTauStabilizationDispLhsBc::OpTauStabilizationDispLhsBc ( boost::shared_ptr< std::vector< BrokenBaseSideData > >  broken_base_side_data, boost::shared_ptr< BcDispVec > &  bc_disp_ptr, ScalarFun  tau_coeff, boost::shared_ptr< Range >  ents_ptr = nullptr  )

inline

Definition at line 1324 of file EshelbianOperators.hpp.

      : OP(broken_base_side_data, tau_coeff, ents_ptr), bcDispPtr(bc_disp_ptr) {
  }

Member Function Documentation

iNtegrate()

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

protected

Examples

/home/lk58p/mofem_install/vanilla_dev_release/mofem-cephas/mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp.

Definition at line 1621 of file EshelbianOperators.cpp.

                                                                         {
  MoFEMFunctionBegin;
  // get entity of face
  EntityHandle fe_ent = OP::getFEEntityHandle();
  // iterate over all boundary data
  for (auto &bc : (*bcDispPtr)) {
    // check if finite element entity is part of boundary condition
    if (bc.faces.find(fe_ent) != bc.faces.end()) {
      int nb_dofs = row_data.getIndices().size();
 
      int nb_integration_pts = OP::getGaussPts().size2();
      auto t_w = OP::getFTensor0IntegrationWeight();
      int nb_base_functions = row_data.getN().size2();
      auto t_row_base_fun = row_data.getFTensor0N();
 
      FTENSOR_INDEX(SPACE_DIM, i);
      FTENSOR_INDEX(SPACE_DIM, j);
 
      FTensor::Tensor1<double, SPACE_DIM> t_bc_mask;
      for (int ii = 0; ii != SPACE_DIM; ++ii)
        t_bc_mask(ii) = bc.flags[ii] ? 1. : 0.;
 
      auto get_t_vec = [&](const int rr) {
        std::array<double *, SPACE_DIM> ptrs;
        for (auto i = 0; i != SPACE_DIM; ++i)
          ptrs[i] = &OP::locMat(rr + i, i);
        return FTensor::Tensor1<FTensor::PackPtr<double *, SPACE_DIM>,
                                SPACE_DIM>(ptrs);
      };
 
      auto area = getMeasure();
      auto t_coords = OP::getFTensor1CoordsAtGaussPts();
      for (int gg = 0; gg != nb_integration_pts; ++gg) {
        auto tau_scale =
            area * t_w * OP::betaCoeff(t_coords(0), t_coords(1), t_coords(2));
        int rr = 0;
        for (; rr != nb_dofs / SPACE_DIM; ++rr) {
          auto t_mat = get_t_vec(SPACE_DIM * rr);
          auto t_col_base_fun = col_data.getFTensor0N(gg, 0);
          for (int cc = 0; cc != nb_dofs / SPACE_DIM; ++cc) {
            t_mat(i) +=
                tau_scale * (t_row_base_fun * t_col_base_fun) * t_bc_mask(i);
            ++t_col_base_fun;
            ++t_mat;
          }
          ++t_row_base_fun;
        }
        for (; rr != nb_base_functions; ++rr)
          ++t_row_base_fun;
 
        ++t_w;
        ++t_coords;
      }
    }
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

bcDispPtr

boost::shared_ptr<BcDispVec> OpTauStabilizationDispLhsBc::bcDispPtr

protected

Definition at line 1333 of file EshelbianOperators.hpp.

---

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

• users_modules/eshelbian_plasticity/src/EshelbianOperators.hpp
• users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp