projects / nihav.git / blame
| Commit | Line | Data |
|---|---|---|
| afe29743 | 1 | use formats; |
| 88c58b1b | 2 | use super::super::*; |
| afe29743 KS |
3 | use io::byteio::*; |
| 4 | use std::io::SeekFrom; | |
| 5 | use std::mem; | |
| 6 | ||
| 7 | struct IviDeltaCB { | |
| 8 | quad_radix: u8, | |
| 9 | data: &'static [i8], | |
| 10 | } | |
| 11 | ||
| 12 | #[derive(Clone, Copy)] | |
| 13 | struct MV { | |
| 14 | x: i8, | |
| 15 | y: i8 | |
| 16 | } | |
| 17 | ||
| 18 | struct Buffers { | |
| 19 | width: usize, | |
| 20 | height: usize, | |
| 21 | buf1: Vec<u8>, | |
| 22 | buf2: Vec<u8>, | |
| 23 | fbuf: bool, | |
| 24 | } | |
| 25 | ||
| 26 | const DEFAULT_PIXEL: u8 = 0x40; | |
| 27 | ||
| 28 | impl Buffers { | |
| 29 | fn new() -> Self { Buffers { width: 0, height: 0, buf1: Vec::new(), buf2: Vec::new(), fbuf: true } } | |
| 30 | fn reset(&mut self) { | |
| 31 | self.width = 0; | |
| 32 | self.height = 0; | |
| 33 | self.buf1.truncate(0); | |
| 34 | self.buf2.truncate(0); | |
| 35 | } | |
| 36 | fn alloc(&mut self, w: usize, h: usize) { | |
| 37 | self.width = w; | |
| 38 | self.height = h; | |
| 39 | self.buf1.resize(w * h + (w >> 2) * (h >> 2) * 2, DEFAULT_PIXEL); | |
| 40 | self.buf2.resize(w * h + (w >> 2) * (h >> 2) * 2, DEFAULT_PIXEL); | |
| 41 | } | |
| 42 | fn flip(&mut self) { self.fbuf = !self.fbuf; } | |
| 43 | fn get_stride(&mut self, planeno: usize) -> usize { | |
| 44 | if planeno == 0 { self.width } else { self.width >> 2 } | |
| 45 | } | |
| 46 | fn get_offset(&mut self, planeno: usize) -> usize { | |
| 47 | match planeno { | |
| 48 | 1 => self.width * self.height, | |
| 49 | 2 => self.width * self.height + (self.width >> 2) * (self.height >> 2), | |
| 50 | _ => 0, | |
| 51 | } | |
| 52 | } | |
| 53 | fn fill_framebuf(&mut self, fbuf: &mut NAVideoBuffer<u8>) { | |
| 54 | for planeno in 0..3 { | |
| 55 | let mut soff = self.get_offset(planeno); | |
| 56 | let mut doff = fbuf.get_offset(planeno); | |
| 57 | let sstride = self.get_stride(planeno); | |
| 58 | let dstride = fbuf.get_stride(planeno); | |
| 59 | let width = if planeno == 0 { self.width } else { self.width >> 2 }; | |
| 60 | let height = if planeno == 0 { self.height } else { self.height >> 2 }; | |
| 61 | let src = if self.fbuf { &self.buf1[0..] } else { &self.buf2[0..] }; | |
| 62 | let mut dst = fbuf.get_data_mut(); | |
| 63 | for _ in 0..height { | |
| 64 | for x in 0..width { | |
| 65 | dst[doff + x] = src[soff + x] * 2; | |
| 66 | } | |
| 67 | soff += sstride; | |
| 68 | doff += dstride; | |
| 69 | } | |
| 70 | } | |
| 71 | } | |
| 72 | fn copy_block(&mut self, doff: usize, soff: usize, stride: usize, w: usize, h: usize) { | |
| 73 | let mut sidx = soff; | |
| 74 | let mut didx = doff; | |
| 75 | if self.fbuf { | |
| 76 | for _ in 0..h { | |
| 77 | for i in 0..w { self.buf1[didx + i] = self.buf2[sidx + i]; } | |
| 78 | sidx += stride; | |
| 79 | didx += stride; | |
| 80 | } | |
| 81 | } else { | |
| 82 | for _ in 0..h { | |
| 83 | for i in 0..w { self.buf2[didx + i] = self.buf1[sidx + i]; } | |
| 84 | sidx += stride; | |
| 85 | didx += stride; | |
| 86 | } | |
| 87 | } | |
| 88 | } | |
| 89 | fn fill_block(&mut self, doff: usize, stride: usize, w: usize, h: usize, topline: bool) { | |
| 90 | let mut didx = doff; | |
| 91 | let mut buf: [u8; 8] = [0; 8]; | |
| 92 | if topline { | |
| 93 | if self.fbuf { | |
| 94 | for _ in 0..h { | |
| 95 | for i in 0..w { self.buf1[didx + i] = DEFAULT_PIXEL; } | |
| 96 | didx += stride; | |
| 97 | } | |
| 98 | } else { | |
| 99 | for _ in 0..h { | |
| 100 | for i in 0..w { self.buf2[didx + i] = DEFAULT_PIXEL; } | |
| 101 | didx += stride; | |
| 102 | } | |
| 103 | } | |
| 104 | } else { | |
| 105 | if self.fbuf { | |
| 106 | for i in 0..w { buf[i] = self.buf1[didx - stride + i]; } | |
| 107 | for _ in 0..h { | |
| 108 | for i in 0..w { self.buf1[didx + i] = buf[i]; } | |
| 109 | didx += stride; | |
| 110 | } | |
| 111 | } else { | |
| 112 | for i in 0..w { buf[i] = self.buf2[didx - stride + i]; } | |
| 113 | for _ in 0..h { | |
| 114 | for i in 0..w { self.buf2[didx + i] = buf[i]; } | |
| 115 | didx += stride; | |
| 116 | } | |
| 117 | } | |
| 118 | } | |
| 119 | } | |
| 120 | } | |
| 121 | ||
| 122 | #[allow(unused_variables)] | |
| 123 | fn apply_delta4x4(bufs: &mut Buffers, off: usize, stride: usize, | |
| 124 | deltas: &[u8], topline: bool, first_line: bool) { | |
| 125 | let mut dst = if bufs.fbuf { &mut bufs.buf1[off..(off + 4)] } | |
| 126 | else { &mut bufs.buf2[off..(off + 4)] }; | |
| 127 | for i in 0..4 { dst[i] = dst[i].wrapping_add(deltas[i]) & 0x7F; } | |
| 128 | } | |
| 129 | ||
| 130 | #[allow(unused_variables)] | |
| 131 | fn apply_delta4x8(bufs: &mut Buffers, off: usize, stride: usize, | |
| 132 | deltas: &[u8], topline: bool, first_line: bool) { | |
| 133 | let mut dst = if bufs.fbuf { &mut bufs.buf1[off..(off + 4 + stride)] } | |
| 134 | else { &mut bufs.buf2[off..(off + 4 + stride)] }; | |
| 135 | for i in 0..4 { dst[i + stride] = dst[i].wrapping_add(deltas[i]) & 0x7F; } | |
| 136 | if !topline { | |
| 137 | for i in 0..4 { dst[i] = (dst[i + stride] + dst[i]) >> 1; } | |
| 138 | } else { | |
| 139 | for i in 0..4 { dst[i] = dst[i + stride]; } | |
| 140 | } | |
| 141 | } | |
| 142 | ||
| 143 | #[allow(unused_variables)] | |
| 144 | fn apply_delta4x8m11(bufs: &mut Buffers, off: usize, stride: usize, | |
| 145 | deltas: &[u8], topline: bool, first_line: bool) { | |
| 146 | let mut dst = if bufs.fbuf { &mut bufs.buf1[off..(off + 4 + stride)] } | |
| 147 | else { &mut bufs.buf2[off..(off + 4 + stride)] }; | |
| 148 | for i in 0..4 { dst[i] = dst[i] .wrapping_add(deltas[i]) & 0x7F; } | |
| 149 | for i in 0..4 { dst[i + stride] = dst[i + stride].wrapping_add(deltas[i]) & 0x7F; } | |
| 150 | } | |
| 151 | ||
| 152 | #[allow(unused_variables)] | |
| 153 | fn apply_delta8x8p(bufs: &mut Buffers, off: usize, stride: usize, | |
| 154 | deltas: &[u8], topline: bool, first_line: bool) { | |
| 155 | let mut dst = if bufs.fbuf { &mut bufs.buf1[off..(off + 8 + stride)] } | |
| 156 | else { &mut bufs.buf2[off..(off + 8 + stride)] }; | |
| 157 | for i in 0..8 { dst[i] = dst[i] .wrapping_add(deltas[i >> 1]) & 0x7F; } | |
| 158 | for i in 0..8 { dst[i + stride] = dst[i + stride].wrapping_add(deltas[i >> 1]) & 0x7F; } | |
| 159 | } | |
| 160 | ||
| 161 | fn apply_delta8x8i(bufs: &mut Buffers, off: usize, stride: usize, | |
| 162 | deltas: &[u8], topline: bool, firstline: bool) { | |
| 163 | let mut dst = if bufs.fbuf { &mut bufs.buf1[off..(off + 8 + stride)] } | |
| 164 | else { &mut bufs.buf2[off..(off + 8 + stride)] }; | |
| 165 | if !firstline { | |
| 166 | for i in 0..8 { dst[i + stride] = dst[i ].wrapping_add(deltas[i >> 1]) & 0x7F; } | |
| 167 | } else { | |
| 168 | for i in 0..8 { dst[i + stride] = dst[i & !1].wrapping_add(deltas[i >> 1]) & 0x7F; } | |
| 169 | } | |
| 170 | if !topline { | |
| 171 | for i in 0..8 { dst[i] = (dst[i + stride] + dst[i]) >> 1; } | |
| 172 | } else { | |
| 173 | for i in 0..8 { dst[i] = dst[i + stride]; } | |
| 174 | } | |
| 175 | } | |
| 176 | ||
| 177 | fn copy_line_top(bufs: &mut Buffers, off: usize, stride: usize, bw: usize, topline: bool) { | |
| 178 | let mut buf: [u8; 8] = [0; 8]; | |
| 179 | if !topline { | |
| 180 | let src = if bufs.fbuf { &bufs.buf1[(off - stride)..(off - stride + bw)] } | |
| 181 | else { &bufs.buf2[(off - stride)..(off - stride + bw)] }; | |
| 182 | for i in 0..bw { buf[i] = src[i]; } | |
| 183 | } else { | |
| 184 | for i in 0..bw { buf[i] = DEFAULT_PIXEL; } | |
| 185 | } | |
| 186 | let mut dst = if bufs.fbuf { &mut bufs.buf1[off..(off + bw)] } | |
| 187 | else { &mut bufs.buf2[off..(off + bw)] }; | |
| 188 | for i in 0..bw { dst[i] = buf[i]; } | |
| 189 | } | |
| 190 | ||
| 191 | fn copy_line_top4x4(bufs: &mut Buffers, off: usize, stride: usize, topline: bool) { | |
| 192 | copy_line_top(bufs, off, stride, 4, topline); | |
| 193 | } | |
| 194 | ||
| 195 | fn copy_line_top4x8(bufs: &mut Buffers, off: usize, stride: usize, topline: bool) { | |
| 196 | copy_line_top(bufs, off, stride, 4, topline); | |
| 197 | copy_line_top(bufs, off + stride, stride, 4, false); | |
| 198 | } | |
| 199 | ||
| 200 | fn copy_line_top8x8(bufs: &mut Buffers, off: usize, stride: usize, topline: bool) { | |
| 201 | let mut buf: [u8; 8] = [0; 8]; | |
| 202 | if !topline { | |
| 203 | let src = if bufs.fbuf { &bufs.buf1[(off - stride)..(off - stride + 8)] } | |
| 204 | else { &bufs.buf2[(off - stride)..(off - stride + 8)] }; | |
| 205 | for i in 0..8 { buf[i] = src[i & !1]; } | |
| 206 | } else { | |
| 207 | for i in 0..8 { buf[i] = DEFAULT_PIXEL; } | |
| 208 | } | |
| 209 | let mut dst = if bufs.fbuf { &mut bufs.buf1[off..(off + 8)] } | |
| 210 | else { &mut bufs.buf2[off..(off + 8)] }; | |
| 211 | for i in 0..8 {dst[i] = buf[i]; } | |
| 212 | } | |
| 213 | ||
| 214 | fn fill_block8x8(bufs: &mut Buffers, doff: usize, stride: usize, h: usize, topline: bool, firstline: bool) { | |
| 215 | let mut didx = doff; | |
| 216 | let mut buf: [u8; 8] = [0; 8]; | |
| 217 | if firstline { | |
| 218 | for i in 0..8 { buf[i] = DEFAULT_PIXEL; } | |
| 219 | } else if bufs.fbuf { | |
| 220 | for i in 0..8 { buf[i] = bufs.buf1[doff - stride + i]; } | |
| 221 | } else { | |
| 222 | for i in 0..8 { buf[i] = bufs.buf1[doff - stride + i]; } | |
| 223 | } | |
| 224 | if topline && !firstline { | |
| 225 | for i in 0..4 { buf[i * 2 + 1] = buf[i * 2]; } | |
| 226 | if bufs.fbuf { | |
| 227 | for i in 0..8 { bufs.buf1[doff + i] = (bufs.buf1[doff - stride + i] + buf[i]) >> 1; } | |
| 228 | } else { | |
| 229 | for i in 0..8 { bufs.buf2[doff + i] = (bufs.buf2[doff - stride + i] + buf[i]) >> 1; } | |
| 230 | } | |
| 231 | } | |
| 232 | ||
| 233 | let start = if !topline { 0 } else { 1 }; | |
| 234 | if bufs.fbuf { | |
| 235 | for _ in start..h { | |
| 236 | for i in 0..8 { bufs.buf1[didx + i] = buf[i]; } | |
| 237 | didx += stride; | |
| 238 | } | |
| 239 | } else { | |
| 240 | for _ in start..h { | |
| 241 | for i in 0..8 { bufs.buf2[didx + i] = buf[i]; } | |
| 242 | didx += stride; | |
| 243 | } | |
| 244 | } | |
| 245 | } | |
| 246 | ||
| 247 | struct Indeo3Decoder { | |
| 248 | info: Rc<NACodecInfo>, | |
| 249 | bpos: u8, | |
| 250 | bbuf: u8, | |
| 251 | width: u16, | |
| 252 | height: u16, | |
| 253 | mvs: Vec<MV>, | |
| 254 | altquant: [u8; 16], | |
| 255 | vq_offset: u8, | |
| 256 | bufs: Buffers, | |
| 257 | } | |
| 258 | ||
| 259 | #[derive(Clone,Copy)] | |
| 260 | struct IV3Cell { | |
| 261 | x: u16, | |
| 262 | y: u16, | |
| 263 | w: u16, | |
| 264 | h: u16, | |
| 265 | d: u8, | |
| 266 | vqt: bool, | |
| 267 | mv: Option<MV>, | |
| 268 | } | |
| 269 | ||
| 270 | impl IV3Cell { | |
| 271 | fn new(w: u16, h: u16) -> Self { | |
| 272 | IV3Cell { x: 0, y: 0, w: w, h: h, d: 20, vqt: false, mv: None } | |
| 273 | } | |
| 274 | fn split_h(&self) -> (Self, Self) { | |
| 275 | let h1 = if self.h > 2 { ((self.h + 2) >> 2) << 1 } else { 1 }; | |
| 276 | let h2 = self.h - h1; | |
| 277 | let mut cell1 = *self; | |
| 278 | cell1.h = h1; | |
| 279 | cell1.d -= 1; | |
| 280 | let mut cell2 = *self; | |
| 281 | cell2.y += h1; | |
| 282 | cell2.h = h2; | |
| 283 | cell2.d -= 1; | |
| 284 | (cell1, cell2) | |
| 285 | } | |
| 286 | fn split_w(&self, stripw: u16) -> (Self, Self) { | |
| 287 | let w1 = if self.w > stripw { | |
| 288 | if self.w > stripw * 2 { stripw * 2 } else { stripw } | |
| 289 | } else { | |
| 290 | if self.w > 2 { ((self.w + 2) >> 2) << 1 } else { 1 } | |
| 291 | }; | |
| 292 | let w2 = self.w - w1; | |
| 293 | let mut cell1 = *self; | |
| 294 | cell1.w = w1; | |
| 295 | cell1.d -= 1; | |
| 296 | let mut cell2 = *self; | |
| 297 | cell2.x += w1; | |
| 298 | cell2.w = w2; | |
| 299 | cell2.d -= 1; | |
| 300 | (cell1, cell2) | |
| 301 | } | |
| 302 | fn no_mv(&self) -> bool { match self.mv { None => true, Some(_) => false } } | |
| 303 | } | |
| 304 | ||
| 305 | struct CellDecParams { | |
| 306 | tab: [usize; 2], | |
| 307 | bw: u16, | |
| 308 | bh: u16, | |
| 309 | swap_q: [bool; 2], | |
| 310 | hq: bool, | |
| 311 | apply_delta: fn (&mut Buffers, usize, usize, &[u8], bool, bool), | |
| 312 | copy_line_top: fn (&mut Buffers, usize, usize, bool), | |
| 313 | } | |
| 314 | ||
| 315 | const FRMH_TAG: u32 = ((b'F' as u32) << 24) | ((b'R' as u32) << 16) | |
| 316 | | ((b'M' as u32) << 8) | (b'H' as u32); | |
| 317 | ||
| 318 | const H_SPLIT: u8 = 0; | |
| 319 | const V_SPLIT: u8 = 1; | |
| 320 | const SKIP_OR_TREE: u8 = 2; | |
| 321 | ||
| 322 | impl Indeo3Decoder { | |
| 323 | fn new() -> Self { | |
| 324 | let dummy_info = Rc::new(DUMMY_CODEC_INFO); | |
| 325 | Indeo3Decoder { info: dummy_info, bpos: 0, bbuf: 0, width: 0, height: 0, | |
| 326 | mvs: Vec::new(), altquant: [0; 16], | |
| 327 | vq_offset: 0, bufs: Buffers::new() } | |
| 328 | } | |
| 329 | ||
| 330 | fn br_reset(&mut self) { | |
| 331 | self.bpos = 0; | |
| 332 | self.bbuf = 0; | |
| 333 | } | |
| 334 | ||
| 335 | fn get_2bits(&mut self, br: &mut ByteReader) -> DecoderResult<u8> { | |
| 336 | if self.bpos == 0 { | |
| 337 | self.bbuf = br.read_byte()?; | |
| 338 | self.bpos = 8; | |
| 339 | } | |
| 340 | self.bpos -= 2; | |
| 341 | Ok((self.bbuf >> self.bpos) & 0x3) | |
| 342 | } | |
| 343 | ||
| 344 | fn decode_cell_data(&mut self, br: &mut ByteReader, cell: IV3Cell, | |
| 345 | off: usize, stride: usize, params: CellDecParams) -> DecoderResult<()> { | |
| 346 | let blk_w = cell.w * 4 / params.bw; | |
| 347 | let blk_h = cell.h * 4 / params.bh; | |
| 348 | let scale: usize = if params.bh == 4 { 1 } else { 2 }; | |
| 349 | ||
| 350 | validate!((((cell.w * 4) % params.bw) == 0) && (((cell.h * 4) % params.bh) == 0)); | |
| 351 | ||
| 352 | let mut run_blocks = 0; | |
| 353 | let mut run_skip = false; | |
| 354 | ||
| 355 | let mut didx: usize = ((cell.x*4) as usize) + ((cell.y * 4) as usize) * stride + off; | |
| 356 | let mut sidx: usize; | |
| 357 | if cell.no_mv() { | |
| 358 | sidx = 0; | |
| 359 | } else { | |
| 360 | let mv = cell.mv.unwrap(); | |
| 361 | let mx = mv.x as i16; | |
| 362 | let my = mv.y as i16; | |
| 363 | let l = (cell.x as i16) * 4 + mx; | |
| 364 | let t = (cell.y as i16) * 4 + my; | |
| 365 | let r = ((cell.x + cell.w) as i16) * 4 + mx; | |
| 366 | let b = ((cell.y + cell.h) as i16) * 4 + my; | |
| 367 | validate!(l >= 0); | |
| 368 | validate!(t >= 0); | |
| 369 | validate!(r <= (self.width as i16)); | |
| 370 | validate!(b <= (self.height as i16)); | |
| 371 | sidx = (l as usize) + (t as usize) * stride + off; | |
| 372 | } | |
| 373 | for y in 0..blk_h { | |
| 374 | let mut xoff: usize = 0; | |
| 375 | for _ in 0..blk_w { | |
| 376 | if run_blocks > 0 { | |
| 377 | if !run_skip || !cell.no_mv() { | |
| 378 | if !(params.bw == 8 && cell.no_mv()) { | |
| 379 | if !cell.no_mv() { | |
| 380 | self.bufs.copy_block(didx + xoff, sidx + xoff, stride, | |
| 381 | params.bw as usize, params.bh as usize); | |
| 382 | } else { | |
| 383 | self.bufs.fill_block(didx + xoff, stride, | |
| 384 | params.bw as usize, params.bh as usize, | |
| 385 | (cell.y == 0) && (y == 0)); | |
| 386 | } | |
| 387 | } else { | |
| 388 | fill_block8x8(&mut self.bufs, | |
| 389 | didx + xoff, stride, 8, | |
| 390 | (y == 0), (cell.y == 0) && (y == 0)); | |
| 391 | } | |
| 392 | } | |
| 393 | run_blocks -= 1; | |
| 394 | } else { | |
| 395 | let mut line: usize = 0; | |
| 396 | while line < 4 { | |
| 397 | let c = br.read_byte()?; | |
| 398 | if c < 0xF8 { | |
| 399 | let delta_tab = if params.hq { | |
| 400 | IVI3_DELTA_CBS[params.tab[line & 1]] | |
| 401 | } else { | |
| 402 | IVI3_DELTA_CBS[params.tab[1]] | |
| 403 | }; | |
| 404 | let mut idx1; | |
| 405 | let mut idx2; | |
| 406 | if (c as usize) < delta_tab.data.len()/2 { | |
| 407 | idx1 = br.read_byte()? as usize; | |
| 408 | validate!(idx1 < delta_tab.data.len()); | |
| 409 | idx2 = c as usize; | |
| 410 | } else { | |
| 411 | let tmp = (c as usize) - delta_tab.data.len()/2; | |
| 412 | idx1 = tmp / (delta_tab.quad_radix as usize); | |
| 413 | idx2 = tmp % (delta_tab.quad_radix as usize); | |
| 414 | if params.swap_q[line & 1] { | |
| 415 | mem::swap(&mut idx1, &mut idx2); | |
| 416 | } | |
| 417 | } | |
| 418 | let deltas: [u8; 4] = [delta_tab.data[idx1 * 2] as u8, | |
| 419 | delta_tab.data[idx1 * 2 + 1] as u8, | |
| 420 | delta_tab.data[idx2 * 2 + 0] as u8, | |
| 421 | delta_tab.data[idx2 * 2 + 1] as u8]; | |
| 422 | let topline = (cell.y == 0) && (y == 0) && (line == 0); | |
| 423 | let first_line = (y == 0) && (line == 0); | |
| 424 | if cell.no_mv() { | |
| 425 | (params.copy_line_top)(&mut self.bufs, | |
| 426 | didx + xoff + line * scale * stride, | |
| 427 | stride, topline); | |
| 428 | } else { | |
| 429 | self.bufs.copy_block(didx + xoff + line * scale * stride, | |
| 430 | sidx + xoff + line * scale * stride, | |
| 431 | stride, params.bw as usize, scale); | |
| 432 | } | |
| 433 | (params.apply_delta)(&mut self.bufs, | |
| 434 | didx + xoff + line * scale * stride, | |
| 435 | stride, &deltas, topline, first_line); | |
| 436 | line += 1; | |
| 437 | } else { | |
| 438 | let mut tocopy: usize = 0; | |
| 439 | let mut do_copy = true; | |
| 440 | if c == 0xF8 { return Err(DecoderError::InvalidData); } | |
| 441 | if c == 0xF9 { | |
| 442 | run_blocks = 1; | |
| 443 | run_skip = true; | |
| 444 | validate!(line == 0); | |
| 445 | tocopy = 4; | |
| 446 | do_copy = !cell.no_mv(); | |
| 447 | } | |
| 448 | if c == 0xFA { | |
| 449 | validate!(line == 0); | |
| 450 | tocopy = 4; | |
| 451 | do_copy = !cell.no_mv(); | |
| 452 | } | |
| 453 | if c == 0xFB { | |
| 454 | let c = br.read_byte()?; | |
| 455 | validate!((c < 64) && ((c & 0x1F) != 0)); | |
| 456 | run_blocks = (c & 0x1F) - 1; | |
| 457 | run_skip = (c & 0x20) != 0; | |
| 458 | tocopy = 4 - line; | |
| 459 | if params.bw == 4 && cell.no_mv() && run_skip { | |
| 460 | do_copy = false; | |
| 461 | } | |
| 462 | } | |
| 463 | if c == 0xFC { | |
| 464 | run_skip = false; | |
| 465 | run_blocks = 1; | |
| 466 | tocopy = 4 - line; | |
| 467 | } | |
| 468 | if c >= 0xFD { | |
| 469 | let nl = 257 - (c as i16) - (line as i16); | |
| 470 | validate!(nl > 0); | |
| 471 | tocopy = nl as usize; | |
| 472 | } | |
| 473 | if do_copy { | |
| 474 | if !(params.bw == 8 && cell.no_mv()) { | |
| 475 | if !cell.no_mv() { | |
| 476 | self.bufs.copy_block(didx + xoff + line * scale * stride, | |
| 477 | sidx + xoff + line * scale * stride, | |
| 478 | stride, params.bw as usize, | |
| 479 | tocopy * scale); | |
| 480 | } else { | |
| 481 | self.bufs.fill_block(didx + xoff + line * scale * stride, | |
| 482 | stride, params.bw as usize, | |
| 483 | tocopy * scale, | |
| 484 | (cell.y == 0) && (y == 0) && (line == 0)); | |
| 485 | } | |
| 486 | } else { | |
| 487 | fill_block8x8(&mut self.bufs, | |
| 488 | didx + xoff + line * 2 * stride, | |
| 489 | stride, tocopy * 2, | |
| 490 | (y == 0) && (line == 0), | |
| 491 | (cell.y == 0) && (y == 0) && (line == 0)); | |
| 492 | } | |
| 493 | } | |
| 494 | line += tocopy; | |
| 495 | } | |
| 496 | } | |
| 497 | } | |
| 498 | xoff += params.bw as usize; | |
| 499 | } | |
| 500 | didx += stride * (params.bh as usize); | |
| 501 | sidx += stride * (params.bh as usize); | |
| 502 | } | |
| 503 | Ok(()) | |
| 504 | } | |
| 505 | ||
| 506 | fn copy_cell(&mut self, cell: IV3Cell, off: usize, stride: usize) -> DecoderResult<()> { | |
| 507 | if cell.no_mv() { return Err(DecoderError::InvalidData); } | |
| 508 | let mv = cell.mv.unwrap(); | |
| 509 | let mx = mv.x as i16; | |
| 510 | let my = mv.y as i16; | |
| 511 | let l = (cell.x as i16) * 4 + mx; | |
| 512 | let t = (cell.y as i16) * 4 + my; | |
| 513 | let r = ((cell.x + cell.w) as i16) * 4 + mx; | |
| 514 | let b = ((cell.y + cell.h) as i16) * 4 + my; | |
| 515 | validate!(l >= 0); | |
| 516 | validate!(t >= 0); | |
| 517 | validate!(r <= (self.width as i16)); | |
| 518 | validate!(b <= (self.height as i16)); | |
| 519 | let sidx: usize = off + (l as usize) + (t as usize) * stride; | |
| 520 | let didx: usize = off + ((cell.x * 4) as usize) + ((cell.y * 4) as usize) * stride; | |
| 521 | self.bufs.copy_block(didx, sidx, stride, (cell.w * 4) as usize, (cell.h * 4) as usize); | |
| 522 | Ok(()) | |
| 523 | } | |
| 524 | ||
| 525 | fn decode_cell(&mut self, br: &mut ByteReader, cell: IV3Cell, off: usize, | |
| 526 | stride: usize, intra: bool) -> DecoderResult<()> { | |
| 527 | let code = br.read_byte()?; | |
| 528 | let mode = code >> 4; | |
| 529 | let vq_idx = code & 0xF; | |
| 530 | ||
| 531 | let mut idx1: usize = vq_idx as usize; | |
| 532 | let mut idx2: usize = vq_idx as usize; | |
| 533 | if (mode == 1) || (mode == 4) { | |
| 534 | let c = self.altquant[vq_idx as usize]; | |
| 535 | idx1 = (c >> 4) as usize; | |
| 536 | idx2 = (c & 0xF) as usize; | |
| 537 | } | |
| 538 | ||
| 539 | idx1 += self.vq_offset as usize; | |
| 540 | idx2 += self.vq_offset as usize; | |
| 541 | validate!((idx1 < 24) && (idx2 < 24)); | |
| 542 | ||
| 543 | let mut cp = CellDecParams { | |
| 544 | tab: [idx2, idx1], | |
| 545 | bw: 0, bh: 0, | |
| 546 | swap_q: [idx2 >= 16, idx1 >= 16], | |
| 547 | hq: false, | |
| 548 | apply_delta: apply_delta4x4, | |
| 549 | copy_line_top: copy_line_top4x4, | |
| 550 | }; | |
| 551 | if (mode == 0) || (mode == 1) { | |
| 552 | cp.bw = 4; | |
| 553 | cp.bh = 4; | |
| 554 | cp.hq = true; | |
| 555 | } else if (mode == 3) || (mode == 4) { | |
| 556 | if !cell.no_mv() { return Err(DecoderError::InvalidData); } | |
| 557 | cp.bw = 4; | |
| 558 | cp.bh = 8; | |
| 559 | cp.hq = true; | |
| 560 | cp.apply_delta = apply_delta4x8; | |
| 561 | cp.copy_line_top = copy_line_top4x8; | |
| 562 | } else if mode == 10 { | |
| 563 | if !cell.no_mv() { | |
| 564 | validate!(!intra); | |
| 565 | cp.apply_delta = apply_delta8x8p; | |
| 566 | } else { | |
| 567 | cp.apply_delta = apply_delta8x8i; | |
| 568 | } | |
| 569 | cp.bw = 8; | |
| 570 | cp.bh = 8; | |
| 571 | cp.copy_line_top = copy_line_top8x8; | |
| 572 | } else if mode == 11 { | |
| 573 | if cell.no_mv() { return Err(DecoderError::InvalidData); } | |
| 574 | validate!(!intra); | |
| 575 | cp.bw = 4; | |
| 576 | cp.bh = 8; | |
| 577 | cp.apply_delta = apply_delta4x8m11; | |
| 578 | cp.copy_line_top = copy_line_top4x8; | |
| 579 | } else { | |
| 580 | return Err(DecoderError::InvalidData); | |
| 581 | } | |
| 582 | self.decode_cell_data(br, cell, off, stride, cp) | |
| 583 | } | |
| 584 | ||
| 585 | fn parse_tree(&mut self, br: &mut ByteReader, cell: IV3Cell, off: usize, | |
| 586 | stride: usize, stripw: u16, intra: bool) -> DecoderResult<()> { | |
| 587 | let op = self.get_2bits(br)?; | |
| 588 | if op == H_SPLIT { | |
| 589 | validate!(cell.h > 1); | |
| 590 | validate!(cell.d > 0); | |
| 591 | let (cell1, cell2) = cell.split_h(); | |
| 592 | self.parse_tree(br, cell1, off, stride, stripw, intra)?; | |
| 593 | self.parse_tree(br, cell2, off, stride, stripw, intra)?; | |
| 594 | Ok(()) | |
| 595 | } else if op == V_SPLIT { | |
| 596 | validate!(cell.w > 1); | |
| 597 | validate!(cell.d > 0); | |
| 598 | let (cell1, cell2) = cell.split_w(stripw); | |
| 599 | self.parse_tree(br, cell1, off, stride, stripw, intra)?; | |
| 600 | self.parse_tree(br, cell2, off, stride, stripw, intra)?; | |
| 601 | Ok(()) | |
| 602 | } else if op == SKIP_OR_TREE { | |
| 603 | if !cell.vqt { | |
| 604 | let mut newcell = cell; | |
| 605 | newcell.vqt = true; | |
| 606 | newcell.d -= 1; | |
| 607 | self.parse_tree(br, newcell, off, stride, stripw, intra) | |
| 608 | } else { | |
| 609 | validate!(!intra); | |
| 610 | let code = self.get_2bits(br)?; | |
| 611 | validate!(code < 2); | |
| 612 | if code == 1 { return Err(DecoderError::NotImplemented); } | |
| 613 | self.copy_cell(cell, off, stride) | |
| 614 | } | |
| 615 | } else { | |
| 616 | if !cell.vqt { | |
| 617 | let mut newcell = cell; | |
| 618 | newcell.vqt = true; | |
| 619 | newcell.d -= 1; | |
| 620 | let mv_idx = br.read_byte()? as usize; | |
| 621 | validate!(mv_idx < self.mvs.len()); | |
| 622 | newcell.mv = Some(self.mvs[mv_idx]); | |
| 623 | self.parse_tree(br, newcell, off, stride, stripw, intra) | |
| 624 | } else { | |
| 625 | self.decode_cell(br, cell, off, stride, intra) | |
| 626 | } | |
| 627 | } | |
| 628 | } | |
| 629 | ||
| 630 | fn decode_plane_intra(&mut self, br: &mut ByteReader, planeno: usize, | |
| 631 | start: u64, end: u64) -> DecoderResult<()> { | |
| 632 | let offs = self.bufs.get_offset(planeno); | |
| 633 | let stride = self.bufs.get_stride(planeno); | |
| 634 | br.seek(SeekFrom::Start(start))?; | |
| 635 | ||
| 636 | let nvec = br.read_u32le()?; | |
| 637 | validate!(nvec == 0); // for intra there should be no mc_vecs | |
| 638 | self.mvs.truncate(0); | |
| 639 | for _ in 0..nvec { | |
| 640 | let x = br.read_byte()? as i8; | |
| 641 | let y = br.read_byte()? as i8; | |
| 642 | self.mvs.push(MV{ x: x, y: y }); | |
| 643 | } | |
| 644 | ||
| 645 | let shift = if planeno == 0 { 2 } else { 4 }; | |
| 646 | let cell = IV3Cell::new((self.bufs.width >> shift) as u16, | |
| 647 | (self.bufs.height >> shift) as u16); | |
| 648 | self.br_reset(); | |
| 649 | self.parse_tree(br, cell, offs, stride, if planeno > 0 { 10 } else { 40 }, true)?; | |
| 650 | validate!(br.tell() <= end); | |
| 651 | Ok(()) | |
| 652 | } | |
| 653 | ||
| 654 | fn decode_plane_inter(&mut self, br: &mut ByteReader, planeno: usize, | |
| 655 | start: u64, end: u64) -> DecoderResult<()> { | |
| 656 | let offs = self.bufs.get_offset(planeno); | |
| 657 | let stride = self.bufs.get_stride(planeno); | |
| 658 | br.seek(SeekFrom::Start(start))?; | |
| 659 | ||
| 660 | let nvec = br.read_u32le()?; | |
| 661 | validate!(nvec <= 256); // for intra there should be no mc_vecs | |
| 662 | self.mvs.truncate(0); | |
| 663 | for _ in 0..nvec { | |
| 664 | let y = br.read_byte()? as i8; | |
| 665 | let x = br.read_byte()? as i8; | |
| 666 | self.mvs.push(MV{ x: x, y: y }); | |
| 667 | } | |
| 668 | ||
| 669 | let shift = if planeno == 0 { 2 } else { 4 }; | |
| 670 | let cell = IV3Cell::new((self.bufs.width >> shift) as u16, | |
| 671 | (self.bufs.height >> shift) as u16); | |
| 672 | self.br_reset(); | |
| 673 | self.parse_tree(br, cell, offs, stride, if planeno > 0 { 10 } else { 40 }, false)?; | |
| 674 | validate!(br.tell() <= end); | |
| 675 | Ok(()) | |
| 676 | } | |
| 677 | } | |
| 678 | ||
| 679 | const FLAG_KEYFRAME: u16 = 1 << 2; | |
| 680 | const FLAG_NONREF: u16 = 1 << 8; | |
| 681 | ||
| 682 | impl NADecoder for Indeo3Decoder { | |
| 683 | fn init(&mut self, info: Rc<NACodecInfo>) -> DecoderResult<()> { | |
| 684 | if let NACodecTypeInfo::Video(vinfo) = info.get_properties() { | |
| 685 | let w = vinfo.get_width(); | |
| 686 | let h = vinfo.get_height(); | |
| 687 | let fmt = formats::YUV410_FORMAT; | |
| 688 | let myinfo = NACodecTypeInfo::Video(NAVideoInfo::new(w, h, false, fmt)); | |
| 689 | self.info = Rc::new(NACodecInfo::new_ref(info.get_name(), myinfo, info.get_extradata())); | |
| 690 | self.bufs.reset(); | |
| 691 | Ok(()) | |
| 692 | } else { | |
| 693 | Err(DecoderError::InvalidData) | |
| 694 | } | |
| 695 | } | |
| 696 | fn decode(&mut self, pkt: &NAPacket) -> DecoderResult<NAFrameRef> { | |
| 697 | let src = pkt.get_buffer(); | |
| 698 | let mut mr = MemoryReader::new_read(&src); | |
| 699 | let mut br = ByteReader::new(&mut mr); | |
| 700 | let frameno = br.read_u32le()?; | |
| 701 | let hdr_2 = br.read_u32le()?; | |
| 702 | let check = br.read_u32le()?; | |
| 703 | let size = br.read_u32le()?; | |
| 704 | ||
| 705 | let data_start = br.tell(); | |
| 706 | ||
| 707 | if (frameno ^ hdr_2 ^ size ^ FRMH_TAG) != check { | |
| 708 | return Err(DecoderError::InvalidData); | |
| 709 | } | |
| 710 | if (size as i64) > br.left() { return Err(DecoderError::InvalidData); } | |
| 711 | let ver = br.read_u16le()?; | |
| 712 | if ver != 32 { return Err(DecoderError::NotImplemented); } | |
| 713 | let flags = br.read_u16le()?; | |
| 714 | let size2 = br.read_u32le()?; | |
| 715 | validate!(((size2 + 7) >> 3) <= size); | |
| 716 | let cb = br.read_byte()?; | |
| 717 | self.vq_offset = cb; | |
| 718 | br.read_skip(3)?; | |
| 719 | let height = br.read_u16le()?; | |
| 720 | let width = br.read_u16le()?; | |
| 721 | validate!((width >= 16) && (width <= 640)); | |
| 722 | validate!((height >= 16) && (height <= 640)); | |
| 723 | validate!(((width & 3) == 0) && ((height & 3) == 0)); | |
| 724 | if (self.bufs.width != (width as usize)) || (self.bufs.height != (height as usize)) { | |
| 725 | self.bufs.alloc(width as usize, height as usize); | |
| 726 | } | |
| 727 | self.width = width; | |
| 728 | self.height = height; | |
| 729 | ||
| 730 | let yoff = br.read_u32le()?; | |
| 731 | let uoff = br.read_u32le()?; | |
| 732 | let voff = br.read_u32le()?; | |
| 733 | if yoff > size { return Err(DecoderError::InvalidData); } | |
| 734 | if uoff > size { return Err(DecoderError::InvalidData); } | |
| 735 | if voff > size { return Err(DecoderError::InvalidData); } | |
| 736 | ||
| 737 | br.read_skip(4)?; | |
| 738 | br.read_buf(&mut self.altquant)?; | |
| 739 | ||
| 740 | let mut yend = src.len() as u32;//size; | |
| 741 | if (uoff < yend) && (uoff > yoff) { yend = uoff; } | |
| 742 | if (voff < yend) && (voff > yoff) { yend = voff; } | |
| 743 | let mut uend = size; | |
| 744 | if (yoff < uend) && (yoff > uoff) { uend = yoff; } | |
| 745 | if (voff < uend) && (voff > uoff) { uend = voff; } | |
| 746 | let mut vend = size; | |
| 747 | if (yoff < vend) && (yoff > voff) { vend = yoff; } | |
| 748 | if (uoff < vend) && (uoff > voff) { vend = uoff; } | |
| 749 | ||
| 5658bdcc | 750 | let intraframe = (flags & FLAG_KEYFRAME) != 0; |
| afe29743 KS |
751 | let vinfo = self.info.get_properties().get_video_info().unwrap(); |
| 752 | let bufret = alloc_video_buffer(vinfo, 2); | |
| 753 | if let Err(_) = bufret { return Err(DecoderError::InvalidData); } | |
| 754 | let mut bufinfo = bufret.unwrap(); | |
| 755 | let mut buf = bufinfo.get_vbuf().unwrap(); | |
| 756 | let ystart = data_start + (yoff as u64); | |
| 757 | let ustart = data_start + (uoff as u64); | |
| 758 | let vstart = data_start + (voff as u64); | |
| 759 | let yendpos = data_start + (yend as u64); | |
| 760 | let uendpos = data_start + (uend as u64); | |
| 761 | let vendpos = data_start + (vend as u64); | |
| 5658bdcc | 762 | if intraframe { |
| afe29743 KS |
763 | self.decode_plane_intra(&mut br, 0, ystart, yendpos)?; |
| 764 | self.decode_plane_intra(&mut br, 1, ustart, uendpos)?; | |
| 765 | self.decode_plane_intra(&mut br, 2, vstart, vendpos)?; | |
| 766 | } else { | |
| 767 | self.decode_plane_inter(&mut br, 0, ystart, yendpos)?; | |
| 768 | self.decode_plane_inter(&mut br, 1, ustart, uendpos)?; | |
| 769 | self.decode_plane_inter(&mut br, 2, vstart, vendpos)?; | |
| 770 | } | |
| 771 | self.bufs.fill_framebuf(&mut buf); | |
| 772 | if (flags & FLAG_NONREF) == 0 { self.bufs.flip(); } | |
| 773 | let mut frm = NAFrame::new_from_pkt(pkt, self.info.clone(), bufinfo); | |
| 5658bdcc KS |
774 | frm.set_keyframe(intraframe); |
| 775 | frm.set_frame_type(if intraframe { FrameType::I } else { FrameType::P }); | |
| afe29743 KS |
776 | Ok(Rc::new(RefCell::new(frm))) |
| 777 | } | |
| 778 | } | |
| 779 | ||
| 780 | #[cfg(test)] | |
| 781 | mod test { | |
| 782 | use codecs::*; | |
| 783 | use demuxers::*; | |
| 784 | use io::byteio::*; | |
| 785 | use std::fs::File; | |
| 786 | ||
| 787 | #[test] | |
| 788 | fn test_indeo3() { | |
| 789 | let avi_dmx = find_demuxer("avi").unwrap(); | |
| 790 | let mut file = File::open("assets/iv32_example.avi").unwrap(); | |
| 791 | let mut fr = FileReader::new_read(&mut file); | |
| 792 | let mut br = ByteReader::new(&mut fr); | |
| 793 | let mut dmx = avi_dmx.new_demuxer(&mut br); | |
| 794 | dmx.open().unwrap(); | |
| 795 | ||
| 796 | let mut decs: Vec<Option<Box<NADecoder>>> = Vec::new(); | |
| 797 | for i in 0..dmx.get_num_streams() { | |
| 798 | let s = dmx.get_stream(i).unwrap(); | |
| 799 | let info = s.get_info(); | |
| 800 | let decfunc = find_decoder(info.get_name()); | |
| 801 | if let Some(df) = decfunc { | |
| 802 | let mut dec = (df)(); | |
| 803 | dec.init(info).unwrap(); | |
| 804 | decs.push(Some(dec)); | |
| 805 | } else { | |
| 806 | decs.push(None); | |
| 807 | panic!("decoder {} not found", info.get_name()); | |
| 808 | } | |
| 809 | } | |
| 810 | ||
| 811 | loop { | |
| 812 | let pktres = dmx.get_frame(); | |
| 813 | if let Err(e) = pktres { | |
| 814 | if e == DemuxerError::EOF { break; } | |
| 815 | } | |
| 816 | let pkt = pktres.unwrap(); | |
| 817 | let streamno = pkt.get_stream().get_id() as usize; | |
| 818 | if let Some(ref mut dec) = decs[streamno] { | |
| df14fb81 KS |
819 | // let frm = |
| 820 | dec.decode(&pkt).unwrap(); | |
| 821 | // if pkt.get_stream().get_info().is_video() { | |
| 822 | // write_pgmyuv("iv3", streamno, pkt.get_pts().unwrap(), frm); | |
| 823 | // } else { | |
| 824 | // write_sound("iv3", streamno, frm, pkt.get_pts().unwrap() == 0); | |
| 825 | // } | |
| afe29743 KS |
826 | } |
| 827 | if pkt.get_pts().unwrap() > 10 { break; } | |
| 828 | } | |
| 829 | } | |
| 830 | } | |
| 831 | pub fn get_decoder() -> Box<NADecoder> { | |
| 832 | Box::new(Indeo3Decoder::new()) | |
| 833 | } | |
| 834 | ||
| 835 | const DT_1_1: IviDeltaCB = IviDeltaCB{ quad_radix: 7, data: &[ | |
| 836 | 0, 0, 2, 2, -2, -2, -1, 3, | |
| 837 | 1, -3, 3, -1, -3, 1, 4, 4, | |
| 838 | -4, -4, 1, 5, -1, -5, 5, 1, | |
| 839 | -5, -1, -4, 4, 4, -4, -2, 6, | |
| 840 | 2, -6, 6, -2, -6, 2, 4, 9, | |
| 841 | -4, -9, 9, 4, -9, -4, 9, 9, | |
| 842 | -9, -9, 1, 10, -1, -10, 10, 1, | |
| 843 | -10, -1, -5, 8, 5, -8, 8, -5, | |
| 844 | -8, 5, 9, 15, -9, -15, 15, 9, | |
| 845 | -15, -9, -3, 12, 3, -12, 12, -3, | |
| 846 | -12, 3, 4, 16, -4, -16, 16, 4, | |
| 847 | -16, -4, 16, 16, -16, -16, 0, 18, | |
| 848 | 0, -18, 18, 0, -18, 0, -12, 12, | |
| 849 | 12, -12, -9, 16, 9, -16, 16, -9, | |
| 850 | -16, 9, 11, 27, -11, -27, 27, 11, | |
| 851 | -27, -11, 19, 28, -19, -28, 28, 19, | |
| 852 | -28, -19, -6, 22, 6, -22, 22, -6, | |
| 853 | -22, 6, 4, 29, -4, -29, 29, 4, | |
| 854 | -29, -4, 30, 30, -30, -30, -2, 33, | |
| 855 | 2, -33, 33, -2, -33, 2, -18, 23, | |
| 856 | 18, -23, 23, -18, -23, 18, -15, 30, | |
| 857 | 15, -30, 30, -15, -30, 15, 22, 46, | |
| 858 | -22, -46, 46, 22, -46, -22, 13, 47, | |
| 859 | -13, -47, 47, 13, -47, -13, 35, 49, | |
| 860 | -35, -49, 49, 35, -49, -35, -11, 41, | |
| 861 | 11, -41, 41, -11, -41, 11, 4, 51, | |
| 862 | -4, -51, 51, 4, -51, -4, 54, 54, | |
| 863 | -54, -54, -34, 34, 34, -34, -29, 42, | |
| 864 | 29, -42, 42, -29, -42, 29, -6, 60, | |
| 865 | 6, -60, 60, -6, -60, 6, 27, 76, | |
| 866 | -27, -76, 76, 27, -76, -27, 43, 77, | |
| 867 | -43, -77, 77, 43, -77, -43, -24, 55, | |
| 868 | 24, -55, 55, -24, -55, 24, 14, 79, | |
| 869 | -14, -79, 79, 14, -79, -14, 63, 83, | |
| 870 | -63, -83, 83, 63, -83, -63, -20, 74, | |
| 871 | 20, -74, 74, -20, -74, 20, 2, 88, | |
| 872 | -2, -88, 88, 2, -88, -2, 93, 93, | |
| 873 | -93, -93, -52, 61, 52, -61, 61, -52, | |
| 874 | -61, 52, 52, 120, -52, -120, 120, 52, | |
| 875 | -120, -52, -45, 75, 45, -75, 75, -45, | |
| 876 | -75, 45, 75, 125, -75, -125, 125, 75, | |
| 877 | -125, -75, 33, 122, -33, -122, 122, 33, | |
| 878 | -122, -33, -13, 103, 13, -103, 103, -13, | |
| 879 | -103, 13, -40, 96, 40, -96, 96, -40, | |
| 880 | -96, 40, -34, 127, 34, -127, 127, -34, | |
| 881 | -127, 34, -89, 89, 89, -89, -78, 105, | |
| 882 | 78, -105, 105, -78, -105, 78, 12, 12, | |
| 883 | -12, -12, 23, 23, -23, -23, 42, 42, | |
| 884 | -42, -42, 73, 73, -73, -73, | |
| 885 | ]}; | |
| 886 | ||
| 887 | const DT_1_2: IviDeltaCB = IviDeltaCB{ quad_radix: 9, data: &[ | |
| 888 | 0, 0, 3, 3, -3, -3, -1, 4, | |
| 889 | 1, -4, 4, -1, -4, 1, 7, 7, | |
| 890 | -7, -7, 2, 8, -2, -8, 8, 2, | |
| 891 | -8, -2, -2, 9, 2, -9, 9, -2, | |
| 892 | -9, 2, -6, 6, 6, -6, 6, 13, | |
| 893 | -6, -13, 13, 6, -13, -6, 13, 13, | |
| 894 | -13, -13, 1, 14, -1, -14, 14, 1, | |
| 895 | -14, -1, -8, 12, 8, -12, 12, -8, | |
| 896 | -12, 8, 14, 23, -14, -23, 23, 14, | |
| 897 | -23, -14, -5, 18, 5, -18, 18, -5, | |
| 898 | -18, 5, 6, 24, -6, -24, 24, 6, | |
| 899 | -24, -6, 24, 24, -24, -24, -1, 27, | |
| 900 | 1, -27, 27, -1, -27, 1, -17, 17, | |
| 901 | 17, -17, -13, 23, 13, -23, 23, -13, | |
| 902 | -23, 13, 16, 40, -16, -40, 40, 16, | |
| 903 | -40, -16, 28, 41, -28, -41, 41, 28, | |
| 904 | -41, -28, -9, 33, 9, -33, 33, -9, | |
| 905 | -33, 9, 6, 43, -6, -43, 43, 6, | |
| 906 | -43, -6, 46, 46, -46, -46, -4, 50, | |
| 907 | 4, -50, 50, -4, -50, 4, -27, 34, | |
| 908 | 27, -34, 34, -27, -34, 27, -22, 45, | |
| 909 | 22, -45, 45, -22, -45, 22, 34, 69, | |
| 910 | -34, -69, 69, 34, -69, -34, 19, 70, | |
| 911 | -19, -70, 70, 19, -70, -19, 53, 73, | |
| 912 | -53, -73, 73, 53, -73, -53, -17, 62, | |
| 913 | 17, -62, 62, -17, -62, 17, 5, 77, | |
| 914 | -5, -77, 77, 5, -77, -5, 82, 82, | |
| 915 | -82, -82, -51, 51, 51, -51, -43, 64, | |
| 916 | 43, -64, 64, -43, -64, 43, -10, 90, | |
| 917 | 10, -90, 90, -10, -90, 10, 41, 114, | |
| 918 | -41, -114, 114, 41, -114, -41, 64, 116, | |
| 919 | -64, -116, 116, 64, -116, -64, -37, 82, | |
| 920 | 37, -82, 82, -37, -82, 37, 22, 119, | |
| 921 | -22, -119, 119, 22, -119, -22, 95, 124, | |
| 922 | -95, -124, 124, 95, -124, -95, -30, 111, | |
| 923 | 30, -111, 111, -30, -111, 30, -78, 92, | |
| 924 | 78, -92, 92, -78, -92, 78, -68, 113, | |
| 925 | 68, -113, 113, -68, -113, 68, 18, 18, | |
| 926 | -18, -18, 34, 34, -34, -34, 63, 63, | |
| 927 | -63, -63, 109, 109, -109, -109, | |
| 928 | ]}; | |
| 929 | ||
| 930 | const DT_1_3: IviDeltaCB = IviDeltaCB{ quad_radix: 10, data: &[ | |
| 931 | 0, 0, 4, 4, -4, -4, -1, 5, | |
| 932 | 1, -5, 5, -1, -5, 1, 3, 10, | |
| 933 | -3, -10, 10, 3, -10, -3, 9, 9, | |
| 934 | -9, -9, -7, 7, 7, -7, -3, 12, | |
| 935 | 3, -12, 12, -3, -12, 3, 8, 17, | |
| 936 | -8, -17, 17, 8, -17, -8, 17, 17, | |
| 937 | -17, -17, 1, 19, -1, -19, 19, 1, | |
| 938 | -19, -1, -11, 16, 11, -16, 16, -11, | |
| 939 | -16, 11, -6, 23, 6, -23, 23, -6, | |
| 940 | -23, 6, 18, 31, -18, -31, 31, 18, | |
| 941 | -31, -18, 8, 32, -8, -32, 32, 8, | |
| 942 | -32, -8, 33, 33, -33, -33, -1, 36, | |
| 943 | 1, -36, 36, -1, -36, 1, -23, 23, | |
| 944 | 23, -23, -17, 31, 17, -31, 31, -17, | |
| 945 | -31, 17, 21, 54, -21, -54, 54, 21, | |
| 946 | -54, -21, 37, 55, -37, -55, 55, 37, | |
| 947 | -55, -37, -12, 44, 12, -44, 44, -12, | |
| 948 | -44, 12, 8, 57, -8, -57, 57, 8, | |
| 949 | -57, -8, 61, 61, -61, -61, -5, 66, | |
| 950 | 5, -66, 66, -5, -66, 5, -36, 45, | |
| 951 | 36, -45, 45, -36, -45, 36, -29, 60, | |
| 952 | 29, -60, 60, -29, -60, 29, 45, 92, | |
| 953 | -45, -92, 92, 45, -92, -45, 25, 93, | |
| 954 | -25, -93, 93, 25, -93, -25, 71, 97, | |
| 955 | -71, -97, 97, 71, -97, -71, -22, 83, | |
| 956 | 22, -83, 83, -22, -83, 22, 7, 102, | |
| 957 | -7, -102, 102, 7, -102, -7, 109, 109, | |
| 958 | -109, -109, -68, 68, 68, -68, -57, 85, | |
| 959 | 57, -85, 85, -57, -85, 57, -13, 120, | |
| 960 | 13, -120, 120, -13, -120, 13, -49, 110, | |
| 961 | 49, -110, 110, -49, -110, 49, -104, 123, | |
| 962 | 104, -123, 123, -104, -123, 104, 24, 24, | |
| 963 | -24, -24, 46, 46, -46, -46, 84, 84, | |
| 964 | -84, -84, | |
| 965 | ]}; | |
| 966 | ||
| 967 | const DT_1_4: IviDeltaCB = IviDeltaCB{ quad_radix: 11, data: &[ | |
| 968 | 0, 0, 5, 5, -5, -5, -2, 7, | |
| 969 | 2, -7, 7, -2, -7, 2, 11, 11, | |
| 970 | -11, -11, 3, 13, -3, -13, 13, 3, | |
| 971 | -13, -3, -9, 9, 9, -9, -4, 15, | |
| 972 | 4, -15, 15, -4, -15, 4, 11, 22, | |
| 973 | -11, -22, 22, 11, -22, -11, 21, 21, | |
| 974 | -21, -21, 2, 24, -2, -24, 24, 2, | |
| 975 | -24, -2, -14, 20, 14, -20, 20, -14, | |
| 976 | -20, 14, 23, 38, -23, -38, 38, 23, | |
| 977 | -38, -23, -8, 29, 8, -29, 29, -8, | |
| 978 | -29, 8, 11, 39, -11, -39, 39, 11, | |
| 979 | -39, -11, 41, 41, -41, -41, -1, 45, | |
| 980 | 1, -45, 45, -1, -45, 1, -29, 29, | |
| 981 | 29, -29, -22, 39, 22, -39, 39, -22, | |
| 982 | -39, 22, 27, 67, -27, -67, 67, 27, | |
| 983 | -67, -27, 47, 69, -47, -69, 69, 47, | |
| 984 | -69, -47, -15, 56, 15, -56, 56, -15, | |
| 985 | -56, 15, 11, 71, -11, -71, 71, 11, | |
| 986 | -71, -11, 76, 76, -76, -76, -6, 83, | |
| 987 | 6, -83, 83, -6, -83, 6, -45, 57, | |
| 988 | 45, -57, 57, -45, -57, 45, -36, 75, | |
| 989 | 36, -75, 75, -36, -75, 36, 56, 115, | |
| 990 | -56, -115, 115, 56, -115, -56, 31, 117, | |
| 991 | -31, -117, 117, 31, -117, -31, 88, 122, | |
| 992 | -88, -122, 122, 88, -122, -88, -28, 104, | |
| 993 | 28, -104, 104, -28, -104, 28, -85, 85, | |
| 994 | 85, -85, -72, 106, 72, -106, 106, -72, | |
| 995 | -106, 72, 30, 30, -30, -30, 58, 58, | |
| 996 | -58, -58, 105, 105, -105, -105, | |
| 997 | ]}; | |
| 998 | ||
| 999 | const DT_1_5: IviDeltaCB = IviDeltaCB{ quad_radix: 12, data: &[ | |
| 1000 | 0, 0, 6, 6, -6, -6, -2, 8, | |
| 1001 | 2, -8, 8, -2, -8, 2, 13, 13, | |
| 1002 | -13, -13, 4, 15, -4, -15, 15, 4, | |
| 1003 | -15, -4, -11, 11, 11, -11, -5, 18, | |
| 1004 | 5, -18, 18, -5, -18, 5, 13, 26, | |
| 1005 | -13, -26, 26, 13, -26, -13, 26, 26, | |
| 1006 | -26, -26, 2, 29, -2, -29, 29, 2, | |
| 1007 | -29, -2, -16, 24, 16, -24, 24, -16, | |
| 1008 | -24, 16, 28, 46, -28, -46, 46, 28, | |
| 1009 | -46, -28, -9, 35, 9, -35, 35, -9, | |
| 1010 | -35, 9, 13, 47, -13, -47, 47, 13, | |
| 1011 | -47, -13, 49, 49, -49, -49, -1, 54, | |
| 1012 | 1, -54, 54, -1, -54, 1, -35, 35, | |
| 1013 | 35, -35, -26, 47, 26, -47, 47, -26, | |
| 1014 | -47, 26, 32, 81, -32, -81, 81, 32, | |
| 1015 | -81, -32, 56, 83, -56, -83, 83, 56, | |
| 1016 | -83, -56, -18, 67, 18, -67, 67, -18, | |
| 1017 | -67, 18, 13, 86, -13, -86, 86, 13, | |
| 1018 | -86, -13, 91, 91, -91, -91, -7, 99, | |
| 1019 | 7, -99, 99, -7, -99, 7, -54, 68, | |
| 1020 | 54, -68, 68, -54, -68, 54, -44, 90, | |
| 1021 | 44, -90, 90, -44, -90, 44, -33, 124, | |
| 1022 | 33, -124, 124, -33, -124, 33, -103, 103, | |
| 1023 | 103, -103, -86, 127, 86, -127, 127, -86, | |
| 1024 | -127, 86, 37, 37, -37, -37, 69, 69, | |
| 1025 | -69, -69, | |
| 1026 | ]}; | |
| 1027 | ||
| 1028 | const DT_1_6: IviDeltaCB = IviDeltaCB{ quad_radix: 12, data: &[ | |
| 1029 | 0, 0, 7, 7, -7, -7, -3, 10, | |
| 1030 | 3, -10, 10, -3, -10, 3, 16, 16, | |
| 1031 | -16, -16, 5, 18, -5, -18, 18, 5, | |
| 1032 | -18, -5, -13, 13, 13, -13, -6, 21, | |
| 1033 | 6, -21, 21, -6, -21, 6, 15, 30, | |
| 1034 | -15, -30, 30, 15, -30, -15, 30, 30, | |
| 1035 | -30, -30, 2, 34, -2, -34, 34, 2, | |
| 1036 | -34, -2, -19, 28, 19, -28, 28, -19, | |
| 1037 | -28, 19, 32, 54, -32, -54, 54, 32, | |
| 1038 | -54, -32, -11, 41, 11, -41, 41, -11, | |
| 1039 | -41, 11, 15, 55, -15, -55, 55, 15, | |
| 1040 | -55, -15, 57, 57, -57, -57, -1, 63, | |
| 1041 | 1, -63, 63, -1, -63, 1, -40, 40, | |
| 1042 | 40, -40, -30, 55, 30, -55, 55, -30, | |
| 1043 | -55, 30, 37, 94, -37, -94, 94, 37, | |
| 1044 | -94, -37, 65, 96, -65, -96, 96, 65, | |
| 1045 | -96, -65, -21, 78, 21, -78, 78, -21, | |
| 1046 | -78, 21, 15, 100, -15, -100, 100, 15, | |
| 1047 | -100, -15, 106, 106, -106, -106, -8, 116, | |
| 1048 | 8, -116, 116, -8, -116, 8, -63, 79, | |
| 1049 | 63, -79, 79, -63, -79, 63, -51, 105, | |
| 1050 | 51, -105, 105, -51, -105, 51, -120, 120, | |
| 1051 | 120, -120, 43, 43, -43, -43, 80, 80, | |
| 1052 | -80, -80, | |
| 1053 | ]}; | |
| 1054 | ||
| 1055 | const DT_1_7: IviDeltaCB = IviDeltaCB{ quad_radix: 12, data: &[ | |
| 1056 | 0, 0, 8, 8, -8, -8, -3, 11, | |
| 1057 | 3, -11, 11, -3, -11, 3, 18, 18, | |
| 1058 | -18, -18, 5, 20, -5, -20, 20, 5, | |
| 1059 | -20, -5, -15, 15, 15, -15, -7, 24, | |
| 1060 | 7, -24, 24, -7, -24, 7, 17, 35, | |
| 1061 | -17, -35, 35, 17, -35, -17, 34, 34, | |
| 1062 | -34, -34, 3, 38, -3, -38, 38, 3, | |
| 1063 | -38, -3, -22, 32, 22, -32, 32, -22, | |
| 1064 | -32, 22, 37, 61, -37, -61, 61, 37, | |
| 1065 | -61, -37, -13, 47, 13, -47, 47, -13, | |
| 1066 | -47, 13, 17, 63, -17, -63, 63, 17, | |
| 1067 | -63, -17, 65, 65, -65, -65, -1, 72, | |
| 1068 | 1, -72, 72, -1, -72, 1, -46, 46, | |
| 1069 | 46, -46, -35, 63, 35, -63, 63, -35, | |
| 1070 | -63, 35, 43, 107, -43, -107, 107, 43, | |
| 1071 | -107, -43, 75, 110, -75, -110, 110, 75, | |
| 1072 | -110, -75, -24, 89, 24, -89, 89, -24, | |
| 1073 | -89, 24, 17, 114, -17, -114, 114, 17, | |
| 1074 | -114, -17, 121, 121, -121, -121, -72, 91, | |
| 1075 | 72, -91, 91, -72, -91, 72, -58, 120, | |
| 1076 | 58, -120, 120, -58, -120, 58, 49, 49, | |
| 1077 | -49, -49, 92, 92, -92, -92, | |
| 1078 | ]}; | |
| 1079 | ||
| 1080 | const DT_1_8: IviDeltaCB = IviDeltaCB{ quad_radix: 13, data: &[ | |
| 1081 | 0, 0, 9, 9, -9, -9, -3, 12, | |
| 1082 | 3, -12, 12, -3, -12, 3, 20, 20, | |
| 1083 | -20, -20, 6, 23, -6, -23, 23, 6, | |
| 1084 | -23, -6, -17, 17, 17, -17, -7, 27, | |
| 1085 | 7, -27, 27, -7, -27, 7, 19, 39, | |
| 1086 | -19, -39, 39, 19, -39, -19, 39, 39, | |
| 1087 | -39, -39, 3, 43, -3, -43, 43, 3, | |
| 1088 | -43, -3, -24, 36, 24, -36, 36, -24, | |
| 1089 | -36, 24, 42, 69, -42, -69, 69, 42, | |
| 1090 | -69, -42, -14, 53, 14, -53, 53, -14, | |
| 1091 | -53, 14, 19, 71, -19, -71, 71, 19, | |
| 1092 | -71, -19, 73, 73, -73, -73, -2, 80, | |
| 1093 | 2, -80, 80, -2, -80, 2, -52, 52, | |
| 1094 | 52, -52, -39, 70, 39, -70, 70, -39, | |
| 1095 | -70, 39, 48, 121, -48, -121, 121, 48, | |
| 1096 | -121, -48, 84, 124, -84, -124, 124, 84, | |
| 1097 | -124, -84, -27, 100, 27, -100, 100, -27, | |
| 1098 | -100, 27, -81, 102, 81, -102, 102, -81, | |
| 1099 | -102, 81, 55, 55, -55, -55, 104, 104, | |
| 1100 | -104, -104, | |
| 1101 | ]}; | |
| 1102 | ||
| 1103 | const DT_2_1: IviDeltaCB = IviDeltaCB{ quad_radix: 7, data: &[ | |
| 1104 | 0, 0, 2, 2, -2, -2, 0, 2, | |
| 1105 | 0, -2, 2, 0, -2, 0, 4, 4, | |
| 1106 | -4, -4, 0, 4, 0, -4, 4, 0, | |
| 1107 | -4, 0, -4, 4, 4, -4, -2, 6, | |
| 1108 | 2, -6, 6, -2, -6, 2, 4, 8, | |
| 1109 | -4, -8, 8, 4, -8, -4, 8, 8, | |
| 1110 | -8, -8, 0, 10, 0, -10, 10, 0, | |
| 1111 | -10, 0, -4, 8, 4, -8, 8, -4, | |
| 1112 | -8, 4, 8, 14, -8, -14, 14, 8, | |
| 1113 | -14, -8, -2, 12, 2, -12, 12, -2, | |
| 1114 | -12, 2, 4, 16, -4, -16, 16, 4, | |
| 1115 | -16, -4, 16, 16, -16, -16, 0, 18, | |
| 1116 | 0, -18, 18, 0, -18, 0, -12, 12, | |
| 1117 | 12, -12, -8, 16, 8, -16, 16, -8, | |
| 1118 | -16, 8, 10, 26, -10, -26, 26, 10, | |
| 1119 | -26, -10, 18, 28, -18, -28, 28, 18, | |
| 1120 | -28, -18, -6, 22, 6, -22, 22, -6, | |
| 1121 | -22, 6, 4, 28, -4, -28, 28, 4, | |
| 1122 | -28, -4, 30, 30, -30, -30, -2, 32, | |
| 1123 | 2, -32, 32, -2, -32, 2, -18, 22, | |
| 1124 | 18, -22, 22, -18, -22, 18, -14, 30, | |
| 1125 | 14, -30, 30, -14, -30, 14, 22, 46, | |
| 1126 | -22, -46, 46, 22, -46, -22, 12, 46, | |
| 1127 | -12, -46, 46, 12, -46, -12, 34, 48, | |
| 1128 | -34, -48, 48, 34, -48, -34, -10, 40, | |
| 1129 | 10, -40, 40, -10, -40, 10, 4, 50, | |
| 1130 | -4, -50, 50, 4, -50, -4, 54, 54, | |
| 1131 | -54, -54, -34, 34, 34, -34, -28, 42, | |
| 1132 | 28, -42, 42, -28, -42, 28, -6, 60, | |
| 1133 | 6, -60, 60, -6, -60, 6, 26, 76, | |
| 1134 | -26, -76, 76, 26, -76, -26, 42, 76, | |
| 1135 | -42, -76, 76, 42, -76, -42, -24, 54, | |
| 1136 | 24, -54, 54, -24, -54, 24, 14, 78, | |
| 1137 | -14, -78, 78, 14, -78, -14, 62, 82, | |
| 1138 | -62, -82, 82, 62, -82, -62, -20, 74, | |
| 1139 | 20, -74, 74, -20, -74, 20, 2, 88, | |
| 1140 | -2, -88, 88, 2, -88, -2, 92, 92, | |
| 1141 | -92, -92, -52, 60, 52, -60, 60, -52, | |
| 1142 | -60, 52, 52, 118, -52, -118, 118, 52, | |
| 1143 | -118, -52, -44, 74, 44, -74, 74, -44, | |
| 1144 | -74, 44, 74, 118, -74, -118, 118, 74, | |
| 1145 | -118, -74, 32, 118, -32, -118, 118, 32, | |
| 1146 | -118, -32, -12, 102, 12, -102, 102, -12, | |
| 1147 | -102, 12, -40, 96, 40, -96, 96, -40, | |
| 1148 | -96, 40, -34, 118, 34, -118, 118, -34, | |
| 1149 | -118, 34, -88, 88, 88, -88, -78, 104, | |
| 1150 | 78, -104, 104, -78, -104, 78, 12, 12, | |
| 1151 | -12, -12, 22, 22, -22, -22, 42, 42, | |
| 1152 | -42, -42, 72, 72, -72, -72, | |
| 1153 | ]}; | |
| 1154 | ||
| 1155 | const DT_2_2: IviDeltaCB = IviDeltaCB{ quad_radix: 9, data: &[ | |
| 1156 | 0, 0, 3, 3, -3, -3, 0, 3, | |
| 1157 | 0, -3, 3, 0, -3, 0, 6, 6, | |
| 1158 | -6, -6, 3, 9, -3, -9, 9, 3, | |
| 1159 | -9, -3, -3, 9, 3, -9, 9, -3, | |
| 1160 | -9, 3, -6, 6, 6, -6, 6, 12, | |
| 1161 | -6, -12, 12, 6, -12, -6, 12, 12, | |
| 1162 | -12, -12, 0, 15, 0, -15, 15, 0, | |
| 1163 | -15, 0, -9, 12, 9, -12, 12, -9, | |
| 1164 | -12, 9, 15, 24, -15, -24, 24, 15, | |
| 1165 | -24, -15, -6, 18, 6, -18, 18, -6, | |
| 1166 | -18, 6, 6, 24, -6, -24, 24, 6, | |
| 1167 | -24, -6, 24, 24, -24, -24, 0, 27, | |
| 1168 | 0, -27, 27, 0, -27, 0, -18, 18, | |
| 1169 | 18, -18, -12, 24, 12, -24, 24, -12, | |
| 1170 | -24, 12, 15, 39, -15, -39, 39, 15, | |
| 1171 | -39, -15, 27, 42, -27, -42, 42, 27, | |
| 1172 | -42, -27, -9, 33, 9, -33, 33, -9, | |
| 1173 | -33, 9, 6, 42, -6, -42, 42, 6, | |
| 1174 | -42, -6, 45, 45, -45, -45, -3, 51, | |
| 1175 | 3, -51, 51, -3, -51, 3, -27, 33, | |
| 1176 | 27, -33, 33, -27, -33, 27, -21, 45, | |
| 1177 | 21, -45, 45, -21, -45, 21, 33, 69, | |
| 1178 | -33, -69, 69, 33, -69, -33, 18, 69, | |
| 1179 | -18, -69, 69, 18, -69, -18, 54, 72, | |
| 1180 | -54, -72, 72, 54, -72, -54, -18, 63, | |
| 1181 | 18, -63, 63, -18, -63, 18, 6, 78, | |
| 1182 | -6, -78, 78, 6, -78, -6, 81, 81, | |
| 1183 | -81, -81, -51, 51, 51, -51, -42, 63, | |
| 1184 | 42, -63, 63, -42, -63, 42, -9, 90, | |
| 1185 | 9, -90, 90, -9, -90, 9, 42, 114, | |
| 1186 | -42, -114, 114, 42, -114, -42, 63, 117, | |
| 1187 | -63, -117, 117, 63, -117, -63, -36, 81, | |
| 1188 | 36, -81, 81, -36, -81, 36, 21, 120, | |
| 1189 | -21, -120, 120, 21, -120, -21, 96, 123, | |
| 1190 | -96, -123, 123, 96, -123, -96, -30, 111, | |
| 1191 | 30, -111, 111, -30, -111, 30, -78, 93, | |
| 1192 | 78, -93, 93, -78, -93, 78, -69, 114, | |
| 1193 | 69, -114, 114, -69, -114, 69, 18, 18, | |
| 1194 | -18, -18, 33, 33, -33, -33, 63, 63, | |
| 1195 | -63, -63, 108, 108, -108, -108, | |
| 1196 | ]}; | |
| 1197 | ||
| 1198 | const DT_2_3: IviDeltaCB = IviDeltaCB{ quad_radix: 10, data: &[ | |
| 1199 | 0, 0, 4, 4, -4, -4, 0, 4, | |
| 1200 | 0, -4, 4, 0, -4, 0, 4, 8, | |
| 1201 | -4, -8, 8, 4, -8, -4, 8, 8, | |
| 1202 | -8, -8, -8, 8, 8, -8, -4, 12, | |
| 1203 | 4, -12, 12, -4, -12, 4, 8, 16, | |
| 1204 | -8, -16, 16, 8, -16, -8, 16, 16, | |
| 1205 | -16, -16, 0, 20, 0, -20, 20, 0, | |
| 1206 | -20, 0, -12, 16, 12, -16, 16, -12, | |
| 1207 | -16, 12, -4, 24, 4, -24, 24, -4, | |
| 1208 | -24, 4, 16, 32, -16, -32, 32, 16, | |
| 1209 | -32, -16, 8, 32, -8, -32, 32, 8, | |
| 1210 | -32, -8, 32, 32, -32, -32, 0, 36, | |
| 1211 | 0, -36, 36, 0, -36, 0, -24, 24, | |
| 1212 | 24, -24, -16, 32, 16, -32, 32, -16, | |
| 1213 | -32, 16, 20, 52, -20, -52, 52, 20, | |
| 1214 | -52, -20, 36, 56, -36, -56, 56, 36, | |
| 1215 | -56, -36, -12, 44, 12, -44, 44, -12, | |
| 1216 | -44, 12, 8, 56, -8, -56, 56, 8, | |
| 1217 | -56, -8, 60, 60, -60, -60, -4, 64, | |
| 1218 | 4, -64, 64, -4, -64, 4, -36, 44, | |
| 1219 | 36, -44, 44, -36, -44, 36, -28, 60, | |
| 1220 | 28, -60, 60, -28, -60, 28, 44, 92, | |
| 1221 | -44, -92, 92, 44, -92, -44, 24, 92, | |
| 1222 | -24, -92, 92, 24, -92, -24, 72, 96, | |
| 1223 | -72, -96, 96, 72, -96, -72, -20, 84, | |
| 1224 | 20, -84, 84, -20, -84, 20, 8, 100, | |
| 1225 | -8, -100, 100, 8, -100, -8, 108, 108, | |
| 1226 | -108, -108, -68, 68, 68, -68, -56, 84, | |
| 1227 | 56, -84, 84, -56, -84, 56, -12, 120, | |
| 1228 | 12, -120, 120, -12, -120, 12, -48, 108, | |
| 1229 | 48, -108, 108, -48, -108, 48, -104, 124, | |
| 1230 | 104, -124, 124, -104, -124, 104, 24, 24, | |
| 1231 | -24, -24, 44, 44, -44, -44, 84, 84, | |
| 1232 | -84, -84, | |
| 1233 | ]}; | |
| 1234 | ||
| 1235 | const DT_2_4: IviDeltaCB = IviDeltaCB{ quad_radix: 11, data: &[ | |
| 1236 | 0, 0, 5, 5, -5, -5, 0, 5, | |
| 1237 | 0, -5, 5, 0, -5, 0, 10, 10, | |
| 1238 | -10, -10, 5, 15, -5, -15, 15, 5, | |
| 1239 | -15, -5, -10, 10, 10, -10, -5, 15, | |
| 1240 | 5, -15, 15, -5, -15, 5, 10, 20, | |
| 1241 | -10, -20, 20, 10, -20, -10, 20, 20, | |
| 1242 | -20, -20, 0, 25, 0, -25, 25, 0, | |
| 1243 | -25, 0, -15, 20, 15, -20, 20, -15, | |
| 1244 | -20, 15, 25, 40, -25, -40, 40, 25, | |
| 1245 | -40, -25, -10, 30, 10, -30, 30, -10, | |
| 1246 | -30, 10, 10, 40, -10, -40, 40, 10, | |
| 1247 | -40, -10, 40, 40, -40, -40, 0, 45, | |
| 1248 | 0, -45, 45, 0, -45, 0, -30, 30, | |
| 1249 | 30, -30, -20, 40, 20, -40, 40, -20, | |
| 1250 | -40, 20, 25, 65, -25, -65, 65, 25, | |
| 1251 | -65, -25, 45, 70, -45, -70, 70, 45, | |
| 1252 | -70, -45, -15, 55, 15, -55, 55, -15, | |
| 1253 | -55, 15, 10, 70, -10, -70, 70, 10, | |
| 1254 | -70, -10, 75, 75, -75, -75, -5, 85, | |
| 1255 | 5, -85, 85, -5, -85, 5, -45, 55, | |
| 1256 | 45, -55, 55, -45, -55, 45, -35, 75, | |
| 1257 | 35, -75, 75, -35, -75, 35, 55, 115, | |
| 1258 | -55, -115, 115, 55, -115, -55, 30, 115, | |
| 1259 | -30, -115, 115, 30, -115, -30, 90, 120, | |
| 1260 | -90, -120, 120, 90, -120, -90, -30, 105, | |
| 1261 | 30, -105, 105, -30, -105, 30, -85, 85, | |
| 1262 | 85, -85, -70, 105, 70, -105, 105, -70, | |
| 1263 | -105, 70, 30, 30, -30, -30, 60, 60, | |
| 1264 | -60, -60, 105, 105, -105, -105, | |
| 1265 | ]}; | |
| 1266 | ||
| 1267 | const DT_2_5: IviDeltaCB = IviDeltaCB{ quad_radix: 12, data: &[ | |
| 1268 | 0, 0, 6, 6, -6, -6, 0, 6, | |
| 1269 | 0, -6, 6, 0, -6, 0, 12, 12, | |
| 1270 | -12, -12, 6, 12, -6, -12, 12, 6, | |
| 1271 | -12, -6, -12, 12, 12, -12, -6, 18, | |
| 1272 | 6, -18, 18, -6, -18, 6, 12, 24, | |
| 1273 | -12, -24, 24, 12, -24, -12, 24, 24, | |
| 1274 | -24, -24, 0, 30, 0, -30, 30, 0, | |
| 1275 | -30, 0, -18, 24, 18, -24, 24, -18, | |
| 1276 | -24, 18, 30, 48, -30, -48, 48, 30, | |
| 1277 | -48, -30, -6, 36, 6, -36, 36, -6, | |
| 1278 | -36, 6, 12, 48, -12, -48, 48, 12, | |
| 1279 | -48, -12, 48, 48, -48, -48, 0, 54, | |
| 1280 | 0, -54, 54, 0, -54, 0, -36, 36, | |
| 1281 | 36, -36, -24, 48, 24, -48, 48, -24, | |
| 1282 | -48, 24, 30, 78, -30, -78, 78, 30, | |
| 1283 | -78, -30, 54, 84, -54, -84, 84, 54, | |
| 1284 | -84, -54, -18, 66, 18, -66, 66, -18, | |
| 1285 | -66, 18, 12, 84, -12, -84, 84, 12, | |
| 1286 | -84, -12, 90, 90, -90, -90, -6, 96, | |
| 1287 | 6, -96, 96, -6, -96, 6, -54, 66, | |
| 1288 | 54, -66, 66, -54, -66, 54, -42, 90, | |
| 1289 | 42, -90, 90, -42, -90, 42, -30, 126, | |
| 1290 | 30, -126, 126, -30, -126, 30, -102, 102, | |
| 1291 | 102, -102, -84, 126, 84, -126, 126, -84, | |
| 1292 | -126, 84, 36, 36, -36, -36, 66, 66, | |
| 1293 | -66, -66, | |
| 1294 | ]}; | |
| 1295 | ||
| 1296 | const DT_2_6: IviDeltaCB = IviDeltaCB{ quad_radix: 12, data: &[ | |
| 1297 | 0, 0, 7, 7, -7, -7, 0, 7, | |
| 1298 | 0, -7, 7, 0, -7, 0, 14, 14, | |
| 1299 | -14, -14, 7, 21, -7, -21, 21, 7, | |
| 1300 | -21, -7, -14, 14, 14, -14, -7, 21, | |
| 1301 | 7, -21, 21, -7, -21, 7, 14, 28, | |
| 1302 | -14, -28, 28, 14, -28, -14, 28, 28, | |
| 1303 | -28, -28, 0, 35, 0, -35, 35, 0, | |
| 1304 | -35, 0, -21, 28, 21, -28, 28, -21, | |
| 1305 | -28, 21, 35, 56, -35, -56, 56, 35, | |
| 1306 | -56, -35, -14, 42, 14, -42, 42, -14, | |
| 1307 | -42, 14, 14, 56, -14, -56, 56, 14, | |
| 1308 | -56, -14, 56, 56, -56, -56, 0, 63, | |
| 1309 | 0, -63, 63, 0, -63, 0, -42, 42, | |
| 1310 | 42, -42, -28, 56, 28, -56, 56, -28, | |
| 1311 | -56, 28, 35, 91, -35, -91, 91, 35, | |
| 1312 | -91, -35, 63, 98, -63, -98, 98, 63, | |
| 1313 | -98, -63, -21, 77, 21, -77, 77, -21, | |
| 1314 | -77, 21, 14, 98, -14, -98, 98, 14, | |
| 1315 | -98, -14, 105, 105, -105, -105, -7, 119, | |
| 1316 | 7, -119, 119, -7, -119, 7, -63, 77, | |
| 1317 | 63, -77, 77, -63, -77, 63, -49, 105, | |
| 1318 | 49, -105, 105, -49, -105, 49, -119, 119, | |
| 1319 | 119, -119, 42, 42, -42, -42, 77, 77, | |
| 1320 | -77, -77, | |
| 1321 | ]}; | |
| 1322 | ||
| 1323 | const DT_2_7: IviDeltaCB = IviDeltaCB{ quad_radix: 12, data: &[ | |
| 1324 | 0, 0, 8, 8, -8, -8, 0, 8, | |
| 1325 | 0, -8, 8, 0, -8, 0, 16, 16, | |
| 1326 | -16, -16, 8, 16, -8, -16, 16, 8, | |
| 1327 | -16, -8, -16, 16, 16, -16, -8, 24, | |
| 1328 | 8, -24, 24, -8, -24, 8, 16, 32, | |
| 1329 | -16, -32, 32, 16, -32, -16, 32, 32, | |
| 1330 | -32, -32, 0, 40, 0, -40, 40, 0, | |
| 1331 | -40, 0, -24, 32, 24, -32, 32, -24, | |
| 1332 | -32, 24, 40, 64, -40, -64, 64, 40, | |
| 1333 | -64, -40, -16, 48, 16, -48, 48, -16, | |
| 1334 | -48, 16, 16, 64, -16, -64, 64, 16, | |
| 1335 | -64, -16, 64, 64, -64, -64, 0, 72, | |
| 1336 | 0, -72, 72, 0, -72, 0, -48, 48, | |
| 1337 | 48, -48, -32, 64, 32, -64, 64, -32, | |
| 1338 | -64, 32, 40, 104, -40, -104, 104, 40, | |
| 1339 | -104, -40, 72, 112, -72, -112, 112, 72, | |
| 1340 | -112, -72, -24, 88, 24, -88, 88, -24, | |
| 1341 | -88, 24, 16, 112, -16, -112, 112, 16, | |
| 1342 | -112, -16, 120, 120, -120, -120, -72, 88, | |
| 1343 | 72, -88, 88, -72, -88, 72, -56, 120, | |
| 1344 | 56, -120, 120, -56, -120, 56, 48, 48, | |
| 1345 | -48, -48, 88, 88, -88, -88, | |
| 1346 | ]}; | |
| 1347 | ||
| 1348 | const DT_2_8: IviDeltaCB = IviDeltaCB{ quad_radix: 13, data: &[ | |
| 1349 | 0, 0, 9, 9, -9, -9, 0, 9, | |
| 1350 | 0, -9, 9, 0, -9, 0, 18, 18, | |
| 1351 | -18, -18, 9, 27, -9, -27, 27, 9, | |
| 1352 | -27, -9, -18, 18, 18, -18, -9, 27, | |
| 1353 | 9, -27, 27, -9, -27, 9, 18, 36, | |
| 1354 | -18, -36, 36, 18, -36, -18, 36, 36, | |
| 1355 | -36, -36, 0, 45, 0, -45, 45, 0, | |
| 1356 | -45, 0, -27, 36, 27, -36, 36, -27, | |
| 1357 | -36, 27, 45, 72, -45, -72, 72, 45, | |
| 1358 | -72, -45, -18, 54, 18, -54, 54, -18, | |
| 1359 | -54, 18, 18, 72, -18, -72, 72, 18, | |
| 1360 | -72, -18, 72, 72, -72, -72, 0, 81, | |
| 1361 | 0, -81, 81, 0, -81, 0, -54, 54, | |
| 1362 | 54, -54, -36, 72, 36, -72, 72, -36, | |
| 1363 | -72, 36, 45, 117, -45, -117, 117, 45, | |
| 1364 | -117, -45, 81, 126, -81, -126, 126, 81, | |
| 1365 | -126, -81, -27, 99, 27, -99, 99, -27, | |
| 1366 | -99, 27, -81, 99, 81, -99, 99, -81, | |
| 1367 | -99, 81, 54, 54, -54, -54, 108, 108, | |
| 1368 | -108, -108, | |
| 1369 | ]}; | |
| 1370 | ||
| 1371 | const DT_3_1: IviDeltaCB = IviDeltaCB{ quad_radix: 11, data: &[ | |
| 1372 | 0, 0, 2, 2, -2, -2, 0, 3, | |
| 1373 | 0, -3, 3, 0, -3, 0, 6, 6, | |
| 1374 | -6, -6, 0, 7, 0, -7, 7, 0, | |
| 1375 | -7, 0, -5, 5, 5, -5, 5, -5, | |
| 1376 | -5, 5, 6, 11, -6, -11, 11, 6, | |
| 1377 | -11, -6, 0, 8, 0, -8, 8, 0, | |
| 1378 | -8, 0, 11, 11, -11, -11, 0, 12, | |
| 1379 | 0, -12, 12, 0, -12, 0, 12, 17, | |
| 1380 | -12, -17, 17, 12, -17, -12, 17, 17, | |
| 1381 | -17, -17, 6, 18, -6, -18, 18, 6, | |
| 1382 | -18, -6, -8, 11, 8, -11, 11, -8, | |
| 1383 | -11, 8, 0, 15, 0, -15, 15, 0, | |
| 1384 | -15, 0, 0, 20, 0, -20, 20, 0, | |
| 1385 | -20, 0, 18, 25, -18, -25, 25, 18, | |
| 1386 | -25, -18, 11, 25, -11, -25, 25, 11, | |
| 1387 | -25, -11, 25, 25, -25, -25, -14, 14, | |
| 1388 | 14, -14, 14, -14, -14, 14, 0, 26, | |
| 1389 | 0, -26, 26, 0, -26, 0, -11, 18, | |
| 1390 | 11, -18, 18, -11, -18, 11, -7, 22, | |
| 1391 | 7, -22, 22, -7, -22, 7, 26, 34, | |
| 1392 | -26, -34, 34, 26, -34, -26, 18, 34, | |
| 1393 | -18, -34, 34, 18, -34, -18, 34, 34, | |
| 1394 | -34, -34, 11, 35, -11, -35, 35, 11, | |
| 1395 | -35, -11, 0, 29, 0, -29, 29, 0, | |
| 1396 | -29, 0, -19, 22, 19, -22, 22, -19, | |
| 1397 | -22, 19, -15, 26, 15, -26, 26, -15, | |
| 1398 | -26, 15, 0, 37, 0, -37, 37, 0, | |
| 1399 | -37, 0, 27, 44, -27, -44, 44, 27, | |
| 1400 | -44, -27, 36, 44, -36, -44, 44, 36, | |
| 1401 | -44, -36, 18, 44, -18, -44, 44, 18, | |
| 1402 | -44, -18, -10, 33, 10, -33, 33, -10, | |
| 1403 | -33, 10, 45, 45, -45, -45, 0, 0, | |
| 1404 | ]}; | |
| 1405 | ||
| 1406 | const DT_3_2: IviDeltaCB = IviDeltaCB{ quad_radix: 13, data: &[ | |
| 1407 | 0, 0, 0, 2, 0, -2, 2, 0, | |
| 1408 | -2, 0, 2, 2, -2, -2, 6, 6, | |
| 1409 | -6, -6, 0, 6, 0, -6, 6, 0, | |
| 1410 | -6, 0, -4, 4, 4, -4, 10, -6, | |
| 1411 | -10, 6, 0, -12, 0, 12, -6, -12, | |
| 1412 | 6, -12, -6, 12, 6, 12, -14, 0, | |
| 1413 | 14, 0, 12, 12, -12, -12, 0, -18, | |
| 1414 | 0, 18, 14, -12, -14, 12, -18, -6, | |
| 1415 | 18, -6, -18, 6, 18, 6, -10, -18, | |
| 1416 | 10, -18, -10, 18, 10, 18, -22, 0, | |
| 1417 | 22, 0, 0, -24, 0, 24, -22, -12, | |
| 1418 | 22, -12, -22, 12, 22, 12, -8, -24, | |
| 1419 | 8, -24, -8, 24, 8, 24, -26, -6, | |
| 1420 | 26, -6, -26, 6, 26, 6, -28, 0, | |
| 1421 | 28, 0, 20, 20, -20, -20, -14, -26, | |
| 1422 | 14, 26, -30, -12, 30, 12, -10, -32, | |
| 1423 | 10, 32, -18, -32, 18, 32, -26, -26, | |
| 1424 | 26, 26, -34, -20, 34, 20, -38, -12, | |
| 1425 | 38, 12, -32, -32, 32, 32, 32, 32, | |
| 1426 | -22, -40, -34, -34, 34, 34, | |
| 1427 | ]}; | |
| 1428 | ||
| 1429 | const DT_3_3: IviDeltaCB = IviDeltaCB{ quad_radix: 13, data: &[ | |
| 1430 | 0, 0, 0, 2, 0, -2, 2, 0, | |
| 1431 | -2, 0, 4, 4, -4, -4, 10, 10, | |
| 1432 | -10, -10, 0, 10, 0, -10, 10, 0, | |
| 1433 | -10, 0, -6, 6, 6, -6, 14, -8, | |
| 1434 | -14, 8, -18, 0, 18, 0, 10, -16, | |
| 1435 | -10, 16, 0, -24, 0, 24, -24, -8, | |
| 1436 | 24, -8, -24, 8, 24, 8, 18, 18, | |
| 1437 | -18, -18, 20, -16, -20, 16, -14, -26, | |
| 1438 | 14, -26, -14, 26, 14, 26, -30, 0, | |
| 1439 | 30, 0, 0, -34, 0, 34, -34, -8, | |
| 1440 | 34, -8, -34, 8, 34, 8, -30, -18, | |
| 1441 | 30, -18, -30, 18, 30, 18, -10, -34, | |
| 1442 | 10, -34, -10, 34, 10, 34, -20, -34, | |
| 1443 | 20, 34, -40, 0, 40, 0, 30, 30, | |
| 1444 | -30, -30, -40, -18, 40, 18, 0, -44, | |
| 1445 | 0, 44, -16, -44, 16, 44, -36, -36, | |
| 1446 | -36, -36, 36, 36, -26, -44, 26, 44, | |
| 1447 | -46, -26, 46, 26, -52, -18, 52, 18, | |
| 1448 | -20, -54, -44, -44, 44, 44, -32, -54, | |
| 1449 | -46, -46, -46, -46, 46, 46, | |
| 1450 | ]}; | |
| 1451 | ||
| 1452 | const DT_3_4: IviDeltaCB = IviDeltaCB{ quad_radix: 13, data: &[ | |
| 1453 | 0, 0, 0, 4, 0, -4, 4, 0, | |
| 1454 | -4, 0, 4, 4, -4, -4, 12, 12, | |
| 1455 | -12, -12, 0, 12, 0, -12, 12, 0, | |
| 1456 | -12, 0, -8, 8, 8, -8, 8, -16, | |
| 1457 | -8, 16, 0, -24, 0, 24, -24, -8, | |
| 1458 | 24, -8, -24, 8, 24, 8, 20, -16, | |
| 1459 | -20, 16, -28, 0, 28, 0, -16, -24, | |
| 1460 | 16, -24, -16, 24, 16, 24, 0, -32, | |
| 1461 | 0, 32, -28, -16, 28, -16, -28, 16, | |
| 1462 | 28, 16, -8, -32, 8, -32, -32, -8, | |
| 1463 | 32, -8, -32, 8, 32, 8, -8, 32, | |
| 1464 | 8, 32, 24, 24, -24, -24, 24, -24, | |
| 1465 | -24, 24, -20, -32, 20, 32, -40, 0, | |
| 1466 | 40, 0, -40, -16, 40, 16, 0, -44, | |
| 1467 | 0, -44, -44, 0, 44, 0, 0, 44, | |
| 1468 | 0, 44, -32, -32, 32, 32, -16, -44, | |
| 1469 | 16, 44, -24, -44, -44, -24, 44, 24, | |
| 1470 | 24, 44, -48, -16, 48, 16, -36, -36, | |
| 1471 | -36, -36, 36, 36, 36, 36, -20, -52, | |
| 1472 | 40, 40, -40, -40, -32, -52, | |
| 1473 | ]}; | |
| 1474 | ||
| 1475 | const DT_3_5: IviDeltaCB = IviDeltaCB{ quad_radix: 13, data: &[ | |
| 1476 | 0, 0, 2, 2, -2, -2, 6, 6, | |
| 1477 | -6, -6, 12, 12, -12, -12, 20, 20, | |
| 1478 | -20, -20, 32, 32, -32, -32, 46, 46, | |
| 1479 | -46, -46, 0, 0, 0, 0, 0, 0, | |
| 1480 | 0, 0, 0, 0, 0, 0, 0, 0, | |
| 1481 | 0, 0, 0, 0, 0, 0, 0, 0, | |
| 1482 | 0, 0, 0, 0, 0, 0, 0, 0, | |
| 1483 | 0, 0, 0, 0, 0, 0, 0, 0, | |
| 1484 | 0, 0, 0, 0, 0, 0, 0, 0, | |
| 1485 | 0, 0, 0, 0, 0, 0, 0, 0, | |
| 1486 | 0, 0, 0, 0, 0, 0, 0, 0, | |
| 1487 | 0, 0, 0, 0, 0, 0, 0, 0, | |
| 1488 | 0, 0, 0, 0, 0, 0, 0, 0, | |
| 1489 | 0, 0, 0, 0, 0, 0, 0, 0, | |
| 1490 | 0, 0, 0, 0, 0, 0, 0, 0, | |
| 1491 | 0, 0, 0, 0, 0, 0, 0, 0, | |
| 1492 | 0, 0, 0, 0, 0, 0, 0, 0, | |
| 1493 | 0, 0, 0, 0, 0, 0, 0, 0, | |
| 1494 | 0, 0, 0, 0, 0, 0, 0, 0, | |
| 1495 | 0, 0, 0, 0, 0, 0, | |
| 1496 | ]}; | |
| 1497 | ||
| 1498 | const IVI3_DELTA_CBS: [&IviDeltaCB; 24] = [ | |
| 1499 | &DT_1_1, &DT_1_2, &DT_1_3, &DT_1_4, &DT_1_5, &DT_1_6, &DT_1_7, &DT_1_8, | |
| 1500 | &DT_2_1, &DT_2_2, &DT_2_3, &DT_2_4, &DT_2_5, &DT_2_6, &DT_2_7, &DT_2_8, | |
| 1501 | &DT_3_1, &DT_3_2, &DT_3_3, &DT_3_4, &DT_3_5, &DT_3_5, &DT_3_5, &DT_3_5 | |
| 1502 | ]; |