Repository navigation
Expand file tree
/
Copy pathnumeric_test.cpp
More file actions
115 lines (91 loc) · 3.03 KB
/
Copy pathnumeric_test.cpp
File metadata and controls
115 lines (91 loc) · 3.03 KB
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
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
// numeric_test.cpp
/*===========================================================================
*
* 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 official 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.
*
* ===========================================================================
*
* Author: Vyacheslav Brover
*
* File Description:
* Test of numeric.{hpp,cpp}
*
*/
#undef NDEBUG
#include "common.hpp"
#include "numeric.hpp"
using namespace Common_sp;
#include "version.inc"
#include "common.inc"
namespace
{
struct ThisApplication final : Application
{
ThisApplication ()
: Application ("Test numeric.cpp")
{
version = VERSION;
addPositional ("go", "Go");
}
void body () const final
{
#if 0
cout << "sizeof(float) = " << sizeof (float) << endl;
cout << "sizeof(double) = " << sizeof (double) << endl;
#endif
ASSERT_EQ (zeta (1.5), 2.612, 1e-3);
ASSERT_EQ (zeta (1.5, 1), 2.612, 1e-3);
ASSERT_EQ (zeta (2), 1.645, 1e-3);
ASSERT_EQ (zeta (3), 1.202, 1e-3);
ASSERT_EQ (zeta (4), 1.0823, 1e-4);
ASSERT_EQ (zeta (2, 2), 0.645, 1e-3);
{
Sum sum;
const uint n = 1000000;
FOR (uint, i, n)
sum. add (1.0 / n);
ASSERT_EQ (sum. get (), 1, 1e-6);
}
{
SumLn sum;
const uint n = 1000000;
FOR (uint, i, n)
sum. addExp (- log (n));
ASSERT_EQ (exp (sum. getLn ()), 1, 1e-6);
}
#if 0
#define REAL2STR(x) cout << #x << ": " << real2str (x,6) << endl;
REAL2STR (0);
REAL2STR (1);
REAL2STR (0.0001);
REAL2STR (1e5);
REAL2STR (0.000123456);
#undef REAL2STR
#endif
}
};
} // namespace
int main (int argc,
const char* argv[])
{
ThisApplication app;
return app. run (argc, argv);
}