v0.16.0
Loading...
Searching...
No Matches
Public Types | Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
OpTauStabilizationOpAnalyticalDispLhsBc Struct Reference

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

Inheritance diagram for OpTauStabilizationOpAnalyticalDispLhsBc:
[legend]
Collaboration diagram for OpTauStabilizationOpAnalyticalDispLhsBc:
[legend]

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

 OpTauStabilizationOpAnalyticalDispLhsBc (boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_base_side_data, boost::shared_ptr< AnalyticalDisplacementBcVec > &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< AnalyticalDisplacementBcVec > 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 1255 of file EshelbianOperators.hpp.

Member Typedef Documentation

◆ OP

Definition at line 1257 of file EshelbianOperators.hpp.

Constructor & Destructor Documentation

◆ OpTauStabilizationOpAnalyticalDispLhsBc()

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

Definition at line 1259 of file EshelbianOperators.hpp.

1263 : OP(broken_base_side_data, tau_coeff, ents_ptr), bcDispPtr(bc_disp_ptr) {
1264 }
boost::shared_ptr< AnalyticalDisplacementBcVec > bcDispPtr

Member Function Documentation

◆ iNtegrate()

MoFEMErrorCode OpTauStabilizationOpAnalyticalDispLhsBc::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 1487 of file EshelbianOperators.cpp.

1488 {
1490
1491 EntityHandle fe_ent = OP::getFEEntityHandle();
1492 for (auto &bc : (*bcDispPtr)) {
1493 if (bc.faces.find(fe_ent) != bc.faces.end()) {
1494 int nb_dofs = row_data.getIndices().size();
1495 int nb_integration_pts = OP::getGaussPts().size2();
1496 auto t_w = OP::getFTensor0IntegrationWeight();
1497 int nb_base_functions = row_data.getN().size2();
1498 auto t_row_base_fun = row_data.getFTensor0N();
1499
1500 auto get_t_vec = [&](const int rr) {
1501 std::array<double *, SPACE_DIM> ptrs;
1502 for (auto i = 0; i != SPACE_DIM; ++i)
1503 ptrs[i] = &OP::locMat(rr + i, i);
1505 SPACE_DIM>(ptrs);
1506 };
1507
1508 auto area = getMeasure();
1509 auto t_coords = OP::getFTensor1CoordsAtGaussPts();
1510 for (int gg = 0; gg != nb_integration_pts; ++gg) {
1511 auto tau_scale =
1512 area * t_w * OP::betaCoeff(t_coords(0), t_coords(1), t_coords(2));
1513 int rr = 0;
1514 for (; rr != nb_dofs / SPACE_DIM; ++rr) {
1515 auto t_mat = get_t_vec(SPACE_DIM * rr);
1516 auto t_col_base_fun = col_data.getFTensor0N(gg, 0);
1517 for (int cc = 0; cc != nb_dofs / SPACE_DIM; ++cc) {
1518 for (int ii = 0; ii != SPACE_DIM; ++ii) {
1519 if (bc.flags[ii]) {
1520 t_mat(ii) += tau_scale * (t_row_base_fun * t_col_base_fun);
1521 }
1522 }
1523 ++t_col_base_fun;
1524 ++t_mat;
1525 }
1526 ++t_row_base_fun;
1527 }
1528 for (; rr != nb_base_functions; ++rr)
1529 ++t_row_base_fun;
1530
1531 ++t_w;
1532 ++t_coords;
1533 }
1534 }
1535 }
1536
1538}
constexpr int SPACE_DIM
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Tensor0< FTensor::PackPtr< double *, 1 > > getFTensor0N(const FieldApproximationBase base)
Get base function as Tensor0.
MatrixDouble & getN(const FieldApproximationBase base)
get base functions this return matrix (nb. of rows is equal to nb. of Gauss pts, nb....
const VectorInt & getIndices() const
Get global indices of degrees of freedom on entity.

Member Data Documentation

◆ bcDispPtr

boost::shared_ptr<AnalyticalDisplacementBcVec> OpTauStabilizationOpAnalyticalDispLhsBc::bcDispPtr
protected

Definition at line 1268 of file EshelbianOperators.hpp.


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