dgindexset.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_DG_INDEX_SET_HH
#define DUNE_FUFEM_DG_INDEX_SET_HH
 
#warning This header is deprecated and will be removed after 2.11.
 
#include <vector>
 
#include <dune/geometry/referenceelements.hh>
#include <dune/grid/common/mcmgmapper.hh>
 
template <class GridView>
class
[[deprecated("This class is deprecated and will be removed after 2.11.")]]
DGIndexSet {
 
public:
 
    DGIndexSet(const GridView& gridView) : mapper_(gridView,Dune::mcmgLayout(Dune::Codim<1>()))
    {
        setup(gridView);
    }
 
    void setup(const GridView& gridView) {
 
        mapper_.update(gridView);
 
        const int dim = GridView::dimension;
 
        // Compute the number of vertices for each face
        faceOffsets_.resize(gridView.size(1) + 1);
 
        for (const auto& e:  elements(gridView)) {
 
            auto refElement = Dune::ReferenceElements<typename GridView::ctype, dim>::general(e.type());
 
            // Alternative without NeighborIterator
            // Get the number of vertices for each face
            for (int i=0; i<refElement.size(1); i++)
                faceOffsets_[mapper_.subIndex(e, i, 1) + 1] = refElement.size(i,1,dim);
 
        }
 
        // Accumulate
        faceOffsets_[0] = 0;
        for (size_t i=1; i<faceOffsets_.size(); i++)
            faceOffsets_[i] += faceOffsets_[i-1];
    }
 
    int size() const {
        return faceOffsets_[faceOffsets_.size()-1];
    }
 
    int operator()(const typename GridView::template Codim<0>::Entity& element, int subFace) const {
        return faceOffsets_[mapper_.subIndex(element, subFace, 1)];
    }
 
    std::vector<int> faceOffsets_;
 
    Dune::MultipleCodimMultipleGeomTypeMapper<GridView> mapper_;
};
 
#endif