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

Public Member Functions

 CohesiveTAOCtxImpl (EshelbianCore *ep, ForcesAndSourcesCore::GaussHookFun set_integration_at_front_face, SmartPetscObj< TS > time_solver)
 
SmartPetscObj< Vec > duplicateKappaVec () override
 
SmartPetscObj< Vec > duplicateDissipationVec () override
 
SmartPetscObj< Vec > duplicateGradientVec () override
 
CommInterface::EntitiesPetscVector & getKappaVec () override
 
- Public Member Functions inherited from EshelbianPlasticity::CohesiveTAOCtx
 CohesiveTAOCtx ()=default
 

Private Attributes

EshelbianCore * ep_ptr
 
ForcesAndSourcesCore::GaussHookFun setIntegrationAtFrontFace
 
SmartPetscObj< TS > timeSolver
 
SmartPetscObj< Vec > kspSolVec
 
SmartPetscObj< Vec > lambdaVec
 
CommInterface::EntitiesPetscVector kappaVec
 
CommInterface::EntitiesPetscVector dissipationVec
 
CommInterface::EntitiesPetscVector gradDissipationVec
 

Friends

MoFEMErrorCode cohesiveEvaluateObjectiveAndGradient (Tao tao, Vec x, PetscReal *f, Vec g, void *ctx)
 

Detailed Description

Definition at line 1270 of file EshelbianCohesive.cpp.

Constructor & Destructor Documentation

◆ CohesiveTAOCtxImpl()

EshelbianPlasticity::CohesiveTAOCtxImpl::CohesiveTAOCtxImpl ( EshelbianCore *  ep,
ForcesAndSourcesCore::GaussHookFun  set_integration_at_front_face,
SmartPetscObj< TS >  time_solver 
)
inline

Definition at line 1271 of file EshelbianCohesive.cpp.

1275 : ep_ptr(ep), setIntegrationAtFrontFace(set_integration_at_front_face),
1276 timeSolver(time_solver),
1277
1278 kspSolVec(createDMVector(ep->dmElastic)),
1279 lambdaVec(createDMVector(ep->dmElastic)),
1280 kappaVec(CommInterface::createEntitiesPetscVector(
1281 ep->mField.get_comm(), ep->mField.get_moab(),
1282 [&ep](Range r) { return intersect(*(ep->interfaceFaces), r); }, 1,
1283 Sev::noisy)),
1284 dissipationVec(CommInterface::createEntitiesPetscVector(
1285 ep->mField.get_comm(), ep->mField.get_moab(),
1286 [&ep](Range r) { return intersect(*(ep->interfaceFaces), r); }, 1,
1287 Sev::noisy)),
1288 gradDissipationVec(CommInterface::createEntitiesPetscVector(
1289 ep->mField.get_comm(), ep->mField.get_moab(),
1290 [&ep](Range r) { return intersect(*(ep->interfaceFaces), r); }, 1,
1291 Sev::noisy)) {}
MoFEM::createDMVector
auto createDMVector(DM dm)
Get smart vector from DM.
Definition DMMoFEM.hpp:1234
EshelbianCore::mField
MoFEM::Interface & mField
Definition EshelbianCore.hpp:191
EshelbianCore::dmElastic
SmartPetscObj< DM > dmElastic
Elastic problem.
Definition EshelbianCore.hpp:205
EshelbianPlasticity::CohesiveTAOCtxImpl::setIntegrationAtFrontFace
ForcesAndSourcesCore::GaussHookFun setIntegrationAtFrontFace
Definition EshelbianCohesive.cpp:1311
EshelbianPlasticity::CohesiveTAOCtxImpl::gradDissipationVec
CommInterface::EntitiesPetscVector gradDissipationVec
Definition EshelbianCohesive.cpp:1318
EshelbianPlasticity::CohesiveTAOCtxImpl::kappaVec
CommInterface::EntitiesPetscVector kappaVec
Definition EshelbianCohesive.cpp:1316
EshelbianPlasticity::CohesiveTAOCtxImpl::dissipationVec
CommInterface::EntitiesPetscVector dissipationVec
Definition EshelbianCohesive.cpp:1317
MoFEM::CommInterface::createEntitiesPetscVector
static EntitiesPetscVector createEntitiesPetscVector(MPI_Comm comm, moab::Interface &moab, std::function< Range(Range)> get_entities_fun, const int nb_coeffs, Sev sev=Sev::verbose, int root_rank=0, bool get_vertices=true)
Create a ghost vector for exchanging data.
Definition CommInterface.cpp:1460
MoFEM::CoreInterface::get_moab
virtual moab::Interface & get_moab()=0
MoFEM::CoreInterface::get_comm
virtual MPI_Comm & get_comm() const =0

