LinearIntervalInterpolation.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include "BaseLib/Error.h"
 
namespace MathLib {

template <typename NUMERIC_TYPE>
class LinearIntervalInterpolation {
public:
    LinearIntervalInterpolation(NUMERIC_TYPE a, NUMERIC_TYPE b, NUMERIC_TYPE c, NUMERIC_TYPE d);
    inline NUMERIC_TYPE operator() (NUMERIC_TYPE x) const;
 
private:
    NUMERIC_TYPE m_;
    NUMERIC_TYPE n_;
};
 
template <typename NUMERIC_TYPE>
LinearIntervalInterpolation<NUMERIC_TYPE>::LinearIntervalInterpolation(
    NUMERIC_TYPE a, NUMERIC_TYPE b, NUMERIC_TYPE c, NUMERIC_TYPE d)
    : m_(d - c), n_(0.0)
{
    if (b == a) {
        OGS_FATAL("LinearIntervalInterpolation::LinearIntervalInterpolation: a == b, empty interval");
    }
    m_ /= (b - a);
    n_ = c - m_ * a;
}
 
template <typename NUMERIC_TYPE>
inline NUMERIC_TYPE LinearIntervalInterpolation<NUMERIC_TYPE>::operator() (NUMERIC_TYPE x) const
{
    return m_ * x + n_;
}
 
} // end namespace MathLib