v0.14.0
Public Member Functions | List of all members
EshelbianPlasticity::OpHMHH Struct Reference
Inheritance diagram for EshelbianPlasticity::OpHMHH:
[legend]
Collaboration diagram for EshelbianPlasticity::OpHMHH:
[legend]

Public Member Functions

 OpHMHH (const int tag, const bool eval_rhs, const bool eval_lhs, boost::shared_ptr< DataAtIntegrationPts > data_ptr, boost::shared_ptr< PhysicalEquations > physics_ptr)
 
MoFEMErrorCode evaluateRhs (EntData &data)
 
MoFEMErrorCode evaluateLhs (EntData &data)
 
- Public Member Functions inherited from OpJacobian
 OpJacobian (const int tag, const bool eval_rhs, const bool eval_lhs, boost::shared_ptr< DataAtIntegrationPts > data_ptr, boost::shared_ptr< PhysicalEquations > physics_ptr)
 
MoFEMErrorCode doWork (int side, EntityType type, EntData &data)
 

Additional Inherited Members

- Protected Member Functions inherited from OpJacobian
 OpJacobian ()
 
- Protected Attributes inherited from OpJacobian
int tAg = -1
 adol-c tape More...
 
bool evalRhs = false
 
bool evalLhs = false
 
boost::shared_ptr< DataAtIntegrationPts > dataAtPts
 data at integration pts More...
 
boost::shared_ptr< PhysicalEquations > physicsPtr
 material physical equations More...
 

Detailed Description

Definition at line 306 of file EshelbianADOL-C.cpp.

Constructor & Destructor Documentation

◆ OpHMHH()

EshelbianPlasticity::OpHMHH::OpHMHH ( const int  tag,
const bool  eval_rhs,
const bool  eval_lhs,
boost::shared_ptr< DataAtIntegrationPts >  data_ptr,
boost::shared_ptr< PhysicalEquations >  physics_ptr 
)
inline

Definition at line 307 of file EshelbianADOL-C.cpp.

310  : OpJacobian(tag, eval_rhs, eval_lhs, data_ptr, physics_ptr) {}

Member Function Documentation

◆ evaluateLhs()

MoFEMErrorCode EshelbianPlasticity::OpHMHH::evaluateLhs ( EntData &  data)
virtual

Implements OpJacobian.

Definition at line 390 of file EshelbianADOL-C.cpp.

