VtkInterface Struct Reference

Collaboration diagram for VtkInterface:

Public Member Functions
	VtkInterface (MoFEM::Interface &m_field, moab::Interface &moab)
MoFEMErrorCode	readVtk ()

Private Member Functions
MoFEMErrorCode	readMesh ()
MoFEMErrorCode	getBoundaryConditions ()
MoFEMErrorCode	setBoundaryConditions ()
MoFEMErrorCode	getMaterialProperties ()
MoFEMErrorCode	getBlockSetBCFromMesh (std::string)
MoFEMErrorCode	writeOutput ()
MoFEMErrorCode	checkResults ()

Private Attributes
MoFEM::Interface &	mField
moab::Interface &	mOab
Simple *	simpleInterface

Detailed Description

Definition at line 11 of file read_vtk.cpp.

Constructor & Destructor Documentation

VtkInterface()

VtkInterface::VtkInterface ( MoFEM::Interface &  m_field, moab::Interface &  moab  )

inline

Definition at line 13 of file read_vtk.cpp.

      : mField(m_field), mOab(moab) {}

Member Function Documentation

checkResults()

MoFEMErrorCode VtkInterface::checkResults ( )

private

Definition at line 237 of file read_vtk.cpp.

                                          {
  MoFEMFunctionBegin;
  PetscBool atom_test = PETSC_FALSE;
  PetscOptionsGetBool(PETSC_NULL, "", "-atom_test", &atom_test, PETSC_NULL);
  if (atom_test) {
    // check mesh has been read correctly
    Range ents;
    CHKERR mOab.get_entities_by_dimension(0, 3, ents);
    if (ents.size() != 125) {
      SETERRQ3(
          PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
          "atom test %d failed: Wrong number of elements %d, should be 125",
          atom_test, ents.size(), 125);
    }
    // check boundary conditions have been set correctly
    MeshsetsManager *meshsets_manager_ptr;
    CHKERR mField.getInterface(meshsets_manager_ptr);
    // iterate meshsets
    for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(mField, BLOCKSET, it)) {
      Range ents;
      CHKERR mField.get_moab().get_entities_by_dimension(it->meshset, 2, ents,
                                                         true);
      if (it->getMeshsetId() == 1) {
        if (it->getName() != "FIX_ALL") {
          SETERRQ3(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
                   "atom test %d failed: Wrong name of meshset %d, should be "
                   "FIX_ALL is %s",
                   atom_test, it->getMeshsetId(), it->getName());
        }
        if (ents.size() != 25) {
          SETERRQ3(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
                   "atom test %d failed: Wrong number of entities in meshset "
                   "%d with %d, should be 25",
                   atom_test, it->getMeshsetId(), ents.size());
        }
      } else if (it->getMeshsetId() == 2) {
        if (it->getName() != "FIX_X") {
          SETERRQ3(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
                   "atom test %d failed: Wrong name of meshset %d, should be "
                   "FIX_X is %s",
                   atom_test, it->getMeshsetId(), it->getName());
        }
        if (ents.size() != 25) {
          SETERRQ3(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
                   "atom test %d failed: Wrong number of entities in meshset "
                   "%d with %d, should be 25",
                   atom_test, it->getMeshsetId(), ents.size());
        }
      } else if (it->getMeshsetId() == 100) {
        if (it->getName() != "ADOLCMAT") {
          SETERRQ3(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
                   "atom test %d failed: Wrong name of meshset %d, should be "
                   "ADOLCMAT is %s",
                   atom_test, it->getMeshsetId(), it->getName());
        }
      }
    }
  }
  MoFEMFunctionReturn(0);
}

getBlockSetBCFromMesh()

MoFEMErrorCode VtkInterface::getBlockSetBCFromMesh ( std::string  block_name )

private

Definition at line 53 of file read_vtk.cpp.

                                                                       {
  MoFEMFunctionBegin;
 
  MeshsetsManager *meshsets_mng;
  CHKERR mField.getInterface(meshsets_mng);
 
  Tag mat_tag;
  rval = mOab.tag_get_handle(block_name.c_str(), 1, MB_TYPE_INTEGER, mat_tag,
                             MB_TAG_EXCL | MB_TAG_SPARSE);
  if (rval != MB_ALREADY_ALLOCATED) {
    // std::cerr << "Error creating tag " << block_name << std::endl;
    MoFEMFunctionReturnHot(0);
  }
  Range block_ents;
  for (auto mbtype : {MBTET, MBTRI, MBHEX, MBQUAD, MBEDGE, MBVERTEX}) {
    CHKERR mOab.get_entities_by_type_and_tag(0, mbtype, &mat_tag, NULL, 1,
                                            block_ents, moab::Interface::UNION);
  }
 
  if (!block_ents.empty()) {
    std::map<int, Range> mat_map;
    for (auto &ent : block_ents) {
      int mat_id;
      CHKERR mOab.tag_get_data(mat_tag, &ent, 1, &mat_id);
      mat_map[mat_id].insert(ent);
    }
 
    for (auto &[mat_id, mat_ents] : mat_map) {
      std::string meshset_name = block_name + "_" + std::to_string(mat_id);
      CHKERR meshsets_mng->addMeshset(BLOCKSET, mat_id, meshset_name);
      CHKERR meshsets_mng->addEntitiesToMeshset(BLOCKSET, mat_id, mat_ents);
      MOFEM_LOG("WORLD", Sev::inform)
          << "Blockset " << meshset_name << " set added.";
    }
  } 
 
  MoFEMFunctionReturn(0);
}

