v0.15.5
Loading...
Searching...
No Matches
Public Member Functions | Private Attributes | Friends | List of all members
EshelbianPlasticity::TopologicalTAOCtxImpl Struct Reference
Inheritance diagram for EshelbianPlasticity::TopologicalTAOCtxImpl:
[legend]
Collaboration diagram for EshelbianPlasticity::TopologicalTAOCtxImpl:
[legend]

Public Member Functions

 TopologicalTAOCtxImpl (EshelbianCore *ep, ForcesAndSourcesCore::GaussHookFun set_integration_at_interior, ForcesAndSourcesCore::GaussHookFun set_integration_at_face, SmartPetscObj< TS > time_solver)
 
- Public Member Functions inherited from EshelbianPlasticity::TopologicalTAOCtx
 TopologicalTAOCtx ()=default
 
virtual ~TopologicalTAOCtx ()=default
 

Private Attributes

EshelbianCore * ep_ptr
 
ForcesAndSourcesCore::GaussHookFun integrationAtInterior
 
ForcesAndSourcesCore::GaussHookFun integrationAtFace
 
SmartPetscObj< TS > timeSolver
 
SmartPetscObj< Vec > primalProblemVec
 
SmartPetscObj< Vec > adjointProblemVec
 

Friends

MoFEMErrorCode evaluatePrimalProblemTopologicalImpl (TopologicalTAOCtxImpl *ctx_impl_ptr)
 
MoFEMErrorCode evaluateAdjointProblemImpl (TopologicalTAOCtxImpl *ctx_impl_ptr)
 
MoFEMErrorCode evaluateGradientImpl (TopologicalTAOCtxImpl *ctx_impl_ptr, double *f, Vec g)
 
MoFEMErrorCode finiteDifferenceGradientTest (TopologicalTAOCtxImpl *ctx_impl_ptr, Vec sol, double *f, Vec g, double epsilon)
 
MoFEMErrorCode topologicalEvaluateObjectiveAndGradient (Tao tao, Vec sol, PetscReal *f, Vec g, void *ctx)
 
MoFEMErrorCode testTopologicalDerivative (TopologicalTAOCtx *ctx_ptr, Vec sol, PetscReal *f, Vec g)
 

Detailed Description

Definition at line 31 of file EshelbianTopologicalDerivative.cpp.

Constructor & Destructor Documentation

◆ TopologicalTAOCtxImpl()

EshelbianPlasticity::TopologicalTAOCtxImpl::TopologicalTAOCtxImpl ( EshelbianCore *  ep,
ForcesAndSourcesCore::GaussHookFun  set_integration_at_interior,
ForcesAndSourcesCore::GaussHookFun  set_integration_at_face,
SmartPetscObj< TS >  time_solver 
)
inline

Definition at line 32 of file EshelbianTopologicalDerivative.cpp.

37 : ep_ptr(ep), integrationAtInterior(set_integration_at_interior),
38 integrationAtFace(set_integration_at_face), timeSolver(time_solver) {
39
40 primalProblemVec = createDMVector(ep->dmElastic, RowColData::COL);
41 adjointProblemVec = createDMVector(ep->dmElastic, RowColData::COL);
42 }
COL
@ COL
Definition definitions.h:136
MoFEM::createDMVector
auto createDMVector(DM dm, RowColData rc=RowColData::COL)
Get smart vector from DM.
Definition DMMoFEM.hpp:1237
EshelbianCore::dmElastic
SmartPetscObj< DM > dmElastic
Elastic problem.
Definition EshelbianCore.hpp:217

Friends And Related Symbol Documentation

◆ evaluateAdjointProblemImpl

MoFEMErrorCode evaluateAdjointProblemImpl ( TopologicalTAOCtxImpl *  ctx_impl_ptr)
friend

Definition at line 119 of file EshelbianTopologicalDerivative.cpp.