390  {
391  MoFEMFunctionBegin;
392  FTensor::Index<'i', 3> i;
393  FTensor::Index<'j', 3> j;
394  FTensor::Index<'k', 3> k;
395  FTensor::Index<'l', 3> l;
396  FTensor::Index<'m', 3> m;
397  FTensor::Index<'n', 3> n;
398  FTensor::Index<'o', 3> o;
399  FTensor::Index<'p', 3> p;
400 
401  const int number_of_active_variables = physicsPtr->activeVariables.size();
402  const int number_of_dependent_variables =
403  physicsPtr->dependentVariablesPiola.size();
404  std::vector<double *> jac_ptr(number_of_dependent_variables);
405  for (unsigned int n = 0; n != number_of_dependent_variables; ++n) {
406  jac_ptr[n] =
407  &(physicsPtr
408  ->dependentVariablesPiolaDirevatives[n *
409  number_of_active_variables]);
410  }
411 
412  const auto nb_integration_pts = getGaussPts().size2();
413 
414  auto create_data_mat = [nb_integration_pts](auto &m) {
415  m.resize(9, nb_integration_pts, false);
416  };
417 
418  dataAtPts->P_du.resize(81, nb_integration_pts, false);
419 
420  auto iu = getFTensor2SymmetricFromMat<3>(dataAtPts->stretchTensorAtPts);
421  auto t_grad_h1 = getFTensor2FromMat<3, 3>(dataAtPts->wGradH1AtPts);
422  auto r_P_du = getFTensor4FromMat<3, 3, 3, 3>(dataAtPts->P_du);
423 
424  FTensor::Number<0> N0;
425  FTensor::Number<1> N1;
426  FTensor::Number<2> N2;
427 
428  constexpr auto t_kd = FTensor::Kronecker_Delta<int>();
429 
430  for (unsigned int gg = 0; gg != nb_integration_pts; ++gg) {
431 
432  FTensor::Tensor2<double, 3, 3> t_h1;
433 
434  switch (EshelbianCore::gradApproximator) {
435  case NO_H1_CONFIGURATION:
436  t_h1(i, j) = t_kd(i, j);
437  physicsPtr->get_h()(i, j) = iu(i, k) * t_h1(k, j);
438  break;
439  case LARGE_ROT:
440  case MODERATE_ROT:
441  t_h1(i, j) = t_grad_h1(i, j) + t_kd(i, j);
442  physicsPtr->get_h()(i, j) = iu(i, k) * t_h1(k, j);
443  break;
444  case SMALL_ROT:
445  physicsPtr->get_h()(i, j) = iu(i, j);
446  break;
447  };
448 
449  // play recorder for jacobians
450  int r = ::jacobian(tAg, number_of_dependent_variables,
451  number_of_active_variables,
452  &physicsPtr->activeVariables[0], &jac_ptr[0]);
453  if (r < 0) {
454  SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
455  "ADOL-C function evaluation with error");
456  }
457 
458  FTensor::Tensor4<double, 3, 3, 3, 3> t_P_dh_tmp;
459  t_P_dh_tmp(i, j, N0, k) = physicsPtr->get_P_dh0()(i, j, k);
460  t_P_dh_tmp(i, j, N1, k) = physicsPtr->get_P_dh1()(i, j, k);
461  t_P_dh_tmp(i, j, N2, k) = physicsPtr->get_P_dh2()(i, j, k);
462 
463  switch (EshelbianCore::gradApproximator) {
464  case NO_H1_CONFIGURATION: {
465  FTensor::Tensor4<double, 3, 3, 3, 3> t_h1_du;
466  t_h1_du(i, j, m, n) = t_kd(i, m) * t_kd(j, n);
467  r_P_du(k, l, m, n) =
468  (t_P_dh_tmp(i, j, o, p) * t_h1_du(o, p, m, n)) * t_h1_du(i, j, k, l);
469  } break;
470  case LARGE_ROT: {
471  FTensor::Tensor4<double, 3, 3, 3, 3> t_h1_du;
472  t_h1_du(i, j, m, n) = t_kd(i, m) * (t_kd(k, n) * t_h1(k, j));
473  r_P_du(k, l, m, n) =
474  (t_P_dh_tmp(i, j, o, p) * t_h1_du(o, p, m, n)) * t_h1_du(i, j, k, l);
475  } break;
476  case MODERATE_ROT:
477  case SMALL_ROT:
478  r_P_du(i, j, m, n) = t_P_dh_tmp(i, j, m, n);
479  break;
480  };
481 
482  ++iu;
483  ++t_grad_h1;
484  ++r_P_du;
485  }
486  MoFEMFunctionReturn(0);
487 }

◆ evaluateRhs()

MoFEMErrorCode EshelbianPlasticity::OpHMHH::evaluateRhs ( EntData &  data)
virtual

Implements OpJacobian.

Definition at line 316 of file EshelbianADOL-C.cpp.

