packed_spatial_index 0.18.0

Packed static spatial index (Hilbert R-tree) for 2D/3D AABBs — SIMD range, kNN, raycast, and spatial-join queries, with zero-copy and streaming serialization.
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
93
94
95
96
97
98
99
100
101
102
103
104

use super::columns::{F32ColumnRefs2D, F32ColumnRefs3D};
use crate::persistence::{MetaFields, PayloadError, write_index_container};

/// Write 2D f32 SoA columns into a `PSINDEX` container.
pub(crate) fn write_columns2d(
    out: &mut Vec<u8>,
    columns: F32ColumnRefs2D<'_>,
    interleaved: bool,
    payloads: Option<&[&[u8]]>,
    record_stride: Option<u32>,
    meta: &MetaFields<'_>,
) -> Result<(), PayloadError> {
    debug_assert_eq!(columns.min_xs.len(), columns.min_ys.len());
    debug_assert_eq!(columns.min_xs.len(), columns.max_xs.len());
    debug_assert_eq!(columns.min_xs.len(), columns.max_ys.len());
    debug_assert_eq!(columns.min_xs.len(), columns.indices.len());
    write_index_container(
        out,
        2,
        4,
        interleaved,
        columns.num_items,
        columns.min_xs.len(),
        columns.node_size,
        |bytes| {
            #[cfg(feature = "stream")]
            if interleaved {
                bytes.write_interleaved_f32_2d(
                    columns.min_xs,
                    columns.min_ys,
                    columns.max_xs,
                    columns.max_ys,
                    columns.indices,
                );
                return;
            }
            bytes.write_soa_boxes_f32_2d(
                columns.min_xs,
                columns.min_ys,
                columns.max_xs,
                columns.max_ys,
            );
            bytes.write_usize_slice_as_u64(columns.indices);
        },
        payloads,
        record_stride,
        &columns.indices[..columns.num_items],
        meta,
    )
}

/// Write 3D f32 SoA columns into a `PSINDEX` container.
pub(crate) fn write_columns3d(
    out: &mut Vec<u8>,
    columns: F32ColumnRefs3D<'_>,
    interleaved: bool,
    payloads: Option<&[&[u8]]>,
    record_stride: Option<u32>,
    meta: &MetaFields<'_>,
) -> Result<(), PayloadError> {
    debug_assert_eq!(columns.min_xs.len(), columns.min_ys.len());
    debug_assert_eq!(columns.min_xs.len(), columns.min_zs.len());
    debug_assert_eq!(columns.min_xs.len(), columns.max_xs.len());
    debug_assert_eq!(columns.min_xs.len(), columns.max_ys.len());
    debug_assert_eq!(columns.min_xs.len(), columns.max_zs.len());
    debug_assert_eq!(columns.min_xs.len(), columns.indices.len());
    write_index_container(
        out,
        3,
        4,
        interleaved,
        columns.num_items,
        columns.min_xs.len(),
        columns.node_size,
        |bytes| {
            #[cfg(feature = "stream")]
            if interleaved {
                bytes.write_interleaved_f32_3d(
                    columns.min_xs,
                    columns.min_ys,
                    columns.min_zs,
                    columns.max_xs,
                    columns.max_ys,
                    columns.max_zs,
                    columns.indices,
                );
                return;
            }
            bytes.write_soa_boxes_f32_3d(
                columns.min_xs,
                columns.min_ys,
                columns.min_zs,
                columns.max_xs,
                columns.max_ys,
                columns.max_zs,
            );
            bytes.write_usize_slice_as_u64(columns.indices);
        },
        payloads,
        record_stride,
        &columns.indices[..columns.num_items],
        meta,
    )
}