v0.16.0
Loading...
Searching...
No Matches
Public Types | Public Member Functions | Protected Attributes | List of all members
EshelbianTestingMonitor Struct Reference
Inheritance diagram for EshelbianTestingMonitor:
[legend]
Collaboration diagram for EshelbianTestingMonitor:
[legend]

Public Types

using SetPtsData = FieldEvaluatorInterface::SetPtsData
 

Public Member Functions

 EshelbianTestingMonitor (EshelbianCore &ep, boost::shared_ptr< EshelbianMonitor > base_monitor_ptr)
 
MoFEMErrorCode preProcess ()
 
MoFEMErrorCode operator() ()
 
MoFEMErrorCode postProcess ()
 
MoFEMErrorCode checkExternalStrain (std::array< double, 3 > point, std::string str, PetscInt atom_test)
 

Protected Attributes

boost::shared_ptr< EshelbianMonitor > baseMonitorPtr
 
EshelbianCore & eP
 
std::map< std::string, std::array< double, 6 > > & reactionForcesMap
 
boost::shared_ptr< SetPtsData > dataFieldEval
 
std::map< std::string, std::array< double, 3 > > & ptsHashMap
 
boost::shared_ptr< double > gEnergy
 

Detailed Description

Definition at line 6 of file EshelbianTestingMonitor.cpp.

Member Typedef Documentation

◆ SetPtsData

using EshelbianTestingMonitor::SetPtsData = FieldEvaluatorInterface::SetPtsData

Definition at line 7 of file EshelbianTestingMonitor.cpp.

Constructor & Destructor Documentation

◆ EshelbianTestingMonitor()

EshelbianTestingMonitor::EshelbianTestingMonitor ( EshelbianCore &  ep,
boost::shared_ptr< EshelbianMonitor >  base_monitor_ptr 
)
inline

Definition at line 9 of file EshelbianTestingMonitor.cpp.

11 : baseMonitorPtr(base_monitor_ptr), eP(baseMonitorPtr->getEpCore()),
12 ptsHashMap(baseMonitorPtr->getHashMap()),
13 reactionForcesMap(baseMonitorPtr->getReactionMap()),
14 gEnergy(baseMonitorPtr->getEnergy()),
15 dataFieldEval(baseMonitorPtr->getDataField()) {
16
17 PetscBool test_cook_flg = PETSC_FALSE;
18 PetscInt atom_test = 0;
19 CHK_THROW_MESSAGE(PetscOptionsGetBool(PETSC_NULLPTR, "", "-test_cook_pts",
20 &test_cook_flg, PETSC_NULLPTR),
21 "get post proc points");
22 CHKERR PetscOptionsGetInt(PETSC_NULLPTR, "", "-atom_test", &atom_test,
23 PETSC_NULLPTR);
24
25 if (test_cook_flg) {
26 ptsHashMap["Point A"] = {48., 60., 4.999};
27 ptsHashMap["Point B"] = {48. / 2., 44. + (60. - 44.) / 2., 0.};
28 ptsHashMap["Point C"] = {48. / 2., (44. - 0.) / 2., 0.};
29 }
30 if (atom_test == 14) {
31 // Points for atom test 14: check external strain
32 ptsHashMap["Point (2.5, 0., 0.)"] = {2.5, 0, 0.};
33 }
34 }
CHK_THROW_MESSAGE
#define CHK_THROW_MESSAGE(err, msg)
Check and throw MoFEM exception.
Definition definitions.h:612
CHKERR
#define CHKERR
Inline error check.
Definition definitions.h:551
MoFEM::PetscOptionsGetInt
PetscErrorCode PetscOptionsGetInt(PetscOptions *, const char pre[], const char name[], PetscInt *ivalue, PetscBool *set)
Definition DeprecatedPetsc.hpp:142
MoFEM::PetscOptionsGetBool
PetscErrorCode PetscOptionsGetBool(PetscOptions *, const char pre[], const char name[], PetscBool *bval, PetscBool *set)
Definition DeprecatedPetsc.hpp:182
EshelbianTestingMonitor::ptsHashMap
std::map< std::string, std::array< double, 3 > > & ptsHashMap
Definition EshelbianTestingMonitor.cpp:152
EshelbianTestingMonitor::gEnergy
boost::shared_ptr< double > gEnergy
Definition EshelbianTestingMonitor.cpp:153
EshelbianTestingMonitor::baseMonitorPtr
boost::shared_ptr< EshelbianMonitor > baseMonitorPtr
Definition EshelbianTestingMonitor.cpp:148
EshelbianTestingMonitor::reactionForcesMap
std::map< std::string, std::array< double, 6 > > & reactionForcesMap
Definition EshelbianTestingMonitor.cpp:150
EshelbianTestingMonitor::dataFieldEval
boost::shared_ptr< SetPtsData > dataFieldEval
Definition EshelbianTestingMonitor.cpp:151
atom_test
int atom_test
Atom test.
Definition plastic.cpp:122

