discretization/facecentered/staggered/subcontrolvolume.hh

Go to the documentation of this file.

// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
// vi: set et ts=4 sw=4 sts=4:
//
// SPDX-FileCopyrightInfo: Copyright © DuMux Project contributors, see AUTHORS.md in root folder
// SPDX-License-Identifier: GPL-3.0-or-later
//
#ifndef DUMUX_DISCRETIZATION_FACECENTERED_STAGGERED_SUBCONTROLVOLUME_HH
#define DUMUX_DISCRETIZATION_FACECENTERED_STAGGERED_SUBCONTROLVOLUME_HH
 
#include <array>
#include <utility>
 
#include <dune/geometry/type.hh>
#include <dune/geometry/axisalignedcubegeometry.hh>
 
#include <dumux/common/indextraits.hh>
 
namespace Dumux {
 
template<class GridView>
struct FaceCenteredDefaultScvGeometryTraits
{
    using GridIndexType = typename IndexTraits<GridView>::GridIndex;
    using LocalIndexType = typename IndexTraits<GridView>::LocalIndex;
    using Scalar = typename GridView::ctype;
    using Element = typename GridView::template Codim<0>::Entity;
    using GlobalPosition = typename Element::Geometry::GlobalCoordinate;
 
    static constexpr int dim = GridView::Grid::dimension;
    static constexpr int dimWorld = GridView::Grid::dimensionworld;
    using CornerStorage = std::array<GlobalPosition, (1<<(dim))>;
    using Geometry = Dune::AxisAlignedCubeGeometry<Scalar, dim, dimWorld>;
};
 
template<class GridView, class T = FaceCenteredDefaultScvGeometryTraits<GridView>>
class FaceCenteredStaggeredSubControlVolume
{
    using Geometry = typename T::Geometry;
    using CornerStorage = typename T::CornerStorage;
    using Element = typename T::Element;
    using GlobalPosition = typename T::GlobalPosition;
    using Scalar = typename T::Scalar;
    using GridIndexType = typename T::GridIndexType;
    using SmallLocalIndexType = typename IndexTraits<GridView>::SmallLocalIndex;
 
    using ElementGeometry = typename Element::Geometry;
    using IntersectionGeometry = typename GridView::Intersection::Geometry;
 
 
public:
    using Traits = T;
 
    FaceCenteredStaggeredSubControlVolume() = default;
 
    FaceCenteredStaggeredSubControlVolume(const ElementGeometry& elementGeometry,
                                          const IntersectionGeometry& intersectionGeometry,
                                          const GridIndexType globalIndex,
                                          const SmallLocalIndexType indexInElement,
                                          const GridIndexType dofIdx,
                                          const SmallLocalIndexType dofAxis,
                                          const GridIndexType eIdx,
                                          const bool boundary)
    : center_(0.5*(intersectionGeometry.center() + elementGeometry.center()))
    , dofPosition_(intersectionGeometry.center())
    , volume_(elementGeometry.volume()*0.5)
    , globalIndex_(globalIndex)
    , indexInElement_(indexInElement)
    , dofIdx_(dofIdx)
    , dofAxis_(dofAxis)
    , eIdx_(eIdx)
    , boundary_(boundary)
    {
        directionSign_ = (indexInElement % 2) ? 1.0 : -1.0;
    }
 
    const GlobalPosition& center() const
    { return center_; }
 
    const GlobalPosition& dofPosition() const
    { return dofPosition_; }
 
    Scalar volume() const
    { return volume_; }
 
    GridIndexType dofIndex() const
    { return dofIdx_; }
 
    GridIndexType index() const
    { return globalIndex_; }
 
    GridIndexType elementIndex() const
    { return eIdx_; }
 
    SmallLocalIndexType indexInElement() const
    { return indexInElement_; }
 
    SmallLocalIndexType localDofIndex() const
    { return indexInElement_; }
 
    SmallLocalIndexType dofAxis() const
    { return dofAxis_; }
 
    std::int_least8_t directionSign() const
    { return directionSign_; }
 
    bool boundary() const
    { return boundary_; }
 
private:
    GlobalPosition center_;
    GlobalPosition dofPosition_;
    Scalar volume_;
    GridIndexType globalIndex_;
    SmallLocalIndexType indexInElement_;
    GridIndexType dofIdx_;
    SmallLocalIndexType dofAxis_;
    std::int_least8_t directionSign_;
    GridIndexType eIdx_;
    bool boundary_;
};
 
} // end namespace Dumux
 
#endif