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

Public Types

using OP = DomainBaseOp
 
using OP = DomainBaseOp
 

Public Member Functions

 OpStateSensitivity (const std::string field_name, boost::shared_ptr< ObjectiveFunctionData > python_ptr, boost::shared_ptr< HookeOps::CommonData > comm_ptr, boost::shared_ptr< MatrixDouble > u_ptr)
 
MoFEMErrorCode iNtegrate (EntitiesFieldData::EntData &data)
 
 OpStateSensitivity (const std::string field_name, boost::shared_ptr< ObjectiveFunctionData > python_ptr, boost::shared_ptr< HookeOps::CommonData > comm_ptr, boost::shared_ptr< MatrixDouble > u_ptr)
 
MoFEMErrorCode iNtegrate (EntitiesFieldData::EntData &data)
 

Private Attributes

boost::shared_ptr< ObjectiveFunctionData > pythonPtr
 
boost::shared_ptr< HookeOps::CommonData > commPtr
 
boost::shared_ptr< MatrixDouble > uPtr
 

Detailed Description

Examples
mofem/tutorials/vec-7_shape_optimisation/adjoint.cpp.

Definition at line 1461 of file adjoint.cpp.

Member Typedef Documentation

◆ OP [1/2]

using OpStateSensitivity::OP = DomainBaseOp

Definition at line 1462 of file adjoint.cpp.

◆ OP [2/2]

using OpStateSensitivity::OP = DomainBaseOp

Definition at line 706 of file gradient.cpp.

Constructor & Destructor Documentation

◆ OpStateSensitivity() [1/2]

OpStateSensitivity::OpStateSensitivity ( const std::string  field_name,
boost::shared_ptr< ObjectiveFunctionData >  python_ptr,
boost::shared_ptr< HookeOps::CommonData >  comm_ptr,
boost::shared_ptr< MatrixDouble >  u_ptr 
)
inline

Definition at line 1463 of file adjoint.cpp.

1467 : OP(field_name, field_name, DomainEleOp::OPROW), pythonPtr(python_ptr),
1468 commPtr(comm_ptr), uPtr(u_ptr) {}
field_name
constexpr auto field_name
Definition poisson_2d_homogeneous.cpp:13
OpStateSensitivity::commPtr
boost::shared_ptr< HookeOps::CommonData > commPtr
Definition adjoint.cpp:1548
OpStateSensitivity::OP
DomainBaseOp OP
Definition adjoint.cpp:1462
OpStateSensitivity::uPtr
boost::shared_ptr< MatrixDouble > uPtr
Definition adjoint.cpp:1549
OpStateSensitivity::pythonPtr
boost::shared_ptr< ObjectiveFunctionData > pythonPtr
Definition adjoint.cpp:1547

◆ OpStateSensitivity() [2/2]

OpStateSensitivity::OpStateSensitivity ( const std::string  field_name,
boost::shared_ptr< ObjectiveFunctionData >  python_ptr,
boost::shared_ptr< HookeOps::CommonData >  comm_ptr,
boost::shared_ptr< MatrixDouble >  u_ptr 
)
inline

Definition at line 707 of file gradient.cpp.

711 : OP(field_name, field_name, DomainEleOp::OPROW), pythonPtr(python_ptr),
712 commPtr(comm_ptr), uPtr(u_ptr) {}

Member Function Documentation

◆ iNtegrate() [1/2]

MoFEMErrorCode OpStateSensitivity::iNtegrate ( EntitiesFieldData::EntData &  data)
inline
Examples
mofem/tutorials/vec-7_shape_optimisation/adjoint.cpp.

Definition at line 1470 of file adjoint.cpp.

