-
Notifications
You must be signed in to change notification settings - Fork 5
/
snowflake.rs
52 lines (47 loc) · 1.58 KB
/
snowflake.rs
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
use pgfplots::{
axis::{plot::*, *},
Engine, Picture,
};
fn main() {
let mut vertices = {
let mut current: Vec<(f64, f64)> = vec![
(0.0, 1.0),
((3.0f64).sqrt() / 2.0, -0.5),
(-(3.0f64).sqrt() / 2.0, -0.5),
];
for _ in 0..5 {
current = snowflake_iter(¤t[..]);
}
current
};
vertices.push(vertices[0]);
// Set snowflake
let mut plot = Plot2D::new();
plot.coordinates = vertices.into_iter().map(|v| v.into()).collect();
plot.add_key(PlotKey::Custom(String::from("fill=gray!20")));
// Customise axis environment
let mut axis = Axis::from(plot);
axis.set_title("Kloch Snowflake");
axis.add_key(AxisKey::Custom(String::from("hide axis")));
#[cfg(feature = "tectonic")]
Picture::from(axis).show_pdf(Engine::Tectonic).unwrap();
#[cfg(not(feature = "tectonic"))]
Picture::from(axis).show_pdf(Engine::PdfLatex).unwrap();
}
// Stolen from plotters crate example
fn snowflake_iter(points: &[(f64, f64)]) -> Vec<(f64, f64)> {
let mut ret = vec![];
for i in 0..points.len() {
let (start, end) = (points[i], points[(i + 1) % points.len()]);
let t = ((end.0 - start.0) / 3.0, (end.1 - start.1) / 3.0);
let s = (
t.0 * 0.5 - t.1 * (0.75f64).sqrt(),
t.1 * 0.5 + (0.75f64).sqrt() * t.0,
);
ret.push(start);
ret.push((start.0 + t.0, start.1 + t.1));
ret.push((start.0 + t.0 + s.0, start.1 + t.1 + s.1));
ret.push((start.0 + t.0 * 2.0, start.1 + t.1 * 2.0));
}
ret
}