insertionindexpermutation.hh

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright © DUNE-FUFEM Project contributors, see file AUTHORS.md
// SPDX-License-Identifier: LicenseRef-GPL-2.0-only-with-DUNE-exception OR LGPL-3.0-or-later
 
#ifndef DUNE_FUFEM_INSERTIONINDEXPERMUTATION_HH
#define DUNE_FUFEM_INSERTIONINDEXPERMUTATION_HH
 
#include <cstddef>
#include <vector>
#include <algorithm>
#include <numeric>
#include <limits>
 
#include <dune/common/rangeutilities.hh>
 
#include <dune/fufem/domains/boundary.hh>
 
 
 
namespace Dune::Fufem {
 
 
 
template<class Permutation, class DataContainer>
void permuteInplace(const Permutation& permutation, DataContainer& data)
{
  auto processed = std::vector<bool>(data.size(), false);
  for(auto k : Dune::IntegralRange(data.size()))
  {
    if (not processed[k])
    {
      auto pos = k;
      while (std::size_t(permutation[pos]) != k)
      {
        std::swap(data[pos], data[permutation[pos]]);
        processed[pos] = true;
        pos = permutation[pos];
      }
      processed[pos] = true;
    }
  }
}
 
 
 
template<class GridFactory, class Grid>
std::vector<std::size_t> boundarySegmentInsertionIndices(const GridFactory& gridFactory, const Grid& grid)
{
  auto gridView = grid.leafGridView();
  auto defaultValue = std::numeric_limits<std::size_t>::max();
  auto indices = std::vector<std::size_t>(grid.numBoundarySegments(), defaultValue);
 
  for(auto&& intersection : Dune::Fufem::Boundary(gridView))
    if (gridFactory.wasInserted(intersection))
      indices[intersection.boundarySegmentIndex()] = gridFactory.insertionIndex(intersection);
    else
      indices[intersection.boundarySegmentIndex()] = defaultValue;
 
  return indices;
}
 
 
 
template<class GridFactory, class Grid, class Container>
void permuteFromBoundarySegmentInsertionIndices(const GridFactory& gridFactory, const Grid& grid, Container& data)
{
  auto permutation = boundarySegmentInsertionIndices(gridFactory, grid);
  Dune::Fufem::permuteInplace(permutation, data);
}
 
 
 
}
 
#endif // DUNE_FUFEM_INSERTIONINDEXPERMUTATION_HH