v0.15.0
Loading...
Searching...
No Matches
Public Types | Public Member Functions | Public Attributes | Private Attributes | List of all members
OpPostProcDataStructure Struct Reference

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

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

Public Types

using OP = VolumeElementForcesAndSourcesCoreOnSide::UserDataOperator
 

Public Member Functions

 OpPostProcDataStructure (moab::Interface &post_proc_mesh, std::vector< EntityHandle > &map_gauss_pts, boost::shared_ptr< DataAtIntegrationPts > data_ptr, int sense)
 
MoFEMErrorCode doWork (int side, EntityType type, EntData &data)
 

Public Attributes

moab::Interface & postProcMesh
 
std::vector< EntityHandle > & mapGaussPts
 
boost::shared_ptr< DataAtIntegrationPtsdataAtPts
 

Private Attributes

int tagSense
 

Detailed Description

Examples
mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianPlasticity.cpp.

Definition at line 873 of file EshelbianOperators.hpp.

Member Typedef Documentation

◆ OP

using OpPostProcDataStructure::OP = VolumeElementForcesAndSourcesCoreOnSide::UserDataOperator

Definition at line 876 of file EshelbianOperators.hpp.

Constructor & Destructor Documentation

◆ OpPostProcDataStructure()

OpPostProcDataStructure::OpPostProcDataStructure ( moab::Interface &  post_proc_mesh,
std::vector< EntityHandle > &  map_gauss_pts,
boost::shared_ptr< DataAtIntegrationPts data_ptr,
int  sense 
)
inline

Definition at line 882 of file EshelbianOperators.hpp.

886 : OP(NOSPACE, UserDataOperator::OPSPACE), postProcMesh(post_proc_mesh),
887 mapGaussPts(map_gauss_pts), dataAtPts(data_ptr), tagSense(sense) {}
@ NOSPACE
Definition definitions.h:83
std::vector< EntityHandle > & mapGaussPts
boost::shared_ptr< DataAtIntegrationPts > dataAtPts
VolumeElementForcesAndSourcesCoreOnSide::UserDataOperator OP

Member Function Documentation

◆ doWork()

MoFEMErrorCode OpPostProcDataStructure::doWork ( int  side,
EntityType  type,
EntData data 
)
Examples
mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp.

Definition at line 3087 of file EshelbianOperators.cpp.

