fvbaseboundarycontext.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:
/*
  Copyright (C) 2011-2013 by Andreas Lauser

  This file is part of the Open Porous Media project (OPM).

  OPM is free software: you can redistribute it and/or modify
  it under the terms of the GNU General Public License as published by
  the Free Software Foundation, either version 2 of the License, or
  (at your option) any later version.

  OPM is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  GNU General Public License for more details.

  You should have received a copy of the GNU General Public License
  along with OPM.  If not, see <http://www.gnu.org/licenses/>.
*/
#ifndef EWOMS_FV_BASE_BOUNDARY_CONTEXT_HH
#define EWOMS_FV_BASE_BOUNDARY_CONTEXT_HH

#include "fvbaseproperties.hh"

#include <dune/common/fvector.hh>

namespace Ewoms {

template<class TypeTag>
class FvBaseBoundaryContext
{
    typedef typename GET_PROP_TYPE(TypeTag, Scalar) Scalar;
    typedef typename GET_PROP_TYPE(TypeTag, Problem) Problem;
    typedef typename GET_PROP_TYPE(TypeTag, Model) Model;
    typedef typename GET_PROP_TYPE(TypeTag, Stencil) Stencil;
    typedef typename GET_PROP_TYPE(TypeTag, ElementContext) ElementContext;
    typedef typename GET_PROP_TYPE(TypeTag, IntensiveQuantities) IntensiveQuantities;
    typedef typename GET_PROP_TYPE(TypeTag, ExtensiveQuantities) ExtensiveQuantities;
    typedef typename GET_PROP_TYPE(TypeTag, GradientCalculator) GradientCalculator;

    typedef typename GET_PROP_TYPE(TypeTag, GridView) GridView;
    typedef typename GridView::template Codim<0>::Entity Element;
    typedef typename GridView::IntersectionIterator IntersectionIterator;
    typedef typename GridView::Intersection Intersection;

    enum { dim = GridView::dimension };
    enum { dimWorld = GridView::dimensionworld };

    typedef typename GridView::ctype CoordScalar;
    typedef Dune::FieldVector<CoordScalar, dim> GlobalPosition;
    typedef Dune::FieldVector<Scalar, dimWorld> Vector;

public:
    explicit FvBaseBoundaryContext(const ElementContext &elemCtx)
        : elemCtx_(elemCtx)
        , intersectionIt_(gridView().ibegin(element()))
    { }

    const Problem &problem() const
    { return elemCtx_.problem(); }

    const Model &model() const
    { return elemCtx_.model(); }

    const GridView &gridView() const
    { return elemCtx_.gridView(); }

    const Element &element() const
    { return elemCtx_.element(); }

    const ElementContext &elementContext() const
    { return elemCtx_; }

    const GradientCalculator &gradientCalculator() const
    { return elemCtx_.gradientCalculator(); }

    int numDof(int timeIdx) const
    { return elemCtx_.numDof(timeIdx); }

    int numPrimaryDof(int timeIdx) const
    { return elemCtx_.numPrimaryDof(timeIdx); }

    int numInteriorFaces(int timeIdx) const
    { return elemCtx_.numInteriorFaces(timeIdx); }

    int numBoundaryFaces(int timeIdx) const
    { return elemCtx_.stencil(timeIdx).numBoundaryFaces(); }

    const Stencil &stencil(int timeIdx) const
    { return elemCtx_.stencil(timeIdx); }

    Vector normal(int boundaryFaceIdx, int timeIdx) const
    {
        auto tmp = stencil(timeIdx).boundaryFace[boundaryFaceIdx].normal;
        tmp /= tmp.two_norm();
        return tmp;
    }

    Scalar boundarySegmentArea(int boundaryFaceIdx, int timeIdx) const
    { return elemCtx_.stencil(timeIdx).boundaryFace(boundaryFaceIdx).area(); }

    const GlobalPosition &pos(int boundaryFaceIdx, int timeIdx) const
    { return stencil(timeIdx).boundaryFace(boundaryFaceIdx).integrationPos(); }

    const GlobalPosition &cvCenter(int boundaryFaceIdx, int timeIdx) const
    {
        int scvIdx = stencil(timeIdx).boundaryFace(boundaryFaceIdx).interiorIndex();
        return stencil(timeIdx).subControlVolume(scvIdx).globalPos();
    }

    short interiorScvIndex(int boundaryFaceIdx, int timeIdx) const
    {
        return stencil(timeIdx).boundaryFace(boundaryFaceIdx).interiorIndex();
    }

    int globalSpaceIndex(int boundaryFaceIdx, int timeIdx) const
    { return elemCtx_.globalSpaceIndex(interiorScvIndex(boundaryFaceIdx, timeIdx), timeIdx); }

    const IntensiveQuantities &intensiveQuantities(int boundaryFaceIdx, int timeIdx) const
    {
        short interiorScvIdx = this->interiorScvIndex(boundaryFaceIdx, timeIdx);
        return elemCtx_.intensiveQuantities(interiorScvIdx, timeIdx);
    }

    const ExtensiveQuantities &extensiveQuantities(int boundaryFaceIdx, int timeIdx) const
    { return elemCtx_.boundaryExtensiveQuantities(boundaryFaceIdx, timeIdx); }

    const Intersection &intersection(int boundaryFaceIdx) const
    { return *intersectionIt_; }

    IntersectionIterator &intersectionIt()
    { return intersectionIt_; }

protected:
    const ElementContext &elemCtx_;
    IntersectionIterator intersectionIt_;
};

} // namespace Ewoms

#endif