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