Loading...
Searching...
No Matches
UniformTableLinear.hpp
/*
Copyright 2010 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_UNIFORMTABLELINEAR_HEADER_INCLUDED
#define OPM_UNIFORMTABLELINEAR_HEADER_INCLUDED
#include <cmath>
#include <exception>
#include <vector>
#include <utility>
#include <ostream>
#include <opm/common/ErrorMacros.hpp>
namespace Opm {
template<typename T>
class UniformTableLinear
{
public:
UniformTableLinear();
UniformTableLinear(double xmin,
double xmax,
const std::vector<T>& y_values);
UniformTableLinear(double xmin,
double xmax,
const T* y_values,
int num_y_values);
std::pair<double, double> domain();
void rescaleDomain(std::pair<double, double> new_domain);
double operator()(const double x) const;
double derivative(const double x) const;
bool operator==(const UniformTableLinear& other) const;
enum RangePolicy {Throw = 0, ClosestValue = 1, Extrapolate = 2};
void setLeftPolicy(RangePolicy rp);
void setRightPolicy(RangePolicy rp);
protected:
double xmin_{};
double xmax_{};
double xdelta_{};
std::vector<T> y_values_{};
RangePolicy left_{};
RangePolicy right_{};
template <typename U>
friend std::ostream& operator<<(std::ostream& os, const UniformTableLinear<U>& t);
};
// Member implementations.
template<typename T>
inline
UniformTableLinear<T>
::UniformTableLinear()
: left_(ClosestValue), right_(ClosestValue)
{
}
template<typename T>
inline
UniformTableLinear<T>
::UniformTableLinear(double xmin,
double xmax,
const std::vector<T>& y_values)
: xmin_(xmin), xmax_(xmax), y_values_(y_values),
left_(ClosestValue), right_(ClosestValue)
{
assert(xmax > xmin);
assert(y_values.size() > 1);
xdelta_ = (xmax - xmin)/(y_values.size() - 1);
}
template<typename T>
inline
UniformTableLinear<T>
::UniformTableLinear(double xmin,
double xmax,
const T* y_values,
int num_y_values)
: xmin_(xmin), xmax_(xmax),
y_values_(y_values, y_values + num_y_values),
left_(ClosestValue), right_(ClosestValue)
{
assert(xmax > xmin);
assert(y_values_.size() > 1);
xdelta_ = (xmax - xmin)/(y_values_.size() - 1);
}
template<typename T>
inline std::pair<double, double>
UniformTableLinear<T>
::domain()
{
return std::make_pair(xmin_, xmax_);
}
template<typename T>
inline void
UniformTableLinear<T>
::rescaleDomain(std::pair<double, double> new_domain)
{
xmin_ = new_domain.first;
xmax_ = new_domain.second;
xdelta_ = (xmax_ - xmin_)/(y_values_.size() - 1);
}
template<typename T>
inline double
UniformTableLinear<T>
::operator()(const double xparam) const
{
// Implements ClosestValue policy.
double x = std::min(xparam, xmax_);
x = std::max(x, xmin_);
// Lookup is easy since we are uniform in x.
double pos = (x - xmin_)/xdelta_;
double posi = std::floor(pos);
int left = int(posi);
if (left == int(y_values_.size()) - 1) {
// We are at xmax_
return y_values_.back();
}
double w = pos - posi;
return (1.0 - w)*y_values_[left] + w*y_values_[left + 1];
}
template<typename T>
inline double
UniformTableLinear<T>
::derivative(const double xparam) const
{
// Implements derivative consistent
// with ClosestValue policy for function
double value;
if (xparam > xmax_ || xparam < xmin_) {
value = 0.0;
} else {
double x = std::min(xparam, xmax_);
x = std::max(x, xmin_);
// Lookup is easy since we are uniform in x.
double pos = (x - xmin_)/xdelta_;
double posi = std::floor(pos);
int left = int(posi);
if (left == int(y_values_.size()) - 1) {
// We are at xmax_
--left;
}
value = (y_values_[left + 1] - y_values_[left])/xdelta_;
}
return value;
}
template<typename T>
inline bool
UniformTableLinear<T>
::operator==(const UniformTableLinear<T>& other) const
{
return xmin_ == other.xmin_
&& xdelta_ == other.xdelta_
&& y_values_ == other.y_values_
&& left_ == other.left_
&& right_ == other.right_;
}
template<typename T>
inline void
UniformTableLinear<T>
::setLeftPolicy(RangePolicy rp)
{
if (rp != ClosestValue) {
OPM_THROW(std::runtime_error, "Only ClosestValue RangePolicy implemented.");
}
left_ = rp;
}
template<typename T>
inline void
UniformTableLinear<T>
::setRightPolicy(RangePolicy rp)
{
if (rp != ClosestValue) {
OPM_THROW(std::runtime_error, "Only ClosestValue RangePolicy implemented.");
}
right_ = rp;
}
template <typename T>
inline std::ostream& operator<<(std::ostream& os, const UniformTableLinear<T>& t)
{
int n = t.y_values_.size();
for (int i = 0; i < n; ++i) {
double f = double(i)/double(n - 1);
os << (1.0 - f)*t.xmin_ + f*t.xmax_
<< " " << t.y_values_[i] << '\n';
}
return os;
}
namespace utils
{
using Opm::UniformTableLinear;
}
} // namespace Opm
#endif // OPM_UNIFORMTABLELINEAR_HEADER_INCLUDED
Opm::UniformTableLinear
This class uses linear interpolation to compute the value (and its derivative) of a function f sample...
Definition UniformTableLinear.hpp:38
Opm::UniformTableLinear::rescaleDomain
void rescaleDomain(std::pair< double, double > new_domain)
Rescale the domain.
Definition UniformTableLinear.hpp:158
Opm::UniformTableLinear::operator()
double operator()(const double x) const
Evaluate the value at x.
Definition UniformTableLinear.hpp:168
Opm::UniformTableLinear::UniformTableLinear
UniformTableLinear()
Default constructor.
Definition UniformTableLinear.hpp:112
Opm::UniformTableLinear::derivative
double derivative(const double x) const
Evaluate the derivative at x.
Definition UniformTableLinear.hpp:189
Opm::UniformTableLinear::setRightPolicy
void setRightPolicy(RangePolicy rp)
Sets the behavioural policy for evaluation to the right of the domain.
Definition UniformTableLinear.hpp:239
Opm::UniformTableLinear::RangePolicy
RangePolicy
Policies for how to behave when trying to evaluate outside the domain.
Definition UniformTableLinear.hpp:85
Opm::UniformTableLinear::domain
std::pair< double, double > domain()
Get the domain.
Definition UniformTableLinear.hpp:150
Opm::UniformTableLinear::setLeftPolicy
void setLeftPolicy(RangePolicy rp)
Sets the behavioural policy for evaluation to the left of the domain.
Definition UniformTableLinear.hpp:228
Opm::UniformTableLinear::operator==
bool operator==(const UniformTableLinear &other) const
Equality operator.
Definition UniformTableLinear.hpp:216
Opm
This class implements a small container which holds the transmissibility mulitpliers for all the face...
Definition Exceptions.hpp:30