33 template <
int dim,
int spacedim>
35 add_indices_recursively_for_first_child_policy(
42 if (cell->level() > 0 && cell->parent()->child(0) == cell)
43 add_indices_recursively_for_first_child_policy(cell->parent(),
50 template <
int dim,
int spacedim>
55 template <
int dim,
int spacedim>
70#ifndef DEAL_II_WITH_MPI
77 const unsigned int process_has_active_locally_owned_cells =
78 tria->n_locally_owned_active_cells() > 0;
79 const unsigned int n_processes_with_active_locally_owned_cells =
82 if (n_processes_with_active_locally_owned_cells ==
86 unsigned int offset = 0;
88 const int ierr = MPI_Exscan(&process_has_active_locally_owned_cells,
92 process_has_active_locally_owned_cells)>,
98 tria->global_active_cell_index_partitioner().lock());
108 template <
int dim,
int spacedim>
111 : n_coarse_cells(tria_fine.n_global_coarse_cells())
112 , n_global_levels(tria_fine.n_global_levels())
117 "FirstChildPolicy is only working for pure hex meshes at the moment."))
124 if (cell->is_locally_owned())
125 add_indices_recursively_for_first_child_policy(cell,
132 template <
int dim,
int spacedim>
145 if (cell->is_locally_owned())
148 std::vector<unsigned int> owning_ranks_of_coarse_cells(
152 process(is_level_partitions,
155 owning_ranks_of_coarse_cells,
160 std::pair<types::global_cell_index, types::global_cell_index>>,
161 std::vector<unsigned int>>
162 consensus_algorithm(process, communicator);
163 consensus_algorithm.run();
173 tria->global_active_cell_index_partitioner().lock());
177 partition[cell->global_active_cell_index()] =
185 template <
int dim,
int spacedim>
187 const unsigned int n_min_cells)
188 : n_min_cells(n_min_cells)
193 template <
int dim,
int spacedim>
205 tria->global_active_cell_index_partitioner().lock());
209 const unsigned int n_locally_owned_active_cells =
213 [](
const auto &cell) { return cell.is_locally_owned(); });
227 const unsigned int n_partitions =
228 std::max<unsigned int>(1,
229 std::min<types::global_cell_index>(
230 n_global_active_cells / n_min_cells,
233 const unsigned int min_cells = n_global_active_cells / n_partitions;
235 const auto convert = [&](
const unsigned int i) {
240 const unsigned int n_partitions_with_additional_cell =
241 n_global_active_cells - min_cells * n_partitions;
243 const unsigned int rank =
244 (i < (min_cells + 1) * n_partitions_with_additional_cell) ?
245 (i / (min_cells + 1)) :
246 ((i - n_partitions_with_additional_cell) / min_cells);
253 for (
const auto i : partition.locally_owned_elements())
254 partition[i] = convert(i);
261 template <
int dim,
int spacedim>
267 : weighting_function(weighting_function)
272 template <
int dim,
int spacedim>
277#ifndef DEAL_II_WITH_MPI
288 const auto partitioner =
289 tria->global_active_cell_index_partitioner().lock();
291 std::vector<unsigned int> weights(partitioner->local_size());
297 for (
const auto &cell :
299 weights[partitioner->global_to_local(cell->global_active_cell_index())] =
304 std::uint64_t process_local_weight = 0;
305 for (
const auto &weight : weights)
306 process_local_weight += weight;
309 std::uint64_t process_local_weight_offset = 0;
311 int ierr = MPI_Exscan(
312 &process_local_weight,
313 &process_local_weight_offset,
321 std::uint64_t total_weight =
322 process_local_weight_offset + process_local_weight;
325 MPI_Bcast(&total_weight,
335 for (std::uint64_t i = 0, weight = process_local_weight_offset;
336 i < partition.locally_owned_size();
337 weight += weights[i], ++i)
338 partition.local_element(i) =
339 static_cast<double>(weight * n_subdomains / total_weight);
351#include "repartitioning_policy_tools.inst"
size_type index_within_set(const size_type global_index) const
size_type n_elements() const
void set_size(const size_type size)
void add_index(const size_type index)
virtual types::global_cell_index n_global_active_cells() const
virtual MPI_Comm get_communicator() const
bool all_reference_cells_are_hyper_cube() const
cell_iterator end() const
active_cell_iterator begin_active(const unsigned int level=0) const
types::global_cell_index translate(const TriaIterator< Accessor > &cell) const
types::global_cell_index size() const
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
IteratorRange< active_cell_iterator > active_cell_iterators() const
static ::ExceptionBase & ExcNotImplemented()
static ::ExceptionBase & ExcNeedsMPI()
#define Assert(cond, exc)
#define AssertThrowMPI(error_code)
#define AssertIndexRange(index, range)
static ::ExceptionBase & ExcMessage(std::string arg1)
T min(const T &t, const MPI_Comm &mpi_communicator)
T sum(const T &t, const MPI_Comm &mpi_communicator)
const MPI_Datatype mpi_type_id_for_type
unsigned int n_mpi_processes(const MPI_Comm &mpi_communicator)
const ::Triangulation< dim, spacedim > & tria