316  {
317  MoFEMFunctionBegin;
318  FTensor::Index<'i', 3> i;
319  FTensor::Index<'j', 3> j;
320  FTensor::Index<'k', 3> k;
321  FTensor::Index<'l', 3> l;
322  FTensor::Index<'m', 3> m;
323  FTensor::Index<'n', 3> n;
324 
325  const auto nb_integration_pts = getGaussPts().size2();
326  auto iu = getFTensor2SymmetricFromMat<3>(dataAtPts->stretchTensorAtPts);
327  auto t_grad_h1 = getFTensor2FromMat<3, 3>(dataAtPts->wGradH1AtPts);
328 
329  auto create_data_vec = [nb_integration_pts](auto &v) {
330  v.resize(nb_integration_pts, false);
331  };
332  auto create_data_mat = [nb_integration_pts](auto &m) {
333  m.resize(9, nb_integration_pts, false);
334  };
335 
336  create_data_mat(dataAtPts->PAtPts);
337 
338  constexpr auto t_kd = FTensor::Kronecker_Delta<int>();
339  auto r_P = getFTensor2FromMat<3, 3>(dataAtPts->PAtPts);
340  for (unsigned int gg = 0; gg != nb_integration_pts; ++gg) {
341 
342  FTensor::Tensor2<double, 3, 3> t_h1;
343 
344  switch (EshelbianCore::gradApproximator) {
345  case NO_H1_CONFIGURATION:
346  t_h1(i, j) = t_kd(i, j);
347  physicsPtr->get_h()(i, j) = iu(i, j);
348  break;
349  case LARGE_ROT:
350  case MODERATE_ROT:
351  t_h1(i, j) = t_grad_h1(i, j) + t_kd(i, j);
352  physicsPtr->get_h()(i, j) = iu(i, k) * t_h1(k, j);
353  break;
354  case SMALL_ROT:
355  physicsPtr->get_h()(i, j) = iu(i, j);
356  break;
357  };
358 
359  int r = ::function(tAg, physicsPtr->dependentVariablesPiola.size(),
360  physicsPtr->activeVariables.size(),
361  &physicsPtr->activeVariables[0],
362  &physicsPtr->dependentVariablesPiola[0]);
363  if (r < 0) { // function is locally analytic
364  SETERRQ1(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
365  "ADOL-C function evaluation with error r = %d", r);
366  }
367 
368  switch (EshelbianCore::gradApproximator) {
369  case NO_H1_CONFIGURATION:
370  r_P(i, j) = physicsPtr->get_P()(i, j);
371  break;
372  case LARGE_ROT:
373  case MODERATE_ROT: {
374  FTensor::Tensor4<double, 3, 3, 3, 3> t_h1_du;
375  t_h1_du(i, j, m, n) = t_kd(i, m) * (t_kd(k, n) * t_h1(k, j));
376  r_P(k, l) = physicsPtr->get_P()(i, j) * t_h1_du(i, j, k, l);
377  }; break;
378  case SMALL_ROT:
379  r_P(i, j) = physicsPtr->get_P()(i, j);
380  break;
381  };
382 
383  ++iu;
384  ++t_grad_h1;
385  ++r_P;
386  }
387  MoFEMFunctionReturn(0);
388 }
The documentation for this struct was generated from the following file:
OpJacobian::tAg
int tAg
adol-c tape
Definition: EshelbianOperators.hpp:30
FTensor::Kronecker_Delta
Kronecker Delta class.
Definition: Kronecker_Delta.hpp:15
EshelbianPlasticity::NO_H1_CONFIGURATION
@ NO_H1_CONFIGURATION
Definition: EshelbianPlasticity.hpp:45
sdf.r
int r
Definition: sdf.py:8
FTensor::Tensor2< double, 3, 3 >
EshelbianCore::gradApproximator
static enum RotSelector gradApproximator
Definition: EshelbianCore.hpp:16
FTensor::Number< 0 >
OpJacobian::dataAtPts
boost::shared_ptr< DataAtIntegrationPts > dataAtPts
data at integration pts
Definition: EshelbianOperators.hpp:35
OpJacobian::OpJacobian
OpJacobian()
Definition: EshelbianOperators.hpp:28
MOFEM_OPERATION_UNSUCCESSFUL
@ MOFEM_OPERATION_UNSUCCESSFUL
Definition: definitions.h:34
FTensor::Tensor4
Definition: Tensor4_value.hpp:18
i
FTensor::Index< 'i', SPACE_DIM > i
Definition: hcurl_divergence_operator_2d.cpp:27
t_kd
constexpr auto t_kd
Definition: free_surface.cpp:139
FTensor::Index< 'i', 3 >
OpJacobian::physicsPtr
boost::shared_ptr< PhysicalEquations > physicsPtr
material physical equations
Definition: EshelbianOperators.hpp:37
convert.n
n
Definition: convert.py:82
v
const double v
phase velocity of light in medium (cm/ns)
Definition: initial_diffusion.cpp:40
EshelbianPlasticity::SMALL_ROT
@ SMALL_ROT
Definition: EshelbianPlasticity.hpp:45
EshelbianPlasticity::LARGE_ROT
@ LARGE_ROT
Definition: EshelbianPlasticity.hpp:45
j
FTensor::Index< 'j', 3 > j
Definition: matrix_function.cpp:19
m
FTensor::Index< 'm', 3 > m
Definition: shallow_wave.cpp:80
EshelbianPlasticity::MODERATE_ROT
@ MODERATE_ROT
Definition: EshelbianPlasticity.hpp:45
k
FTensor::Index< 'k', 3 > k
Definition: matrix_function.cpp:20
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:429
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:359
l
FTensor::Index< 'l', 3 > l
Definition: matrix_function.cpp:21
Generated by Doxygen 1.8.17 and hosted at University of Glasgow