upscale.hpp

Go to the documentation of this file.

#ifndef OPM_VERTEQ_UPSCALE_HPP_INCLUDED
#define OPM_VERTEQ_UPSCALE_HPP_INCLUDED

// Copyright (C) 2013 Uni Research AS
// This file is licensed under the GNU General Public License v3.0

#include <utility> // pair

#ifndef OPM_VERTEQ_VISIBILITY_HPP_INCLUDED
#include <opm/verteq/visibility.hpp>
#endif /* OPM_VERTEQ_VISIBILITY_HPP_INCLUDED */

#ifndef OPM_VERTEQ_TOPSURF_HPP_INCLUDED
#include <opm/verteq/topsurf.hpp>
#endif /* OPM_VERTEQ_TOPSURF_HPP_INCLUDED */

#ifndef OPM_VERTEQ_RUNLEN_HPP_INCLUDED
#include <opm/verteq/utility/runlen.hpp>
#endif /* OPM_VERTEQ_RUNLEN_HPP_INCLUDED */

namespace Opm {

struct Elevation : private std::pair <int, double> {
    Elevation (int aBlock, double aFraction)
      : std::pair <int, double> (aBlock, aFraction) {
    }

    int block () const { return this->first; }

    double fraction () const { return this->second; }

    bool operator< (const Elevation& rhs) const;
};

struct VertEqUpscaler {
    VertEqUpscaler (const TopSurf& topSurf)
      : ts (topSurf) {
    }

    void gather (int col, double* buf, const double* data, int stride, int offset) const;

    void wgt_dpt (int col, const double* val, rlw_double& res) const;

    double dpt_avg (int col, const double* val) const;

    double sum (int col, const double* val) const;

    int num_rows (int col) const;

    Elevation bottom (int col) const;

    double eval (int col, const rlw_double& dpt, const Elevation zeta) const;

    Elevation find (int col, const double* dpt, const double target) const;

protected:
    const TopSurf& ts;
};

double snapToRange (const double value,
                    const double lo,
                    const double hi);

} // namespace Opm

#endif // OPM_VERTEQ_UPSCALE_HPP_INCLUDED