119 {
120 MoFEMFunctionBegin;
121 auto &ep = *ctx_impl_ptr->ep_ptr;
122 auto &m_field = ep.mField;
123 auto interior_integration_hook = ctx_impl_ptr->integrationAtInterior;
124 auto boundary_integration_hook = ctx_impl_ptr->integrationAtFace;
125
126 auto fe_spatial =
127 boost::make_shared<VolumeElementForcesAndSourcesCore>(m_field);
128 auto dJ_dx_vec = pushTopologicalSpatialOps(
129 ep, fe_spatial, interior_integration_hook, boundary_integration_hook);
130 CHKERR DMoFEMLoopFiniteElements(ep.dmElastic, ep.elementVolumeName,
131 fe_spatial);
132 CHKERR VecAssemblyBegin(dJ_dx_vec);
133 CHKERR VecAssemblyEnd(dJ_dx_vec);
134 CHKERR VecGhostUpdateBegin(dJ_dx_vec, ADD_VALUES, SCATTER_REVERSE);
135 CHKERR VecGhostUpdateEnd(dJ_dx_vec, ADD_VALUES, SCATTER_REVERSE);
136 CHKERR VecGhostUpdateBegin(dJ_dx_vec, INSERT_VALUES, SCATTER_FORWARD);
137 CHKERR VecGhostUpdateEnd(dJ_dx_vec, INSERT_VALUES, SCATTER_FORWARD);
138
139 auto ksp = snesGetKSP(tsGetSNES(ctx_impl_ptr->timeSolver));
140 CHKERR KSPSolveTranspose(ksp, dJ_dx_vec, ctx_impl_ptr->adjointProblemVec);
141 CHKERR VecGhostUpdateBegin(ctx_impl_ptr->adjointProblemVec, INSERT_VALUES,
142 SCATTER_FORWARD);
143 CHKERR VecGhostUpdateEnd(ctx_impl_ptr->adjointProblemVec, INSERT_VALUES,
144 SCATTER_FORWARD);
145
146 MoFEMFunctionReturn(0);
147}
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition definitions.h:359
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition definitions.h:432
CHKERR
#define CHKERR
Inline error check.
Definition definitions.h:551
MoFEM::DMoFEMLoopFiniteElements
PetscErrorCode DMoFEMLoopFiniteElements(DM dm, const char fe_name[], MoFEM::FEMethod *method, CacheTupleWeakPtr cache_ptr=CacheTupleSharedPtr())
Executes FEMethod for finite elements in DM.
Definition DMMoFEM.cpp:576
EshelbianPlasticity::pushTopologicalSpatialOps
SmartPetscObj< Vec > pushTopologicalSpatialOps(EshelbianCore &ep, boost::shared_ptr< VolumeElementForcesAndSourcesCore > fe, ForcesAndSourcesCore::GaussHookFun interior_integration_hook, ForcesAndSourcesCore::GaussHookFun boundary_integration_hook)
Definition EshelbianTopologicalDerivativePush.cpp:195
MoFEM::snesGetKSP
auto snesGetKSP(SNES snes)
Definition PetscSmartObj.hpp:416
MoFEM::tsGetSNES
auto tsGetSNES(TS ts)
Definition PetscSmartObj.hpp:410
EshelbianCore::mField
MoFEM::Interface & mField
Definition EshelbianCore.hpp:203

◆ evaluateGradientImpl

MoFEMErrorCode evaluateGradientImpl ( TopologicalTAOCtxImpl *  ctx_impl_ptr,
double *  f,
Vec  g 
)
friend

Definition at line 149 of file EshelbianTopologicalDerivative.cpp.

