From e89c947f764b579caa0f73a07c8c8060bc8d0bf2 Mon Sep 17 00:00:00 2001 From: Timo Heister Date: Fri, 17 May 2019 14:03:37 -0600 Subject: [PATCH] clarify a comment found while trying to understand the code --- source/numerics/data_out.cc | 15 ++++++++------- 1 file changed, 8 insertions(+), 7 deletions(-) diff --git a/source/numerics/data_out.cc b/source/numerics/data_out.cc index 4b83db88cc..4b3857f942 100644 --- a/source/numerics/data_out.cc +++ b/source/numerics/data_out.cc @@ -457,17 +457,18 @@ DataOut::build_patches( // that maps the cell indices to the patch numbers, as this will be needed // for generation of neighborship information. // Note, there is a confusing mess of different indices here at play: - // patch_index - the index of a patch in all_cells - // cell->index - only unique on each level, used in cell_to_patch_index_map - // active_index - index for a cell when counting from begin_active() using - // ++cell cell_index - unique index of a cell counted using - // next_locally_owned_cell() - // starting from first_locally_owned_cell() + // - patch_index: the index of a patch in all_cells + // - cell->index: only unique on each level, used in cell_to_patch_index_map + // - active_index: index for a cell when counting from begin_active() using + // ++cell (identical to cell->active_cell_index()) + // - cell_index: unique index of a cell counted using + // next_locally_owned_cell() starting from first_locally_owned_cell() // // It turns out that we create one patch for each selected cell, so // patch_index==cell_index. // - // will be cell_to_patch_index_map[cell->level][cell->index] = patch_index + // Now construct the map such that + // cell_to_patch_index_map[cell->level][cell->index] = patch_index std::vector> cell_to_patch_index_map; cell_to_patch_index_map.resize(this->triangulation->n_levels()); for (unsigned int l = 0; l < this->triangulation->n_levels(); ++l) -- 2.39.5