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

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

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

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 1236 of file EshelbianOperators.hpp.

Member Typedef Documentation

◆ OP

Definition at line 1239 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 1241 of file EshelbianOperators.hpp.

1246 : OpBrokenBaseTimesBrokenDisp(broken_base_side_data, tau_coeff, ents_ptr),
1247 bcDispPtr(bc_disp_ptr), scalingMethodsMap(smv) {}
boost::shared_ptr< AnalyticalDisplacementBcVec > bcDispPtr
std::map< std::string, boost::shared_ptr< ScalingMethod > > scalingMethodsMap

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 1416 of file EshelbianOperators.cpp.

1416 {
1418
1419 EntityHandle fe_ent = OP::getFEEntityHandle();
1420 for (auto &bc : (*bcDispPtr)) {
1421 if (bc.faces.find(fe_ent) != bc.faces.end()) {
1423 auto analytical_data = getAnalyticalExpr(this, analytical_expr, bc.blockName);
1424 auto &v_analytical_expr = std::get<1>(analytical_data);
1425
1426 int nb_dofs = data.getIndices().size();
1427 int nb_integration_pts = OP::getGaussPts().size2();
1428 auto t_w = OP::getFTensor0IntegrationWeight();
1429 int nb_base_functions = data.getN().size2();
1430 auto t_row_base_fun = data.getFTensor0N();
1431
1432#ifndef NDEBUG
1433 if (!this->sourceVec) {
1434 SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
1435 "Source vector for OpTauStabilizationOpAnalyticalDispBc is not "
1436 "set");
1437 }
1438 if (data.getN().size1() != nb_integration_pts) {
1439 SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
1440 "Number of integration points in data should be %d but is %d",
1441 nb_integration_pts, (int)data.getN().size1());
1442 }
1443 if (nb_base_functions < nb_dofs / SPACE_DIM) {
1444 SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
1445 "Number of base functions in data should be at least %d but is "
1446 "%d",
1447 nb_dofs / SPACE_DIM, nb_base_functions);
1448 }
1449#endif
1450
1451 auto t_disp_val =
1453 *this->sourceVec, nb_integration_pts)();
1454 auto t_bc_disp = getFTensor1FromMat<3, -1, DL>(v_analytical_expr);
1455
1456 auto area = getMeasure();
1457 auto t_coords = OP::getFTensor1CoordsAtGaussPts();
1458 for (int gg = 0; gg != nb_integration_pts; ++gg) {
1459 auto tau_scale =
1460 area * t_w * OP::betaCoeff(t_coords(0), t_coords(1), t_coords(2));
1461 auto t_nf = getFTensor1FromPtr<3, 3>(OP::locF.data().data());
1462 int bb = 0;
1463 for (; bb != nb_dofs / SPACE_DIM; ++bb) {
1464 for (auto ii = 0; ii != SPACE_DIM; ++ii) {
1465 if (bc.flags[ii]) {
1466 t_nf(ii) +=
1467 (tau_scale * t_row_base_fun) * (t_disp_val(ii) - t_bc_disp(ii));
1468 }
1469 }
1470 ++t_nf;
1471 ++t_row_base_fun;
1472 }
1473 for (; bb != nb_base_functions; ++bb)
1474 ++t_row_base_fun;
1475
1476 ++t_w;
1477 ++t_coords;
1478 ++t_disp_val;
1479 ++t_bc_disp;
1480 }
1481 }
1482 }
1483
1485}
constexpr int SPACE_DIM
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
@ MOFEM_DATA_INCONSISTENCY
Definition definitions.h:31
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
std::tuple< std::string, MatrixDouble > getAnalyticalExpr(OP_PTR op_ptr, MatrixDouble &analytical_expr, const std::string block_name)
DataLayoutTraits< DataLayout::GaussByCoeffs > DL
Definition MatHuHu.hpp:33
auto getFTensor1FromMat(M &data, int rr=0, int cc=0)
Get tensor rank 1 (vector) form data matrix.
decltype(GetFTensor1FromMatImpl< Tensor_Dim, S, DL, M >::get(std::declval< M & >(), 0, 0)) GetFTensor1FromMatType
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> OpTauStabilizationOpAnalyticalDispBc::bcDispPtr
protected

Definition at line 1251 of file EshelbianOperators.hpp.

◆ scalingMethodsMap

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

Definition at line 1252 of file EshelbianOperators.hpp.


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