150 {
151 MoFEMFunctionBegin;
152 auto &ep = *ctx_impl_ptr->ep_ptr;
153 auto &m_field = ep.mField;
154 auto interior_integration_hook = ctx_impl_ptr->integrationAtInterior;
155 auto boundary_integration_hook = ctx_impl_ptr->integrationAtFace;
156
157 auto fe_material =
158 boost::make_shared<VolumeElementForcesAndSourcesCore>(m_field);
159 auto fe_adjoint =
160 boost::make_shared<VolumeElementForcesAndSourcesCore>(m_field);
161
162 boost::shared_ptr<double> J_ptr(fe_material, f);
163 auto dJ_dX_vec = SmartPetscObj<Vec>(g, true);
164 CHKERR VecZeroEntries(dJ_dX_vec);
165 CHKERR VecGhostUpdateBegin(dJ_dX_vec, INSERT_VALUES, SCATTER_FORWARD);
166 CHKERR VecGhostUpdateEnd(dJ_dX_vec, INSERT_VALUES, SCATTER_FORWARD);
167
168 CHKERR pushTopologicalMaterialOps(ep, fe_material, interior_integration_hook,
169 boundary_integration_hook, J_ptr,
170 dJ_dX_vec);
171 auto alpha = ep.alphaW;
172 auto rho = ep.alphaRho;
173 auto alpha_omega = ep.alphaOmega;
174 CHKERR pushTopologicalOps_dJ_adjoint_gradient(
175 ep, fe_adjoint, interior_integration_hook, boundary_integration_hook,
176 ctx_impl_ptr->adjointProblemVec, dJ_dX_vec, alpha, rho, alpha_omega);
177
178 CHKERR DMoFEMLoopFiniteElements(ep.dmMaterial, ep.elementVolumeName,
179 fe_material);
180 CHKERR DMoFEMLoopFiniteElements(ep.dmMaterial, ep.elementVolumeName,
181 fe_adjoint);
182 CHKERR VecAssemblyBegin(dJ_dX_vec);
183 CHKERR VecAssemblyEnd(dJ_dX_vec);
184 CHKERR VecGhostUpdateBegin(dJ_dX_vec, ADD_VALUES, SCATTER_REVERSE);
185 CHKERR VecGhostUpdateEnd(dJ_dX_vec, ADD_VALUES, SCATTER_REVERSE);
186 CHKERR VecGhostUpdateBegin(dJ_dX_vec, INSERT_VALUES, SCATTER_FORWARD);
187 CHKERR VecGhostUpdateEnd(dJ_dX_vec, INSERT_VALUES, SCATTER_FORWARD);
188
189 MoFEMFunctionReturn(0);
190}
EshelbianPlasticity::pushTopologicalOps_dJ_adjoint_gradient
MoFEMErrorCode pushTopologicalOps_dJ_adjoint_gradient(EshelbianCore &ep, boost::shared_ptr< VolumeElementForcesAndSourcesCore > fe, ForcesAndSourcesCore::GaussHookFun interior_integration_hook, ForcesAndSourcesCore::GaussHookFun boundary_integration_hook, SmartPetscObj< Vec > lambda_vec, SmartPetscObj< Vec > dJ_dX_vec, const double alpha, const double rho, const double alpha_omega)
Definition EshelbianTopologicalDerivativePush.cpp:263
EshelbianPlasticity::pushTopologicalMaterialOps
MoFEMErrorCode pushTopologicalMaterialOps(EshelbianCore &ep, boost::shared_ptr< VolumeElementForcesAndSourcesCore > fe, ForcesAndSourcesCore::GaussHookFun interior_integration_hook, ForcesAndSourcesCore::GaussHookFun boundary_integration_hook, boost::shared_ptr< double > J_ptr, SmartPetscObj< Vec > dJ_dX_vec)
Definition EshelbianTopologicalDerivativePush.cpp:229
g
constexpr double g
Definition shallow_wave.cpp:63
MoFEM::SmartPetscObj
intrusive_ptr for managing petsc objects
Definition PetscSmartObj.hpp:85
rho
double rho
Definition plastic.cpp:144

◆ evaluatePrimalProblemTopologicalImpl

MoFEMErrorCode evaluatePrimalProblemTopologicalImpl ( TopologicalTAOCtxImpl *  ctx_impl_ptr)
friend

Definition at line 91 of file EshelbianTopologicalDerivative.cpp.