Member Function Documentation

◆ checkExternalStrain()

MoFEMErrorCode EshelbianTestingMonitor::checkExternalStrain ( std::array< double, 3 >  point,
std::string  str,
PetscInt  atom_test 
)

Definition at line 156 of file EshelbianTestingMonitor.cpp.

157 {
158 MoFEMFunctionBegin;
159
160 dataFieldEval->setEvalPoints(point.data(), point.size() / 3);
161 auto vec = createVectorMPI(eP.mField.get_comm(), 1, 1);
162 auto w_l2_at_pts = eP.dataAtPts->getSmallWL2AtPts();
163 w_l2_at_pts->resize(0, 0, false);
164
165 if (auto fe_ptr = dataFieldEval->feMethodPtr) {
166 CHKERR eP.mField.getInterface<FieldEvaluatorInterface>()
167 ->evalFEAtThePoint<SPACE_DIM>(
168 point.data(), 1e-12, problemPtr->getName(), "EP", dataFieldEval,
169 eP.mField.get_comm_rank(), eP.mField.get_comm_rank(), nullptr,
170 MF_EXIST, QUIET);
171 }
172 if (w_l2_at_pts->size1() == 0 || w_l2_at_pts->size2() == 0) {
173 CHKERR VecSetValue(vec, 0, 0.0, ADD_VALUES);
174 } else if (w_l2_at_pts->size1() == 1 || w_l2_at_pts->size2() == 1) {
175 auto add = [&]() {
176 std::ostringstream s;
177 s << str << " elem " << getFEEntityHandle() << " ";
178 return s.str();
179 };
180 MOFEM_LOG("EPSYNC", Sev::inform)
181 << add() << "comm rank " << eP.mField.get_comm_rank();
182 MOFEM_LOG("EPSYNC", Sev::inform)
183 << add() << "point " << getVectorAdaptor(point.data(), 3);
184 MOFEM_LOG("EPSYNC", Sev::inform)
185 << add() << "w " << *w_l2_at_pts;
186 double disp_at_point = (*w_l2_at_pts)(0, 0);
187 CHKERR VecSetValue(vec, 0, disp_at_point, ADD_VALUES);
188 } else {
189 SETERRQ(PETSC_COMM_WORLD, MOFEM_DATA_INCONSISTENCY,
190 "Unexpected displacement data shape (%ld, %ld).",
191 static_cast<long>(w_l2_at_pts->size1()),
192 static_cast<long>(w_l2_at_pts->size2()));
193 }
194
195 MOFEM_LOG_SEVERITY_SYNC(eP.mField.get_comm(), Sev::inform);
196 if (ts_t > 0.0) {
197 CHKERR VecAssemblyBegin(vec);
198 CHKERR VecAssemblyEnd(vec);
199 double error;
200 PetscInt idx = 0;
201 if (eP.mField.get_comm_rank() == 0) {
202 CHKERR VecGetValues(vec, 1, &idx, &error);
203 }
204 MPI_Bcast(&error, 1, MPI_DOUBLE, 0, PETSC_COMM_WORLD);
205 // Check if the displacement is correct for the applied external
206 // strain For the bar problem, we expect a displacement of 0.25 at
207 // point A
208 if (std::abs(error - 0.25) > 1e-5) {
209 SETERRQ(PETSC_COMM_WORLD, MOFEM_ATOM_TEST_INVALID,
210 "Atom test %d failed: wrong displacement %.12g.", atom_test,
211 error);
212 }
213 }
214 MoFEMFunctionReturn(0);
215}
MOFEM_LOG_SEVERITY_SYNC
#define MOFEM_LOG_SEVERITY_SYNC(comm, severity)
Synchronise "SYNC" on curtain severity level.
Definition LogManager.hpp:360
QUIET
@ QUIET
Definition definitions.h:221
MF_EXIST
@ MF_EXIST
Definition definitions.h:113
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition definitions.h:359
MOFEM_ATOM_TEST_INVALID
@ MOFEM_ATOM_TEST_INVALID
Definition definitions.h:40
MOFEM_DATA_INCONSISTENCY
@ MOFEM_DATA_INCONSISTENCY
Definition definitions.h:31
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition definitions.h:432
MOFEM_LOG
#define MOFEM_LOG(channel, severity)
Log.
Definition LogManager.hpp:316
MoFEM::getVectorAdaptor
auto getVectorAdaptor(T1 ptr, const size_t n)
Get Vector adaptor.
Definition Templates.hpp:49
MoFEM::createVectorMPI
auto createVectorMPI(MPI_Comm comm, PetscInt n, PetscInt N)
Create MPI Vector.
Definition PetscSmartObj.hpp:204
EshelbianCore::mField
MoFEM::Interface & mField
Definition EshelbianCore.hpp:228
EshelbianCore::dataAtPts
boost::shared_ptr< DataAtIntegrationPts > dataAtPts
Definition EshelbianCore.hpp:230
MoFEM::CoreInterface::get_comm
virtual MPI_Comm & get_comm() const =0
MoFEM::CoreInterface::get_comm_rank
virtual int get_comm_rank() const =0
MoFEM::UnknownInterface::getInterface
MoFEMErrorCode getInterface(IFACE *&iface) const
Get interface reference to pointer of interface.
Definition UnknownInterface.hpp:93

