morton 0.3.1

Morton space filling curve functions
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

use rand::{rng, Rng};

use morton::utils::{idx_tile, idx_tile_tuple};
use morton::{deinterleave_morton, interleave_morton};

#[test]
fn interleave() {
    let mut tile_morton = [0; 32 * 32]; // 1024 locations
    let mut tile_normal = [0; 32 * 32]; // 1024 locations
                                        // fill tiles with same random numbers
    for x in 0..32 {
        for y in 0..32 {
            let random = rng().random::<u32>();
            tile_morton[interleave_morton(x as u16, y as u16) as usize] = random;
            tile_normal[idx_tile(x, y, 32)] = random;
        }
    }
    // check that the same random numbers are stored there
    // (morton curve did not override it's own elements)
    for x in 0..32 {
        for y in 0..32 {
            let morton = tile_morton[interleave_morton(x as u16, y as u16) as usize];
            let normal = tile_normal[idx_tile(x, y, 32)];
            assert!(morton == normal);
        }
    }
}

#[test]
fn deinterleave() {
    let mut tile_morton = [0; 32 * 32]; // 1024 locations
    let mut tile_normal = [0; 32 * 32]; // 1024 locations
                                        // fill tiles with same random numbers
    for x in 0..32 {
        for y in 0..32 {
            let random = rng().random::<u32>();
            tile_morton[interleave_morton(x as u16, y as u16) as usize] = random;
            tile_normal[idx_tile(x, y, 32)] = random;
        }
    }
    // check that the same random numbers are stored there
    // (morton curve did not override it's own elements)
    for z in 0..1024 {
        let morton = tile_morton[z];
        let normal = tile_normal[idx_tile_tuple(deinterleave_morton(z as u32), 32)];
        assert!(morton == normal);
    }
}

#[test]
fn deinterleave_interleave() {
    for z in 0..65536 {
        let (x, y) = deinterleave_morton(z);
        let morton = interleave_morton(x, y);
        assert!(morton == z);
    }
}

#[test]
fn interleave_deinterleave() {
    for x in 0..1024 {
        for y in 0..1024 {
            let morton = interleave_morton(x, y);
            let (d_x, d_y) = deinterleave_morton(morton);
            assert!(d_x == x && d_y == y);
        }
    }
}

// tests with random input
#[test]
fn rand_interleave_deinterleave_1000() {
    for _ in 0..1024 {
        let x = rng().random::<u16>();
        let y = rng().random::<u16>();
        let morton = interleave_morton(x, y);
        let (d_x, d_y) = deinterleave_morton(morton);
        assert!(d_x == x && d_y == y);
    }
}

#[test]
fn rand_deinterleave_interleave_1000() {
    for _ in 0..1024 {
        let z = rng().random::<u32>();
        let (x, y) = deinterleave_morton(z);
        let morton = interleave_morton(x, y);
        assert!(morton == z);
    }
}