ElasticTie::OpTiePrismLhs Struct Reference

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

Inheritance diagram for ElasticTie::OpTiePrismLhs:

Collaboration diagram for ElasticTie::OpTiePrismLhs:

Public Member Functions
	OpTiePrismLhs (const std::string &row_field_name, const std::string &col_field_name, const bool assemble_transpose=false)
double	getSlaveFaceArea ()

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

Protected Attributes
std::string	rowFieldName
std::string	colFieldName

Detailed Description

Definition at line 602 of file ElasticTie.hpp.

Constructor & Destructor Documentation

OpTiePrismLhs()

ElasticTie::OpTiePrismLhs::OpTiePrismLhs ( const std::string &  row_field_name, const std::string &  col_field_name, const bool  assemble_transpose = false  )

inline

Definition at line 604 of file ElasticTie.hpp.

      : TiePrismLhsBase(row_field_name, col_field_name, OPROWCOL, nullptr),
        rowFieldName(row_field_name), colFieldName(col_field_name) {
    sYmm = false;
    this->assembleTranspose = assemble_transpose;
  }

Member Function Documentation

getSlaveFaceArea()

double ElasticTie::OpTiePrismLhs::getSlaveFaceArea ( )

inline

Definition at line 613 of file ElasticTie.hpp.

{ return getAreaF3(); }

iNtegrate()

MoFEMErrorCode ElasticTie::OpTiePrismLhs::iNtegrate ( EntitiesFieldData::EntData &  row_data, EntitiesFieldData::EntData &  col_data  )

inlineoverrideprotected

Definition at line 619 of file ElasticTie.hpp.

                                                                        {
    MoFEMFunctionBeginHot;
 
    const auto *fe =
        dynamic_cast<const FlatPrismElementForcesAndSourcesCore *>(ptrFE);
    if (!fe)
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "TIE prism operator is attached to a wrong finite element");
 
    if (nbRows % SPACE_DIM || nbCols % SPACE_DIM)
      SETERRQ(
          PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
          "Unexpected number of TIE prism vector dofs on rows or columns: "
          "rowField=%s colField=%s rowSide=%d colSide=%d rowType=%d colType=%d "
          "nbRows=%d nbCols=%d",
          rowFieldName.c_str(), colFieldName.c_str(), rowSide, colSide,
          rowType, colType, nbRows, nbCols);
 
    if (colFieldName != "LAMBDA")
      MoFEMFunctionReturnHot(0);
 
    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;
    const int nb_row_base = nbRows / SPACE_DIM;
    const int nb_col_base = nbCols / SPACE_DIM;
 
    auto &col_n = col_data.getN();
    if (!col_n.size1() || !col_n.size2())
      MoFEMFunctionReturnHot(0);
 
    auto &row_n = row_data.getN();
    if (!row_n.size1() || !row_n.size2())
      MoFEMFunctionReturnHot(0);
 
    if (row_n.size1() != nb_gauss_pts || col_n.size1() != nb_gauss_pts)
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "Unexpected number of Gauss points on TIE prism blocks");
 
    if (static_cast<size_t>(nb_row_base) > row_n.size2() ||
        static_cast<size_t>(nb_col_base) > col_n.size2())
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "More TIE prism bases than shape functions");
    if (locMat.size1() != static_cast<size_t>(nbRows) ||
        locMat.size2() != static_cast<size_t>(nbCols))
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "Unexpected TIE LHS block size: rowField=%s colField=%s "
              "rowSide=%d colSide=%d rowType=%d colType=%d nbRows=%d "
              "nbCols=%d locMat=(%d,%d) row_n=(%d,%d) col_n=(%d,%d)",
              rowFieldName.c_str(), colFieldName.c_str(), rowSide, colSide,
              rowType, colType, nbRows, nbCols,
              static_cast<int>(locMat.size1()), static_cast<int>(locMat.size2()),
              static_cast<int>(row_n.size1()), static_cast<int>(row_n.size2()),
              static_cast<int>(col_n.size1()), static_cast<int>(col_n.size2()));
 
    if (!isOnPrismTieTrace(rowSide, rowType))
      MoFEMFunctionReturnHot(0);
 
    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 LHS Gauss index out of range: rowField=%s colField=%s "
                  "rowSide=%d rowType=%d rr=%d row_gg=%d row_n.size1()=%d "
                  "nb_slave_gauss=%d",
                  rowFieldName.c_str(), colFieldName.c_str(), 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);
        auto t_col_base = col_data.getFTensor0N(gg, 0);
        for (int cc = 0; cc != nb_col_base; ++cc) {
          const double value = row_sign * slave_alpha * t_row_base * t_col_base;
          for (int dd = 0; dd != SPACE_DIM; ++dd)
            locMat(SPACE_DIM * rr + dd, SPACE_DIM * cc + dd) += value;
          ++t_col_base;
        }
        ++t_row_base;
      }
      ++t_w;
    }
 
    MoFEMFunctionReturnHot(0);
  }

Member Data Documentation

colFieldName

std::string ElasticTie::OpTiePrismLhs::colFieldName

protected

Definition at line 617 of file ElasticTie.hpp.

rowFieldName

std::string ElasticTie::OpTiePrismLhs::rowFieldName

protected

Definition at line 616 of file ElasticTie.hpp.

---

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

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