Member Function Documentation

◆ duplicateDissipationVec()

SmartPetscObj< Vec > EshelbianPlasticity::CohesiveTAOCtxImpl::duplicateDissipationVec ( )
inlineoverridevirtual

Implements EshelbianPlasticity::CohesiveTAOCtx.

Definition at line 1297 of file EshelbianCohesive.cpp.

1297 {
1298 return vectorDuplicate(dissipationVec.second);
1299 }
MoFEM::vectorDuplicate
SmartPetscObj< Vec > vectorDuplicate(Vec vec)
Create duplicate vector of smart vector.
Definition PetscSmartObj.hpp:223

◆ duplicateGradientVec()

SmartPetscObj< Vec > EshelbianPlasticity::CohesiveTAOCtxImpl::duplicateGradientVec ( )
inlineoverridevirtual

Implements EshelbianPlasticity::CohesiveTAOCtx.

Definition at line 1301 of file EshelbianCohesive.cpp.

1301 {
1302 return vectorDuplicate(gradDissipationVec.second);
1303 }

◆ duplicateKappaVec()

SmartPetscObj< Vec > EshelbianPlasticity::CohesiveTAOCtxImpl::duplicateKappaVec ( )
inlineoverridevirtual

Implements EshelbianPlasticity::CohesiveTAOCtx.

Definition at line 1293 of file EshelbianCohesive.cpp.

1293 {
1294 return vectorDuplicate(kappaVec.second);
1295 }

◆ getKappaVec()

CommInterface::EntitiesPetscVector & EshelbianPlasticity::CohesiveTAOCtxImpl::getKappaVec ( )
inlineoverridevirtual

Implements EshelbianPlasticity::CohesiveTAOCtx.

Definition at line 1305 of file EshelbianCohesive.cpp.

1305 {
1306 return kappaVec;
1307 }

Friends And Related Symbol Documentation

◆ cohesiveEvaluateObjectiveAndGradient

MoFEMErrorCode cohesiveEvaluateObjectiveAndGradient ( Tao  tao,
Vec  x,
PetscReal *  f,
Vec  g,
void *  ctx 
)
friend

Definition at line 1332 of file EshelbianCohesive.cpp.

