Compute objective function contributions at element level.
Evaluates Python objective function with current displacement and stress state, and accumulates global objective value and gradients.
1599 {
1601
1602
1607
1610
1612
1614
1615 auto nb_gauss_pts = getGaussPts().size2();
1616 auto objective_ptr = boost::make_shared<VectorDouble>(nb_gauss_pts);
1617 auto objective_dstress =
1618 boost::make_shared<MatrixDouble>(symm_size, nb_gauss_pts);
1619 auto objective_dstrain =
1620 boost::make_shared<MatrixDouble>(symm_size, nb_gauss_pts);
1621 auto objective_du = boost::make_shared<MatrixDouble>(
SPACE_DIM, nb_gauss_pts);
1622
1623
1624 auto evaluate_python = [&]() {
1626 auto &coords = OP::getCoordsAtGaussPts();
1629 objective_ptr);
1632 objective_dstress);
1635 objective_dstrain);
1638 objective_du);
1639
1640 auto t_grad_u =
1641 getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*(
commPtr->matGradPtr));
1642 auto t_D =
1643 getFTensor4DdgFromMat<SPACE_DIM, SPACE_DIM, 0>(*(
commPtr->matDPtr));
1644 auto t_jac = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*(
jacPtr));
1645 auto t_diff_jac = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*(
diffJacPtr));
1647 auto t_d_grad = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*(
dGradPtr));
1648
1650 auto t_obj_dstress =
1651 getFTensor2SymmetricFromMat<SPACE_DIM>(*objective_dstress);
1652 auto t_obj_dstrain =
1653 getFTensor2SymmetricFromMat<SPACE_DIM>(*objective_dstrain);
1654 auto t_obj_du = getFTensor1FromMat<SPACE_DIM>(*objective_du);
1655 auto t_d_u = getFTensor1FromMat<SPACE_DIM>(*
dUPtr);
1656
1657 auto vol = OP::getMeasure();
1658 auto t_w = getFTensor0IntegrationWeight();
1659 for (auto gg = 0; gg != nb_gauss_pts; ++gg) {
1660
1664
1666 t_diff_inv_jac(
i,
j) =
1667 -(t_inv_jac(
i,
I) * t_diff_jac(
I,
J)) * t_inv_jac(
J,
j);
1669 t_diff_grad(
i,
j) = t_grad_u(
i,
k) * t_diff_inv_jac(
k,
j);
1670
1672 t_d_strain(
i,
j) = t_diff_symm(
i,
j,
k,
l) * (
1673
1675
1676 +
1677
1679
1680 );
1681
1682 auto alpha = t_w * vol;
1683
1684 (*globObjectivePtr) += alpha * t_obj;
1685 (*globObjectiveGradPtr) +=
1686 alpha *
1687 (
1688
1689 t_obj_dstress(
i,
j) * (t_D(
i,
j,
k,
l) * t_d_strain(
k,
l))
1690
1691 +
1692
1693 t_obj_dstrain(
i,
j) * t_d_strain(
i,
j)
1694
1695 +
1696
1697 t_obj_du(
i) * t_d_u(
i)
1698
1699 +
1700
1701 t_obj * t_cof
1702
1703 );
1704
1705 ++t_w;
1706 ++t_jac;
1707 ++t_diff_jac;
1708 ++t_cof;
1709
1710 ++t_obj;
1711 ++t_obj_dstress;
1712 ++t_obj_dstrain;
1713 ++t_obj_du;
1714
1715 ++t_grad_u;
1716 ++t_d_grad;
1717 ++t_d_u;
1718 }
1720 };
1721
1722 CHKERR evaluate_python();
1723
1725 }
#define FTENSOR_INDEX(DIM, I)
auto diff_symmetrize(FTensor::Number< DIM >)
constexpr int SPACE_DIM
[Define dimension]
constexpr IntegrationType I
Use Gauss quadrature for integration.
#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.
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.
static auto getFTensor0FromVec(ublas::vector< T, A > &data)
Get tensor rank 0 (scalar) form data vector.
