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