ElasticTie::OpTiePrismRhsU Struct Reference

#include "tutorials/vec-11_elastic_tie_mesh/src/ElasticTie.hpp"

Inheritance diagram for ElasticTie::OpTiePrismRhsU:

Collaboration diagram for ElasticTie::OpTiePrismRhsU:

Public Member Functions
	OpTiePrismRhsU (const std::string &field_name, boost::shared_ptr< MatrixDouble > lambda_ptr)
double	getSlaveFaceArea ()

Protected Member Functions
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data) override

Private Attributes
boost::shared_ptr< MatrixDouble >	lambdaPtr

Detailed Description

Definition at line 651 of file ElasticTie.hpp.

Constructor & Destructor Documentation

OpTiePrismRhsU()

ElasticTie::OpTiePrismRhsU::OpTiePrismRhsU ( const std::string &  field_name, boost::shared_ptr< MatrixDouble >  lambda_ptr  )

inline

Definition at line 653 of file ElasticTie.hpp.

      : TiePrismRhsBase(field_name, field_name, OPROW, nullptr),
        lambdaPtr(lambda_ptr) {
    if (!lambdaPtr)
      THROW_MESSAGE("Pointers are not set");
  }

Member Function Documentation

getSlaveFaceArea()

double ElasticTie::OpTiePrismRhsU::getSlaveFaceArea ( )

inline

Definition at line 661 of file ElasticTie.hpp.

{ return getAreaF3(); }

iNtegrate()

MoFEMErrorCode ElasticTie::OpTiePrismRhsU::iNtegrate ( EntitiesFieldData::EntData &  row_data )

inlineoverrideprotected

Definition at line 664 of file ElasticTie.hpp.

                                                                        {
    MoFEMFunctionBeginHot;
 
    const auto *fe =
        dynamic_cast<const FlatPrismElementForcesAndSourcesCore *>(ptrFE);
    if (!fe)
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "TIE prism RHS operator is attached to a wrong finite element");
 
    const auto nb_gauss_pts = getGaussPts().size2();
    if (!nb_gauss_pts || nb_gauss_pts % 2) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "Expected concatenated [slave|master] TIE Gauss points");
    }
    const auto nb_slave_gauss = nb_gauss_pts / 2;
    if (lambdaPtr->size1() != nb_gauss_pts) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "TIE multiplier values at Gauss points were not prepared "
              "before RHS assembly");
    }
 
    if (nbRows % SPACE_DIM)
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "Unexpected number of TIE prism vector dofs on rows");
 
    if (!isOnPrismTieTrace(rowSide, rowType))
      MoFEMFunctionReturnHot(0);
 
    auto &row_n = row_data.getN();
    if (!row_n.size1() || !row_n.size2())
      MoFEMFunctionReturnHot(0);
    if (row_n.size1() != nb_gauss_pts)
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "Unexpected number of Gauss points on TIE displacement block");
 
    const int nb_row_base = nbRows / SPACE_DIM;
    if (static_cast<size_t>(nb_row_base) > row_n.size2())
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "More TIE displacement bases than shape functions");
    if (lambdaPtr->size2() != SPACE_DIM)
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "Unexpected TIE multiplier matrix size: size2=%d expected=%d",
              static_cast<int>(lambdaPtr->size2()), SPACE_DIM);
    if (locF.size() != static_cast<size_t>(nbRows))
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "Unexpected TIE U RHS block size: rowSide=%d rowType=%d "
              "nbRows=%d locF.size()=%d row_n=(%d,%d) lambda=(%d,%d)",
              rowSide, rowType, nbRows, static_cast<int>(locF.size()),
              static_cast<int>(row_n.size1()), static_cast<int>(row_n.size2()),
              static_cast<int>(lambdaPtr->size1()),
              static_cast<int>(lambdaPtr->size2()));
 
    auto t_w = getFTensor0IntegrationWeight();
    for (size_t gg = 0; gg != nb_slave_gauss; ++gg) {
      const double slave_alpha = getSlaveFaceArea() * t_w;
      for (int rr = 0; rr != nb_row_base; ++rr) {
        const bool base_on_slave = isOnSlaveTieRowBase(rr, rowSide, rowType);
        const double row_sign = base_on_slave ? 1.0 : -1.0;
        const auto row_gg = base_on_slave ? gg : gg + nb_slave_gauss;
        if (row_gg >= row_n.size1())
          SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                  "TIE U RHS Gauss index out of range: rowSide=%d rowType=%d "
                  "rr=%d row_gg=%d row_n.size1()=%d nb_slave_gauss=%d",
                  rowSide, rowType, rr, static_cast<int>(row_gg),
                  static_cast<int>(row_n.size1()),
                  static_cast<int>(nb_slave_gauss));
        auto t_row_base = row_data.getFTensor0N(row_gg, rr);
        const double row_base = t_row_base;
        if (row_base != 0.) {
          auto t_row_rhs = getFTensor1FromPtr<SPACE_DIM>(&locF[SPACE_DIM * rr]);
          for (int dd = 0; dd != SPACE_DIM; ++dd)
            t_row_rhs(dd) +=
                row_sign * row_base * slave_alpha * (*lambdaPtr)(gg, dd);
        }
        ++t_row_base;
      }
 
      ++t_w;
    }
 
    MoFEMFunctionReturnHot(0);
  }

Member Data Documentation

lambdaPtr

boost::shared_ptr<MatrixDouble> ElasticTie::OpTiePrismRhsU::lambdaPtr

private

Definition at line 748 of file ElasticTie.hpp.

---

The documentation for this struct was generated from the following file:

• tutorials/vec-11_elastic_tie_mesh/src/ElasticTie.hpp