MortarContactStructures Struct Reference

#include <users_modules/mortar_contact/src/MakeStructures.hpp>

Static Public Member Functions
static MoFEMErrorCode	masterSlaveTrianglesCreation (MoFEM::Interface &m_field, Range &range_surf_master, Range &range_surf_slave, int step=1)
static MoFEMErrorCode	findContactSurfacePairs (MoFEM::Interface &m_field, std::vector< std::pair< Range, Range >> &contact_surface_pairs)

Detailed Description

Definition at line 18 of file MakeStructures.hpp.

Member Function Documentation

findContactSurfacePairs()

static MoFEMErrorCode MortarContactStructures::findContactSurfacePairs ( MoFEM::Interface &  m_field, std::vector< std::pair< Range, Range >> &  contact_surface_pairs  )

inlinestatic

Definition at line 43 of file MakeStructures.hpp.

                                                               {
    MoFEMFunctionBegin;
 
    std::map<std::string, Range> slave_surface_map, master_surface_map;
 
    auto fill_surface_map = [&](auto &surface_map, std::string surface_name) {
      MoFEMFunctionBegin;
 
      for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
        if (it->getName().compare(0, surface_name.size(), surface_name) == 0) {
          Range tris;
          CHKERR m_field.get_moab().get_entities_by_type(it->meshset, MBTRI,
                                                         tris, true);
          std::string key = it->getName().substr(surface_name.size());
          auto res =
              surface_map.insert(std::pair<std::string, Range>(key, tris));
          if (!res.second) {
            SETERRQ1(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                     "Contact surface name %s used more than once",
                     it->getName().c_str());
          }
        }
      }
 
      MoFEMFunctionReturn(0);
    };
 
    CHKERR fill_surface_map(slave_surface_map, "MORTAR_SLAVE");
    CHKERR fill_surface_map(master_surface_map, "MORTAR_MASTER");
 
    for (const auto &[key, slave_range] : slave_surface_map) {
      auto master_search = master_surface_map.find(key);
      if (master_search != master_surface_map.end()) {
        contact_surface_pairs.push_back(
            std::pair<Range, Range>(slave_range, master_search->second));
        master_surface_map.erase(master_search);
      } else {
        SETERRQ2(
            PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "Cannot find pairing contact surface MORTAR_MASTER%s for existing "
            "surface MORTAR_SLAVE%s",
            key.c_str(), key.c_str());
      }
    }
 
    if (!master_surface_map.empty()) {
      std::string key = master_surface_map.begin()->first;
      SETERRQ2(
          PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
          "Cannot find pairing contact surface MORTAR_SLAVE%s for existing "
          "surface MORTAR_MASTER%s",
          key.c_str(), key.c_str());
    }
 
    MoFEMFunctionReturn(0);
  }

masterSlaveTrianglesCreation()

static MoFEMErrorCode MortarContactStructures::masterSlaveTrianglesCreation ( MoFEM::Interface &  m_field, Range &  range_surf_master, Range &  range_surf_slave, int  step = 1  )

inlinestatic

Definition at line 20 of file MakeStructures.hpp.

                                                                   {
    MoFEMFunctionBegin;
    for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
      if (it->getName().compare(0, 13, "MORTAR_MASTER") == 0) {
        CHKERR m_field.get_moab().get_entities_by_type(it->meshset, MBTRI,
                                                       range_surf_master, true);
      }
    }
 
    for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
      if (it->getName().compare(0, 12, "MORTAR_SLAVE") == 0) {
        rval = m_field.get_moab().get_entities_by_type(it->meshset, MBTRI,
                                                       range_surf_slave, true);
        CHKERRQ_MOAB(rval);
      }
    }
 
    MoFEMFunctionReturn(0);
  }

---

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

• users_modules/mortar_contact/src/MakeStructures.hpp