AdaptiveOrderRef.hpp File Reference

Calculate error indicators. More...

Go to the source code of this file.

Classes
struct	OpCalculateTractionJump

Functions
static MoFEMErrorCode	getTagHandle (MoFEM::Interface &m_field, const char *name, DataType type, Tag &tag_handle)
static MoFEMErrorCode	getDomainEntsPart (MoFEM::Interface &m_field, Range &domain_ents_part)
static MoFEMErrorCode	copyTagOnSkin (MoFEM::Interface &m_field, const char *name, DataType type)

Detailed Description

Calculate error indicators.

Date

2024-04-14

Copyright

Copyright (c) 2025

Definition in file AdaptiveOrderRef.hpp.

Function Documentation

copyTagOnSkin()

static MoFEMErrorCode copyTagOnSkin ( MoFEM::Interface &  m_field, const char *  name, DataType  type  )

static

Definition at line 45 of file AdaptiveOrderRef.hpp.

                                                   {
  MoFEMFunctionBegin;
  auto simple = m_field.getInterface<Simple>();
 
  Range all_boundary_ents;
  CHKERR m_field.get_moab().get_entities_by_handle(simple->getBoundaryMeshSet(),
                                                   all_boundary_ents);
 
  Range domain_ents_part;
  CHKERR getDomainEntsPart(m_field, domain_ents_part);
  Range domain_faces, boundary_faces;
  CHKERR m_field.get_moab().get_adjacencies(domain_ents_part, SPACE_DIM - 1,
                                            true, domain_faces,
                                            moab::Interface::UNION);
  boundary_faces = intersect(domain_faces, all_boundary_ents);
 
  Tag th;
  CHKERR getTagHandle(m_field, name, type, th);
 
  for (auto &face_ent : boundary_faces) {
    Range adj_tet;
    CHKERR m_field.get_moab().get_adjacencies(&face_ent, 1, SPACE_DIM, true,
                                              adj_tet, moab::Interface::UNION);
    if (adj_tet.size() != 1)
      CHK_THROW_MESSAGE(
          MOFEM_DATA_INCONSISTENCY,
          "There should be only one ent adjacent to the boundary ent");
 
    std::variant<int, double> data;
 
    switch (type) {
    case MB_TYPE_INTEGER:
      data = int{};
      break;
    case MB_TYPE_DOUBLE:
      data = double{};
      break;
    default:
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
               "Wrong data type %d for tag", type);
    }
 
    std::visit(
        [&](auto &val) {
          CHKERR m_field.get_moab().tag_get_data(th, adj_tet, &val);
          CHKERR m_field.get_moab().tag_set_data(th, &face_ent, 1, &val);
        },
        data);
  }
 
  MoFEMFunctionReturn(0);
}

getDomainEntsPart()

static MoFEMErrorCode getDomainEntsPart ( MoFEM::Interface &  m_field, Range &  domain_ents_part  )

static

Definition at line 31 of file AdaptiveOrderRef.hpp.

                                                                 {
  MoFEMFunctionBegin;
  auto simple = m_field.getInterface<Simple>();
  Range domain_ents;
  CHKERR m_field.get_moab().get_entities_by_dimension(simple->getMeshset(),
                                                      SPACE_DIM, domain_ents);
  ParallelComm *pcomm =
      ParallelComm::get_pcomm(&m_field.get_moab(), MYPCOMM_INDEX);
  CHKERR pcomm->filter_pstatus(domain_ents, PSTATUS_NOT_OWNED, PSTATUS_NOT, -1,
                               &domain_ents_part);
  MoFEMFunctionReturn(0);
}

getTagHandle()

static MoFEMErrorCode getTagHandle ( MoFEM::Interface &  m_field, const char *  name, DataType  type, Tag &  tag_handle  )

static

Definition at line 10 of file AdaptiveOrderRef.hpp.

                                                                   {
  MoFEMFunctionBegin;
  int int_val = 0;
  double double_val = 0;
  switch (type) {
  case MB_TYPE_INTEGER:
    CHKERR m_field.get_moab().tag_get_handle(
        name, 1, type, tag_handle, MB_TAG_CREAT | MB_TAG_SPARSE, &int_val);
    break;
  case MB_TYPE_DOUBLE:
    CHKERR m_field.get_moab().tag_get_handle(
        name, 1, type, tag_handle, MB_TAG_CREAT | MB_TAG_SPARSE, &double_val);
    break;
  default:
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
             "Wrong data type %d for tag", type);
  }
  MoFEMFunctionReturn(0);
}