91 {
92 MoFEMFunctionBegin;
93 auto &ep = *ctx_impl_ptr->ep_ptr;
94 auto ts = ctx_impl_ptr->timeSolver;
95
96 CHKERR VecGhostUpdateBegin(ctx_impl_ptr->primalProblemVec, INSERT_VALUES,
97 SCATTER_FORWARD);
98 CHKERR VecGhostUpdateEnd(ctx_impl_ptr->primalProblemVec, INSERT_VALUES,
99 SCATTER_FORWARD);
100 CHKERR DMoFEMMeshToLocalVector(ep.dmElastic, ctx_impl_ptr->primalProblemVec,
101 INSERT_VALUES, SCATTER_FORWARD,
102 RowColData::ROW);
103
104 CHKERR TSSetFromOptions(ts);
105 CHKERR TSSetStepNumber(ts, 0);
106 CHKERR TSSetTime(ts, 0);
107 CHKERR TSSetSolution(ts, ctx_impl_ptr->primalProblemVec);
108 double dt;
109 CHKERR TSGetTimeStep(ts, &dt);
110 MOFEM_LOG("EP", Sev::inform)
111 << "evaluatePrimalProblemTopologicalImpl: Time step dt: " << dt;
112 CHKERR TSSolve(ts, PETSC_NULLPTR);
113 CHKERR TSSetTimeStep(ts, dt);
114
115
116 MoFEMFunctionReturn(0);
117}
ROW
@ ROW
Definition definitions.h:136
MoFEM::DMoFEMMeshToLocalVector
PetscErrorCode DMoFEMMeshToLocalVector(DM dm, Vec l, InsertMode mode, ScatterMode scatter_mode, RowColData rc=RowColData::COL)
set local (or ghosted) vector values on mesh for partition only
Definition DMMoFEM.cpp:514
MOFEM_LOG
#define MOFEM_LOG(channel, severity)
Log.
Definition LogManager.hpp:316
dt
double dt
Definition heat_method.cpp:27

◆ finiteDifferenceGradientTest

MoFEMErrorCode finiteDifferenceGradientTest ( TopologicalTAOCtxImpl *  ctx_impl_ptr,
Vec  sol,
double *  f,
Vec  g,
double  epsilon 
)
friend

Definition at line 192 of file EshelbianTopologicalDerivative.cpp.

