|
1 | 1 | use crate::internal_prelude::*; |
2 | | -use hdf5_sys::h5d::H5Dchunk_iter; |
3 | 2 |
|
4 | | -/// Borrowed version of [ChunkInfo](crate::dataset::ChunkInfo) |
5 | | -#[derive(Debug)] |
6 | | -pub struct ChunkInfoBorrowed<'a> { |
7 | | - pub offset: &'a [u64], |
| 3 | +#[cfg(feature = "1.10.5")] |
| 4 | +use hdf5_sys::h5d::{H5Dget_chunk_info, H5Dget_num_chunks}; |
| 5 | + |
| 6 | +#[cfg(feature = "1.10.5")] |
| 7 | +#[derive(Clone, Debug, PartialEq, Eq)] |
| 8 | +pub struct ChunkInfo { |
| 9 | + /// Array with a size equal to the dataset’s rank whose elements contain 0-based |
| 10 | + /// logical positions of the chunk’s first element in each dimension. |
| 11 | + pub offset: Vec<u64>, |
| 12 | + /// Filter mask that indicates which filters were used with the chunk when written. |
| 13 | + /// A zero value indicates that all enabled filters are applied on the chunk. |
| 14 | + /// A filter is skipped if the bit corresponding to the filter’s position in |
| 15 | + /// the pipeline (0 ≤ position < 32) is turned on. |
8 | 16 | pub filter_mask: u32, |
| 17 | + /// Chunk address in the file. |
9 | 18 | pub addr: u64, |
| 19 | + /// Chunk size in bytes. |
10 | 20 | pub size: u64, |
11 | 21 | } |
12 | 22 |
|
13 | | -struct RustCallback<F> { |
14 | | - ndims: usize, |
15 | | - callback: F, |
16 | | -} |
| 23 | +#[cfg(feature = "1.10.5")] |
| 24 | +impl ChunkInfo { |
| 25 | + pub(crate) fn new(ndim: usize) -> Self { |
| 26 | + let offset = vec![0; ndim]; |
| 27 | + Self { offset, filter_mask: 0, addr: 0, size: 0 } |
| 28 | + } |
17 | 29 |
|
18 | | -extern "C" fn chunks_callback<F>( |
19 | | - offset: *const hsize_t, filter_mask: u32, addr: haddr_t, nbytes: u32, op_data: *mut c_void, |
20 | | -) -> herr_t |
21 | | -where |
22 | | - F: FnMut(ChunkInfoBorrowed) -> i32, |
23 | | -{ |
24 | | - unsafe { |
25 | | - std::panic::catch_unwind(|| { |
26 | | - let data: *mut RustCallback<F> = op_data.cast::<RustCallback<F>>(); |
27 | | - let ndims = (*data).ndims; |
28 | | - let callback = &mut (*data).callback; |
29 | | - |
30 | | - let offset = std::slice::from_raw_parts(offset, ndims); |
31 | | - |
32 | | - let info = |
33 | | - ChunkInfoBorrowed { offset, filter_mask, addr: addr as u64, size: nbytes as u64 }; |
34 | | - |
35 | | - callback(info) |
36 | | - }) |
37 | | - .unwrap_or(-1) |
| 30 | + /// Returns positional indices of disabled filters. |
| 31 | + pub fn disabled_filters(&self) -> Vec<usize> { |
| 32 | + (0..32).filter(|i| self.filter_mask & (1 << i) != 0).collect() |
38 | 33 | } |
39 | 34 | } |
40 | 35 |
|
41 | | -pub(crate) fn visit<F>(ds: &Dataset, callback: F) -> Result<()> |
42 | | -where |
43 | | - F: for<'a> FnMut(ChunkInfoBorrowed<'a>) -> i32, |
44 | | -{ |
45 | | - let mut data = RustCallback::<F> { ndims: ds.ndim(), callback }; |
46 | | - |
47 | | - h5try!(H5Dchunk_iter( |
48 | | - ds.id(), |
49 | | - H5P_DEFAULT, |
50 | | - Some(chunks_callback::<F>), |
51 | | - std::ptr::addr_of_mut!(data).cast() |
52 | | - )); |
| 36 | +#[cfg(feature = "1.10.5")] |
| 37 | +pub(crate) fn chunk_info(ds: &Dataset, index: usize) -> Option<ChunkInfo> { |
| 38 | + if !ds.is_chunked() { |
| 39 | + return None; |
| 40 | + } |
| 41 | + h5lock!(ds.space().map_or(None, |s| { |
| 42 | + let mut chunk_info = ChunkInfo::new(ds.ndim()); |
| 43 | + h5check(H5Dget_chunk_info( |
| 44 | + ds.id(), |
| 45 | + s.id(), |
| 46 | + index as _, |
| 47 | + chunk_info.offset.as_mut_ptr(), |
| 48 | + &mut chunk_info.filter_mask, |
| 49 | + &mut chunk_info.addr, |
| 50 | + &mut chunk_info.size, |
| 51 | + )) |
| 52 | + .map(|_| chunk_info) |
| 53 | + .ok() |
| 54 | + })) |
| 55 | +} |
53 | 56 |
|
54 | | - Ok(()) |
| 57 | +#[cfg(feature = "1.10.5")] |
| 58 | +pub(crate) fn get_num_chunks(ds: &Dataset) -> Option<usize> { |
| 59 | + if !ds.is_chunked() { |
| 60 | + return None; |
| 61 | + } |
| 62 | + h5lock!(ds.space().map_or(None, |s| { |
| 63 | + let mut n: hsize_t = 0; |
| 64 | + h5check(H5Dget_num_chunks(ds.id(), s.id(), &mut n)).map(|_| n as _).ok() |
| 65 | + })) |
55 | 66 | } |
56 | 67 |
|
57 | | -#[cfg(test)] |
58 | | -mod test { |
| 68 | +#[cfg(feature = "1.13.0")] |
| 69 | +mod one_thirteen { |
59 | 70 | use super::*; |
| 71 | + use hdf5_sys::h5d::H5Dchunk_iter; |
| 72 | + /// Borrowed version of [ChunkInfo](crate::dataset::ChunkInfo) |
| 73 | + #[derive(Debug)] |
| 74 | + pub struct ChunkInfoBorrowed<'a> { |
| 75 | + pub offset: &'a [u64], |
| 76 | + pub filter_mask: u32, |
| 77 | + pub addr: u64, |
| 78 | + pub size: u64, |
| 79 | + } |
| 80 | + |
| 81 | + struct RustCallback<F> { |
| 82 | + ndims: usize, |
| 83 | + callback: F, |
| 84 | + } |
| 85 | + |
| 86 | + extern "C" fn chunks_callback<F>( |
| 87 | + offset: *const hsize_t, filter_mask: u32, addr: haddr_t, nbytes: u32, op_data: *mut c_void, |
| 88 | + ) -> herr_t |
| 89 | + where |
| 90 | + F: FnMut(ChunkInfoBorrowed) -> i32, |
| 91 | + { |
| 92 | + unsafe { |
| 93 | + std::panic::catch_unwind(|| { |
| 94 | + let data: *mut RustCallback<F> = op_data.cast::<RustCallback<F>>(); |
| 95 | + let ndims = (*data).ndims; |
| 96 | + let callback = &mut (*data).callback; |
| 97 | + |
| 98 | + let offset = std::slice::from_raw_parts(offset, ndims); |
| 99 | + |
| 100 | + let info = ChunkInfoBorrowed { |
| 101 | + offset, |
| 102 | + filter_mask, |
| 103 | + addr: addr as u64, |
| 104 | + size: nbytes as u64, |
| 105 | + }; |
| 106 | + |
| 107 | + callback(info) |
| 108 | + }) |
| 109 | + .unwrap_or(-1) |
| 110 | + } |
| 111 | + } |
| 112 | + |
| 113 | + pub(crate) fn visit<F>(ds: &Dataset, callback: F) -> Result<()> |
| 114 | + where |
| 115 | + F: for<'a> FnMut(ChunkInfoBorrowed<'a>) -> i32, |
| 116 | + { |
| 117 | + let mut data = RustCallback::<F> { ndims: ds.ndim(), callback }; |
| 118 | + |
| 119 | + h5try!(H5Dchunk_iter( |
| 120 | + ds.id(), |
| 121 | + H5P_DEFAULT, |
| 122 | + Some(chunks_callback::<F>), |
| 123 | + std::ptr::addr_of_mut!(data).cast() |
| 124 | + )); |
| 125 | + |
| 126 | + Ok(()) |
| 127 | + } |
| 128 | + |
| 129 | + #[cfg(test)] |
| 130 | + mod test { |
| 131 | + use super::*; |
| 132 | + |
| 133 | + #[test] |
| 134 | + fn chunks_visit() { |
| 135 | + with_tmp_file(|f| { |
| 136 | + let ds = f.new_dataset::<i16>().no_chunk().shape((4, 4)).create("nochunk").unwrap(); |
| 137 | + assert_err_re!(visit(&ds, |_| 0), "not a chunked dataset"); |
| 138 | + |
| 139 | + let ds = |
| 140 | + f.new_dataset::<i16>().shape([3, 2]).chunk([1, 1]).create("chunk").unwrap(); |
| 141 | + ds.write(&ndarray::arr2(&[[1, 2], [3, 4], [5, 6]])).unwrap(); |
60 | 142 |
|
61 | | - #[test] |
62 | | - fn chunks_visit() { |
63 | | - with_tmp_file(|f| { |
64 | | - let ds = f.new_dataset::<i16>().no_chunk().shape((4, 4)).create("nochunk").unwrap(); |
65 | | - assert_err_re!(visit(&ds, |_| 0), "not a chunked dataset"); |
66 | | - |
67 | | - let ds = f.new_dataset::<i16>().shape([3, 2]).chunk([1, 1]).create("chunk").unwrap(); |
68 | | - ds.write(&ndarray::arr2(&[[1, 2], [3, 4], [5, 6]])).unwrap(); |
69 | | - |
70 | | - let mut i = 0; |
71 | | - let f = |c: ChunkInfoBorrowed| { |
72 | | - match i { |
73 | | - 0 => assert_eq!(c.offset, [0, 0]), |
74 | | - 1 => assert_eq!(c.offset, [0, 1]), |
75 | | - 2 => assert_eq!(c.offset, [1, 0]), |
76 | | - 3 => assert_eq!(c.offset, [1, 1]), |
77 | | - 4 => assert_eq!(c.offset, [2, 0]), |
78 | | - 5 => assert_eq!(c.offset, [2, 1]), |
79 | | - _ => unreachable!(), |
80 | | - } |
81 | | - assert_eq!(c.size, std::mem::size_of::<i16>() as u64); |
82 | | - i += 1; |
83 | | - 0 |
84 | | - }; |
85 | | - visit(&ds, f).unwrap(); |
86 | | - assert_eq!(i, 6); |
87 | | - }) |
| 143 | + let mut i = 0; |
| 144 | + let f = |c: ChunkInfoBorrowed| { |
| 145 | + match i { |
| 146 | + 0 => assert_eq!(c.offset, [0, 0]), |
| 147 | + 1 => assert_eq!(c.offset, [0, 1]), |
| 148 | + 2 => assert_eq!(c.offset, [1, 0]), |
| 149 | + 3 => assert_eq!(c.offset, [1, 1]), |
| 150 | + 4 => assert_eq!(c.offset, [2, 0]), |
| 151 | + 5 => assert_eq!(c.offset, [2, 1]), |
| 152 | + _ => unreachable!(), |
| 153 | + } |
| 154 | + assert_eq!(c.size, std::mem::size_of::<i16>() as u64); |
| 155 | + i += 1; |
| 156 | + 0 |
| 157 | + }; |
| 158 | + visit(&ds, f).unwrap(); |
| 159 | + assert_eq!(i, 6); |
| 160 | + }) |
| 161 | + } |
88 | 162 | } |
89 | 163 | } |
| 164 | +#[cfg(feature = "1.13.0")] |
| 165 | +pub use one_thirteen::*; |
0 commit comments