19#include <range/v3/algorithm/transform.hpp>
20#include <range/v3/range/conversion.hpp>
21#include <unordered_map>
48 OGS_FATAL(
"Unsupported MeshItemType {:s}.",
54 std::ostream& os, std::vector<std::size_t>
const& global_node_ids)
const
56 std::vector<MeshLib::IO::NodeData> nodes_buffer;
57 nodes_buffer.reserve(
nodes.size());
59 for (
const auto* node :
nodes)
61 double const* coords = node->data();
62 nodes_buffer.emplace_back(global_node_ids[node->getID()], coords[0],
63 coords[1], coords[2]);
65 return os.write(
reinterpret_cast<const char*
>(nodes_buffer.data()),
74 std::vector<const MeshLib::Element*>
const& elements)
76 return 3 * elements.size() +
77 std::accumulate(begin(elements), end(elements), 0,
78 [](
auto const nnodes,
auto const* e)
79 {
return nnodes + e->getNumberOfNodes(); });
85 static_cast<long>(
nodes.size()),
99 return os.write(
reinterpret_cast<const char*
>(data),
sizeof(data));
103 std::vector<std::size_t>
const& partition_ids,
104 std::vector<std::size_t>
const& node_id_mapping)
106 return partition_ids[node_id_mapping[node_id]];
109std::pair<std::vector<MeshLib::Node const*>, std::vector<MeshLib::Node const*>>
114 std::vector<MeshLib::Node const*> base_nodes;
117 base_nodes.reserve(nodes.size() / 2);
118 std::vector<MeshLib::Node const*> higher_order_nodes;
121 higher_order_nodes.reserve(nodes.size() / 2);
125 begin(nodes), end(nodes), std::back_inserter(base_nodes),
126 std::back_inserter(higher_order_nodes),
130 return {base_nodes, higher_order_nodes};
136std::tuple<std::vector<MeshLib::Node*>, std::vector<MeshLib::Node*>>
138 std::size_t
const part_id,
139 std::vector<MeshLib::Node*>
const& nodes,
140 std::vector<MeshLib::Element const*>
const& ghost_elements,
141 std::vector<std::size_t>
const& partition_ids,
143 std::vector<std::size_t>
const& node_id_mapping)
145 std::vector<MeshLib::Node*> base_ghost_nodes;
146 std::vector<MeshLib::Node*> higher_order_ghost_nodes;
148 std::vector<bool> is_ghost_node(nodes.size(),
false);
149 for (
const auto* ghost_elem : ghost_elements)
151 for (
unsigned i = 0; i < ghost_elem->getNumberOfNodes(); i++)
153 auto const& n = ghost_elem->getNode(i);
154 auto const node_id = n->getID();
155 if (is_ghost_node[node_id])
165 base_ghost_nodes.push_back(nodes[node_id]);
169 higher_order_ghost_nodes.push_back(nodes[node_id]);
171 is_ghost_node[node_id] =
true;
175 return std::tuple<std::vector<MeshLib::Node*>, std::vector<MeshLib::Node*>>{
176 base_ghost_nodes, higher_order_ghost_nodes};
184 std::size_t
const offset,
188 auto const& nodes = p.nodes;
189 auto const nnodes = nodes.size();
191 for (std::size_t i = 0; i < nnodes; ++i)
193 const auto global_id = nodes[i]->getID();
194 std::copy_n(&pv[n_components * global_id], n_components,
195 &partitioned_pv[offset + n_components * i]);
197 return n_components * nnodes;
206 std::size_t
const offset,
210 std::size_t
const n_regular(p.regular_elements.size());
212 for (std::size_t i = 0; i < n_regular; ++i)
214 const auto id = p.regular_elements[i]->getID();
215 std::copy_n(&pv[n_components *
id], n_components,
216 &partitioned_pv[offset + n_components * i]);
219 std::size_t
const n_ghost(p.ghost_elements.size());
220 for (std::size_t i = 0; i < n_ghost; ++i)
222 const auto id = p.ghost_elements[i]->getID();
223 std::copy_n(&pv[n_components *
id], n_components,
224 &partitioned_pv[offset + n_components * (n_regular + i)]);
226 return n_components * (n_regular + n_ghost);
237 std::size_t
const id_offset_partition,
238 std::vector<std::size_t>
const& element_ip_data_offsets,
247 std::copy_n(&pv[0], pv.
size(), &partitioned_pv[id_offset_partition]);
253 std::size_t id_offset = 0;
256 [&](std::vector<const MeshLib::Element*>
const& elements)
261 for (
auto const element : elements)
263 int const number_of_element_field_data =
268 auto const element_id = element->getID();
269 int const begin_pos = element_ip_data_offsets[element_id];
270 int const end_pos = element_ip_data_offsets[element_id + 1];
272 std::copy(pv.begin() + begin_pos, pv.begin() + end_pos,
273 &partitioned_pv[id_offset + id_offset_partition]);
274 id_offset += number_of_element_field_data;
278 copyFieldData(p.regular_elements);
279 copyFieldData(p.ghost_elements);
285 std::vector<Partition>& partitions)
287 for (
auto const& [name, property] : properties)
289 auto const item_type =
property->getMeshItemType();
298 if (property->getPropertyName().find(
"_ip") == std::string::npos)
303 std::string
const property_name =
property->getPropertyName();
304 auto countIntegrationPoints =
305 [&](std::vector<const MeshLib::Element*>
const& elements)
307 auto const ip_meta_data =
309 std::size_t counter = 0;
310 for (
auto const element : elements)
312 int const number_of_integration_points =
314 ip_meta_data, *element);
315 counter += number_of_integration_points;
320 for (
auto& p : partitions)
322 p.number_of_integration_points =
323 countIntegrationPoints(p.regular_elements) +
324 countIntegrationPoints(p.ghost_elements);
332 std::vector<MeshLib::Element*>
const& global_mesh_elements,
335 std::vector<Partition>
const& partitions,
337 std::map<MeshLib::MeshItemType, std::size_t>
const& total_number_of_tuples)
345 std::size_t partitioned_pv_size = total_number_of_tuples.at(item_type) *
348 std::vector<std::size_t> element_ip_data_offsets;
356 partitioned_pv_size = pv->
size() * partitions.size();
359 element_ip_data_offsets =
361 global_mesh_elements, *pv, properties);
367 if (partitioned_pv ==
nullptr)
370 "Could not create partitioned property vector {:s} for {} data "
374 partitioned_pv->resize(partitioned_pv_size);
376 auto copy_property_vector_values =
377 [&](
Partition const& p, std::size_t offset)
382 element_ip_data_offsets,
383 *pv, *partitioned_pv);
398 "Copying of property vector values for mesh item type {:s} is not "
400 toString(item_type));
403 std::size_t position_offset(0);
404 for (
auto p : partitions)
406 position_offset += copy_property_vector_values(p, position_offset);
412 std::vector<Partition>
const& partitions,
413 std::size_t
const total_number_of_cells)
415 auto* vtk_ghost_type =
418 if (vtk_ghost_type ==
nullptr)
420 OGS_FATAL(
"Could not create vtkGhostType cell data array.");
423 vtk_ghost_type->resize(total_number_of_cells);
424 std::size_t offset = 0;
425 for (
auto const& partition : partitions)
427 offset += partition.regular_elements.size();
428 for (std::size_t i = 0; i < partition.ghost_elements.size(); ++i)
430 if (partition.duplicate_ghost_cell[i])
432 (*vtk_ghost_type)[offset + i] |=
433 vtkDataSetAttributes::DUPLICATECELL;
436 offset += partition.ghost_elements.size();
442 std::unique_ptr<MeshLib::Mesh>
const& mesh,
443 std::vector<Partition>& partitions)
453 auto count_tuples = [&](
MeshItemType const mesh_item_type)
455 return std::accumulate(
456 begin(partitions), end(partitions), 0,
457 [&](std::size_t
const sum,
Partition const& p)
458 {
return sum + p.numberOfMeshItems(mesh_item_type); });
461 std::map<MeshItemType, std::size_t>
const total_number_of_tuples = {
462 {MeshItemType::Cell, count_tuples(MeshItemType::Cell)},
463 {MeshItemType::Node, count_tuples(MeshItemType::Node)},
464 {MeshItemType::IntegrationPoint,
465 count_tuples(MeshItemType::IntegrationPoint)}};
468 "total number of tuples after partitioning defined for cells is {:d} "
469 "and for nodes {:d} and for integration points {:d}.",
470 total_number_of_tuples.at(MeshItemType::Cell),
471 total_number_of_tuples.at(MeshItemType::Node),
472 total_number_of_tuples.at(MeshItemType::IntegrationPoint));
479 [&](
auto type,
auto const property)
482 mesh->getElements(), partitioned_properties, properties,
485 total_number_of_tuples);
490 total_number_of_tuples.at(MeshItemType::Cell));
492 return partitioned_properties;
496 std::vector<Partition>& partitions)
498 std::vector<bool> cell_visited(mesh.
getElements().size(),
false);
500 for (
auto& partition : partitions)
502 partition.duplicate_ghost_cell.resize(partition.ghost_elements.size(),
505 for (std::size_t i = 0; i < partition.ghost_elements.size(); i++)
507 const auto& ghost_element = *partition.ghost_elements[i];
508 if (!cell_visited[ghost_element.getID()])
510 cell_visited[ghost_element.getID()] =
true;
511 partition.duplicate_ghost_cell[i] =
false;
518 std::vector<MeshLib::Element*>
const& global_mesh_elements,
521 for (
auto const& [name, property] : properties)
523 auto const item_type =
property->getMeshItemType();
532 if (property->getPropertyName().find(
"_ip") == std::string::npos)
537 std::size_t number_of_total_integration_points = 0;
539 properties, property->getPropertyName());
540 for (
auto const element : global_mesh_elements)
542 int const number_of_integration_points =
545 number_of_total_integration_points += number_of_integration_points;
551 if (pv->size() != number_of_total_integration_points * component_number)
554 "The property vector's size {:d} for integration point data "
555 "{:s} does not match its actual size {:d}. The field data in "
556 "the vtu file are wrong.",
558 number_of_total_integration_points * component_number);
564 std::vector<std::size_t>
const& node_partition_map,
565 std::vector<MeshLib::Element*>
const& elements,
566 std::vector<std::size_t>
const& bulk_node_ids)
568 auto node_partition_ids = ranges::views::transform(
571 auto node_lookup = ranges::views::transform(
572 [&](std::size_t
const i)
573 {
return node_partition_map[bulk_node_ids[i]]; });
576 ranges::to<std::vector>;
579 return elements | node_partition_ids | ranges::to<std::vector>;
583 std::vector<Partition>& partitions,
584 std::vector<std::size_t>
const& nodes_partition_ids,
585 std::vector<MeshLib::Node*>
const& nodes,
586 std::vector<std::size_t>
const& bulk_node_ids)
588 for (
auto const*
const node : nodes)
590 partitions[nodes_partition_ids[bulk_node_ids[node->getID()]]]
591 .nodes.push_back(node);
598 std::vector<MeshLib::Node const*> higher_order_nodes;
601 std::tie(partition.
nodes, higher_order_nodes) =
604 std::copy(begin(higher_order_nodes), end(higher_order_nodes),
605 std::back_inserter(partition.
nodes));
610 std::vector<Partition>& partitions,
MeshLib::Mesh const& mesh)
612 for (
auto& partition : partitions)
623 for (
auto& partition : partitions)
625 partition.number_of_regular_nodes = partition.nodes.size();
626 partition.number_of_mesh_base_nodes = number_of_mesh_base_nodes;
627 partition.number_of_mesh_all_nodes = number_of_mesh_all_nodes;
632 std::vector<Partition>& partitions,
634 std::vector<std::vector<std::size_t>>
const& partition_ids_per_element)
638 auto const element_id = element->getID();
639 auto node_partition_ids = partition_ids_per_element[element_id];
641 std::sort(node_partition_ids.begin(), node_partition_ids.end());
643 std::unique(node_partition_ids.begin(), node_partition_ids.end());
644 node_partition_ids.erase(last, node_partition_ids.end());
647 if (node_partition_ids.size() == 1)
649 partitions[node_partition_ids[0]].regular_elements.push_back(
654 for (
auto const partition_id : node_partition_ids)
656 partitions[partition_id].ghost_elements.push_back(element);
666 std::vector<Partition>& partitions,
MeshLib::Mesh const& mesh,
667 std::vector<std::size_t>
const& nodes_partition_ids,
668 std::vector<std::size_t>
const& node_id_mapping)
670 for (std::size_t part_id = 0; part_id < partitions.size(); part_id++)
672 auto& partition = partitions[part_id];
673 std::vector<MeshLib::Node*> base_ghost_nodes;
674 std::vector<MeshLib::Node*> higher_order_ghost_nodes;
675 std::tie(base_ghost_nodes, higher_order_ghost_nodes) =
677 part_id, mesh.
getNodes(), partition.ghost_elements,
678 nodes_partition_ids, mesh, node_id_mapping);
680 std::copy(begin(base_ghost_nodes), end(base_ghost_nodes),
681 std::back_inserter(partition.nodes));
683 partition.number_of_base_nodes =
684 partition.number_of_regular_base_nodes + base_ghost_nodes.size();
686 std::copy(begin(higher_order_ghost_nodes),
687 end(higher_order_ghost_nodes),
688 std::back_inserter(partition.nodes));
694 std::vector<std::size_t>
const& nodes_partition_ids,
695 std::vector<std::size_t>
const& bulk_node_ids)
700 nodes_partition_ids, mesh.
getElements(), bulk_node_ids);
701 INFO(
"partitionMesh(): Partition IDs per element computed in {:g} s",
707 INFO(
"partitionMesh(): distribute nodes to partitions took {:g} s",
713 "partitionMesh(): sorting [base nodes | higher order nodes] took {:g} "
720 "partitionMesh(): setting number of nodes and of all mesh base nodes "
726 partition_ids_per_element);
727 INFO(
"partitionMesh(): distribute elements into partitions took {:g} s",
732 partitions, mesh, nodes_partition_ids, bulk_node_ids);
733 INFO(
"partitionMesh(): determine / append ghost nodes took {:g} s",
738 INFO(
"partitionMesh(): markDuplicateGhostCells took {:g} s",
744 std::vector<std::size_t> bulk_node_ids(
_mesh->getNumberOfNodes());
745 std::iota(bulk_node_ids.begin(), bulk_node_ids.end(), 0);
761 std::vector<Partition>
const& partitions,
764 auto const bulk_node_ids_string =
773 auto const bulk_element_ids_string =
776 static_cast<std::string
>(bulk_element_ids_string),
788 std::vector<Partition>
const& local_partitions)
const
790 if (bulk_node_ids_pv ==
nullptr)
795 auto& bulk_node_ids = *bulk_node_ids_pv;
797 std::size_t offset = 0;
799 assert(
_partitions.size() == local_partitions.size());
800 int const n_partitions =
static_cast<int>(
_partitions.size());
801 for (
int partition_id = 0; partition_id < n_partitions; ++partition_id)
803 auto const& bulk_partition =
_partitions[partition_id];
804 auto const& local_partition = local_partitions[partition_id];
807 auto const& bulk_nodes = bulk_partition.nodes;
808 auto const n_bulk_nodes = bulk_nodes.size();
809 std::map<std::size_t, std::size_t> global_to_local;
810 for (std::size_t local_node_id = 0; local_node_id < n_bulk_nodes;
813 global_to_local[bulk_nodes[local_node_id]->getID()] = local_node_id;
816 auto const& local_nodes = local_partition.nodes;
817 auto const n_local_nodes = local_nodes.size();
818 for (std::size_t local_node_id = 0; local_node_id < n_local_nodes;
821 bulk_node_ids[offset + local_node_id] =
822 global_to_local[bulk_node_ids[offset + local_node_id]];
824 offset += n_local_nodes;
830 std::vector<Partition>
const& local_partitions)
const
832 if (bulk_element_ids_pv ==
nullptr)
837 auto& bulk_element_ids = *bulk_element_ids_pv;
839 std::size_t offset = 0;
841 assert(
_partitions.size() == local_partitions.size());
842 int const n_partitions =
static_cast<int>(
_partitions.size());
843 for (
int partition_id = 0; partition_id < n_partitions; ++partition_id)
845 auto const& bulk_partition =
_partitions[partition_id];
846 auto const& local_partition = local_partitions[partition_id];
849 std::map<std::size_t, std::size_t> global_to_local;
852 std::vector<MeshLib::Element const*>
const& elements,
853 std::size_t
const offset)
855 auto const n_elements = elements.size();
856 for (std::size_t e = 0; e < n_elements; ++e)
858 global_to_local[elements[e]->getID()] = offset + e;
862 map_elements(bulk_partition.regular_elements, 0);
863 map_elements(bulk_partition.ghost_elements,
864 bulk_partition.regular_elements.size());
867 auto renumber_elements =
868 [&bulk_element_ids, &global_to_local](
869 std::vector<MeshLib::Element const*>
const& elements,
870 std::size_t
const offset)
872 auto const n_elements = elements.size();
873 for (std::size_t e = 0; e < n_elements; ++e)
875 bulk_element_ids[offset + e] =
876 global_to_local[bulk_element_ids[offset + e]];
881 offset += renumber_elements(local_partition.regular_elements, offset);
882 offset += renumber_elements(local_partition.ghost_elements, offset);
889 auto const bulk_node_ids_string =
891 auto const& bulk_node_ids =
895 std::vector<Partition> partitions(
_partitions.size());
904 std::size_t node_global_id_offset = 0;
908 for (std::size_t i = 0; i < partition.number_of_regular_nodes; i++)
911 node_global_id_offset++;
920 os.write(
reinterpret_cast<const char*
>(pv.data()), pv.
size() *
sizeof(T));
926 std::ostream& out_val, std::ostream& out_meta)
951 std::vector<Partition>
const& partitions,
954 auto const number_of_properties =
955 partitioned_properties.
size(mesh_item_type);
956 if (number_of_properties == 0)
961 auto const file_name_infix = toString(mesh_item_type);
963 auto const file_name_cfg = file_name_base +
"_partitioned_" +
964 file_name_infix +
"_properties_cfg" +
965 std::to_string(partitions.size()) +
".bin";
966 std::ofstream out(file_name_cfg, std::ios::binary);
969 OGS_FATAL(
"Could not open file '{:s}' for output.", file_name_cfg);
972 auto const file_name_val = file_name_base +
"_partitioned_" +
973 file_name_infix +
"_properties_val" +
974 std::to_string(partitions.size()) +
".bin";
975 std::ofstream out_val(file_name_val, std::ios::binary);
978 OGS_FATAL(
"Could not open file '{:s}' for output.", file_name_val);
984 partitioned_properties,
985 [&](
auto type,
auto const& property)
990 mesh_item_type, out_val, out);
993 unsigned long offset = 0;
994 for (
const auto& partition : partitions)
997 offset,
static_cast<unsigned long>(
998 partition.numberOfMeshItems(mesh_item_type))};
1000 "Write meta data for node-based PropertyVector: global offset "
1001 "{:d}, number of tuples {:d}",
1002 pvpmd.offset, pvpmd.number_of_tuples);
1004 offset += pvpmd.number_of_tuples;
1014 std::ostream&
writeConfig(std::ostream& os)
const;
1019 os.write(
reinterpret_cast<const char*
>(
this),
sizeof(
ConfigOffsets));
1021 static long reserved = 0;
1023 return os.write(
reinterpret_cast<const char*
>(&reserved),
sizeof(
long));
1035 return {
static_cast<long>(partition.
nodes.size()),
1067 const std::string& file_name_base, std::vector<Partition>
const& partitions)
1069 auto const file_name_cfg = file_name_base +
"_partitioned_msh_cfg" +
1070 std::to_string(partitions.size()) +
".bin";
1071 std::ofstream of_bin_cfg(file_name_cfg, std::ios::binary);
1074 OGS_FATAL(
"Could not open file '{:s}' for output.", file_name_cfg);
1077 std::vector<long> partitions_element_offsets;
1078 partitions_element_offsets.reserve(partitions.size());
1079 std::vector<long> partitions_ghost_element_offsets;
1080 partitions_ghost_element_offsets.reserve(partitions.size());
1083 for (
const auto& partition : partitions)
1092 partitions_element_offsets.push_back(
1093 partition_offsets.regular_elements);
1094 partitions_ghost_element_offsets.push_back(
1095 partition_offsets.ghost_elements);
1098 return std::make_tuple(partitions_element_offsets,
1099 partitions_ghost_element_offsets);
1111 const std::unordered_map<std::size_t, long>& local_node_ids,
1112 std::vector<long>& elem_info,
1115 constexpr unsigned mat_id =
1118 elem_info[counter++] = mat_id;
1119 elem_info[counter++] =
static_cast<long>(elem.
getCellType());
1120 elem_info[counter++] = nn;
1122 for (
long i = 0; i < nn; i++)
1124 auto const& n = *elem.
getNode(i);
1125 elem_info[counter++] = local_node_ids.at(n.getID());
1131 std::vector<MeshLib::Node const*>
const& nodes)
1133 std::unordered_map<std::size_t, long> local_ids;
1134 local_ids.reserve(nodes.size());
1136 long local_node_id = 0;
1137 for (
const auto* node : nodes)
1139 local_ids[node->getID()] = local_node_id++;
1151 std::vector<Partition>
const& partitions,
1152 std::vector<long>
const& regular_element_offsets,
1153 std::vector<long>
const& ghost_element_offsets)
1155 const std::string npartitions_str = std::to_string(partitions.size());
1157 auto const file_name_ele =
1158 file_name_base +
"_partitioned_msh_ele" + npartitions_str +
".bin";
1159 std::ofstream element_info_os(file_name_ele, std::ios::binary);
1160 if (!element_info_os)
1162 OGS_FATAL(
"Could not open file '{:s}' for output.", file_name_ele);
1165 auto const file_name_ele_g =
1166 file_name_base +
"_partitioned_msh_ele_g" + npartitions_str +
".bin";
1167 std::ofstream ghost_element_info_os(file_name_ele_g, std::ios::binary);
1168 if (!ghost_element_info_os)
1170 OGS_FATAL(
"Could not open file '{:s}' for output.", file_name_ele_g);
1173 for (std::size_t i = 0; i < partitions.size(); i++)
1175 const auto& partition = partitions[i];
1180 std::vector<long> ele_info(regular_element_offsets[i]);
1182 auto writeElementData =
1184 std::vector<MeshLib::Element const*>
const& elements,
1185 long const element_offsets,
1186 std::ofstream& output_stream)
1188 long counter = elements.size();
1189 std::vector<long> ele_info(element_offsets);
1191 for (std::size_t j = 0; j < elements.size(); j++)
1193 const auto* elem = elements[j];
1194 ele_info[j] = counter;
1199 output_stream.write(
reinterpret_cast<const char*
>(ele_info.data()),
1200 ele_info.size() *
sizeof(
long));
1204 writeElementData(partition.regular_elements, regular_element_offsets[i],
1207 writeElementData(partition.ghost_elements, ghost_element_offsets[i],
1208 ghost_element_info_os);
1217 std::vector<Partition>
const& partitions,
1218 std::vector<std::size_t>
const& global_node_ids)
1220 auto const file_name = file_name_base +
"_partitioned_msh_nod" +
1221 std::to_string(partitions.size()) +
".bin";
1222 std::ofstream os(file_name, std::ios::binary);
1225 OGS_FATAL(
"Could not open file '{:s}' for output.", file_name);
1228 for (
const auto& partition : partitions)
1230 partition.writeNodes(os, global_node_ids);
1245 const std::vector<IntegerType>& regular_element_offsets =
1246 std::get<0>(elements_offsets);
1247 const std::vector<IntegerType>& ghost_element_offsets =
1248 std::get<1>(elements_offsets);
1250 ghost_element_offsets);
1256 std::string
const& output_filename_base,
1257 std::vector<Partition>
const& partitions,
1262 const auto elem_integers =
1265 const std::vector<IntegerType>& num_elem_integers =
1266 std::get<0>(elem_integers);
1267 const std::vector<IntegerType>& num_g_elem_integers =
1268 std::get<1>(elem_integers);
1269 writeElements(output_filename_base, partitions, num_elem_integers,
1270 num_g_elem_integers);
1272 writeProperties(output_filename_base, partitioned_properties, partitions,
1274 writeProperties(output_filename_base, partitioned_properties, partitions,
void INFO(fmt::format_string< Args... > fmt, Args &&... args)
void DBUG(fmt::format_string< Args... > fmt, Args &&... args)
Definition of mesh-related Enumerations.
Declare a class to perform node wise mesh partitioning.
void applyToPropertyVectors(Properties const &properties, Function f)
Definition of the RunTime class.
Implementation of the VtuInterface class.
void renumberBulkElementIdsProperty(MeshLib::PropertyVector< std::size_t > *const bulk_element_ids_pv, std::vector< Partition > const &local_partitions) const
std::vector< Partition > partitionOtherMesh(MeshLib::Mesh const &mesh) const
MeshLib::Mesh const & mesh() const
void write(const std::string &file_name_base)
void renumberBulkIdsProperty(std::vector< Partition > const &partitions, MeshLib::Properties &partitioned_properties)
void renumberBulkNodeIdsProperty(MeshLib::PropertyVector< std::size_t > *const bulk_node_ids, std::vector< Partition > const &local_partitions) const
void renumberNodeIndices()
std::unique_ptr< MeshLib::Mesh > _mesh
Pointer to a mesh object.
MeshLib::Properties _partitioned_properties
Properties where values at ghost nodes and extra nodes are inserted.
std::vector< std::size_t > _nodes_global_ids
Global IDs of all nodes after partitioning.
void partitionByMETIS()
Partition by node.
std::vector< Partition > _partitions
Data for all partitions.
std::vector< std::size_t > _nodes_partition_ids
Partition IDs of each nodes.
void writeOtherMesh(std::string const &output_filename_base, std::vector< Partition > const &partitions, MeshLib::Properties const &partitioned_properties) const
double elapsed() const
Get the elapsed time in seconds.
void start()
Start the timer.
virtual CellType getCellType() const =0
virtual unsigned getNumberOfNodes() const =0
virtual const Node * getNode(unsigned idx) const =0
constexpr std::span< Node *const > nodes() const
Span of element's nodes, their pointers actually.
std::vector< Node * > const & getNodes() const
Get the nodes-vector for the mesh.
std::vector< Element * > const & getElements() const
Get the element-vector for the mesh.
std::size_t computeNumberOfBaseNodes() const
Get the number of base nodes.
Properties & getProperties()
std::size_t getNumberOfNodes() const
Get the number of nodes.
std::vector< Element const * > const & getElementsConnectedToNode(std::size_t node_id) const
Property manager on mesh items. Class Properties manages scalar, vector or matrix properties....
bool hasPropertyVector(std::string_view name) const
std::map< std::string, PropertyVectorBase * >::size_type size() const
PropertyVector< T > * createNewPropertyVector(std::string_view name, MeshItemType mesh_item_type, std::size_t n_components=1)
PropertyVector< T > const * getPropertyVector(std::string_view name) const
MeshItemType getMeshItemType() const
int getNumberOfGlobalComponents() const
std::string const & getPropertyName() const
std::size_t getNumberOfTuples() const
std::tuple< std::vector< MeshLib::Node * >, std::vector< MeshLib::Node * > > findGhostNodesInPartition(std::size_t const part_id, std::vector< MeshLib::Node * > const &nodes, std::vector< MeshLib::Element const * > const &ghost_elements, std::vector< std::size_t > const &partition_ids, MeshLib::Mesh const &mesh, std::vector< std::size_t > const &node_id_mapping)
std::size_t copyCellPropertyVectorValues(Partition const &p, std::size_t const offset, MeshLib::PropertyVector< T > const &pv, MeshLib::PropertyVector< T > &partitioned_pv)
std::unordered_map< std::size_t, long > enumerateLocalNodeIds(std::vector< MeshLib::Node const * > const &nodes)
Generates a mapping of given node ids to a new local (renumbered) node ids.
void addVtkGhostTypeProperty(MeshLib::Properties &partitioned_properties, std::vector< Partition > const &partitions, std::size_t const total_number_of_cells)
bool writePropertyVector(MeshLib::PropertyVector< T > const *const pv, MeshLib::MeshItemType const mesh_item_type, std::ostream &out_val, std::ostream &out_meta)
void checkFieldPropertyVectorSize(std::vector< MeshLib::Element * > const &global_mesh_elements, MeshLib::Properties const &properties)
void writeNodes(const std::string &file_name_base, std::vector< Partition > const &partitions, std::vector< std::size_t > const &global_node_ids)
ConfigOffsets incrementConfigOffsets(ConfigOffsets const &oldConfig, PartitionOffsets const &offsets)
void setIntegrationPointNumberOfPartition(MeshLib::Properties const &properties, std::vector< Partition > &partitions)
void distributeNodesToPartitions(std::vector< Partition > &partitions, std::vector< std::size_t > const &nodes_partition_ids, std::vector< MeshLib::Node * > const &nodes, std::vector< std::size_t > const &bulk_node_ids)
void reorderNodesIntoBaseAndHigherOrderNodes(Partition &partition, MeshLib::Mesh const &mesh)
void reorderNodesIntoBaseAndHigherOrderNodesPerPartition(std::vector< Partition > &partitions, MeshLib::Mesh const &mesh)
MeshLib::Properties partitionProperties(std::unique_ptr< MeshLib::Mesh > const &mesh, std::vector< Partition > &partitions)
Partition existing properties and add vtkGhostType cell data array property.
void setNumberOfNodesInPartitions(std::vector< Partition > &partitions, MeshLib::Mesh const &mesh)
std::pair< std::vector< MeshLib::Node const * >, std::vector< MeshLib::Node const * > > splitIntoBaseAndHigherOrderNodes(std::vector< MeshLib::Node const * > const &nodes, MeshLib::Mesh const &mesh)
PartitionOffsets computePartitionOffsets(Partition const &partition)
NodeWiseMeshPartitioner::IntegerType getNumberOfIntegerVariablesOfElements(std::vector< const MeshLib::Element * > const &elements)
void determineAndAppendGhostNodesToPartitions(std::vector< Partition > &partitions, MeshLib::Mesh const &mesh, std::vector< std::size_t > const &nodes_partition_ids, std::vector< std::size_t > const &node_id_mapping)
void partitionMesh(std::vector< Partition > &partitions, MeshLib::Mesh const &mesh, std::vector< std::size_t > const &nodes_partition_ids, std::vector< std::size_t > const &bulk_node_ids)
std::tuple< std::vector< long >, std::vector< long > > writeConfigData(const std::string &file_name_base, std::vector< Partition > const &partitions)
std::size_t copyFieldPropertyDataToPartitions(MeshLib::Properties const &properties, Partition const &p, std::size_t const id_offset_partition, std::vector< std::size_t > const &element_ip_data_offsets, MeshLib::PropertyVector< T > const &pv, MeshLib::PropertyVector< T > &partitioned_pv)
void writePropertyVectorValues(std::ostream &os, MeshLib::PropertyVector< T > const &pv)
void getElementIntegerVariables(const MeshLib::Element &elem, const std::unordered_map< std::size_t, long > &local_node_ids, std::vector< long > &elem_info, long &counter)
void markDuplicateGhostCells(MeshLib::Mesh const &mesh, std::vector< Partition > &partitions)
void distributeElementsIntoPartitions(std::vector< Partition > &partitions, MeshLib::Mesh const &mesh, std::vector< std::vector< std::size_t > > const &partition_ids_per_element)
std::size_t copyNodePropertyVectorValues(Partition const &p, std::size_t const offset, MeshLib::PropertyVector< T > const &pv, MeshLib::PropertyVector< T > &partitioned_pv)
std::size_t partitionLookup(std::size_t const &node_id, std::vector< std::size_t > const &partition_ids, std::vector< std::size_t > const &node_id_mapping)
std::vector< std::vector< std::size_t > > computePartitionIDPerElement(std::vector< std::size_t > const &node_partition_map, std::vector< MeshLib::Element * > const &elements, std::vector< std::size_t > const &bulk_node_ids)
void writeProperties(const std::string &file_name_base, MeshLib::Properties const &partitioned_properties, std::vector< Partition > const &partitions, MeshLib::MeshItemType const mesh_item_type)
void writeElements(std::string const &file_name_base, std::vector< Partition > const &partitions, std::vector< long > const ®ular_element_offsets, std::vector< long > const &ghost_element_offsets)
bool copyPropertyVector(std::vector< MeshLib::Element * > const &global_mesh_elements, MeshLib::Properties &partitioned_properties, MeshLib::Properties const &properties, std::vector< Partition > const &partitions, MeshLib::PropertyVector< T > const *const pv, std::map< MeshLib::MeshItemType, std::size_t > const &total_number_of_tuples)
void writeValueBinary(std::ostream &out, T const &val)
write value as binary into the given output stream
void writePropertyVectorPartitionMetaData(std::ostream &os, PropertyVectorPartitionMetaData const &pvpmd)
void writePropertyVectorMetaData(std::ostream &os, PropertyVectorMetaData const &pvmd)
constexpr ranges::views::view_closure ids
For an element of a range view return its id.
IntegrationPointMetaData getIntegrationPointMetaData(MeshLib::Properties const &properties, std::string const &name)
constexpr std::string_view getBulkIDString(MeshItemType mesh_item_type)
static constexpr char const * toString(const MeshItemType t)
Returns a char array for a specific MeshItemType.
long ghost_element_rank_offset
std::ostream & writeConfig(std::ostream &os) const
std::size_t number_of_mesh_base_nodes
std::size_t number_of_base_nodes
std::size_t number_of_regular_nodes
std::ostream & writeConfig(std::ostream &os) const
std::size_t number_of_mesh_all_nodes
std::size_t numberOfMeshItems(MeshLib::MeshItemType const item_type) const
std::size_t number_of_regular_base_nodes
std::vector< const MeshLib::Element * > regular_elements
Non ghost elements.
std::vector< const MeshLib::Element * > ghost_elements
std::ostream & writeNodes(std::ostream &os, std::vector< std::size_t > const &global_node_ids) const
std::vector< MeshLib::Node const * > nodes
nodes.
std::size_t number_of_integration_points
struct NodeData used for parallel reading and also partitioning