merger.hh

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
// SPDX-FileCopyrightInfo: Copyright © DUNE Project contributors, see file LICENSE.md in module root
// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-GPL-2.0-only-with-dune-grid-glue-exception
#ifndef DUNE_GRIDGLUE_MERGING_MERGER_HH
#define DUNE_GRIDGLUE_MERGING_MERGER_HH
 
#include <vector>
 
#include <dune/common/fvector.hh>
#include <dune/geometry/type.hh>
 
#include <dune/grid-glue/merging/intersectionlist.hh>
 
namespace Dune {
namespace GridGlue {
 
template <class ctype, int grid1Dim, int grid2Dim, int dimworld>
class Merger
{
public:
 
  typedef Dune::FieldVector<ctype, grid1Dim>  Grid1Coords;
 
  typedef Dune::FieldVector<ctype, grid2Dim>  Grid2Coords;
 
  typedef Dune::FieldVector<ctype, dimworld>  WorldCoords;
 
  using IntersectionList = Dune::GridGlue::IntersectionList<Grid1Coords, Grid2Coords>;
 
  virtual void build(const std::vector<Dune::FieldVector<ctype,dimworld> >& grid1_coords,
                     const std::vector<unsigned int>& grid1_elements,
                     const std::vector<Dune::GeometryType>& grid1_element_types,
                     const std::vector<Dune::FieldVector<ctype,dimworld> >& grid2_coords,
                     const std::vector<unsigned int>& grid2_elements,
                     const std::vector<Dune::GeometryType>& grid2_element_types) = 0;
 
  unsigned int nSimplices() const
    { return intersectionList()->size(); }
 
  virtual void clear() = 0;
 
  virtual std::shared_ptr<IntersectionList> intersectionList() const = 0;
 
  template<int n>
  unsigned int parents(unsigned int idx) const {
    return intersectionList()->template parents<n>(idx);
  }
 
  template<int n>
  unsigned int parent(unsigned int idx, unsigned int parId = 0) const
  {
    return intersectionList()->template parent<n>(idx, parId);
  }
 
  template<int n>
  auto parentLocal(unsigned int idx, unsigned int corner, unsigned int parId = 0) const
  {
    return intersectionList()->template corner<n>(idx, corner, parId);
  }
 
  unsigned int counter;
};
 
} /* namespace GridGlue */
} /* namespace Dune */
 
#endif