v0.14.0
Public Member Functions | Private Attributes | List of all members
HenckyOps::OpCalculateHenckyStressImpl< DIM, GAUSS, DomainEleOp, S > Struct Template Reference

#include <tutorials/vec-2/src/HenckyOps.hpp>

Inheritance diagram for HenckyOps::OpCalculateHenckyStressImpl< DIM, GAUSS, DomainEleOp, S >:
[legend]
Collaboration diagram for HenckyOps::OpCalculateHenckyStressImpl< DIM, GAUSS, DomainEleOp, S >:
[legend]

Public Member Functions

 OpCalculateHenckyStressImpl (const std::string field_name, boost::shared_ptr< CommonData > common_data)
 
MoFEMErrorCode doWork (int side, EntityType type, EntData &data)
 

Private Attributes

boost::shared_ptr< CommonDatacommonDataPtr
 

Detailed Description

template<int DIM, typename DomainEleOp, int S>
struct HenckyOps::OpCalculateHenckyStressImpl< DIM, GAUSS, DomainEleOp, S >

Definition at line 291 of file HenckyOps.hpp.

Constructor & Destructor Documentation

◆ OpCalculateHenckyStressImpl()

template<int DIM, typename DomainEleOp , int S>
HenckyOps::OpCalculateHenckyStressImpl< DIM, GAUSS, DomainEleOp, S >::OpCalculateHenckyStressImpl ( const std::string  field_name,
boost::shared_ptr< CommonData common_data 
)
inline

Definition at line 294 of file HenckyOps.hpp.

296  : DomainEleOp(field_name, DomainEleOp::OPROW),
297  commonDataPtr(common_data) {
298  std::fill(&DomainEleOp::doEntities[MBEDGE],
299  &DomainEleOp::doEntities[MBMAXTYPE], false);
300  }

Member Function Documentation

◆ doWork()

template<int DIM, typename DomainEleOp , int S>
MoFEMErrorCode HenckyOps::OpCalculateHenckyStressImpl< DIM, GAUSS, DomainEleOp, S >::doWork ( int  side,
EntityType  type,
EntData data 
)
inline

Definition at line 302 of file HenckyOps.hpp.

302  {
304 
309 
310  // const size_t nb_gauss_pts = matGradPtr->size2();
311  const size_t nb_gauss_pts = DomainEleOp::getGaussPts().size2();
312  auto t_D = getFTensor4DdgFromMat<DIM, DIM, S>(*commonDataPtr->matDPtr);
313  auto t_logC = getFTensor2SymmetricFromMat<DIM>(commonDataPtr->matLogC);
314  constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
315  commonDataPtr->matHenckyStress.resize(size_symm, nb_gauss_pts, false);
316  auto t_T = getFTensor2SymmetricFromMat<DIM>(commonDataPtr->matHenckyStress);
317 
318  for (size_t gg = 0; gg != nb_gauss_pts; ++gg) {
319  t_T(i, j) = t_D(i, j, k, l) * t_logC(k, l);
320  ++t_logC;
321  ++t_T;
322  ++t_D;
323  }
324 
326  }

Member Data Documentation

◆ commonDataPtr

template<int DIM, typename DomainEleOp , int S>
boost::shared_ptr<CommonData> HenckyOps::OpCalculateHenckyStressImpl< DIM, GAUSS, DomainEleOp, S >::commonDataPtr
private

Definition at line 329 of file HenckyOps.hpp.


The documentation for this struct was generated from the following file:
DomainEleOp
DomainEle::UserDataOperator DomainEleOp
Finire element operator type.
Definition: child_and_parent.cpp:36
size_symm
constexpr auto size_symm
Definition: plastic.cpp:42
i
FTensor::Index< 'i', SPACE_DIM > i
Definition: hcurl_divergence_operator_2d.cpp:27
HenckyOps::OpCalculateHenckyStressImpl< DIM, GAUSS, DomainEleOp, S >::commonDataPtr
boost::shared_ptr< CommonData > commonDataPtr
Definition: HenckyOps.hpp:329
field_name
constexpr auto field_name
Definition: poisson_2d_homogeneous.cpp:13
FTensor::Index< 'i', DIM >
j
FTensor::Index< 'j', 3 > j
Definition: matrix_function.cpp:19
k
FTensor::Index< 'k', 3 > k
Definition: matrix_function.cpp:20
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:429
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:359
l
FTensor::Index< 'l', 3 > l
Definition: matrix_function.cpp:21