194 {
195 MoFEMFunctionBegin;
196
197 auto &ep = *ctx_impl_ptr->ep_ptr;
198 auto &m_field = ep.mField;
199 auto interior_integration_hook = ctx_impl_ptr->integrationAtInterior;
200 auto boundary_integration_hook = ctx_impl_ptr->integrationAtFace;
201
202 auto fe_material =
203 boost::make_shared<VolumeElementForcesAndSourcesCore>(m_field);
204 auto J_ptr = boost::make_shared<double>(0);
205 CHKERR pushTopologicalMaterialOps(ep, fe_material, interior_integration_hook,
206 boundary_integration_hook, J_ptr,
207 SmartPetscObj<Vec>());
208
209 auto opt = ep.mField.getInterface<OperatorsTester>();
210 constexpr double eps = 1e-6;
211 auto random_vec = opt->setRandomFields(
212 ep.dmMaterial, {{ep.materialH1Positions, {-eps, eps}}}, nullptr,
213 RowColData::ROW);
214 auto a_vec = vectorDuplicate(random_vec);
215 CHKERR VecCopy(sol, a_vec);
216 CHKERR VecAXPY(a_vec, 1, random_vec);
217 auto b_vec = vectorDuplicate(random_vec);
218 CHKERR VecCopy(sol, b_vec);
219 CHKERR VecAXPY(b_vec, -1, random_vec);
220
221 *J_ptr = 0;
222 CHKERR VecGhostUpdateBegin(a_vec, INSERT_VALUES, SCATTER_FORWARD);
223 CHKERR VecGhostUpdateEnd(a_vec, INSERT_VALUES, SCATTER_FORWARD);
224 CHKERR DMoFEMMeshToLocalVector(ep.dmMaterial, a_vec, INSERT_VALUES,
225 SCATTER_REVERSE, RowColData::ROW);
226 CHKERR DMoFEMLoopFiniteElements(ep.dmMaterial, ep.elementVolumeName,
227 fe_material);
228 double J_plus = *J_ptr;
229 *J_ptr = 0;
230 CHKERR VecGhostUpdateBegin(b_vec, INSERT_VALUES, SCATTER_FORWARD);
231 CHKERR VecGhostUpdateEnd(b_vec, INSERT_VALUES, SCATTER_FORWARD);
232 CHKERR DMoFEMMeshToLocalVector(ep.dmMaterial, b_vec, INSERT_VALUES,
233 SCATTER_REVERSE, RowColData::ROW);
234 CHKERR DMoFEMLoopFiniteElements(ep.dmMaterial, ep.elementVolumeName,
235 fe_material);
236 double J_minus = *J_ptr;
237 double dJ_da = (J_plus - J_minus) / (2 * eps);
238
239 double exact_dJ;
240 CHKERR VecDot(g, random_vec, &exact_dJ);
241 double error = std::abs(dJ_da - exact_dJ);
242
243 MOFEM_LOG("EP", Sev::inform)
244 << "J = " << *f << ", dJ/da = " << dJ_da << ", exact dJ/da = " << exact_dJ
245 << ", error = " << error;
246
247 CHKERR VecGhostUpdateBegin(sol, INSERT_VALUES, SCATTER_FORWARD);
248 CHKERR VecGhostUpdateEnd(sol, INSERT_VALUES, SCATTER_FORWARD);
249 CHKERR DMoFEMMeshToLocalVector(ep.dmMaterial, sol, INSERT_VALUES,
250 SCATTER_REVERSE, RowColData::ROW);
251
252 MoFEMFunctionReturn(0);
253}
eps
static const double eps
Definition check_base_functions_derivatives_on_tet.cpp:11
MoFEM::vectorDuplicate
SmartPetscObj< Vec > vectorDuplicate(Vec vec)
Create duplicate vector of smart vector.
Definition PetscSmartObj.hpp:223
MoFEM::OperatorsTester
Calculate directional derivative of the right hand side and compare it with tangent matrix derivative...
Definition OperatorsTester.hpp:21

◆ testTopologicalDerivative

MoFEMErrorCode testTopologicalDerivative ( TopologicalTAOCtx *  ctx_ptr,
Vec  sol,
PetscReal *  f,
Vec  g 
)
friend

Definition at line 255 of file EshelbianTopologicalDerivative.cpp.

256 {
257
258 MoFEMFunctionBegin;
259 auto *ctx_impl_ptr = dynamic_cast<TopologicalTAOCtxImpl *>(ctx_ptr);
260 if (!ctx_impl_ptr)
261 SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
262 "Invalid context pointer type");
263 auto &ep = *ctx_impl_ptr->ep_ptr;
264
265 CHKERR VecGhostUpdateBegin(sol, INSERT_VALUES, SCATTER_FORWARD);
266 CHKERR VecGhostUpdateEnd(sol, INSERT_VALUES, SCATTER_FORWARD);
267 CHKERR DMoFEMMeshToLocalVector(ep.dmMaterial, sol, INSERT_VALUES,
268 SCATTER_REVERSE, RowColData::ROW);
269
270 CHKERR VecZeroEntries(g);
271 *f = 0;
272
273 CHKERR evaluatePrimalProblemTopologicalImpl(ctx_impl_ptr);
274 double norm2_x;
275 CHKERR VecNorm(ctx_impl_ptr->primalProblemVec, NORM_2, &norm2_x);
276 MOFEM_LOG("EP", Sev::inform)
277 << "evaluatePrimalProblemTopologicalImpl: Norm of displacement vector: "
278 << norm2_x;
279 if (norm2_x < 1e-12) {
280 MoFEMFunctionReturnHot(0);
281 }
282
283 CHKERR evaluateAdjointProblemImpl(ctx_impl_ptr);
284 CHKERR evaluateGradientImpl(ctx_impl_ptr, f, g);
285
286 CHKERR finiteDifferenceGradientTest(ctx_impl_ptr, sol, f, g, 1e-6);
287
288 MoFEMFunctionReturn(0);
289}
MoFEMFunctionReturnHot
#define MoFEMFunctionReturnHot(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition definitions.h:463
MOFEM_DATA_INCONSISTENCY
@ MOFEM_DATA_INCONSISTENCY
Definition definitions.h:31
EshelbianPlasticity::TopologicalTAOCtxImpl::evaluatePrimalProblemTopologicalImpl
friend MoFEMErrorCode evaluatePrimalProblemTopologicalImpl(TopologicalTAOCtxImpl *ctx_impl_ptr)
Definition EshelbianTopologicalDerivative.cpp:91
EshelbianPlasticity::TopologicalTAOCtxImpl::evaluateGradientImpl
friend MoFEMErrorCode evaluateGradientImpl(TopologicalTAOCtxImpl *ctx_impl_ptr, double *f, Vec g)
Definition EshelbianTopologicalDerivative.cpp:149
EshelbianPlasticity::TopologicalTAOCtxImpl::finiteDifferenceGradientTest
friend MoFEMErrorCode finiteDifferenceGradientTest(TopologicalTAOCtxImpl *ctx_impl_ptr, Vec sol, double *f, Vec g, double epsilon)
Definition EshelbianTopologicalDerivative.cpp:192
EshelbianPlasticity::TopologicalTAOCtxImpl::evaluateAdjointProblemImpl
friend MoFEMErrorCode evaluateAdjointProblemImpl(TopologicalTAOCtxImpl *ctx_impl_ptr)
Definition EshelbianTopologicalDerivative.cpp:119

