1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
|
#ifndef INCLUDED_NJN_APPROX
#define INCLUDED_NJN_APPROX
/* $Id: $
* ===========================================================================
*
* PUBLIC DOMAIN NOTICE
* National Center for Biotechnology Information
*
* This software/database is a "United States Government Work" under the
* terms of the United States Copyright Act. It was written as part of
* the author's offical duties as a United States Government employee and
* thus cannot be copyrighted. This software/database is freely available
* to the public for use. The National Library of Medicine and the U.S.
* Government have not placed any restriction on its use or reproduction.
*
* Although all reasonable efforts have been taken to ensure the accuracy
* and reliability of the software and data, the NLM and the U.S.
* Government do not and cannot warrant the performance or results that
* may be obtained by using this software or data. The NLM and the U.S.
* Government disclaim all warranties, express or implied, including
* warranties of performance, merchantability or fitness for any particular
* purpose.
*
* Please cite the author in any work or product based on this material.
*
* ===========================================================================*/
/*****************************************************************************
File name: njn_approx.hpp
Author: John Spouge
Contents:
******************************************************************************/
#include <float.h>
#include <cassert>
#include <math.h>
namespace Njn {
namespace Approx{
// Approximation
const float FLT_THRESHOLD = 10.0F; // Rounding threshold
const float FLT_ROUND = FLT_THRESHOLD * FLT_EPSILON; // Rounding error threshold
const double DBL_THRESHOLD = 100.0; // Rounding threshold
const double DBL_ROUND = DBL_THRESHOLD * DBL_EPSILON; // Rounding error threshold
template <typename T> inline bool approx (T x_, T y_, T eps_);
template <typename T> inline bool relApprox (T x_, T y_, T eps_);
template <typename T> inline bool absRelApprox (T x_, T y_, T tol_, T rtol_);
// Rounding
template <typename T> inline bool eq (T x_, T y_, T round_);
template <typename T> inline bool ge (T x_, T y_, T round_);
template <typename T> inline bool gt (T x_, T y_, T round_);
template <typename T> inline bool ne (T x_, T y_, T round_);
template <typename T> inline bool lt (T x_, T y_, T round_);
template <typename T> inline bool le (T x_, T y_, T round_);
// Approximation
template <typename T> bool approx (T x_, T y_, T eps_) {return fabs (x_ - y_) <= fabs (eps_);}
template <typename T> bool relApprox (T x_, T y_, T eps_) {return approx <T> (x_, y_, eps_ * y_);}
template <typename T> bool absRelApprox (T x_, T y_, T tol_, T rtol_) {return approx <T> (x_, y_, tol_) || relApprox <T> (x_, y_, rtol_);}
// Rounding
template <typename T> bool eq (T x_, T y_, T round_) {return relApprox <T> (x_, y_, round_);}
template <typename T> bool ge (T x_, T y_, T round_) {return (x_ - y_) >= -(fabs (y_)) * round_;}
template <typename T> bool gt (T x_, T y_, T round_) {return (x_ - y_) > (fabs (y_)) * round_;}
template <typename T> bool ne (T x_, T y_, T round_) {return ! eq (x_, y_, round_);}
template <typename T> bool lt (T x_, T y_, T round_) {return ! ge (x_, y_, round_);}
template <typename T> bool le (T x_, T y_, T round_) {return ! gt (x_, y_, round_);}
}
}
#endif //! INCLUDED
|