1470 {
1471 MoFEMFunctionBegin;
1472 FTENSOR_INDEX(SPACE_DIM, i);
1473 FTENSOR_INDEX(SPACE_DIM, j);
1474 FTENSOR_INDEX(SPACE_DIM, k);
1475 FTENSOR_INDEX(SPACE_DIM, l);
1476
1477 constexpr auto symm_size = (SPACE_DIM * (SPACE_DIM + 1)) / 2;
1478 auto nb_gauss_pts = getGaussPts().size2();
1479
1480 auto objective_dstress =
1481 boost::make_shared<MatrixDouble>(nb_gauss_pts, symm_size);
1482 auto objective_dstrain =
1483 boost::make_shared<MatrixDouble>(nb_gauss_pts, symm_size);
1484 auto objective_du =
1485 boost::make_shared<MatrixDouble>(nb_gauss_pts, SPACE_DIM);
1486
1487 auto evaluate_python = [&]() {
1488 MoFEMFunctionBegin;
1489 auto &coords = OP::getCoordsAtGaussPts();
1490 CHKERR pythonPtr->evalInteriorObjectiveGradientStress(
1491 coords, uPtr, commPtr->getMatCauchyStress(), commPtr->getMatStrain(),
1492 objective_dstress);
1493 CHKERR pythonPtr->evalInteriorObjectiveGradientStrain(
1494 coords, uPtr, commPtr->getMatCauchyStress(), commPtr->getMatStrain(),
1495 objective_dstrain);
1496 CHKERR pythonPtr->evalInteriorObjectiveGradientU(
1497 coords, uPtr, commPtr->getMatCauchyStress(), commPtr->getMatStrain(),
1498 objective_du);
1499
1500 auto vol = OP::getMeasure();
1501 auto t_w = OP::getFTensor0IntegrationWeight();
1502
1503 auto t_D =
1504 getFTensor4DdgFromMat<SPACE_DIM, SPACE_DIM, 0>(*(commPtr->matDPtr));
1505 auto t_row_grad = data.getFTensor1DiffN<SPACE_DIM>();
1506 auto t_row_base = data.getFTensor0N();
1507
1508 auto t_obj_dstress =
1509 getFTensor2SymmetricFromMat<SPACE_DIM>(*objective_dstress);
1510 auto t_obj_dstrain =
1511 getFTensor2SymmetricFromMat<SPACE_DIM>(*objective_dstrain);
1512 auto t_obj_du = getFTensor1FromMat<SPACE_DIM>(*objective_du);
1513
1514 for (int gg = 0; gg != nb_gauss_pts; ++gg) {
1515 const double alpha = t_w * vol;
1516 FTensor::Tensor2<double, SPACE_DIM, SPACE_DIM> t_adjoint_stress;
1517 t_adjoint_stress(i, j) =
1518 t_D(i, j, k, l) * t_obj_dstress(k, l) + t_obj_dstrain(i, j);
1519
1520 auto t_nf = OP::template getNf<SPACE_DIM>();
1521 int rr = 0;
1522 for (; rr != OP::nbRows / SPACE_DIM; rr++) {
1523 t_nf(j) += alpha * t_row_grad(i) * t_adjoint_stress(i, j);
1524 t_nf(j) += alpha * t_row_base * t_obj_du(j);
1525
1526 ++t_row_grad;
1527 ++t_row_base;
1528 ++t_nf;
1529 }
1530
1531 for (; rr < OP::nbRowBaseFunctions; ++rr) {
1532 ++t_row_grad;
1533 ++t_row_base;
1534 }
1535 ++t_obj_dstrain;
1536 ++t_obj_dstress;
1537 ++t_obj_du;
1538 ++t_w;
1539 }
1540 MoFEMFunctionReturn(0);
1541 };
1542 CHKERR evaluate_python();
1543 MoFEMFunctionReturn(0);
1544 }
FTENSOR_INDEX
#define FTENSOR_INDEX(DIM, I)
Definition Templates.hpp:2673
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
SPACE_DIM
constexpr int SPACE_DIM
[Define dimension]
Definition gradient.cpp:19
i
FTensor::Index< 'i', SPACE_DIM > i
Definition hcurl_divergence_operator_2d.cpp:27
l
FTensor::Index< 'l', 3 > l
Definition matrix_function.cpp:23
j
FTensor::Index< 'j', 3 > j
Definition matrix_function.cpp:21
k
FTensor::Index< 'k', 3 > k
Definition matrix_function.cpp:22
MoFEM::EntitiesFieldData::EntData::getFTensor0N
FTensor::Tensor0< FTensor::PackPtr< double *, 1 > > getFTensor0N(const FieldApproximationBase base)
Get base function as Tensor0.
Definition EntitiesFieldData.hpp:1806
MoFEM::EntitiesFieldData::EntData::getFTensor1DiffN
auto getFTensor1DiffN(const FieldApproximationBase base)
Get derivatives of base functions.
Definition EntitiesFieldData.hpp:802

◆ iNtegrate() [2/2]

MoFEMErrorCode OpStateSensitivity::iNtegrate ( EntitiesFieldData::EntData &  data)
inline

Definition at line 714 of file gradient.cpp.

