boundarydata.hh

Go to the documentation of this file.

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
 
// 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_GRID_BOUNDARYDATA_HH
#define DUNE_FUFEM_GRID_BOUNDARYDATA_HH
 
 
#include <cstddef>
#include <limits>
#include <utility>
#include <vector>
 
#include <dune/common/exceptions.hh>
 
#include <dune/grid/common/gridfactory.hh>
#include <dune/grid/io/file/gmshreader.hh>
 
 
namespace Dune::Fufem {
 
 
 
template<class T=int>
class BoundaryData
{
public:
 
  BoundaryData(const std::vector<T>& data)
    : data_(data)
  {}
 
  BoundaryData(std::vector<T>&& data)
    : data_(std::move(data))
  {}
 
  template<class Grid>
  BoundaryData(const std::vector<T>& insertedBoundaryData, const Dune::GridFactory<Grid>& factory, const Grid& grid)
    : data_()
  {
    int defaultValue = std::numeric_limits<int>::max();
    auto coarseGridView = grid.leafGridView();
    data_.resize(grid.numBoundarySegments());
    for(const auto& element : elements(coarseGridView))
      for(const auto& intersection : intersections(coarseGridView, element))
        if (intersection.boundary())
        {
          if (factory.wasInserted(intersection))
            data_[intersection.boundarySegmentIndex()] = insertedBoundaryData[factory.insertionIndex(intersection)];
          else
            data_[intersection.boundarySegmentIndex()] = defaultValue;
        }
  }
 
  template<class Grid>
  BoundaryData(const Dune::GmshReader<Grid>& reader, const Dune::GridFactory<Grid>& factory, const Grid& grid)
    : data_()
  {
    if (not reader.hasBoundaryData())
      DUNE_THROW(Dune::InvalidStateException, "The GmshReader passed to BoundaryData has no stored boundary data.");
    (*this) = BoundaryData(reader.boundaryData(), factory, grid);
  }
 
  template<class Intersection>
  const T& operator()(const Intersection& intersection) const
  {
    return data_[intersection.boundarySegmentIndex()];
  }
 
  const T& operator[](std::size_t boundarySegmentIndex) const
  {
    return data_[boundarySegmentIndex];
  }
 
  std::size_t size() const
  {
    return data_.size();
  }
 
  const std::vector<T>& data() const
  {
    return data_;
  }
 
  const std::vector<T>& coarseData() const
  {
    return data_;
  }
 
private:
  std::vector<T> data_;
};
 
template<class Grid>
BoundaryData(const Dune::GmshReader<Grid>&, const Dune::GridFactory<Grid>&, const Grid&)
  -> BoundaryData<int>;
 
 
 
} // namespace Dune::Fufem
 
 
#endif // DUNE_FUFEM_GRID_BOUNDARYDATA_HH