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
|
// SPDX-License-Identifier: LGPL-3.0-or-later
// Author: Kristian Lytje
#include <utility/Axis.h>
#include <utility/Limit.h>
#include <ostream>
#include <cmath>
using namespace ausaxs;
Axis::Axis() noexcept : bins(0), min(0), max(0) {}
Axis::Axis(const Limit& limits, int bins) noexcept : bins(bins), min(limits.min), max(limits.max) {}
Axis& Axis::operator=(std::initializer_list<double> list) noexcept {
std::vector<double> d = list;
bins = std::round(d[0]);
min = d[1];
max = d[2];
return *this;
}
std::string Axis::to_string() const noexcept {
return "Axis: (" + std::to_string(min) + ", " + std::to_string(max) + ") with " + std::to_string(bins) + " bins";
}
bool Axis::operator==(const Axis& rhs) const noexcept = default;
void Axis::resize(unsigned int bins) noexcept {
auto w = width();
this->bins = bins;
this->max = min + bins*w;
}
bool Axis::empty() const noexcept {return bins==0;}
Limit Axis::limits() const noexcept {return Limit(min, max);}
unsigned int Axis::get_bin(double value) const noexcept {
if (bins == 0) [[unlikely]] {return 0;}
if (value <= min) {return 0;}
if (value >= max) {return bins;}
return std::floor((value+1e-6-min)/width()); // +1e-6 to avoid flooring floating point errors, and we will likely never have bins this small anyway
}
double Axis::get_bin_value(unsigned int bin) const noexcept {
if (bins == 0) [[unlikely]] {return 0;}
return min + bin*width();
}
Axis Axis::sub_axis(double vmin, double vmax) const noexcept {
unsigned int min_bin = get_bin(vmin);
unsigned int max_bin = get_bin(vmax);
double new_min = get_bin_value(min_bin);
double new_max = get_bin_value(max_bin);
return Axis(new_min, new_max, max_bin - min_bin);
}
std::ostream& operator<<(std::ostream& os, const Axis& axis) noexcept {os << axis.to_string(); return os;}
|