◆ topologicalEvaluateObjectiveAndGradient

MoFEMErrorCode topologicalEvaluateObjectiveAndGradient ( Tao  tao,
Vec  sol,
PetscReal *  f,
Vec  g,
void *  ctx 
)
friend

Definition at line 291 of file EshelbianTopologicalDerivative.cpp.

293 {
294
295 MoFEMFunctionBegin;
296// auto *ctx_impl_ptr = static_cast<TopologicalTAOCtxImpl *>(ctx);
297// if (!ctx_impl_ptr)
298// SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
299// "Invalid context pointer type");
300// auto &ep = *ctx_impl_ptr->ep_ptr;
301// auto &m_field = ep.mField;
302// auto interior_integration_hook = ctx_impl_ptr->integrationAtInterior;
303// auto boundary_integration_hook = ctx_impl_ptr->integrationAtFace;
304// auto alpha = ep.alphaW;
305// auto rho = ep.alphaRho;
306// auto alpha_omega = ep.alphaOmega;
307
308
309 MoFEMFunctionReturn(0);
310}

Member Data Documentation

◆ adjointProblemVec

SmartPetscObj<Vec> EshelbianPlasticity::TopologicalTAOCtxImpl::adjointProblemVec
private

Definition at line 72 of file EshelbianTopologicalDerivative.cpp.

◆ ep_ptr

EshelbianCore* EshelbianPlasticity::TopologicalTAOCtxImpl::ep_ptr
private

Definition at line 45 of file EshelbianTopologicalDerivative.cpp.

◆ integrationAtFace

ForcesAndSourcesCore::GaussHookFun EshelbianPlasticity::TopologicalTAOCtxImpl::integrationAtFace
private

Definition at line 47 of file EshelbianTopologicalDerivative.cpp.

◆ integrationAtInterior

ForcesAndSourcesCore::GaussHookFun EshelbianPlasticity::TopologicalTAOCtxImpl::integrationAtInterior
private

Definition at line 46 of file EshelbianTopologicalDerivative.cpp.

◆ primalProblemVec

SmartPetscObj<Vec> EshelbianPlasticity::TopologicalTAOCtxImpl::primalProblemVec
private

Definition at line 71 of file EshelbianTopologicalDerivative.cpp.

◆ timeSolver

SmartPetscObj<TS> EshelbianPlasticity::TopologicalTAOCtxImpl::timeSolver
private

Definition at line 48 of file EshelbianTopologicalDerivative.cpp.

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