v0.15.0
Loading...
Searching...
No Matches
MoFEM::OpCalculateHVecTensorTrace< Tensor_Dim, OpBase > Struct Template Reference

Calculate trace of vector (Hdiv/Hcurl) space. More...

#include "src/finite_elements/UserDataOperators.hpp"

Inheritance diagram for MoFEM::OpCalculateHVecTensorTrace< Tensor_Dim, OpBase >:
[legend]
Collaboration diagram for MoFEM::OpCalculateHVecTensorTrace< Tensor_Dim, OpBase >:
[legend]

Public Member Functions

 OpCalculateHVecTensorTrace (const std::string field_name, boost::shared_ptr< MatrixDouble > data_ptr, boost::shared_ptr< double > scale_ptr, const EntityType zero_type=MBEDGE, const int zero_side=0)
 
 OpCalculateHVecTensorTrace (const std::string field_name, boost::shared_ptr< MatrixDouble > data_ptr, const EntityType zero_type=MBEDGE, const int zero_side=0)
 
MoFEMErrorCode doWork (int side, EntityType type, EntitiesFieldData::EntData &data)
 
- Public Member Functions inherited from MoFEM::OpBaseImpl< A, EleOp >
 OpBaseImpl (const std::string row_field_name, const std::string col_field_name, const OpType type, boost::shared_ptr< Range > ents_ptr=nullptr)
 
MoFEMErrorCode doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntData &row_data, EntData &col_data)
 Do calculations for the left hand side.
 
MoFEMErrorCode doWork (int row_side, EntityType row_type, EntData &row_data)
 Do calculations for the right hand side.
 

Private Attributes

boost::shared_ptr< MatrixDoubledataPtr
 
boost::shared_ptr< doublescalePtr
 
const EntityHandle zeroType
 
const int zeroSide
 
FTensor::Index< 'i', Tensor_Dim > i
 
FTensor::Index< 'j', Tensor_Dim > j
 

Additional Inherited Members

- Public Types inherited from MoFEM::OpBaseImpl< A, EleOp >
using OpType = typename EleOp::OpType
 
using EntData = EntitiesFieldData::EntData
 
using MatSetValuesHook
 
- Public Attributes inherited from MoFEM::OpBaseImpl< A, EleOp >
TimeFun timeScalingFun
 assumes that time variable is set
 
FEFun feScalingFun
 assumes that time variable is set
 
boost::shared_ptr< RangeentsPtr
 Entities on which element is run.
 
- Static Public Attributes inherited from MoFEM::OpBaseImpl< A, EleOp >
static MatSetValuesHook matSetValuesHook
 
- Protected Member Functions inherited from MoFEM::OpBaseImpl< A, EleOp >
template<int DIM>
FTensor::Tensor1< FTensor::PackPtr< double *, DIM >, DIM > getNf ()
 
template<int DIM>
FTensor::Tensor2< FTensor::PackPtr< double *, DIM >, DIM, DIM > getLocMat (const int rr)
 
virtual MoFEMErrorCode iNtegrate (EntData &row_data, EntData &col_data)
 Integrate grad-grad operator.
 
virtual MoFEMErrorCode aSsemble (EntData &row_data, EntData &col_data, const bool trans)
 
virtual MoFEMErrorCode iNtegrate (EntData &data)
 Class dedicated to integrate operator.
 
virtual MoFEMErrorCode aSsemble (EntData &data)
 
virtual size_t getNbOfBaseFunctions (EntitiesFieldData::EntData &data)
 Get number of base functions.
 
- Protected Attributes inherited from MoFEM::OpBaseImpl< A, EleOp >
int nbRows
 number of dofs on rows
 
int nbCols
 number if dof on column
 
int nbIntegrationPts
 number of integration points
 
int nbRowBaseFunctions
 number or row base functions
 
int rowSide
 row side number
 
int colSide
 column side number
 
EntityType rowType
 row type
 
EntityType colType
 column type
 
bool assembleTranspose
 
bool onlyTranspose
 
MatrixDouble locMat
 local entity block matrix
 
MatrixDouble locMatTranspose
 local entity block matrix
 
VectorDouble locF
 local entity vector
 

Detailed Description

template<int Tensor_Dim, typename OpBase>
struct MoFEM::OpCalculateHVecTensorTrace< Tensor_Dim, OpBase >

Calculate trace of vector (Hdiv/Hcurl) space.

Template Parameters
Tensor_Dim
OpBase
Examples
EshelbianPlasticity.cpp, and tensor_divergence_operator.cpp.

Definition at line 3213 of file UserDataOperators.hpp.

Constructor & Destructor Documentation

◆ OpCalculateHVecTensorTrace() [1/2]

template<int Tensor_Dim, typename OpBase >
MoFEM::OpCalculateHVecTensorTrace< Tensor_Dim, OpBase >::OpCalculateHVecTensorTrace ( const std::string field_name,
boost::shared_ptr< MatrixDouble > data_ptr,
boost::shared_ptr< double > scale_ptr,
const EntityType zero_type = MBEDGE,
const int zero_side = 0 )
inline

Definition at line 3215 of file UserDataOperators.hpp.

