Compute objective function contributions at element level.
Evaluates Python objective function with current displacement and stress state, and accumulates global objective value and gradients.
1577 {
1579
1580
1585
1588
1590 2;
1591
1592 auto t_diff_symm = FTensor::diff_symmetrize<double>();
1593
1594
1595 auto nb_gauss_pts =
1596 getGaussPts().size2();
1597 auto objective_ptr = boost::make_shared<MatrixDouble>(
1598 1, nb_gauss_pts);
1599 auto objective_dstress = boost::make_shared<MatrixDouble>(
1600 nb_gauss_pts,
1601 symm_size);
1602 auto objective_dstrain = boost::make_shared<MatrixDouble>(
1603 nb_gauss_pts,
1604 symm_size);
1605 auto objective_du =
1606 boost::make_shared<MatrixDouble>(nb_gauss_pts,
SPACE_DIM);
1607
1608
1609 auto evaluate_python = [&]() {
1611 auto &coords = OP::getCoordsAtGaussPts();
1614 objective_ptr);
1617 objective_dstress);
1620 objective_dstrain);
1623 objective_du);
1624
1625 auto t_grad_u =
1626 getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*(
commPtr->matGradPtr));
1627 auto t_D =
1628 getFTensor4DdgFromMat<SPACE_DIM, SPACE_DIM, 0>(*(
commPtr->matDPtr));
1629 auto t_jac = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*(
jacPtr));
1630 auto t_diff_jac = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*(
diffJacPtr));
1632 auto t_d_grad = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*(
dGradPtr));
1633
1635 auto t_obj_dstress =
1636 getFTensor2SymmetricFromMat<SPACE_DIM>(*objective_dstress);
1637 auto t_obj_dstrain =
1638 getFTensor2SymmetricFromMat<SPACE_DIM>(*objective_dstrain);
1639 auto t_obj_du = getFTensor1FromMat<SPACE_DIM>(*objective_du);
1640 auto t_d_u = getFTensor1FromMat<SPACE_DIM>(*
dUPtr);
1641
1642 auto vol = OP::getMeasure();
1643 auto t_w = getFTensor0IntegrationWeight();
1644 for (auto gg = 0; gg != nb_gauss_pts; ++gg) {
1645
1649
1651 t_diff_inv_jac(
i,
j) =
1652 -(t_inv_jac(
i,
I) * t_diff_jac(
I,
J)) * t_inv_jac(
J,
j);
1654 t_diff_grad(
i,
j) = t_grad_u(
i,
k) * t_diff_inv_jac(
k,
j);
1655
1657 t_d_strain(
i,
j) = t_diff_symm(
i,
j,
k,
l) * (
1658
1660
1661 +
1662
1664
1665 );
1666
1667 auto alpha = t_w * vol;
1668
1669 (*globObjectivePtr) += alpha * t_obj;
1670 (*globObjectiveGradPtr) +=
1671 alpha *
1672 (
1673
1674 t_obj_dstress(
i,
j) * (t_D(
i,
j,
k,
l) * t_d_strain(
k,
l))
1675
1676 +
1677
1678 t_obj_dstrain(
i,
j) * t_d_strain(
i,
j)
1679
1680 +
1681
1682 t_obj_du(
i) * t_d_u(
i)
1683
1684 +
1685
1686 t_obj * t_cof
1687
1688 );
1689
1690 ++t_w;
1691 ++t_jac;
1692 ++t_diff_jac;
1693 ++t_cof;
1694
1695 ++t_obj;
1696 ++t_obj_dstress;
1697 ++t_obj_dstrain;
1698 ++t_obj_du;
1699
1700 ++t_grad_u;
1701 ++t_d_grad;
1702 ++t_d_u;
1703 }
1705 };
1706
1707 CHKERR evaluate_python();
1708
1710 }
#define FTENSOR_INDEX(DIM, I)
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
#define CHKERR
Inline error check.
constexpr int SPACE_DIM
[Define dimension]
constexpr IntegrationType I
Use Gauss quadrature for integration.
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Index< 'J', DIM1 > J
FTensor::Index< 'l', 3 > l
FTensor::Index< 'j', 3 > j
FTensor::Index< 'k', 3 > k
static MoFEMErrorCode invertTensor(FTensor::Tensor2< T1, DIM, DIM > &t, T2 &det, FTensor::Tensor2< T3, DIM, DIM > &inv_t)
static auto determinantTensor(FTensor::Tensor2< T, DIM, DIM > &t)
Calculate the determinant of a tensor of rank DIM.
auto getFTensor0FromMat(M &data)
Get tensor rank 0 (scalar) form data vector.
static auto getFTensor0FromVec(V &data)
Get tensor rank 0 (scalar) form data vector.