3088 {
3090
3091 if (tagSense != getSkeletonSense())
3093
3094 auto get_tag = [&](auto name) {
3095 auto &mob = getPtrFE()->mField.get_moab();
3096 Tag tag;
3097 CHK_MOAB_THROW(mob.tag_get_handle(name, tag), "get tag");
3098 return tag;
3099 };
3100
3101 auto get_tag_value = [&](auto &&tag, int dim) {
3102 auto &mob = getPtrFE()->mField.get_moab();
3103 auto face = getSidePtrFE()->getFEEntityHandle();
3104 std::vector<double> value(dim);
3105 CHK_MOAB_THROW(mob.tag_get_data(tag, &face, 1, value.data()), "set tag");
3106 return value;
3107 };
3108
3109 auto create_tag = [this](const std::string tag_name, const int size) {
3110 double def_VAL[] = {0, 0, 0, 0, 0, 0, 0, 0, 0};
3111 Tag th;
3112 CHKERR postProcMesh.tag_get_handle(tag_name.c_str(), size, MB_TYPE_DOUBLE,
3113 th, MB_TAG_CREAT | MB_TAG_SPARSE,
3114 def_VAL);
3115 return th;
3116 };
3117
3118 Tag th_cauchy_streess = create_tag("CauchyStress", 9);
3119 Tag th_detF = create_tag("detF", 1);
3120 Tag th_traction = create_tag("traction", 3);
3121 Tag th_disp_error = create_tag("DisplacementError", 1);
3122
3123 Tag th_energy = create_tag("Energy", 1);
3124
3125 auto t_w = getFTensor1FromMat<3>(dataAtPts->wL2AtPts);
3126 auto t_h = getFTensor2FromMat<3, 3>(dataAtPts->hAtPts);
3127 auto t_approx_P = getFTensor2FromMat<3, 3>(dataAtPts->approxPAtPts);
3128
3129 auto t_normal = getFTensor1NormalsAtGaussPts();
3130 auto t_disp = getFTensor1FromMat<3>(dataAtPts->wH1AtPts);
3131
3132 auto next = [&]() {
3133 ++t_w;
3134 ++t_h;
3135 ++t_approx_P;
3136 ++t_normal;
3137 ++t_disp;
3138 };
3139
3140 if (dataAtPts->energyAtPts.size() == 0) {
3141 // that is for case that energy is not calculated
3142 dataAtPts->energyAtPts.resize(getGaussPts().size2());
3143 dataAtPts->energyAtPts.clear();
3144 }
3145 auto t_energy = getFTensor0FromVec(dataAtPts->energyAtPts);
3146
3147 FTensor::Index<'i', 3> i;
3148 FTensor::Index<'j', 3> j;
3149 FTensor::Index<'k', 3> k;
3150 FTensor::Index<'l', 3> l;
3151
3152 auto set_float_precision = [](const double x) {
3153 if (std::abs(x) < std::numeric_limits<float>::epsilon())
3154 return 0.;
3155 else
3156 return x;
3157 };
3158
3159 // scalars
3160 auto save_scal_tag = [&](auto &th, auto v, const int gg) {
3162 v = set_float_precision(v);
3163 CHKERR postProcMesh.tag_set_data(th, &mapGaussPts[gg], 1, &v);
3165 };
3166
3167 // vectors
3168 VectorDouble3 v(3);
3169 FTensor::Tensor1<FTensor::PackPtr<double *, 0>, 3> t_v(&v[0], &v[1], &v[2]);
3170 auto save_vec_tag = [&](auto &th, auto &t_d, const int gg) {
3172 t_v(i) = t_d(i);
3173 for (auto &a : v.data())
3174 a = set_float_precision(a);
3175 CHKERR postProcMesh.tag_set_data(th, &mapGaussPts[gg], 1,
3176 &*v.data().begin());
3178 };
3179
3180 // tensors
3181
3182 MatrixDouble3by3 m(3, 3);
3184 &m(0, 0), &m(0, 1), &m(0, 2),
3185
3186 &m(1, 0), &m(1, 1), &m(1, 2),
3187
3188 &m(2, 0), &m(2, 1), &m(2, 2));
3189
3190 auto save_mat_tag = [&](auto &th, auto &t_d, const int gg) {
3192 t_m(i, j) = t_d(i, j);
3193 for (auto &v : m.data())
3194 v = set_float_precision(v);
3195 CHKERR postProcMesh.tag_set_data(th, &mapGaussPts[gg], 1,
3196 &*m.data().begin());
3198 };
3199
3200 const auto nb_gauss_pts = getGaussPts().size2();
3201 for (auto gg = 0; gg != nb_gauss_pts; ++gg) {
3202
3203 FTensor::Tensor1<double, 3> t_traction;
3204 t_traction(i) = t_approx_P(i, j) * t_normal(j) / t_normal.l2();
3205 // vectors
3206 CHKERR save_vec_tag(th_traction, t_traction, gg);
3207
3208 double u_error = sqrt((t_disp(i) - t_w(i)) * (t_disp(i) - t_w(i)));
3209 CHKERR save_scal_tag(th_disp_error, u_error, gg);
3210 CHKERR save_scal_tag(th_energy, t_energy, gg);
3211
3212 const double jac = determinantTensor3by3(t_h);
3214 t_cauchy(i, j) = (1. / jac) * (t_approx_P(i, k) * t_h(j, k));
3215 CHKERR save_mat_tag(th_cauchy_streess, t_cauchy, gg);
3216 CHKERR postProcMesh.tag_set_data(th_detF, &mapGaussPts[gg], 1, &jac);
3217
3218 next();
3219 }
3220
3222}
constexpr double a
#define MoFEMFunctionReturnHot(a)
Last executable line of each PETSc function used for error handling. Replaces return()
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
#define CHK_MOAB_THROW(err, msg)
Check error code of MoAB function and throw MoFEM exception.
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
#define CHKERR
Inline error check.
FTensor::Index< 'i', SPACE_DIM > i
const double v
phase velocity of light in medium (cm/ns)
FTensor::Index< 'l', 3 > l
FTensor::Index< 'j', 3 > j
FTensor::Index< 'k', 3 > k
VectorBoundedArray< double, 3 > VectorDouble3
Definition Types.hpp:92
static auto getFTensor0FromVec(ublas::vector< T, A > &data)
Get tensor rank 0 (scalar) form data vector.
static auto determinantTensor3by3(T &t)
Calculate the determinant of a 3x3 matrix or a tensor of rank 2.
FTensor::Index< 'm', 3 > m

Member Data Documentation

◆ dataAtPts

boost::shared_ptr<DataAtIntegrationPts> OpPostProcDataStructure::dataAtPts

Definition at line 880 of file EshelbianOperators.hpp.

◆ mapGaussPts

std::vector<EntityHandle>& OpPostProcDataStructure::mapGaussPts

Definition at line 879 of file EshelbianOperators.hpp.

◆ postProcMesh

moab::Interface& OpPostProcDataStructure::postProcMesh

Definition at line 878 of file EshelbianOperators.hpp.

◆ tagSense

int OpPostProcDataStructure::tagSense
private

Definition at line 892 of file EshelbianOperators.hpp.


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