3220 : OpBase(field_name, OpBase::OPROW), dataPtr(data_ptr),
3221 scalePtr(scale_ptr), zeroType(zero_type), zeroSide(zero_side) {
3222 if (!dataPtr)
3223 THROW_MESSAGE("Pointer is not set");
3224 }
#define THROW_MESSAGE(msg)
Throw MoFEM exception.
constexpr auto field_name
OpBaseImpl< PETSC, EdgeEleOp > OpBase
Definition radiation.cpp:29
boost::shared_ptr< MatrixDouble > dataPtr
boost::shared_ptr< double > scalePtr

◆ OpCalculateHVecTensorTrace() [2/2]

template<int Tensor_Dim, typename OpBase >
MoFEM::OpCalculateHVecTensorTrace< Tensor_Dim, OpBase >::OpCalculateHVecTensorTrace ( const std::string field_name,
boost::shared_ptr< MatrixDouble > data_ptr,
const EntityType zero_type = MBEDGE,
const int zero_side = 0 )
inline

Definition at line 3226 of file UserDataOperators.hpp.

3230 : OpCalculateHVecTensorTrace(field_name, data_ptr, nullptr, zero_type,
3231 zero_side) {}
OpCalculateHVecTensorTrace(const std::string field_name, boost::shared_ptr< MatrixDouble > data_ptr, boost::shared_ptr< double > scale_ptr, const EntityType zero_type=MBEDGE, const int zero_side=0)

Member Function Documentation

◆ doWork()

template<int Tensor_Dim, typename OpBase >
MoFEMErrorCode MoFEM::OpCalculateHVecTensorTrace< Tensor_Dim, OpBase >::doWork ( int side,
EntityType type,
EntitiesFieldData::EntData & data )
inline

Definition at line 3233 of file UserDataOperators.hpp.

3234 {
3236 const size_t nb_integration_points = OpBase::getGaussPts().size2();
3237 if (type == zeroType && side == 0) {
3238 dataPtr->resize(Tensor_Dim, nb_integration_points, false);
3239 dataPtr->clear();
3240 }
3241 const size_t nb_dofs = data.getFieldData().size();
3242 if (nb_dofs) {
3243 double scale_val = (scalePtr) ? *scalePtr : 1.0;
3244 auto t_normal = OpBase::getFTensor1NormalsAtGaussPts();
3245 const size_t nb_base_functions = data.getN().size2() / 3;
3246 auto t_base = data.getFTensor1N<3>();
3248 for (size_t gg = 0; gg != nb_integration_points; ++gg) {
3249 FTensor::Tensor1<double, Tensor_Dim> t_normalized_normal;
3250 t_normalized_normal(j) = t_normal(j);
3251 t_normalized_normal.normalize();
3252 auto t_dof = data.getFTensor1FieldData<Tensor_Dim>();
3253 size_t bb = 0;
3254 for (; bb != nb_dofs / Tensor_Dim; ++bb) {
3255 t_data(i) +=
3256 (scale_val * t_dof(i)) * (t_base(j) * t_normalized_normal(j));
3257 ++t_base;
3258 ++t_dof;
3259 }
3260 for (; bb < nb_base_functions; ++bb) {
3261 ++t_base;
3262 }
3263 ++t_data;
3264 ++t_normal;
3265 }
3266 }
3268 }
Tensor1< T, Tensor_Dim > normalize()
#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::Tensor1< FTensor::PackPtr< T *, S >, Tensor_Dim > getFTensor1FromMat(ublas::matrix< T, L, A > &data)
Get tensor rank 1 (vector) form data matrix.
FTensor::Index< 'j', Tensor_Dim > j
FTensor::Index< 'i', Tensor_Dim > i

Member Data Documentation

◆ dataPtr

template<int Tensor_Dim, typename OpBase >
boost::shared_ptr<MatrixDouble> MoFEM::OpCalculateHVecTensorTrace< Tensor_Dim, OpBase >::dataPtr
private

Definition at line 3271 of file UserDataOperators.hpp.

◆ i

template<int Tensor_Dim, typename OpBase >
FTensor::Index<'i', Tensor_Dim> MoFEM::OpCalculateHVecTensorTrace< Tensor_Dim, OpBase >::i
private

Definition at line 3275 of file UserDataOperators.hpp.

◆ j

template<int Tensor_Dim, typename OpBase >
FTensor::Index<'j', Tensor_Dim> MoFEM::OpCalculateHVecTensorTrace< Tensor_Dim, OpBase >::j
private

Definition at line 3276 of file UserDataOperators.hpp.

◆ scalePtr

template<int Tensor_Dim, typename OpBase >
boost::shared_ptr<double> MoFEM::OpCalculateHVecTensorTrace< Tensor_Dim, OpBase >::scalePtr
private

Definition at line 3272 of file UserDataOperators.hpp.

◆ zeroSide

template<int Tensor_Dim, typename OpBase >
const int MoFEM::OpCalculateHVecTensorTrace< Tensor_Dim, OpBase >::zeroSide
private

Definition at line 3274 of file UserDataOperators.hpp.

◆ zeroType

template<int Tensor_Dim, typename OpBase >
const EntityHandle MoFEM::OpCalculateHVecTensorTrace< Tensor_Dim, OpBase >::zeroType
private

Definition at line 3273 of file UserDataOperators.hpp.


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