getBoundaryConditions()

MoFEMErrorCode VtkInterface::getBoundaryConditions ( )

private

Definition at line 92 of file read_vtk.cpp.

                                                   {
  MoFEMFunctionBegin;
  MeshsetsManager *meshsets_manager_ptr;
  CHKERR mField.getInterface(meshsets_manager_ptr);
 
  Tag bc_tag;
  CHKERR mOab.tag_get_handle("BOUNDARY_CONDITIONS", 1, MB_TYPE_INTEGER, bc_tag,
                             MB_TAG_CREAT | MB_TAG_SPARSE);
 
  std::vector<std::tuple<int, std::string, std::string>> conditions = {
      {1, "FIX_ALL", "FIX_ALL"},
      {2, "FIX_X", "FIX_X"},
      {3, "FIX_Y", "FIX_Y"},
      {4, "FIX_Z", "FIX_Z"},
      {5, "FIX_X_2", "DISP_X"},
      {6, "FIX_Y_2", "DISP_Y"},
      {7, "FIX_Z_2", "DISP_Z"},
      {8, "FORCE_X", "FORCE_X"},
      {9, "FORCE_Y", "FORCE_Y"},
      {10, "FORCE_Z", "FORCE_Z"},
      {11, "VEL_X", "VEL_X"},
      {12, "VEL_Y", "VEL_Y"},
      {13, "VEL_Z", "VEL_Z"},
      {14, "ACCL_X", "ACCL_X"},
      {15, "ACCL_Y", "ACCL_Y"},
      {16, "ACCL_Z", "ACCL_Z"},
      {17, "TEMP", "TEMP"},
      {18, "PRESSURE", "PRESSURE"},
      {19, "HEAT_FLUX", "HEAT_FLUX"},
      {20, "CONTACT", "CONTACT"},
      {21, "DISPLACEMENT", "DISPLACEMENT"},
      {22, "ROTATE_ALL", "ROTATE_ALL"},
      {23, "ROTATE_X", "ROTATE_X"},
      {24, "ROTATE_Y", "ROTATE_Y"},
      {25, "ROTATE_Z", "ROTATE_Z"},
      {26, "TEMPERATURE", "TEMPERATURE"},
      {50, "TIE_MATRIX", "TIE_MATRIX"}};
 
  for (auto &condition : conditions) {
    int desired_val = std::get<0>(condition);
    const void *desired_val_ptr = &desired_val;
    Range bc;
    CHKERR mOab.get_entities_by_type_and_tag(0, MBVERTEX, &bc_tag,
                                             &desired_val_ptr, 1, bc);
 
    if (!bc.empty()) {
      Range faces;
      CHKERR mOab.get_adjacencies(bc, 2, true, faces, moab::Interface::UNION);
      Range edges;
      CHKERR mOab.get_adjacencies(bc, 1, true, edges, moab::Interface::UNION);
 
      Range adj_nodes;
      CHKERR mOab.get_adjacencies(faces, 0, false, adj_nodes,
                                  moab::Interface::UNION);
 
      Range non_BC_nodes;
      non_BC_nodes = subtract(adj_nodes, bc);
      Range non_BC_faces;
      CHKERR mOab.get_adjacencies(non_BC_nodes, 2, false, non_BC_faces,
                                  moab::Interface::UNION);
      Range non_BC_edges;
      CHKERR mOab.get_adjacencies(non_BC_nodes, 1, false, non_BC_edges,
                                  moab::Interface::UNION);
 
      faces = subtract(faces, non_BC_faces);
      edges = subtract(edges, non_BC_edges);
 
      MOFEM_LOG_C("WORLD", Sev::inform, "Number of nodes assigned to %s = %d",
                  std::get<2>(condition).c_str(), bc.size());
      MOFEM_LOG_C("WORLD", Sev::inform, "Number of edges assigned to %s = %d",
                  std::get<2>(condition).c_str(), edges.size());
      MOFEM_LOG_C("WORLD", Sev::inform, "Number of faces assigned to %s = %d",
                  std::get<2>(condition).c_str(), faces.size());
 
      CHKERR meshsets_manager_ptr->addMeshset(BLOCKSET, desired_val,
                                              std::get<1>(condition));
      CHKERR meshsets_manager_ptr->addEntitiesToMeshset(BLOCKSET, desired_val,
                                                        bc);
      CHKERR meshsets_manager_ptr->addEntitiesToMeshset(BLOCKSET, desired_val,
                                                        edges);
      CHKERR meshsets_manager_ptr->addEntitiesToMeshset(BLOCKSET, desired_val,
                                                        faces);
    } else {
      CHKERR getBlockSetBCFromMesh(std::get<1>(condition));
    }
  }
 
  MoFEMFunctionReturn(0);
}

