DGBasis.hpp

Go to the documentation of this file.

/*
  Copyright 2013 SINTEF ICT, Applied Mathematics.

  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 3 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 OPM_DGBASIS_HEADER_INCLUDED
#define OPM_DGBASIS_HEADER_INCLUDED

#include <vector>

struct UnstructuredGrid;

namespace Opm
{

    class DGBasisInterface
    {
    public:
        virtual ~DGBasisInterface();

        virtual int numBasisFunc() const = 0;

        virtual int dimensions() const = 0;

        virtual int degree() const = 0;

        virtual void eval(const int cell,
                          const double* x,
                          double* f_x) const = 0;

        virtual void evalGrad(const int cell,
                              const double* x,
                              double* grad_f_x) const = 0;

        virtual void addConstant(const double increment,
                                 double* coefficients) const = 0;

        virtual void multiplyGradient(const double factor,
                                      double* coefficients) const = 0;

        double evalFunc(const int cell,
                        const double* coefficients,
                        const double* x) const;

        virtual double functionAverage(const double* coefficients) const = 0;

    private:
        mutable std::vector<double> bvals_; // For evalFunc().
    };





    class DGBasisBoundedTotalDegree : public DGBasisInterface
    {
    public:
        DGBasisBoundedTotalDegree(const UnstructuredGrid& grid, const int degree);

        virtual ~DGBasisBoundedTotalDegree();

        virtual int numBasisFunc() const;

        virtual int dimensions() const;

        virtual int degree() const;

        virtual void eval(const int cell,
                          const double* x,
                          double* f_x) const;

        virtual void evalGrad(const int cell,
                              const double* x,
                              double* grad_f_x) const;

        virtual void addConstant(const double increment,
                                 double* coefficients) const;

        virtual void multiplyGradient(const double factor,
                                      double* coefficients) const;

        virtual double functionAverage(const double* coefficients) const;

    private:
        const UnstructuredGrid& grid_;
        const int degree_;
    };




    class DGBasisMultilin : public DGBasisInterface
    {
    public:
        DGBasisMultilin(const UnstructuredGrid& grid, const int degree);

        virtual ~DGBasisMultilin();

        virtual int numBasisFunc() const;

        virtual int dimensions() const;

        virtual int degree() const;

        virtual void eval(const int cell,
                          const double* x,
                          double* f_x) const;

        virtual void evalGrad(const int cell,
                              const double* x,
                              double* grad_f_x) const;

        virtual void addConstant(const double increment,
                                 double* coefficients) const;

        virtual void multiplyGradient(const double factor,
                                      double* coefficients) const;

        virtual double functionAverage(const double* coefficients) const;

    private:
        const UnstructuredGrid& grid_;
        const int degree_;

    };




} // namespace Opm


#endif // OPM_DGBASIS_HEADER_INCLUDED