34 ParallelComm *pcomm = ParallelComm::get_pcomm(&moab,
MYPCOMM_INDEX);
36 pcomm =
new ParallelComm(&moab, PETSC_COMM_WORLD);
40 char mesh_out_file[255] =
"out.h5m";
41 PetscBool flg_file = PETSC_FALSE;
42 PetscBool flg_n_part = PETSC_FALSE;
43 PetscBool flg_part = PETSC_FALSE;
44 PetscBool only_tags = PETSC_FALSE;
45 PetscInt n_partas = 1;
46 PetscBool create_lower_dim_ents = PETSC_TRUE;
50 ierr = PetscOptionsBegin(PETSC_COMM_WORLD,
"",
"none",
"none");
53 CHKERR PetscOptionsString(
"-my_file",
"mesh file name",
"",
"mesh.h5m",
55 if (flg_file != PETSC_TRUE)
56 CHKERR PetscOptionsString(
"-file_name",
"mesh file name",
"",
"mesh.h5m",
58 if (flg_file != PETSC_TRUE)
59 SETERRQ(PETSC_COMM_SELF, 1,
60 "*** ERROR -my_file (-file_name) (MESH FILE NEEDED)");
66 CHKERR PetscOptionsString(
"-output_file",
"output mesh file name",
"",
67 "out.h5m", mesh_out_file, 255, PETSC_NULL);
68 CHKERR PetscOptionsInt(
"-my_nparts",
"number of parts",
"", n_partas,
69 &n_partas, &flg_n_part);
70 CHKERR PetscOptionsInt(
"-nparts",
"number of parts",
"", n_partas,
71 &n_partas, &flg_part);
73 if (!flg_n_part && !flg_part)
74 SETERRQ(PETSC_COMM_SELF, 1,
75 "*** ERROR partitioning number not given (-nparts)");
77 auto get_nb_ents_by_dim = [&](
const int dim) {
79 CHKERR moab.get_number_entities_by_dimension(0,
dim, nb);
83 if (get_nb_ents_by_dim(
dim))
89 "Dimension of entities to partition not found");
91 CHKERR PetscOptionsInt(
"-dim",
"entities dim",
"",
dim, &
dim, PETSC_NULL);
93 CHKERR PetscOptionsInt(
"-adj_dim",
"adjacency dim",
"", adj_dim, &adj_dim,
96 "-my_create_lower_dim_ents",
"if tru create lower dimension entireties",
97 "", create_lower_dim_ents, &create_lower_dim_ents, PETSC_NULL);
98 CHKERR PetscOptionsBool(
"-block_tags",
"obly block and meshsests tags",
"",
99 only_tags, &only_tags, PETSC_NULL);
101 ierr = PetscOptionsEnd();
109 CHKERR meshsets_interface_ptr->setMeshsetFromFile();
114 meshsets_interface_ptr->getBegin();
115 cit != meshsets_interface_ptr->getEnd(); cit++)
121 CHKERR moab.get_entities_by_dimension(0,
dim, ents_dim,
false);
122 if (create_lower_dim_ents) {
125 CHKERR moab.get_adjacencies(ents_dim, 2,
true, faces,
126 moab::Interface::UNION);
130 CHKERR moab.get_adjacencies(ents_dim, 1,
true, edges,
131 moab::Interface::UNION);
136 ents_dim,
dim, adj_dim, n_partas,
nullptr,
nullptr,
nullptr,
VERBOSE);
139 auto get_tag_list = [&]() {
140 std::vector<Tag> tags_list;
143 for (
auto &
m : list) {
144 auto meshset =
m.getMeshset();
145 std::vector<Tag> tmp_tags_list;
148 for (
auto t : tmp_tags_list) {
149 tags_list.push_back(
t);
159 auto tags_list = get_tag_list();
160 std::sort(tags_list.begin(), tags_list.end());
161 auto new_end = std::unique(tags_list.begin(), tags_list.end());
162 tags_list.resize(std::distance(tags_list.begin(), new_end));
163 tags_list.push_back(pcomm->part_tag());
165 CHKERR moab.write_file(mesh_out_file,
"MOAB",
"", &root_mesh, 1,
166 &*tags_list.begin(), tags_list.size());
168 CHKERR moab.write_file(mesh_out_file,
"MOAB",
"");
#define CATCH_ERRORS
Catch errors.
#define MYPCOMM_INDEX
default communicator number PCOMM
@ MOFEM_DATA_INCONSISTENCY
#define CHKERR
Inline error check.
#define MOFEM_LOG(channel, severity)
Log.
SeverityLevel
Severity levels.
#define MOFEM_LOG_TAG(channel, tag)
Tag channel.
#define MOFEM_LOG_CHANNEL(channel)
Set and reset channel.
MoFEMErrorCode partitionMesh(const Range &ents, const int dim, const int adj_dim, const int n_parts, Tag *th_vertex_weights=nullptr, Tag *th_edge_weights=nullptr, Tag *th_part_weights=nullptr, int verb=VERBOSE, const bool debug=false)
Set partition tag to each finite element in the problem.
static MoFEMErrorCodeGeneric< PetscErrorCode > ierr
multi_index_container< CubitMeshSets, indexed_by< hashed_unique< tag< Meshset_mi_tag >, member< CubitMeshSets, EntityHandle, &CubitMeshSets::meshset > >, ordered_non_unique< tag< CubitMeshsetType_mi_tag >, const_mem_fun< CubitMeshSets, unsigned long int, &CubitMeshSets::getBcTypeULong > >, ordered_non_unique< tag< CubitMeshsetMaskedType_mi_tag >, const_mem_fun< CubitMeshSets, unsigned long int, &CubitMeshSets::getMaksedBcTypeULong > >, ordered_non_unique< tag< CubitMeshsets_name >, const_mem_fun< CubitMeshSets, std::string, &CubitMeshSets::getName > >, hashed_unique< tag< Composite_Cubit_msId_And_MeshsetType_mi_tag >, composite_key< CubitMeshSets, const_mem_fun< CubitMeshSets, int, &CubitMeshSets::getMeshsetId >, const_mem_fun< CubitMeshSets, unsigned long int, &CubitMeshSets::getMaksedBcTypeULong > > > > > CubitMeshSet_multiIndex
Stores data about meshsets (see CubitMeshSets) storing data about boundary conditions,...
DeprecatedCoreInterface Interface
constexpr double t
plate stiffness
FTensor::Index< 'm', 3 > m
virtual moab::Interface & get_moab()=0
virtual int get_comm_rank() const =0
static MoFEMErrorCode Initialize(int *argc, char ***args, const char file[], const char help[])
Initializes the MoFEM database PETSc, MOAB and MPI.
static MoFEMErrorCode Finalize()
Checks for options to be called at the conclusion of the program.
Deprecated interface functions.
Interface for managing meshsets containing materials and boundary conditions.
CubitMeshSet_multiIndex & getMeshsetsMultindex()
MoFEMErrorCode getInterface(IFACE *&iface) const
Get interface refernce to pointer of interface.