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
114
115
116
117
118
119
120
121
122
123
124
125
126
|
#include <imgen.h>
#include <cmath>
#include <thread>
#include <vector>
#include <algorithm>
std::ptrdiff_t minkowski_metric(
const double p,
const std::ptrdiff_t refX, const std::ptrdiff_t refY,
const std::ptrdiff_t x, const std::ptrdiff_t y
) {
return static_cast<std::ptrdiff_t>(
std::nearbyint(
std::pow(
std::pow(std::abs(refX - x), p) + std::pow(std::abs(refY - y), p),
1.0 / p
)
)
);
}
std::ptrdiff_t manhattan_metric(
const std::ptrdiff_t refX, const std::ptrdiff_t refY,
const std::ptrdiff_t x, const std::ptrdiff_t y
) {
return minkowski_metric(1, refX, refY, x, y);
}
std::ptrdiff_t euclidean_metric(
const std::ptrdiff_t refX, const std::ptrdiff_t refY,
const std::ptrdiff_t x, const std::ptrdiff_t y
) {
return minkowski_metric(2, refX, refY, x, y);
}
void generate_minkowski_voronoi(const double p) {
constexpr std::array<imgen::colored_vector, 9> ref{
imgen::colored_vector{ 0, 0, imgen::red() },
imgen::colored_vector{ 240, 200, imgen::color{220, 220, 220} },
imgen::colored_vector{-100, 230, imgen::color{ 94, 113, 106} },
imgen::colored_vector{ 120, -100, imgen::color{140, 146, 172} },
imgen::colored_vector{ -42, -200, imgen::color{128, 128, 128} },
imgen::colored_vector{ 120, 40, imgen::color{ 16, 20, 22} },
imgen::colored_vector{-150, 50, imgen::color{192, 192, 192} },
imgen::colored_vector{ 60, -128, imgen::color{178, 190, 181} },
imgen::colored_vector{-240, -20, imgen::color{ 54, 69, 79} }
};
imgen::write_ppm(
"voronoi_" + std::to_string(p) + ".ppm",
512,
512,
[&ref, p](std::ptrdiff_t x, std::ptrdiff_t y) -> imgen::color {
std::array<std::ptrdiff_t, 9> distances;
std::transform(
ref.begin(),
ref.end(),
distances.begin(),
[x, y, p](const imgen::colored_vector& pos) {
return minkowski_metric(p, std::get<0>(pos), std::get<1>(pos), x, y);
});
const auto& minimal_distance = std::min_element(
distances.begin(),
distances.end()
);
const imgen::colored_vector& nearest = ref[
std::distance(distances.begin(), minimal_distance)
];
if ( *minimal_distance <= 5 ) {
return imgen::black();
} else {
return std::get<2>(nearest);
}
}
);
}
void generate_minkowski_voronoi(
const double lower,
const double upper,
const double epsilon
) {
const auto n{ std::nearbyint( ( upper - lower ) / epsilon) };
for ( std::size_t i = 0; i < n; ++i ) {
generate_minkowski_voronoi(lower + i * epsilon);
}
}
void generate_parallel_minkowski_voronoi(
const unsigned int thread_count,
const double lower,
const double upper,
const double epsilon
) {
std::vector<std::thread> threads;
const double step = ( upper - lower ) / thread_count;
double offset = lower;
while ( threads.size() < thread_count ) {
threads.emplace_back([offset, step, epsilon]{
generate_minkowski_voronoi(
offset,
offset + step,
epsilon
);
});
offset += step;
}
generate_minkowski_voronoi(upper);
for ( auto& thread : threads ) {
thread.join();
}
}
int main(int, char*[]) {
generate_parallel_minkowski_voronoi(4, 1.0, 2.0, 0.025);
}
|