OpTauStabilizationOpAnalyticalDispBc Struct Reference

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

Inheritance diagram for OpTauStabilizationOpAnalyticalDispBc:

Collaboration diagram for OpTauStabilizationOpAnalyticalDispBc:

Public Types
using	OP = OpBrokenBaseTimesBrokenDisp
Public Types inherited from OpBrokenBaseTimesBrokenDisp
using	OP = OpStabBrokenBaseImpl< OpBaseTimesVectorFace >
Public Types inherited from OpStabBrokenBaseImpl< OpBaseTimesVectorFace >
using	BASE = OpBrokenBaseImpl< OpBaseTimesVectorFace >

Public Member Functions
	OpTauStabilizationOpAnalyticalDispBc (boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_base_side_data, boost::shared_ptr< AnalyticalDisplacementBcVec > &bc_disp_ptr, std::map< std::string, boost::shared_ptr< ScalingMethod > > smv, ScalarFun tau_coeff, boost::shared_ptr< Range > ents_ptr=nullptr)
Public Member Functions inherited from OpBrokenBaseTimesBrokenDisp
	OpBrokenBaseTimesBrokenDisp (boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_base_side_data, ScalarFun beta_coeff=[](double, double, double) constexpr { return 1;}, boost::shared_ptr< Range > ents_ptr=nullptr)
Public Member Functions inherited from OpStabBrokenBaseImpl< OpBaseTimesVectorFace >
MoFEMErrorCode	doWork (int row_side, EntityType row_type, EntitiesFieldData::EntData &row_data)

Protected Member Functions
MoFEMErrorCode	iNtegrate (EntData &data)
Protected Member Functions inherited from OpBrokenBaseTimesBrokenDisp
MoFEMErrorCode	doWork (int row_side, EntityType row_type, EntitiesFieldData::EntData &row_data)

Protected Attributes
boost::shared_ptr< AnalyticalDisplacementBcVec >	bcDispPtr
std::map< std::string, boost::shared_ptr< ScalingMethod > >	scalingMethodsMap
Protected Attributes inherited from OpStabBrokenBaseImpl< OpBaseTimesVectorFace >
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 1280 of file EshelbianOperators.hpp.

Member Typedef Documentation

OP

using OpTauStabilizationOpAnalyticalDispBc::OP = OpBrokenBaseTimesBrokenDisp

Definition at line 1283 of file EshelbianOperators.hpp.

Constructor & Destructor Documentation

OpTauStabilizationOpAnalyticalDispBc()

OpTauStabilizationOpAnalyticalDispBc::OpTauStabilizationOpAnalyticalDispBc ( boost::shared_ptr< std::vector< BrokenBaseSideData > >  broken_base_side_data, boost::shared_ptr< AnalyticalDisplacementBcVec > &  bc_disp_ptr, std::map< std::string, boost::shared_ptr< ScalingMethod > >  smv, ScalarFun  tau_coeff, boost::shared_ptr< Range >  ents_ptr = nullptr  )

inline

Definition at line 1285 of file EshelbianOperators.hpp.

      : OpBrokenBaseTimesBrokenDisp(broken_base_side_data, tau_coeff, ents_ptr),
        bcDispPtr(bc_disp_ptr), scalingMethodsMap(smv) {}

Member Function Documentation

iNtegrate()

MoFEMErrorCode OpTauStabilizationOpAnalyticalDispBc::iNtegrate ( EntData &  data )

protected

Examples

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

Definition at line 1760 of file EshelbianOperators.cpp.

                                                                            {
  MoFEMFunctionBegin;
 
  EntityHandle fe_ent = OP::getFEEntityHandle();
  for (auto &bc : (*bcDispPtr)) {
    if (bc.faces.find(fe_ent) != bc.faces.end()) {
      MatrixDouble analytical_expr;
      auto analytical_data = getAnalyticalExpr(this, analytical_expr, bc.blockName);
      auto &v_analytical_expr = std::get<1>(analytical_data);
 
      int nb_dofs = data.getIndices().size();
      int nb_integration_pts = OP::getGaussPts().size2();
      auto t_w = OP::getFTensor0IntegrationWeight();
      int nb_base_functions = data.getN().size2();
      auto t_row_base_fun = data.getFTensor0N();
 
#ifndef NDEBUG
      if (!this->sourceVec) {
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                "Source vector for OpTauStabilizationOpAnalyticalDispBc is not "
                "set");
      }
      if (data.getN().size1() != nb_integration_pts) {
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                "Number of integration points in data should be %d but is %d",
                nb_integration_pts, (int)data.getN().size1());
      }
      if (nb_base_functions < nb_dofs / SPACE_DIM) {
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                "Number of base functions in data should be at least %d but is "
                "%d",
                nb_dofs / SPACE_DIM, nb_base_functions);
      }
#endif
 
      auto t_disp_val =
          MatrixSizeHelper<GetFTensor1FromMatType<3, -1, DL>, DL>::get(
              *this->sourceVec, nb_integration_pts)();
      auto t_bc_disp = getFTensor1FromMat<3, -1, DL>(v_analytical_expr);
 
      FTENSOR_INDEX(SPACE_DIM, i);
 
      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));
        FTensor::Tensor1<double, SPACE_DIM> t_bc_residual;
        for (int ii = 0; ii != SPACE_DIM; ++ii)
          t_bc_residual(ii) =
              bc.flags[ii] ? t_disp_val(ii) - t_bc_disp(ii) : 0.;
        auto t_nf = getFTensor1FromPtr<3, 3>(OP::locF.data().data());
        int bb = 0;
        for (; bb != nb_dofs / SPACE_DIM; ++bb) {
          t_nf(i) +=
              (tau_scale * t_row_base_fun) * t_bc_residual(i);
          ++t_nf;
          ++t_row_base_fun;
        }
        for (; bb != nb_base_functions; ++bb)
          ++t_row_base_fun;
 
        ++t_w;
        ++t_coords;
        ++t_disp_val;
        ++t_bc_disp;
      }
    }
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

bcDispPtr

boost::shared_ptr<AnalyticalDisplacementBcVec> OpTauStabilizationOpAnalyticalDispBc::bcDispPtr

protected

Definition at line 1295 of file EshelbianOperators.hpp.

scalingMethodsMap

std::map<std::string, boost::shared_ptr<ScalingMethod> > OpTauStabilizationOpAnalyticalDispBc::scalingMethodsMap

protected

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