Fullerene 0.2.1

Fullerene on the ODE (Open Dynamics Engine) for Rust
Documentation 
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

//! icosahedron.rs
//!
//! Icosahedron
//! vertices (3 x 4) = 12
//! 3 x 20 faces (vertices 60 / 5 = 12)
//! edges 20 x 3 / 2 = 30
//!

// use num::{Float, zero, one};
use num::Float;

use ph_faces::Polyhedron;

/// Icosahedron
#[derive(Debug)]
pub struct Icosahedron<F: Float> {
  /// polyhedron tri: Vec 20 of Vec 1 indexed triangles
  pub ph: Polyhedron<F>,
  /// edges (duplex)
  pub edges: Vec<(u8, [u8; 5])>
}

/// Icosahedron
impl<F: Float> Icosahedron<F> {
  /// construct
  pub fn new(r: F) -> Self {
    let o = <F>::from(0).unwrap(); // = zero(); // = zero::<F>();
    let i = <F>::from(1).unwrap() * r; // = one(); // = one::<F>();
/*
    let d = <F>::from(2).unwrap();
    let r5 = <F>::from(5).unwrap().sqrt() * r;
    let p = (r5 + i) / d; // golden ratio
*/
    let p = <F>::from((5.0.sqrt() + 1.0) / 2.0).unwrap() * r; // golden ratio
    let ph = Polyhedron{
    vtx: vec![
      [o, -i, -p], [o, i, -p], [o, i, p], [o, -i, p],
      [-p, o, -i], [-p, o, i], [p, o, i], [p, o, -i],
      [-i, -p, o], [i, -p, o], [i, p, o], [-i, p, o]
    ],
    tri: vec![
/*
      [0, 8, 4], // triple x 1
      [0, 9, 8], [4, 1, 0], [8, 5, 4], // double x 3
      [0, 1, 7], [0, 7, 9], // single x 2 (x = 0)
      [4, 5, 11], [4, 11, 1], // single x 2 (y = 0)
      [8, 9, 3], [8, 3, 5], // single x 2 (z = 0)
      [6, 9, 7], [6, 3, 9], // single x 2 (y = 0)
      [10, 1, 11], [10, 7, 1], // single x 2 (z = 0)
      [2, 5, 3], [2, 11, 5], // single x 2 (x = 0)
      [2, 3, 6], [6, 7, 10], [10, 11, 2], // double x 3
      [2, 6, 10] // triple x 1
*/
      vec![[0, 8, 4]],
      vec![[0, 9, 8]],
      vec![[0, 7, 9]],
      vec![[0, 1, 7]],
      vec![[4, 1, 0]],
      vec![[4, 11, 1]],
      vec![[4, 5, 11]],
      vec![[8, 5, 4]],
      vec![[8, 3, 5]],
      vec![[8, 9, 3]],
      vec![[6, 3, 9]],
      vec![[6, 9, 7]],
      vec![[6, 7, 10]],
      vec![[10, 7, 1]],
      vec![[10, 1, 11]],
      vec![[10, 11, 2]],
      vec![[2, 11, 5]],
      vec![[2, 5, 3]],
      vec![[2, 3, 6]],
      vec![[2, 6, 10]]
    ]
    };
    let edges = vec![ // (duplex)
      (0, [1, 7, 9, 8, 4]), // 6 2 11 8 9 (x = 0) 6 + 8 9 9 1 9
      (4, [5, 11, 1, 0, 8]), // (y = 0)
      (8, [9, 3, 5, 4, 0]), // (z = 0)
      (9, [8, 0, 7, 6, 3]), // (z = 0)
      (7, [6, 9, 0, 1, 10]), // (y = 0)
      (1, [0, 4, 11, 10, 7]), // 4 7 11 9 5 (x = 0) 4 + 3 4 10 8 11
      (11, [10, 1, 4, 5, 2]), // (z = 0)
      (5, [4, 8, 3, 2, 11]), // (y = 0)
      (3, [2, 5, 8, 9, 6]), // 3 3 1 9 8 (x = 0) 3 + 0 10 8 11 7
      (6, [7, 10, 2, 3, 9]), // (y = 0)
      (10, [11, 2, 6, 7, 1]), // (z = 0)
      (2, [3, 6, 10, 11, 5]) // 3 4 1 6 10 (x = 0) 3 + 1 9 5 4 5
    ];
    Icosahedron{ph, edges}
  }
}