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
|
#include "gtest/gtest.h"
#include "bessel.hpp"
#include "mathutil.hpp"
#include <iostream>
using namespace libecpint;
class BesselTest : public testing::Test {
protected:
BesselFunction bessie;
void SetUp() override {
initFactorials();
bessie.init(2, 1600, 200, 1e-15);
}
};
TEST_F(BesselTest, CalculateBigL) {
std::vector<double> values(4, 0.0);
// check it catches maxL > lMax
bessie.calculate(SMALL, 3, values);
}
TEST_F(BesselTest, CalculateSmallZ) {
// check negative
std::vector<double> values(3, 0.0);
bessie.calculate(-1.0, 2, values);
EXPECT_DOUBLE_EQ(values[0], 1.0);
values.assign(3, 0.0);
bessie.calculate(SMALL/2, 2, values);
EXPECT_NEAR(values[0], 1.0, 1e-7);
EXPECT_NEAR(values[1], 1.66667e-8, 1e-12);
double v0 = bessie.calculate(SMALL/2, 0);
double v1 = bessie.calculate(SMALL/2, 1);
double v2 = bessie.calculate(SMALL/2, 2);
EXPECT_NEAR(values[0], v0, 1e-12);
EXPECT_NEAR(values[1], v1, 1e-12);
EXPECT_NEAR(values[2], v2, 1e-12);
}
TEST_F(BesselTest, CalculateBigZ) {
std::vector<double> values(2, 0.0);
bessie.calculate(17.0, 1, values);
EXPECT_NEAR(values[0], 0.0294118, 1e-7);
EXPECT_NEAR(values[1], 0.0276817, 1e-7);
double v0 = bessie.calculate(17.0, 0);
double v1 = bessie.calculate(17.0, 1);
EXPECT_NEAR(values[0], v0, 1e-12);
EXPECT_NEAR(values[1], v1, 1e-12);
}
TEST_F(BesselTest, CalculateMidZ) {
std::vector<double> values(2, 0.0);
bessie.calculate(5.0, 1, values);
EXPECT_NEAR(values[0], 0.0999955, 1e-7);
EXPECT_NEAR(values[1], 0.0800054, 1e-7);
double v0 = bessie.calculate(5.0, 0);
double v1 = bessie.calculate(5.0, 1);
EXPECT_NEAR(values[0], v0, 1e-12);
EXPECT_NEAR(values[1], v1, 1e-12);
}
TEST_F(BesselTest, UpperBoundSmallZ) {
std::vector<double> values(3, 0.0);
bessie.calculate(SMALL/2, 2, values);
double ub0 = bessie.upper_bound(SMALL/2, 0);
double ub1 = bessie.upper_bound(SMALL/2, 1);
double ub2 = bessie.upper_bound(SMALL/2, 2);
EXPECT_TRUE(ub0 >= values[0]);
EXPECT_TRUE(ub1 >= values[1]);
EXPECT_TRUE(ub2 >= values[2]);
}
TEST_F(BesselTest, UpperBoundBigZ) {
std::vector<double> values(3, 0.0);
bessie.calculate(17.0, 2, values);
double ub0 = bessie.upper_bound(17.0, 0);
double ub1 = bessie.upper_bound(17.0, 1);
double ub2 = bessie.upper_bound(17.0, 2);
EXPECT_TRUE(ub0 >= values[0]);
EXPECT_TRUE(ub1 >= values[1]);
EXPECT_TRUE(ub2 >= values[2]);
}
TEST_F(BesselTest, UpperBoundMidZ) {
std::vector<double> values(3, 0.0);
bessie.calculate(5.054, 2, values);
double ub0 = bessie.upper_bound(5.054, 0);
double ub1 = bessie.upper_bound(5.054, 1);
double ub2 = bessie.upper_bound(5.054, 2);
EXPECT_TRUE(ub0 >= values[0]);
EXPECT_TRUE(ub1 >= values[1]);
EXPECT_TRUE(ub2 >= values[2]);
}
int main(int argc, char** argv)
{
::testing::InitGoogleTest(&argc, argv);
return RUN_ALL_TESTS();
}
|
