forked from npshub/mantid
-
Notifications
You must be signed in to change notification settings - Fork 0
/
XValidationTest.h
112 lines (95 loc) · 3.72 KB
/
XValidationTest.h
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
// Mantid Repository : https://github.com/mantidproject/mantid
//
// Copyright © 2018 ISIS Rutherford Appleton Laboratory UKRI,
// NScD Oak Ridge National Laboratory, European Spallation Source,
// Institut Laue - Langevin & CSNS, Institute of High Energy Physics, CAS
// SPDX - License - Identifier: GPL - 3.0 +
#pragma once
#include <cxxtest/TestSuite.h>
#include "MantidHistogramData/BinEdges.h"
#include "MantidHistogramData/Points.h"
#include "MantidHistogramData/XValidation.h"
#include <numeric>
using Mantid::HistogramData::BinEdges;
using Mantid::HistogramData::HistogramX;
using Mantid::HistogramData::isValid;
using Mantid::HistogramData::Points;
class XValidationTest : public CxxTest::TestSuite {
public:
// This pair of boilerplate methods prevent the suite being created statically
// This means the constructor isn't called when running other tests
static XValidationTest *createSuite() { return new XValidationTest(); }
static void destroySuite(XValidationTest *suite) { delete suite; }
// FixedLengthVector contains validation, so we need to take a detour for
// creating a potentially invalid HistogramX.
HistogramX makeX(std::initializer_list<double> list) {
std::vector<double> helper(list.size());
std::iota(helper.begin(), helper.end(), 0.0);
HistogramX x(helper);
std::copy(list.begin(), list.end(), x.begin());
return x;
}
void test_works_for_HistogramX() {
HistogramX data{1.0, 2.0};
TS_ASSERT(isValid(data));
data[1] = data[0];
TS_ASSERT(!isValid(data));
}
void test_works_for_BinEdges() {
BinEdges data{1.0, 2.0};
TS_ASSERT(isValid(data));
data.mutableRawData()[1] = data[0];
TS_ASSERT(!isValid(data));
}
void test_works_for_Points() {
Points data{1.0, 2.0};
TS_ASSERT(isValid(data));
data.mutableRawData()[1] = data[0];
TS_ASSERT(!isValid(data));
}
void test_detects_zero_width() { TS_ASSERT(!isValid(makeX({1.0, 2.0, 2.0, 3.0}))); }
void test_detects_non_increasing() { TS_ASSERT(!isValid(makeX({1.0, 3.0, 2.0, 4.0}))); }
void test_accepts_nan() {
TS_ASSERT(isValid(makeX({NAN})));
TS_ASSERT(isValid(makeX({NAN, 1.0})));
TS_ASSERT(isValid(makeX({NAN, -1.0})));
TS_ASSERT(isValid(makeX({1.0, NAN})));
TS_ASSERT(isValid(makeX({-1.0, NAN})));
}
void test_detects_non_boundary_nan() {
TS_ASSERT(!isValid(makeX({-1.0, NAN, 1.0})));
TS_ASSERT(isValid(makeX({NAN, -1.0, 0.0, 1.0})));
TS_ASSERT(!isValid(makeX({NAN, -1.0, NAN, 1.0})));
TS_ASSERT(isValid(makeX({-1.0, 0.0, 1.0, NAN})));
TS_ASSERT(!isValid(makeX({-1.0, NAN, 1.0, NAN})));
}
void test_accepts_inf() {
TS_ASSERT(isValid(makeX({INFINITY})));
TS_ASSERT(isValid(makeX({-INFINITY})));
TS_ASSERT(isValid(makeX({-INFINITY, 0.0})));
TS_ASSERT(isValid(makeX({0.0, INFINITY})));
TS_ASSERT(isValid(makeX({-INFINITY, INFINITY})));
TS_ASSERT(isValid(makeX({-DBL_MAX / 2.0, DBL_MAX / 2.0})));
TS_ASSERT(isValid(makeX({-DBL_MAX, DBL_MAX})));
}
void test_detects_non_increasing_inf() {
// INF is ok, but order must be correct
TS_ASSERT(!isValid(makeX({0.0, -INFINITY})));
TS_ASSERT(!isValid(makeX({INFINITY, 0.0})));
TS_ASSERT(!isValid(makeX({INFINITY, -INFINITY})));
TS_ASSERT(!isValid(makeX({DBL_MAX, -DBL_MAX})));
}
void test_denormal() {
// Denormal values are ok
TS_ASSERT(isValid(makeX({0.0})));
TS_ASSERT(isValid(makeX({DBL_MIN / 2.0})));
TS_ASSERT(isValid(makeX({DBL_MIN / 2.0, 1.0})));
TS_ASSERT(isValid(makeX({-1.0, DBL_MIN / 2.0})));
}
void test_detects_denormal() {
// Denormal differences are not ok
TS_ASSERT(isValid(makeX({0.0, DBL_MIN})));
TS_ASSERT(!isValid(makeX({0.0, DBL_MIN / 2.0})));
TS_ASSERT(!isValid(makeX({DBL_MIN / 2.0, DBL_MIN})));
}
};