◆ operator()()

MoFEMErrorCode EshelbianTestingMonitor::operator() ( )
inline

Definition at line 38 of file EshelbianTestingMonitor.cpp.

38{ return 0; }

◆ postProcess()

MoFEMErrorCode EshelbianTestingMonitor::postProcess ( )
inline

Definition at line 40 of file EshelbianTestingMonitor.cpp.

40 {
41 MoFEMFunctionBegin;
42
43 MOFEM_LOG("EP", Sev::inform) << "Testing Monitor postProcess";
44
45 PetscBool test_cook_flg = PETSC_FALSE;
46 PetscInt atom_test = 0;
47 char reaction_block_name[255] = "FIX_ALL";
48 CHK_THROW_MESSAGE(PetscOptionsGetBool(PETSC_NULLPTR, "", "-test_cook_pts",
49 &test_cook_flg, PETSC_NULLPTR),
50 "get post proc points");
51 CHKERR PetscOptionsGetInt(PETSC_NULLPTR, "", "-atom_test", &atom_test,
52 PETSC_NULLPTR);
53 CHKERR PetscOptionsGetString(PETSC_NULLPTR, "", "-atom_test_reaction_block",
54 reaction_block_name, 255, PETSC_NULLPTR);
55
56 if (test_cook_flg) {
57 // Points for Cook beam
58 for (auto &pts : ptsHashMap) {
59 CHKERR baseMonitorPtr->postProcAtPoints(pts.second, pts.first);
60 MOFEM_LOG_SEVERITY_SYNC(eP.mField.get_comm(), Sev::inform);
61 }
62 }
63
64 switch (atom_test) {
65 case 14:
66 // Points for external strain
67 for (auto &pts : ptsHashMap) {
68 CHKERR checkExternalStrain(pts.second, pts.first, atom_test);
69 }
70 break;
71 case 15:
72 if (ts_t > 0.0) {
73 // L = 5, ty = tz = 1, A = 1, UDL Top face 0.1,
74 // Moment at x = 0 should be 0.25.
75 if (std::abs(reactionForcesMap[reaction_block_name][5] - 0.25) > 1e-3) {
76 MOFEM_LOG_SYNCHRONISE(eP.mField.get_comm());
77 SETERRQ(
78 PETSC_COMM_WORLD, MOFEM_ATOM_TEST_INVALID,
79 "Atom test 15 failed: reaction moment does not match expected "
80 "value. Got [%3.6e], expected [0.25].",
81 reactionForcesMap[reaction_block_name][5]);
82 }
83 }
84 break;
85 case 16:
86 if (ts_t > 0.0) {
87 // L = 5, ty = tz = 1, A = 1, E = 1000 * exp(0.1*x), u = 0.1, F = 1000 *
88 // 0.1/(1-exp(-0.5)) = 25.4149, so reaction force should
89 // be -25.4149
90 if (std::abs(reactionForcesMap[reaction_block_name][0] - (-25.4149)) >
91 1e-4) {
92 MOFEM_LOG_SYNCHRONISE(eP.mField.get_comm());
93 SETERRQ(PETSC_COMM_WORLD, MOFEM_ATOM_TEST_INVALID,
94 "Atom test 16 failed: reaction force does not match expected "
95 "value. Got [%3.6e, %3.6e, %3.6e], expected [-25.4149, 0.0, "
96 "0.0].",
97 reactionForcesMap[reaction_block_name][0],
98 reactionForcesMap[reaction_block_name][1],
99 reactionForcesMap[reaction_block_name][2]);
100 }
101 }
102 break;
103 case 17:
104 if (ts_t > 0.0) {
105 if (std::abs(*gEnergy - 1.27096) > 1e-5) {
106 MOFEM_LOG_SYNCHRONISE(eP.mField.get_comm());
107 SETERRQ(PETSC_COMM_WORLD, MOFEM_ATOM_TEST_INVALID,
108 "Atom test 17 failed: strain energy does not match expected "
109 "value. Got %3.6e, expected 1.27096.",
110 *gEnergy);
111 }
112 }
113 break;
114 case 18:
115 if (ts_step == 4) {
116 if (std::abs(eP.loadFactor - 12.5325) > 1e-5) {
117 MOFEM_LOG_SYNCHRONISE(eP.mField.get_comm());
118 SETERRQ(PETSC_COMM_WORLD, MOFEM_ATOM_TEST_INVALID,
119 "Atom test 18 failed: load factor does not match expected "
120 "value. Got %3.6e, expected 12.5325.",
121 eP.loadFactor);
122 }
123 }
124 break;
125 case 19:
126 if (ts_t > 0.0) {
127 if (std::abs(reactionForcesMap["SPRING_BC"][2] - 2.016e1) > 1e-3) {
128 MOFEM_LOG_SYNCHRONISE(eP.mField.get_comm());
129 SETERRQ(
130 PETSC_COMM_WORLD, MOFEM_ATOM_TEST_INVALID,
131 "Atom test 19 failed: reaction force does not match expected "
132 "value. Got [%3.6e], expected [2.016e1].",
133 reactionForcesMap["SPRING_BC"][2]);
134 }
135 }
136 break;
137 default:
138 break;
139 }
140
141 MoFEMFunctionReturn(0);
142 }
MOFEM_LOG_SYNCHRONISE
#define MOFEM_LOG_SYNCHRONISE(comm)
Synchronise "SYNC" channel.
Definition LogManager.hpp:353
MoFEM::PetscOptionsGetString
PetscErrorCode PetscOptionsGetString(PetscOptions *, const char pre[], const char name[], char str[], size_t size, PetscBool *set)
Definition DeprecatedPetsc.hpp:172
EshelbianTestingMonitor::checkExternalStrain
MoFEMErrorCode checkExternalStrain(std::array< double, 3 > point, std::string str, PetscInt atom_test)
Definition EshelbianTestingMonitor.cpp:156