1334 {
1335
1336 MoFEMFunctionBegin;
1337 auto cohesive_ctx = static_cast<CohesiveTAOCtxImpl *>(ctx);
1338 auto &ep = *(cohesive_ctx->ep_ptr);
1339
1340#ifndef NDEBUG
1341 CHKERR VecView(delta_kappa, PETSC_VIEWER_STDOUT_WORLD);
1342#endif
1343 // Set delta_kappa values to the tag
1344 CHKERR VecCopy(delta_kappa, cohesive_ctx->kappaVec.second);
1345 CHKERR VecGhostUpdateBegin(cohesive_ctx->kappaVec.second, INSERT_VALUES,
1346 SCATTER_FORWARD);
1347 CHKERR VecGhostUpdateEnd(cohesive_ctx->kappaVec.second, INSERT_VALUES,
1348 SCATTER_FORWARD);
1349 CHKERR CommInterface::setTagFromVector(
1350 ep.mField.get_moab(), cohesive_ctx->kappaVec,
1351 get_delta_kappa_tag(ep.mField.get_moab()));
1352
1353 // solve primal problem
1354 CHKERR evaluatePrimalProblemCohesiveImpl(ep, cohesive_ctx->timeSolver,
1355 cohesive_ctx->kspSolVec);
1356 // solve adjoint problem to get lambda
1357 auto ksp = snesGetKSP(tsGetSNES(cohesive_ctx->timeSolver));
1358 CHKERR evaluateCohesiveLambdaImpl(ep, cohesive_ctx->setIntegrationAtFrontFace,
1359 ksp, cohesive_ctx->lambdaVec);
1360 // evaluate dissipation and its gradient
1361 CHKERR evaluateDissipationAndGradImpl(
1362 ep, cohesive_ctx->setIntegrationAtFrontFace, cohesive_ctx->lambdaVec,
1363 cohesive_ctx->dissipationVec, cohesive_ctx->gradDissipationVec);
1364
1365 CHKERR VecSum(cohesive_ctx->dissipationVec.second, f);
1366 CHKERR VecCopy(cohesive_ctx->gradDissipationVec.second, g);
1367
1368 MOFEM_LOG("EP", Sev::inform)
1369 << "Cohesive objective function (negative total dissipation): " << *f;
1370
1371 MoFEMFunctionReturn(0);
1372}
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_LOG
#define MOFEM_LOG(channel, severity)
Log.
Definition LogManager.hpp:316
EshelbianPlasticity::get_delta_kappa_tag
Tag get_delta_kappa_tag(moab::Interface &moab)
Definition EshelbianCohesive.cpp:278
EshelbianPlasticity::evaluateDissipationAndGradImpl
static MoFEMErrorCode evaluateDissipationAndGradImpl(EshelbianCore &ep, ForcesAndSourcesCore::GaussHookFun set_integration_at_front_face, SmartPetscObj< Vec > lambda_vec, CommInterface::EntitiesPetscVector &dissipation_vec, CommInterface::EntitiesPetscVector &grad_dissipation_vec)
Definition EshelbianCohesive.cpp:1120
EshelbianPlasticity::evaluatePrimalProblemCohesiveImpl
static MoFEMErrorCode evaluatePrimalProblemCohesiveImpl(EshelbianCore &ep, SmartPetscObj< TS > ts, SmartPetscObj< Vec > x)
Definition EshelbianCohesive.cpp:1240
EshelbianPlasticity::evaluateCohesiveLambdaImpl
static MoFEMErrorCode evaluateCohesiveLambdaImpl(EshelbianCore &ep, ForcesAndSourcesCore::GaussHookFun set_integration_at_front_face, SmartPetscObj< KSP > ksp, SmartPetscObj< Vec > lambda_vec)
Definition EshelbianCohesive.cpp:1175
MoFEM::snesGetKSP
auto snesGetKSP(SNES snes)
Definition PetscSmartObj.hpp:416
MoFEM::tsGetSNES
auto tsGetSNES(TS ts)
Definition PetscSmartObj.hpp:410
g
constexpr double g
Definition shallow_wave.cpp:63
EshelbianPlasticity::CohesiveTAOCtxImpl::CohesiveTAOCtxImpl
CohesiveTAOCtxImpl(EshelbianCore *ep, ForcesAndSourcesCore::GaussHookFun set_integration_at_front_face, SmartPetscObj< TS > time_solver)
Definition EshelbianCohesive.cpp:1271
MoFEM::CommInterface::setTagFromVector
static MoFEMErrorCode setTagFromVector(moab::Interface &moab, EntitiesPetscVector &vec, Tag tag)
Definition CommInterface.cpp:1694

Member Data Documentation

◆ dissipationVec

CommInterface::EntitiesPetscVector EshelbianPlasticity::CohesiveTAOCtxImpl::dissipationVec
private

Definition at line 1317 of file EshelbianCohesive.cpp.

◆ ep_ptr

EshelbianCore* EshelbianPlasticity::CohesiveTAOCtxImpl::ep_ptr
private

Definition at line 1310 of file EshelbianCohesive.cpp.

◆ gradDissipationVec

CommInterface::EntitiesPetscVector EshelbianPlasticity::CohesiveTAOCtxImpl::gradDissipationVec
private

Definition at line 1318 of file EshelbianCohesive.cpp.

◆ kappaVec

CommInterface::EntitiesPetscVector EshelbianPlasticity::CohesiveTAOCtxImpl::kappaVec
private

Definition at line 1316 of file EshelbianCohesive.cpp.

◆ kspSolVec

SmartPetscObj<Vec> EshelbianPlasticity::CohesiveTAOCtxImpl::kspSolVec
private

Definition at line 1314 of file EshelbianCohesive.cpp.

◆ lambdaVec

SmartPetscObj<Vec> EshelbianPlasticity::CohesiveTAOCtxImpl::lambdaVec
private

Definition at line 1315 of file EshelbianCohesive.cpp.

◆ setIntegrationAtFrontFace

ForcesAndSourcesCore::GaussHookFun EshelbianPlasticity::CohesiveTAOCtxImpl::setIntegrationAtFrontFace
private

Definition at line 1311 of file EshelbianCohesive.cpp.

◆ timeSolver

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

Definition at line 1312 of file EshelbianCohesive.cpp.

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