714 {
715 MoFEMFunctionBegin;
716 FTENSOR_INDEX(SPACE_DIM, i);
717 FTENSOR_INDEX(SPACE_DIM, j);
718 FTENSOR_INDEX(SPACE_DIM, k);
719 FTENSOR_INDEX(SPACE_DIM, l);
720
721 constexpr auto symm_size = (SPACE_DIM * (SPACE_DIM + 1)) / 2;
722 auto nb_gauss_pts = getGaussPts().size2();
723
724 auto objective_dstress =
725 boost::make_shared<MatrixDouble>(nb_gauss_pts, symm_size);
726 auto objective_dstrain =
727 boost::make_shared<MatrixDouble>(nb_gauss_pts, symm_size);
728 auto objective_du =
729 boost::make_shared<MatrixDouble>(nb_gauss_pts, SPACE_DIM);
730
731 auto evaluate_python = [&]() {
732 MoFEMFunctionBegin;
733 auto &coords = OP::getCoordsAtGaussPts();
734 CHKERR pythonPtr->evalInteriorObjectiveGradientStress(
735 coords, uPtr, commPtr->getMatCauchyStress(), commPtr->getMatStrain(),
736 objective_dstress);
737 CHKERR pythonPtr->evalInteriorObjectiveGradientStrain(
738 coords, uPtr, commPtr->getMatCauchyStress(), commPtr->getMatStrain(),
739 objective_dstrain);
740 CHKERR pythonPtr->evalInteriorObjectiveGradientU(
741 coords, uPtr, commPtr->getMatCauchyStress(), commPtr->getMatStrain(),
742 objective_du);
743
744 auto vol = OP::getMeasure();
745 auto t_w = OP::getFTensor0IntegrationWeight();
746
747 auto t_D =
748 getFTensor4DdgFromMat<SPACE_DIM, SPACE_DIM, 0>(*(commPtr->matDPtr));
749 auto t_row_grad = data.getFTensor1DiffN<SPACE_DIM>();
750 auto t_row_base = data.getFTensor0N();
751
752 auto t_obj_dstress =
753 getFTensor2SymmetricFromMat<SPACE_DIM>(*objective_dstress);
754 auto t_obj_dstrain =
755 getFTensor2SymmetricFromMat<SPACE_DIM>(*objective_dstrain);
756 auto t_obj_du = getFTensor1FromMat<SPACE_DIM>(*objective_du);
757
758 for (int gg = 0; gg != nb_gauss_pts; ++gg) {
759 const double alpha = t_w * vol;
760 FTensor::Tensor2<double, SPACE_DIM, SPACE_DIM> t_adjoint_stress;
761 t_adjoint_stress(i, j) =
762 t_D(i, j, k, l) * t_obj_dstress(k, l) + t_obj_dstrain(i, j);
763
764 auto t_nf = OP::template getNf<SPACE_DIM>();
765 int rr = 0;
766 for (; rr != OP::nbRows / SPACE_DIM; rr++) {
767 t_nf(j) += alpha * t_row_grad(i) * t_adjoint_stress(i, j);
768 t_nf(j) += alpha * t_row_base * t_obj_du(j);
769
770 ++t_row_grad;
771 ++t_row_base;
772 ++t_nf;
773 }
774
775 for (; rr < OP::nbRowBaseFunctions; ++rr) {
776 ++t_row_grad;
777 ++t_row_base;
778 }
779 ++t_obj_dstrain;
780 ++t_obj_dstress;
781 ++t_obj_du;
782 ++t_w;
783 }
784 MoFEMFunctionReturn(0);
785 };
786 CHKERR evaluate_python();
787 MoFEMFunctionReturn(0);
788 }

Member Data Documentation

◆ commPtr

boost::shared_ptr< HookeOps::CommonData > OpStateSensitivity::commPtr
private
Examples
mofem/tutorials/vec-7_shape_optimisation/adjoint.cpp.

Definition at line 1548 of file adjoint.cpp.

◆ pythonPtr

boost::shared_ptr< ObjectiveFunctionData > OpStateSensitivity::pythonPtr
private
Examples
mofem/tutorials/vec-7_shape_optimisation/adjoint.cpp.

Definition at line 1547 of file adjoint.cpp.

◆ uPtr

boost::shared_ptr< MatrixDouble > OpStateSensitivity::uPtr
private
Examples
mofem/tutorials/vec-7_shape_optimisation/adjoint.cpp.

Definition at line 1549 of file adjoint.cpp.

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