◆ preProcess()

MoFEMErrorCode EshelbianTestingMonitor::preProcess ( )
inline

Definition at line 36 of file EshelbianTestingMonitor.cpp.

36{ return 0; }

Member Data Documentation

◆ baseMonitorPtr

boost::shared_ptr<EshelbianMonitor> EshelbianTestingMonitor::baseMonitorPtr
protected

Definition at line 148 of file EshelbianTestingMonitor.cpp.

◆ dataFieldEval

boost::shared_ptr<SetPtsData> EshelbianTestingMonitor::dataFieldEval
protected

Definition at line 151 of file EshelbianTestingMonitor.cpp.

◆ eP

EshelbianCore& EshelbianTestingMonitor::eP
protected

Definition at line 149 of file EshelbianTestingMonitor.cpp.

◆ gEnergy

boost::shared_ptr<double> EshelbianTestingMonitor::gEnergy
protected

Definition at line 153 of file EshelbianTestingMonitor.cpp.

◆ ptsHashMap

std::map<std::string, std::array<double, 3> >& EshelbianTestingMonitor::ptsHashMap
protected

Definition at line 152 of file EshelbianTestingMonitor.cpp.

◆ reactionForcesMap

std::map<std::string, std::array<double, 6> >& EshelbianTestingMonitor::reactionForcesMap
protected

Definition at line 150 of file EshelbianTestingMonitor.cpp.

The documentation for this struct was generated from the following file:
Generated by Doxygen 1.9.8 and hosted at University of Glasgow