geom3d 0.2.0

Data structures and algorithms for 3D geometric modeling.
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

use crate::Float;
// use cached::{proc_macro::cached, SizedCache};

// fn make_key(n: usize, u: Float) -> u64 {
//     u64::from_le_bytes(u.to_le_bytes()) | (n as u64)
// }
//
// #[cached(
//     type = "SizedCache<u64, Vec<Float>>",
//     create = "{ SizedCache::with_size(1000) }",
//     convert = r#"{ make_key(n, u) }"#
// )]

/// Compute values of (n-1)th-degree Bernstein polynomial
pub fn bernstein(n: usize, u: Float) -> Vec<Float> {
    let mut values = vec![0.0; n];
    values[0] = 1.0;
    let u1 = 1.0 - u;

    for j in 1..n {
        let mut saved = 0.0;
        for k in 0..j {
            let temp = values[k];
            values[k] = saved + u1 * temp;
            saved = u * temp;
        }
        values[j] = saved;
    }

    values
}

#[test]
fn test_bernstein() {
    assert_eq!(
        bernstein(4, 0.25),
        vec![0.421875, 0.421875, 0.140625, 0.015625,]
    );
}