getMaterialProperties()

MoFEMErrorCode VtkInterface::getMaterialProperties ( )

private

Definition at line 201 of file read_vtk.cpp.

                                                   {
  MoFEMFunctionBegin;
  MeshsetsManager *meshsets_manager_ptr;
  CHKERR mField.getInterface(meshsets_manager_ptr);
 
  Range ents;
  rval = mOab.get_entities_by_dimension(0, 3, ents, true);
  if (rval != MB_SUCCESS) {
    MOFEM_LOG("WORLD", Sev::warning)
        << "No hexes/tets in the mesh, no material block set. Not Implemented";
  } else {
    CHKERR getBlockSetBCFromMesh("MAT_ELASTIC");
    // for backward compatibility
    CHKERR meshsets_manager_ptr->addMeshset(BLOCKSET, 100, "ADOLCMAT");
    CHKERR meshsets_manager_ptr->addEntitiesToMeshset(BLOCKSET, 100, ents);
    MOFEM_LOG("WORLD", Sev::inform) << "Material block ADOLCMAT set added.";
  }
  MoFEMFunctionReturn(0);
}

readMesh()

MoFEMErrorCode VtkInterface::readMesh ( )

private

Definition at line 43 of file read_vtk.cpp.

                                      {
  MoFEMFunctionBegin;
 
  CHKERR mField.getInterface(simpleInterface);
  CHKERR simpleInterface->getOptions();
  CHKERR simpleInterface->loadFile();
 
  MoFEMFunctionReturn(0);
}

readVtk()

MoFEMErrorCode VtkInterface::readVtk

(

)

Definition at line 32 of file read_vtk.cpp.

                                     {
  MoFEMFunctionBegin;
  CHKERR readMesh();
  CHKERR getBoundaryConditions();
  CHKERR getMaterialProperties();
  CHKERR setBoundaryConditions();
  CHKERR writeOutput();
  CHKERR checkResults();
  MoFEMFunctionReturn(0);
}

setBoundaryConditions()

MoFEMErrorCode VtkInterface::setBoundaryConditions ( )

private

Definition at line 182 of file read_vtk.cpp.

                                                   {
  MoFEMFunctionBegin;
  // Add meshsets if config file provided
  MeshsetsManager *meshsets_interface_ptr;
  CHKERR mField.getInterface(meshsets_interface_ptr);
  CHKERR meshsets_interface_ptr->setMeshsetFromFile();
 
  MOFEM_LOG_CHANNEL("WORLD");
  MOFEM_LOG_TAG("WORLD", "read_vtk")
  MOFEM_LOG("WORLD", Sev::inform)
      << "Print all meshsets (old and added from meshsets "
         "configurational file)";
  for (auto cit = meshsets_interface_ptr->getBegin();
       cit != meshsets_interface_ptr->getEnd(); cit++)
    MOFEM_LOG("WORLD", Sev::inform) << *cit;
 
  MoFEMFunctionReturn(0);
}

writeOutput()

MoFEMErrorCode VtkInterface::writeOutput ( )

private

Definition at line 221 of file read_vtk.cpp.

                                         {
  MoFEMFunctionBegin;
 
  char mesh_out_file[255] = "out_vtk.h5m";
 
  CHKERR PetscOptionsBegin(mField.get_comm(), "", "Read vtk tool", "none");
  CHKERR PetscOptionsString("-output_file", "output mesh file name", "",
                            "out.h5m", mesh_out_file, 255, PETSC_NULL);
  ierr = PetscOptionsEnd();
  CHKERRQ(ierr);
 
  CHKERR mOab.write_file(mesh_out_file);
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

mField

MoFEM::Interface& VtkInterface::mField

private

Definition at line 19 of file read_vtk.cpp.

mOab

moab::Interface& VtkInterface::mOab

private

Definition at line 20 of file read_vtk.cpp.

simpleInterface

Simple* VtkInterface::simpleInterface

private

Definition at line 21 of file read_vtk.cpp.

---

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

• tools/read_vtk.cpp