You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

1280 lines
55KB

  1. /*
  2. * common functions for Indeo Video Interactive codecs (Indeo4 and Indeo5)
  3. *
  4. * Copyright (c) 2009 Maxim Poliakovski
  5. *
  6. * This file is part of Libav.
  7. *
  8. * Libav is free software; you can redistribute it and/or
  9. * modify it under the terms of the GNU Lesser General Public
  10. * License as published by the Free Software Foundation; either
  11. * version 2.1 of the License, or (at your option) any later version.
  12. *
  13. * Libav is distributed in the hope that it will be useful,
  14. * but WITHOUT ANY WARRANTY; without even the implied warranty of
  15. * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
  16. * Lesser General Public License for more details.
  17. *
  18. * You should have received a copy of the GNU Lesser General Public
  19. * License along with Libav; if not, write to the Free Software
  20. * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
  21. */
  22. /**
  23. * @file
  24. * This file contains functions and data shared by both Indeo4 and
  25. * Indeo5 decoders.
  26. */
  27. #define BITSTREAM_READER_LE
  28. #include "libavutil/attributes.h"
  29. #include "avcodec.h"
  30. #include "get_bits.h"
  31. #include "internal.h"
  32. #include "mathops.h"
  33. #include "ivi_common.h"
  34. #include "ivi_dsp.h"
  35. extern const IVIHuffDesc ff_ivi_mb_huff_desc[8]; ///< static macroblock huffman tables
  36. extern const IVIHuffDesc ff_ivi_blk_huff_desc[8]; ///< static block huffman tables
  37. VLC ff_ivi_mb_vlc_tabs [8];
  38. VLC ff_ivi_blk_vlc_tabs[8];
  39. /**
  40. * Reverse "nbits" bits of the value "val" and return the result
  41. * in the least significant bits.
  42. */
  43. static uint16_t inv_bits(uint16_t val, int nbits)
  44. {
  45. uint16_t res;
  46. if (nbits <= 8) {
  47. res = ff_reverse[val] >> (8-nbits);
  48. } else
  49. res = ((ff_reverse[val & 0xFF] << 8) + (ff_reverse[val >> 8])) >> (16-nbits);
  50. return res;
  51. }
  52. int ff_ivi_create_huff_from_desc(const IVIHuffDesc *cb, VLC *vlc, int flag)
  53. {
  54. int pos, i, j, codes_per_row, prefix, not_last_row;
  55. uint16_t codewords[256]; /* FIXME: move this temporal storage out? */
  56. uint8_t bits[256];
  57. pos = 0; /* current position = 0 */
  58. for (i = 0; i < cb->num_rows; i++) {
  59. codes_per_row = 1 << cb->xbits[i];
  60. not_last_row = (i != cb->num_rows - 1);
  61. prefix = ((1 << i) - 1) << (cb->xbits[i] + not_last_row);
  62. for (j = 0; j < codes_per_row; j++) {
  63. if (pos >= 256) /* Some Indeo5 codebooks can have more than 256 */
  64. break; /* elements, but only 256 codes are allowed! */
  65. bits[pos] = i + cb->xbits[i] + not_last_row;
  66. if (bits[pos] > IVI_VLC_BITS)
  67. return -1; /* invalid descriptor */
  68. codewords[pos] = inv_bits((prefix | j), bits[pos]);
  69. if (!bits[pos])
  70. bits[pos] = 1;
  71. pos++;
  72. }//for j
  73. }//for i
  74. /* number of codewords = pos */
  75. return init_vlc(vlc, IVI_VLC_BITS, pos, bits, 1, 1, codewords, 2, 2,
  76. (flag ? INIT_VLC_USE_NEW_STATIC : 0) | INIT_VLC_LE);
  77. }
  78. void ff_ivi_init_static_vlc(void)
  79. {
  80. int i;
  81. static VLC_TYPE table_data[8192 * 16][2];
  82. static int initialized_vlcs = 0;
  83. if (initialized_vlcs)
  84. return;
  85. for (i = 0; i < 8; i++) {
  86. ff_ivi_mb_vlc_tabs[i].table = table_data + i * 2 * 8192;
  87. ff_ivi_mb_vlc_tabs[i].table_allocated = 8192;
  88. ff_ivi_create_huff_from_desc(&ff_ivi_mb_huff_desc[i], &ff_ivi_mb_vlc_tabs[i], 1);
  89. ff_ivi_blk_vlc_tabs[i].table = table_data + (i * 2 + 1) * 8192;
  90. ff_ivi_blk_vlc_tabs[i].table_allocated = 8192;
  91. ff_ivi_create_huff_from_desc(&ff_ivi_blk_huff_desc[i], &ff_ivi_blk_vlc_tabs[i], 1);
  92. }
  93. initialized_vlcs = 1;
  94. }
  95. int ff_ivi_dec_huff_desc(GetBitContext *gb, int desc_coded, int which_tab,
  96. IVIHuffTab *huff_tab, AVCodecContext *avctx)
  97. {
  98. int i, result;
  99. IVIHuffDesc new_huff;
  100. if (!desc_coded) {
  101. /* select default table */
  102. huff_tab->tab = (which_tab) ? &ff_ivi_blk_vlc_tabs[7]
  103. : &ff_ivi_mb_vlc_tabs [7];
  104. } else {
  105. huff_tab->tab_sel = get_bits(gb, 3);
  106. if (huff_tab->tab_sel == 7) {
  107. /* custom huffman table (explicitly encoded) */
  108. new_huff.num_rows = get_bits(gb, 4);
  109. if (!new_huff.num_rows) {
  110. av_log(avctx, AV_LOG_ERROR, "Empty custom Huffman table!\n");
  111. return AVERROR_INVALIDDATA;
  112. }
  113. for (i = 0; i < new_huff.num_rows; i++)
  114. new_huff.xbits[i] = get_bits(gb, 4);
  115. /* Have we got the same custom table? Rebuild if not. */
  116. if (ff_ivi_huff_desc_cmp(&new_huff, &huff_tab->cust_desc)) {
  117. ff_ivi_huff_desc_copy(&huff_tab->cust_desc, &new_huff);
  118. if (huff_tab->cust_tab.table)
  119. ff_free_vlc(&huff_tab->cust_tab);
  120. result = ff_ivi_create_huff_from_desc(&huff_tab->cust_desc,
  121. &huff_tab->cust_tab, 0);
  122. if (result) {
  123. huff_tab->cust_desc.num_rows = 0; // reset faulty description
  124. av_log(avctx, AV_LOG_ERROR,
  125. "Error while initializing custom vlc table!\n");
  126. return result;
  127. }
  128. }
  129. huff_tab->tab = &huff_tab->cust_tab;
  130. } else {
  131. /* select one of predefined tables */
  132. huff_tab->tab = (which_tab) ? &ff_ivi_blk_vlc_tabs[huff_tab->tab_sel]
  133. : &ff_ivi_mb_vlc_tabs [huff_tab->tab_sel];
  134. }
  135. }
  136. return 0;
  137. }
  138. int ff_ivi_huff_desc_cmp(const IVIHuffDesc *desc1, const IVIHuffDesc *desc2)
  139. {
  140. return desc1->num_rows != desc2->num_rows
  141. || memcmp(desc1->xbits, desc2->xbits, desc1->num_rows);
  142. }
  143. void ff_ivi_huff_desc_copy(IVIHuffDesc *dst, const IVIHuffDesc *src)
  144. {
  145. dst->num_rows = src->num_rows;
  146. memcpy(dst->xbits, src->xbits, src->num_rows);
  147. }
  148. av_cold int ff_ivi_init_planes(IVIPlaneDesc *planes, const IVIPicConfig *cfg)
  149. {
  150. int p, b;
  151. uint32_t b_width, b_height, align_fac, width_aligned, height_aligned, buf_size;
  152. IVIBandDesc *band;
  153. ff_ivi_free_buffers(planes);
  154. /* fill in the descriptor of the luminance plane */
  155. planes[0].width = cfg->pic_width;
  156. planes[0].height = cfg->pic_height;
  157. planes[0].num_bands = cfg->luma_bands;
  158. /* fill in the descriptors of the chrominance planes */
  159. planes[1].width = planes[2].width = (cfg->pic_width + 3) >> 2;
  160. planes[1].height = planes[2].height = (cfg->pic_height + 3) >> 2;
  161. planes[1].num_bands = planes[2].num_bands = cfg->chroma_bands;
  162. for (p = 0; p < 3; p++) {
  163. planes[p].bands = av_mallocz(planes[p].num_bands * sizeof(IVIBandDesc));
  164. if (!planes[p].bands)
  165. return AVERROR(ENOMEM);
  166. /* select band dimensions: if there is only one band then it
  167. * has the full size, if there are several bands each of them
  168. * has only half size */
  169. b_width = planes[p].num_bands == 1 ? planes[p].width : (planes[p].width + 1) >> 1;
  170. b_height = planes[p].num_bands == 1 ? planes[p].height : (planes[p].height + 1) >> 1;
  171. /* luma band buffers will be aligned on 16x16 (max macroblock size) */
  172. /* chroma band buffers will be aligned on 8x8 (max macroblock size) */
  173. align_fac = p ? 8 : 16;
  174. width_aligned = FFALIGN(b_width , align_fac);
  175. height_aligned = FFALIGN(b_height, align_fac);
  176. buf_size = width_aligned * height_aligned * sizeof(int16_t);
  177. for (b = 0; b < planes[p].num_bands; b++) {
  178. band = &planes[p].bands[b]; /* select appropriate plane/band */
  179. band->plane = p;
  180. band->band_num = b;
  181. band->width = b_width;
  182. band->height = b_height;
  183. band->pitch = width_aligned;
  184. band->aheight = height_aligned;
  185. band->bufs[0] = av_mallocz(buf_size);
  186. band->bufs[1] = av_mallocz(buf_size);
  187. if (!band->bufs[0] || !band->bufs[1])
  188. return AVERROR(ENOMEM);
  189. /* allocate the 3rd band buffer for scalability mode */
  190. if (cfg->luma_bands > 1) {
  191. band->bufs[2] = av_mallocz(buf_size);
  192. if (!band->bufs[2])
  193. return AVERROR(ENOMEM);
  194. }
  195. planes[p].bands[0].blk_vlc.cust_desc.num_rows = 0; /* reset custom vlc */
  196. }
  197. }
  198. return 0;
  199. }
  200. av_cold void ff_ivi_free_buffers(IVIPlaneDesc *planes)
  201. {
  202. int p, b, t;
  203. for (p = 0; p < 3; p++) {
  204. for (b = 0; b < planes[p].num_bands; b++) {
  205. av_freep(&planes[p].bands[b].bufs[0]);
  206. av_freep(&planes[p].bands[b].bufs[1]);
  207. av_freep(&planes[p].bands[b].bufs[2]);
  208. if (planes[p].bands[b].blk_vlc.cust_tab.table)
  209. ff_free_vlc(&planes[p].bands[b].blk_vlc.cust_tab);
  210. for (t = 0; t < planes[p].bands[b].num_tiles; t++)
  211. av_freep(&planes[p].bands[b].tiles[t].mbs);
  212. av_freep(&planes[p].bands[b].tiles);
  213. }
  214. av_freep(&planes[p].bands);
  215. }
  216. }
  217. av_cold int ff_ivi_init_tiles(IVIPlaneDesc *planes, int tile_width, int tile_height)
  218. {
  219. int p, b, x, y, x_tiles, y_tiles, t_width, t_height;
  220. IVIBandDesc *band;
  221. IVITile *tile, *ref_tile;
  222. for (p = 0; p < 3; p++) {
  223. t_width = !p ? tile_width : (tile_width + 3) >> 2;
  224. t_height = !p ? tile_height : (tile_height + 3) >> 2;
  225. if (!p && planes[0].num_bands == 4) {
  226. t_width >>= 1;
  227. t_height >>= 1;
  228. }
  229. for (b = 0; b < planes[p].num_bands; b++) {
  230. band = &planes[p].bands[b];
  231. x_tiles = IVI_NUM_TILES(band->width, t_width);
  232. y_tiles = IVI_NUM_TILES(band->height, t_height);
  233. band->num_tiles = x_tiles * y_tiles;
  234. av_freep(&band->tiles);
  235. band->tiles = av_mallocz(band->num_tiles * sizeof(IVITile));
  236. if (!band->tiles)
  237. return AVERROR(ENOMEM);
  238. tile = band->tiles;
  239. /* use the first luma band as reference for motion vectors
  240. * and quant */
  241. ref_tile = planes[0].bands[0].tiles;
  242. for (y = 0; y < band->height; y += t_height) {
  243. for (x = 0; x < band->width; x += t_width) {
  244. tile->xpos = x;
  245. tile->ypos = y;
  246. tile->mb_size = band->mb_size;
  247. tile->width = FFMIN(band->width - x, t_width);
  248. tile->height = FFMIN(band->height - y, t_height);
  249. tile->is_empty = tile->data_size = 0;
  250. /* calculate number of macroblocks */
  251. tile->num_MBs = IVI_MBs_PER_TILE(tile->width, tile->height,
  252. band->mb_size);
  253. av_freep(&tile->mbs);
  254. tile->mbs = av_malloc(tile->num_MBs * sizeof(IVIMbInfo));
  255. if (!tile->mbs)
  256. return AVERROR(ENOMEM);
  257. tile->ref_mbs = 0;
  258. if (p || b) {
  259. tile->ref_mbs = ref_tile->mbs;
  260. ref_tile++;
  261. }
  262. tile++;
  263. }
  264. }
  265. }// for b
  266. }// for p
  267. return 0;
  268. }
  269. int ff_ivi_dec_tile_data_size(GetBitContext *gb)
  270. {
  271. int len;
  272. len = 0;
  273. if (get_bits1(gb)) {
  274. len = get_bits(gb, 8);
  275. if (len == 255)
  276. len = get_bits_long(gb, 24);
  277. }
  278. /* align the bitstream reader on the byte boundary */
  279. align_get_bits(gb);
  280. return len;
  281. }
  282. int ff_ivi_decode_blocks(GetBitContext *gb, IVIBandDesc *band, IVITile *tile)
  283. {
  284. int mbn, blk, num_blocks, num_coeffs, blk_size, scan_pos, run, val,
  285. pos, is_intra, mc_type = 0, mv_x, mv_y, col_mask;
  286. uint8_t col_flags[8];
  287. int32_t prev_dc, trvec[64];
  288. uint32_t cbp, sym, lo, hi, quant, buf_offs, q;
  289. IVIMbInfo *mb;
  290. RVMapDesc *rvmap = band->rv_map;
  291. void (*mc_with_delta_func)(int16_t *buf, const int16_t *ref_buf, uint32_t pitch, int mc_type);
  292. void (*mc_no_delta_func) (int16_t *buf, const int16_t *ref_buf, uint32_t pitch, int mc_type);
  293. const uint16_t *base_tab;
  294. const uint8_t *scale_tab;
  295. prev_dc = 0; /* init intra prediction for the DC coefficient */
  296. blk_size = band->blk_size;
  297. col_mask = blk_size - 1; /* column mask for tracking non-zero coeffs */
  298. num_blocks = (band->mb_size != blk_size) ? 4 : 1; /* number of blocks per mb */
  299. num_coeffs = blk_size * blk_size;
  300. if (blk_size == 8) {
  301. mc_with_delta_func = ff_ivi_mc_8x8_delta;
  302. mc_no_delta_func = ff_ivi_mc_8x8_no_delta;
  303. } else {
  304. mc_with_delta_func = ff_ivi_mc_4x4_delta;
  305. mc_no_delta_func = ff_ivi_mc_4x4_no_delta;
  306. }
  307. for (mbn = 0, mb = tile->mbs; mbn < tile->num_MBs; mb++, mbn++) {
  308. is_intra = !mb->type;
  309. cbp = mb->cbp;
  310. buf_offs = mb->buf_offs;
  311. quant = av_clip(band->glob_quant + mb->q_delta, 0, 23);
  312. base_tab = is_intra ? band->intra_base : band->inter_base;
  313. scale_tab = is_intra ? band->intra_scale : band->inter_scale;
  314. if (scale_tab)
  315. quant = scale_tab[quant];
  316. if (!is_intra) {
  317. mv_x = mb->mv_x;
  318. mv_y = mb->mv_y;
  319. if (band->is_halfpel) {
  320. mc_type = ((mv_y & 1) << 1) | (mv_x & 1);
  321. mv_x >>= 1;
  322. mv_y >>= 1; /* convert halfpel vectors into fullpel ones */
  323. }
  324. if (mb->type) {
  325. int dmv_x, dmv_y, cx, cy;
  326. dmv_x = mb->mv_x >> band->is_halfpel;
  327. dmv_y = mb->mv_y >> band->is_halfpel;
  328. cx = mb->mv_x & band->is_halfpel;
  329. cy = mb->mv_y & band->is_halfpel;
  330. if ( mb->xpos + dmv_x < 0
  331. || mb->xpos + dmv_x + band->mb_size + cx > band->pitch
  332. || mb->ypos + dmv_y < 0
  333. || mb->ypos + dmv_y + band->mb_size + cy > band->aheight) {
  334. return AVERROR_INVALIDDATA;
  335. }
  336. }
  337. }
  338. for (blk = 0; blk < num_blocks; blk++) {
  339. /* adjust block position in the buffer according to its number */
  340. if (blk & 1) {
  341. buf_offs += blk_size;
  342. } else if (blk == 2) {
  343. buf_offs -= blk_size;
  344. buf_offs += blk_size * band->pitch;
  345. }
  346. if (cbp & 1) { /* block coded ? */
  347. scan_pos = -1;
  348. memset(trvec, 0, num_coeffs*sizeof(trvec[0])); /* zero transform vector */
  349. memset(col_flags, 0, sizeof(col_flags)); /* zero column flags */
  350. while (scan_pos <= num_coeffs) {
  351. sym = get_vlc2(gb, band->blk_vlc.tab->table, IVI_VLC_BITS, 1);
  352. if (sym == rvmap->eob_sym)
  353. break; /* End of block */
  354. if (sym == rvmap->esc_sym) { /* Escape - run/val explicitly coded using 3 vlc codes */
  355. run = get_vlc2(gb, band->blk_vlc.tab->table, IVI_VLC_BITS, 1) + 1;
  356. lo = get_vlc2(gb, band->blk_vlc.tab->table, IVI_VLC_BITS, 1);
  357. hi = get_vlc2(gb, band->blk_vlc.tab->table, IVI_VLC_BITS, 1);
  358. val = IVI_TOSIGNED((hi << 6) | lo); /* merge them and convert into signed val */
  359. } else {
  360. if (sym >= 256U) {
  361. av_log(NULL, AV_LOG_ERROR, "Invalid sym encountered: %d.\n", sym);
  362. return -1;
  363. }
  364. run = rvmap->runtab[sym];
  365. val = rvmap->valtab[sym];
  366. }
  367. /* de-zigzag and dequantize */
  368. scan_pos += run;
  369. if (scan_pos >= num_coeffs)
  370. break;
  371. pos = band->scan[scan_pos];
  372. if (!val)
  373. av_dlog(NULL, "Val = 0 encountered!\n");
  374. q = (base_tab[pos] * quant) >> 9;
  375. if (q > 1)
  376. val = val * q + FFSIGN(val) * (((q ^ 1) - 1) >> 1);
  377. trvec[pos] = val;
  378. col_flags[pos & col_mask] |= !!val; /* track columns containing non-zero coeffs */
  379. }// while
  380. if (scan_pos >= num_coeffs && sym != rvmap->eob_sym)
  381. return -1; /* corrupt block data */
  382. /* undoing DC coeff prediction for intra-blocks */
  383. if (is_intra && band->is_2d_trans) {
  384. prev_dc += trvec[0];
  385. trvec[0] = prev_dc;
  386. col_flags[0] |= !!prev_dc;
  387. }
  388. /* apply inverse transform */
  389. band->inv_transform(trvec, band->buf + buf_offs,
  390. band->pitch, col_flags);
  391. /* apply motion compensation */
  392. if (!is_intra)
  393. mc_with_delta_func(band->buf + buf_offs,
  394. band->ref_buf + buf_offs + mv_y * band->pitch + mv_x,
  395. band->pitch, mc_type);
  396. } else {
  397. /* block not coded */
  398. /* for intra blocks apply the dc slant transform */
  399. /* for inter - perform the motion compensation without delta */
  400. if (is_intra && band->dc_transform) {
  401. band->dc_transform(&prev_dc, band->buf + buf_offs,
  402. band->pitch, blk_size);
  403. } else
  404. mc_no_delta_func(band->buf + buf_offs,
  405. band->ref_buf + buf_offs + mv_y * band->pitch + mv_x,
  406. band->pitch, mc_type);
  407. }
  408. cbp >>= 1;
  409. }// for blk
  410. }// for mbn
  411. align_get_bits(gb);
  412. return 0;
  413. }
  414. /**
  415. * Handle empty tiles by performing data copying and motion
  416. * compensation respectively.
  417. *
  418. * @param[in] avctx ptr to the AVCodecContext
  419. * @param[in] band pointer to the band descriptor
  420. * @param[in] tile pointer to the tile descriptor
  421. * @param[in] mv_scale scaling factor for motion vectors
  422. */
  423. static int ivi_process_empty_tile(AVCodecContext *avctx, IVIBandDesc *band,
  424. IVITile *tile, int32_t mv_scale)
  425. {
  426. int x, y, need_mc, mbn, blk, num_blocks, mv_x, mv_y, mc_type;
  427. int offs, mb_offset, row_offset;
  428. IVIMbInfo *mb, *ref_mb;
  429. const int16_t *src;
  430. int16_t *dst;
  431. void (*mc_no_delta_func)(int16_t *buf, const int16_t *ref_buf, uint32_t pitch,
  432. int mc_type);
  433. if (tile->num_MBs != IVI_MBs_PER_TILE(tile->width, tile->height, band->mb_size)) {
  434. av_log(avctx, AV_LOG_ERROR, "Allocated tile size %d mismatches "
  435. "parameters %d in ivi_process_empty_tile()\n",
  436. tile->num_MBs, IVI_MBs_PER_TILE(tile->width, tile->height, band->mb_size));
  437. return AVERROR_INVALIDDATA;
  438. }
  439. offs = tile->ypos * band->pitch + tile->xpos;
  440. mb = tile->mbs;
  441. ref_mb = tile->ref_mbs;
  442. row_offset = band->mb_size * band->pitch;
  443. need_mc = 0; /* reset the mc tracking flag */
  444. for (y = tile->ypos; y < (tile->ypos + tile->height); y += band->mb_size) {
  445. mb_offset = offs;
  446. for (x = tile->xpos; x < (tile->xpos + tile->width); x += band->mb_size) {
  447. mb->xpos = x;
  448. mb->ypos = y;
  449. mb->buf_offs = mb_offset;
  450. mb->type = 1; /* set the macroblocks type = INTER */
  451. mb->cbp = 0; /* all blocks are empty */
  452. if (!band->qdelta_present && !band->plane && !band->band_num) {
  453. mb->q_delta = band->glob_quant;
  454. mb->mv_x = 0;
  455. mb->mv_y = 0;
  456. }
  457. if (band->inherit_qdelta && ref_mb)
  458. mb->q_delta = ref_mb->q_delta;
  459. if (band->inherit_mv) {
  460. /* motion vector inheritance */
  461. if (mv_scale) {
  462. mb->mv_x = ivi_scale_mv(ref_mb->mv_x, mv_scale);
  463. mb->mv_y = ivi_scale_mv(ref_mb->mv_y, mv_scale);
  464. } else {
  465. mb->mv_x = ref_mb->mv_x;
  466. mb->mv_y = ref_mb->mv_y;
  467. }
  468. need_mc |= mb->mv_x || mb->mv_y; /* tracking non-zero motion vectors */
  469. }
  470. mb++;
  471. if (ref_mb)
  472. ref_mb++;
  473. mb_offset += band->mb_size;
  474. } // for x
  475. offs += row_offset;
  476. } // for y
  477. if (band->inherit_mv && need_mc) { /* apply motion compensation if there is at least one non-zero motion vector */
  478. num_blocks = (band->mb_size != band->blk_size) ? 4 : 1; /* number of blocks per mb */
  479. mc_no_delta_func = (band->blk_size == 8) ? ff_ivi_mc_8x8_no_delta
  480. : ff_ivi_mc_4x4_no_delta;
  481. for (mbn = 0, mb = tile->mbs; mbn < tile->num_MBs; mb++, mbn++) {
  482. mv_x = mb->mv_x;
  483. mv_y = mb->mv_y;
  484. if (!band->is_halfpel) {
  485. mc_type = 0; /* we have only fullpel vectors */
  486. } else {
  487. mc_type = ((mv_y & 1) << 1) | (mv_x & 1);
  488. mv_x >>= 1;
  489. mv_y >>= 1; /* convert halfpel vectors into fullpel ones */
  490. }
  491. for (blk = 0; blk < num_blocks; blk++) {
  492. /* adjust block position in the buffer according with its number */
  493. offs = mb->buf_offs + band->blk_size * ((blk & 1) + !!(blk & 2) * band->pitch);
  494. mc_no_delta_func(band->buf + offs,
  495. band->ref_buf + offs + mv_y * band->pitch + mv_x,
  496. band->pitch, mc_type);
  497. }
  498. }
  499. } else {
  500. /* copy data from the reference tile into the current one */
  501. src = band->ref_buf + tile->ypos * band->pitch + tile->xpos;
  502. dst = band->buf + tile->ypos * band->pitch + tile->xpos;
  503. for (y = 0; y < tile->height; y++) {
  504. memcpy(dst, src, tile->width*sizeof(band->buf[0]));
  505. src += band->pitch;
  506. dst += band->pitch;
  507. }
  508. }
  509. return 0;
  510. }
  511. #ifdef DEBUG
  512. static uint16_t ivi_calc_band_checksum(IVIBandDesc *band)
  513. {
  514. int x, y;
  515. int16_t *src, checksum;
  516. src = band->buf;
  517. checksum = 0;
  518. for (y = 0; y < band->height; src += band->pitch, y++)
  519. for (x = 0; x < band->width; x++)
  520. checksum += src[x];
  521. return checksum;
  522. }
  523. #endif
  524. void ff_ivi_output_plane(IVIPlaneDesc *plane, uint8_t *dst, int dst_pitch)
  525. {
  526. int x, y;
  527. const int16_t *src = plane->bands[0].buf;
  528. uint32_t pitch = plane->bands[0].pitch;
  529. if (!src)
  530. return;
  531. for (y = 0; y < plane->height; y++) {
  532. for (x = 0; x < plane->width; x++)
  533. dst[x] = av_clip_uint8(src[x] + 128);
  534. src += pitch;
  535. dst += dst_pitch;
  536. }
  537. }
  538. /**
  539. * Decode an Indeo 4 or 5 band.
  540. *
  541. * @param[in,out] ctx ptr to the decoder context
  542. * @param[in,out] band ptr to the band descriptor
  543. * @param[in] avctx ptr to the AVCodecContext
  544. * @return result code: 0 = OK, -1 = error
  545. */
  546. static int decode_band(IVI45DecContext *ctx,
  547. IVIBandDesc *band, AVCodecContext *avctx)
  548. {
  549. int result, i, t, idx1, idx2, pos;
  550. IVITile *tile;
  551. band->buf = band->bufs[ctx->dst_buf];
  552. if (!band->buf) {
  553. av_log(avctx, AV_LOG_ERROR, "Band buffer points to no data!\n");
  554. return AVERROR_INVALIDDATA;
  555. }
  556. band->ref_buf = band->bufs[ctx->ref_buf];
  557. band->data_ptr = ctx->frame_data + (get_bits_count(&ctx->gb) >> 3);
  558. result = ctx->decode_band_hdr(ctx, band, avctx);
  559. if (result) {
  560. av_log(avctx, AV_LOG_ERROR, "Error while decoding band header: %d\n",
  561. result);
  562. return result;
  563. }
  564. if (band->is_empty) {
  565. av_log(avctx, AV_LOG_ERROR, "Empty band encountered!\n");
  566. return AVERROR_INVALIDDATA;
  567. }
  568. band->rv_map = &ctx->rvmap_tabs[band->rvmap_sel];
  569. /* apply corrections to the selected rvmap table if present */
  570. for (i = 0; i < band->num_corr; i++) {
  571. idx1 = band->corr[i * 2];
  572. idx2 = band->corr[i * 2 + 1];
  573. FFSWAP(uint8_t, band->rv_map->runtab[idx1], band->rv_map->runtab[idx2]);
  574. FFSWAP(int16_t, band->rv_map->valtab[idx1], band->rv_map->valtab[idx2]);
  575. }
  576. pos = get_bits_count(&ctx->gb);
  577. for (t = 0; t < band->num_tiles; t++) {
  578. tile = &band->tiles[t];
  579. if (tile->mb_size != band->mb_size) {
  580. av_log(avctx, AV_LOG_ERROR, "MB sizes mismatch: %d vs. %d\n",
  581. band->mb_size, tile->mb_size);
  582. return AVERROR_INVALIDDATA;
  583. }
  584. tile->is_empty = get_bits1(&ctx->gb);
  585. if (tile->is_empty) {
  586. result = ivi_process_empty_tile(avctx, band, tile,
  587. (ctx->planes[0].bands[0].mb_size >> 3) - (band->mb_size >> 3));
  588. if (result < 0)
  589. break;
  590. av_dlog(avctx, "Empty tile encountered!\n");
  591. } else {
  592. tile->data_size = ff_ivi_dec_tile_data_size(&ctx->gb);
  593. if (!tile->data_size) {
  594. av_log(avctx, AV_LOG_ERROR, "Tile data size is zero!\n");
  595. return AVERROR_INVALIDDATA;
  596. }
  597. result = ctx->decode_mb_info(ctx, band, tile, avctx);
  598. if (result < 0)
  599. break;
  600. result = ff_ivi_decode_blocks(&ctx->gb, band, tile);
  601. if (result < 0 || ((get_bits_count(&ctx->gb) - pos) >> 3) != tile->data_size) {
  602. av_log(avctx, AV_LOG_ERROR, "Corrupted tile data encountered!\n");
  603. break;
  604. }
  605. pos += tile->data_size << 3; // skip to next tile
  606. }
  607. }
  608. /* restore the selected rvmap table by applying its corrections in reverse order */
  609. for (i = band->num_corr-1; i >= 0; i--) {
  610. idx1 = band->corr[i*2];
  611. idx2 = band->corr[i*2+1];
  612. FFSWAP(uint8_t, band->rv_map->runtab[idx1], band->rv_map->runtab[idx2]);
  613. FFSWAP(int16_t, band->rv_map->valtab[idx1], band->rv_map->valtab[idx2]);
  614. }
  615. #ifdef DEBUG
  616. if (band->checksum_present) {
  617. uint16_t chksum = ivi_calc_band_checksum(band);
  618. if (chksum != band->checksum) {
  619. av_log(avctx, AV_LOG_ERROR,
  620. "Band checksum mismatch! Plane %d, band %d, received: %x, calculated: %x\n",
  621. band->plane, band->band_num, band->checksum, chksum);
  622. }
  623. }
  624. #endif
  625. align_get_bits(&ctx->gb);
  626. return result;
  627. }
  628. int ff_ivi_decode_frame(AVCodecContext *avctx, void *data, int *got_frame,
  629. AVPacket *avpkt)
  630. {
  631. IVI45DecContext *ctx = avctx->priv_data;
  632. const uint8_t *buf = avpkt->data;
  633. int buf_size = avpkt->size;
  634. int result, p, b;
  635. init_get_bits(&ctx->gb, buf, buf_size * 8);
  636. ctx->frame_data = buf;
  637. ctx->frame_size = buf_size;
  638. result = ctx->decode_pic_hdr(ctx, avctx);
  639. if (result) {
  640. av_log(avctx, AV_LOG_ERROR,
  641. "Error while decoding picture header: %d\n", result);
  642. return -1;
  643. }
  644. if (ctx->gop_invalid)
  645. return AVERROR_INVALIDDATA;
  646. if (ctx->gop_flags & IVI5_IS_PROTECTED) {
  647. av_log(avctx, AV_LOG_ERROR, "Password-protected clip!\n");
  648. return -1;
  649. }
  650. ctx->switch_buffers(ctx);
  651. //{ START_TIMER;
  652. if (ctx->is_nonnull_frame(ctx)) {
  653. for (p = 0; p < 3; p++) {
  654. for (b = 0; b < ctx->planes[p].num_bands; b++) {
  655. result = decode_band(ctx, &ctx->planes[p].bands[b], avctx);
  656. if (result) {
  657. av_log(avctx, AV_LOG_ERROR,
  658. "Error while decoding band: %d, plane: %d\n", b, p);
  659. return -1;
  660. }
  661. }
  662. }
  663. }
  664. //STOP_TIMER("decode_planes"); }
  665. /* If the bidirectional mode is enabled, next I and the following P frame will */
  666. /* be sent together. Unfortunately the approach below seems to be the only way */
  667. /* to handle the B-frames mode. That's exactly the same Intel decoders do. */
  668. if (avctx->codec_id == AV_CODEC_ID_INDEO4 && ctx->frame_type == 0/*FRAMETYPE_INTRA*/) {
  669. while (get_bits(&ctx->gb, 8)); // skip version string
  670. skip_bits_long(&ctx->gb, 64); // skip padding, TODO: implement correct 8-bytes alignment
  671. if (get_bits_left(&ctx->gb) > 18 && show_bits(&ctx->gb, 18) == 0x3FFF8)
  672. av_log(avctx, AV_LOG_ERROR, "Buffer contains IP frames!\n");
  673. }
  674. if (ctx->frame.data[0])
  675. avctx->release_buffer(avctx, &ctx->frame);
  676. ctx->frame.reference = 0;
  677. avcodec_set_dimensions(avctx, ctx->planes[0].width, ctx->planes[0].height);
  678. if ((result = ff_get_buffer(avctx, &ctx->frame)) < 0) {
  679. av_log(avctx, AV_LOG_ERROR, "get_buffer() failed\n");
  680. return result;
  681. }
  682. if (ctx->is_scalable) {
  683. if (avctx->codec_id == AV_CODEC_ID_INDEO4)
  684. ff_ivi_recompose_haar(&ctx->planes[0], ctx->frame.data[0], ctx->frame.linesize[0]);
  685. else
  686. ff_ivi_recompose53 (&ctx->planes[0], ctx->frame.data[0], ctx->frame.linesize[0]);
  687. } else {
  688. ff_ivi_output_plane(&ctx->planes[0], ctx->frame.data[0], ctx->frame.linesize[0]);
  689. }
  690. ff_ivi_output_plane(&ctx->planes[2], ctx->frame.data[1], ctx->frame.linesize[1]);
  691. ff_ivi_output_plane(&ctx->planes[1], ctx->frame.data[2], ctx->frame.linesize[2]);
  692. *got_frame = 1;
  693. *(AVFrame*)data = ctx->frame;
  694. return buf_size;
  695. }
  696. /**
  697. * Close Indeo5 decoder and clean up its context.
  698. */
  699. av_cold int ff_ivi_decode_close(AVCodecContext *avctx)
  700. {
  701. IVI45DecContext *ctx = avctx->priv_data;
  702. ff_ivi_free_buffers(&ctx->planes[0]);
  703. if (ctx->mb_vlc.cust_tab.table)
  704. ff_free_vlc(&ctx->mb_vlc.cust_tab);
  705. if (ctx->frame.data[0])
  706. avctx->release_buffer(avctx, &ctx->frame);
  707. #if IVI4_STREAM_ANALYSER
  708. if (avctx->codec_id == AV_CODEC_ID_INDEO4) {
  709. if (ctx->is_scalable)
  710. av_log(avctx, AV_LOG_ERROR, "This video uses scalability mode!\n");
  711. if (ctx->uses_tiling)
  712. av_log(avctx, AV_LOG_ERROR, "This video uses local decoding!\n");
  713. if (ctx->has_b_frames)
  714. av_log(avctx, AV_LOG_ERROR, "This video contains B-frames!\n");
  715. if (ctx->has_transp)
  716. av_log(avctx, AV_LOG_ERROR, "Transparency mode is enabled!\n");
  717. if (ctx->uses_haar)
  718. av_log(avctx, AV_LOG_ERROR, "This video uses Haar transform!\n");
  719. if (ctx->uses_fullpel)
  720. av_log(avctx, AV_LOG_ERROR, "This video uses fullpel motion vectors!\n");
  721. }
  722. #endif
  723. return 0;
  724. }
  725. /**
  726. * These are 2x8 predefined Huffman codebooks for coding macroblock/block
  727. * signals. They are specified using "huffman descriptors" in order to
  728. * avoid huge static tables. The decoding tables will be generated at
  729. * startup from these descriptors.
  730. */
  731. const IVIHuffDesc ff_ivi_mb_huff_desc[8] = {
  732. {8, {0, 4, 5, 4, 4, 4, 6, 6}},
  733. {12, {0, 2, 2, 3, 3, 3, 3, 5, 3, 2, 2, 2}},
  734. {12, {0, 2, 3, 4, 3, 3, 3, 3, 4, 3, 2, 2}},
  735. {12, {0, 3, 4, 4, 3, 3, 3, 3, 3, 2, 2, 2}},
  736. {13, {0, 4, 4, 3, 3, 3, 3, 2, 3, 3, 2, 1, 1}},
  737. {9, {0, 4, 4, 4, 4, 3, 3, 3, 2}},
  738. {10, {0, 4, 4, 4, 4, 3, 3, 2, 2, 2}},
  739. {12, {0, 4, 4, 4, 3, 3, 2, 3, 2, 2, 2, 2}}
  740. };
  741. const IVIHuffDesc ff_ivi_blk_huff_desc[8] = {
  742. {10, {1, 2, 3, 4, 4, 7, 5, 5, 4, 1}},
  743. {11, {2, 3, 4, 4, 4, 7, 5, 4, 3, 3, 2}},
  744. {12, {2, 4, 5, 5, 5, 5, 6, 4, 4, 3, 1, 1}},
  745. {13, {3, 3, 4, 4, 5, 6, 6, 4, 4, 3, 2, 1, 1}},
  746. {11, {3, 4, 4, 5, 5, 5, 6, 5, 4, 2, 2}},
  747. {13, {3, 4, 5, 5, 5, 5, 6, 4, 3, 3, 2, 1, 1}},
  748. {13, {3, 4, 5, 5, 5, 6, 5, 4, 3, 3, 2, 1, 1}},
  749. {9, {3, 4, 4, 5, 5, 5, 6, 5, 5}}
  750. };
  751. /**
  752. * Scan patterns shared between indeo4 and indeo5
  753. */
  754. const uint8_t ff_ivi_vertical_scan_8x8[64] = {
  755. 0, 8, 16, 24, 32, 40, 48, 56,
  756. 1, 9, 17, 25, 33, 41, 49, 57,
  757. 2, 10, 18, 26, 34, 42, 50, 58,
  758. 3, 11, 19, 27, 35, 43, 51, 59,
  759. 4, 12, 20, 28, 36, 44, 52, 60,
  760. 5, 13, 21, 29, 37, 45, 53, 61,
  761. 6, 14, 22, 30, 38, 46, 54, 62,
  762. 7, 15, 23, 31, 39, 47, 55, 63
  763. };
  764. const uint8_t ff_ivi_horizontal_scan_8x8[64] = {
  765. 0, 1, 2, 3, 4, 5, 6, 7,
  766. 8, 9, 10, 11, 12, 13, 14, 15,
  767. 16, 17, 18, 19, 20, 21, 22, 23,
  768. 24, 25, 26, 27, 28, 29, 30, 31,
  769. 32, 33, 34, 35, 36, 37, 38, 39,
  770. 40, 41, 42, 43, 44, 45, 46, 47,
  771. 48, 49, 50, 51, 52, 53, 54, 55,
  772. 56, 57, 58, 59, 60, 61, 62, 63
  773. };
  774. const uint8_t ff_ivi_direct_scan_4x4[16] = {
  775. 0, 1, 4, 8, 5, 2, 3, 6, 9, 12, 13, 10, 7, 11, 14, 15
  776. };
  777. /**
  778. * Run-value (RLE) tables.
  779. */
  780. const RVMapDesc ff_ivi_rvmap_tabs[9] = {
  781. { /* MapTab0 */
  782. 5, /* eob_sym */
  783. 2, /* esc_sym */
  784. /* run table */
  785. {1, 1, 0, 1, 1, 0, 1, 1, 2, 2, 1, 1, 1, 1, 3, 3,
  786. 1, 1, 2, 2, 1, 1, 4, 4, 1, 1, 1, 1, 2, 2, 5, 5,
  787. 1, 1, 3, 3, 1, 1, 6, 6, 1, 2, 1, 2, 7, 7, 1, 1,
  788. 8, 8, 1, 1, 4, 2, 1, 4, 2, 1, 3, 3, 1, 1, 1, 9,
  789. 9, 1, 2, 1, 2, 1, 5, 5, 1, 1, 10, 10, 1, 1, 3, 3,
  790. 2, 2, 1, 1, 11, 11, 6, 4, 4, 1, 6, 1, 2, 1, 2, 12,
  791. 8, 1, 12, 7, 8, 7, 1, 16, 1, 16, 1, 3, 3, 13, 1, 13,
  792. 2, 2, 1, 15, 1, 5, 14, 15, 1, 5, 14, 1, 17, 8, 17, 8,
  793. 1, 4, 4, 2, 2, 1, 25, 25, 24, 24, 1, 3, 1, 3, 1, 8,
  794. 6, 7, 6, 1, 18, 8, 18, 1, 7, 23, 2, 2, 23, 1, 1, 21,
  795. 22, 9, 9, 22, 19, 1, 21, 5, 19, 5, 1, 33, 20, 33, 20, 8,
  796. 4, 4, 1, 32, 2, 2, 8, 3, 32, 26, 3, 1, 7, 7, 26, 6,
  797. 1, 6, 1, 1, 16, 1, 10, 1, 10, 2, 16, 29, 28, 2, 29, 28,
  798. 1, 27, 5, 8, 5, 27, 1, 8, 3, 7, 3, 31, 41, 31, 1, 41,
  799. 6, 1, 6, 7, 4, 4, 1, 1, 2, 1, 2, 11, 34, 30, 11, 1,
  800. 30, 15, 15, 34, 36, 40, 36, 40, 35, 35, 37, 37, 39, 39, 38, 38},
  801. /* value table */
  802. { 1, -1, 0, 2, -2, 0, 3, -3, 1, -1, 4, -4, 5, -5, 1, -1,
  803. 6, -6, 2, -2, 7, -7, 1, -1, 8, -8, 9, -9, 3, -3, 1, -1,
  804. 10, -10, 2, -2, 11, -11, 1, -1, 12, 4, -12, -4, 1, -1, 13, -13,
  805. 1, -1, 14, -14, 2, 5, 15, -2, -5, -15, -3, 3, 16, -16, 17, 1,
  806. -1, -17, 6, 18, -6, -18, 2, -2, 19, -19, 1, -1, 20, -20, 4, -4,
  807. 7, -7, 21, -21, 1, -1, 2, 3, -3, 22, -2, -22, 8, 23, -8, 1,
  808. 2, -23, -1, 2, -2, -2, 24, 1, -24, -1, 25, 5, -5, 1, -25, -1,
  809. 9, -9, 26, 1, -26, 3, 1, -1, 27, -3, -1, -27, 1, 3, -1, -3,
  810. 28, -4, 4, 10, -10, -28, 1, -1, 1, -1, 29, 6, -29, -6, 30, -4,
  811. 3, 3, -3, -30, 1, 4, -1, 31, -3, 1, 11, -11, -1, -31, 32, -1,
  812. -1, 2, -2, 1, 1, -32, 1, 4, -1, -4, 33, -1, 1, 1, -1, 5,
  813. 5, -5, -33, -1, -12, 12, -5, -7, 1, 1, 7, 34, 4, -4, -1, 4,
  814. -34, -4, 35, 36, -2, -35, -2, -36, 2, 13, 2, -1, 1, -13, 1, -1,
  815. 37, 1, -5, 6, 5, -1, 38, -6, -8, 5, 8, -1, 1, 1, -37, -1,
  816. 5, 39, -5, -5, 6, -6, -38, -39, -14, 40, 14, 2, 1, 1, -2, -40,
  817. -1, -2, 2, -1, -1, -1, 1, 1, 1, -1, 1, -1, 1, -1, 1, -1}
  818. },{
  819. /* MapTab1 */
  820. 0, /* eob_sym */
  821. 38, /* esc_sym */
  822. /* run table */
  823. {0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 8, 6, 8, 7,
  824. 7, 9, 9, 10, 10, 11, 11, 1, 12, 1, 12, 13, 13, 16, 14, 16,
  825. 14, 15, 15, 17, 17, 18, 0, 18, 19, 20, 21, 19, 22, 21, 20, 22,
  826. 25, 24, 2, 25, 24, 23, 23, 2, 26, 28, 26, 28, 29, 27, 29, 27,
  827. 33, 33, 1, 32, 1, 3, 32, 30, 36, 3, 36, 30, 31, 31, 35, 34,
  828. 37, 41, 34, 35, 37, 4, 41, 4, 49, 8, 8, 49, 40, 38, 5, 38,
  829. 40, 39, 5, 39, 42, 43, 42, 7, 57, 6, 43, 44, 6, 50, 7, 44,
  830. 57, 48, 50, 48, 45, 45, 46, 47, 51, 46, 47, 58, 1, 51, 58, 1,
  831. 52, 59, 53, 9, 52, 55, 55, 59, 53, 56, 54, 56, 54, 9, 64, 64,
  832. 60, 63, 60, 63, 61, 62, 61, 62, 2, 10, 2, 10, 11, 1, 11, 13,
  833. 12, 1, 12, 13, 16, 16, 8, 8, 14, 3, 3, 15, 14, 15, 4, 4,
  834. 1, 17, 17, 5, 1, 7, 7, 5, 6, 1, 2, 2, 6, 22, 1, 25,
  835. 21, 22, 8, 24, 1, 21, 25, 24, 8, 18, 18, 23, 9, 20, 23, 33,
  836. 29, 33, 20, 1, 19, 1, 29, 36, 9, 36, 19, 41, 28, 57, 32, 3,
  837. 28, 3, 1, 27, 49, 49, 1, 32, 26, 26, 2, 4, 4, 7, 57, 41,
  838. 2, 7, 10, 5, 37, 16, 10, 27, 8, 8, 13, 16, 37, 13, 1, 5},
  839. /* value table */
  840. {0, 1, -1, 1, -1, 1, -1, 1, -1, 1, -1, 1, 1, -1, -1, 1,
  841. -1, 1, -1, 1, -1, 1, -1, 2, 1, -2, -1, 1, -1, 1, 1, -1,
  842. -1, 1, -1, 1, -1, 1, 0, -1, 1, 1, 1, -1, 1, -1, -1, -1,
  843. 1, 1, 2, -1, -1, 1, -1, -2, 1, 1, -1, -1, 1, 1, -1, -1,
  844. 1, -1, 3, 1, -3, 2, -1, 1, 1, -2, -1, -1, -1, 1, 1, 1,
  845. 1, 1, -1, -1, -1, 2, -1, -2, 1, 2, -2, -1, 1, 1, 2, -1,
  846. -1, 1, -2, -1, 1, 1, -1, 2, 1, 2, -1, 1, -2, -1, -2, -1,
  847. -1, 1, 1, -1, 1, -1, 1, 1, 1, -1, -1, 1, 4, -1, -1, -4,
  848. 1, 1, 1, 2, -1, -1, 1, -1, -1, 1, -1, -1, 1, -2, 1, -1,
  849. 1, 1, -1, -1, 1, 1, -1, -1, 3, 2, -3, -2, 2, 5, -2, 2,
  850. 2, -5, -2, -2, -2, 2, -3, 3, 2, 3, -3, 2, -2, -2, 3, -3,
  851. 6, 2, -2, 3, -6, 3, -3, -3, 3, 7, -4, 4, -3, 2, -7, 2,
  852. 2, -2, -4, 2, 8, -2, -2, -2, 4, 2, -2, 2, 3, 2, -2, -2,
  853. 2, 2, -2, -8, -2, 9, -2, 2, -3, -2, 2, -2, 2, 2, 2, 4,
  854. -2, -4, 10, 2, 2, -2, -9, -2, 2, -2, 5, 4, -4, 4, -2, 2,
  855. -5, -4, -3, 4, 2, -3, 3, -2, -5, 5, 3, 3, -2, -3, -10, -4}
  856. },{
  857. /* MapTab2 */
  858. 2, /* eob_sym */
  859. 11, /* esc_sym */
  860. /* run table */
  861. {1, 1, 0, 2, 2, 1, 1, 3, 3, 4, 4, 0, 1, 1, 5, 5,
  862. 2, 2, 6, 6, 7, 7, 1, 8, 1, 8, 3, 3, 9, 9, 1, 2,
  863. 2, 1, 4, 10, 4, 10, 11, 11, 1, 5, 12, 12, 1, 5, 13, 13,
  864. 3, 3, 6, 6, 2, 2, 14, 14, 16, 16, 15, 7, 15, 8, 8, 7,
  865. 1, 1, 17, 17, 4, 4, 1, 1, 18, 18, 2, 2, 5, 5, 25, 3,
  866. 9, 3, 25, 9, 19, 24, 19, 24, 1, 21, 20, 1, 21, 22, 20, 22,
  867. 23, 23, 8, 6, 33, 6, 8, 33, 7, 7, 26, 26, 1, 32, 1, 32,
  868. 28, 4, 28, 10, 29, 27, 27, 10, 41, 4, 29, 2, 2, 41, 36, 31,
  869. 49, 31, 34, 30, 34, 36, 30, 35, 1, 49, 11, 5, 35, 11, 1, 3,
  870. 3, 5, 37, 37, 8, 40, 8, 40, 12, 12, 42, 42, 1, 38, 16, 57,
  871. 1, 6, 16, 39, 38, 6, 7, 7, 13, 13, 39, 43, 2, 43, 57, 2,
  872. 50, 9, 44, 9, 50, 4, 15, 48, 44, 4, 1, 15, 48, 14, 14, 1,
  873. 45, 45, 8, 3, 5, 8, 51, 47, 3, 46, 46, 47, 5, 51, 1, 17,
  874. 17, 58, 1, 58, 2, 52, 52, 2, 53, 7, 59, 6, 6, 56, 53, 55,
  875. 7, 55, 1, 54, 59, 56, 54, 10, 1, 10, 4, 60, 1, 60, 8, 4,
  876. 8, 64, 64, 61, 1, 63, 3, 63, 62, 61, 5, 11, 5, 3, 11, 62},
  877. /* value table */
  878. { 1, -1, 0, 1, -1, 2, -2, 1, -1, 1, -1, 0, 3, -3, 1, -1,
  879. 2, -2, 1, -1, 1, -1, 4, 1, -4, -1, 2, -2, 1, -1, 5, 3,
  880. -3, -5, 2, 1, -2, -1, 1, -1, 6, 2, 1, -1, -6, -2, 1, -1,
  881. 3, -3, 2, -2, 4, -4, 1, -1, 1, -1, 1, 2, -1, 2, -2, -2,
  882. 7, -7, 1, -1, 3, -3, 8, -8, 1, -1, 5, -5, 3, -3, 1, 4,
  883. 2, -4, -1, -2, 1, 1, -1, -1, 9, 1, 1, -9, -1, 1, -1, -1,
  884. 1, -1, 3, -3, 1, 3, -3, -1, 3, -3, 1, -1, 10, 1, -10, -1,
  885. 1, 4, -1, 2, 1, -1, 1, -2, 1, -4, -1, 6, -6, -1, 1, 1,
  886. 1, -1, 1, 1, -1, -1, -1, 1, 11, -1, -2, 4, -1, 2, -11, 5,
  887. -5, -4, -1, 1, 4, 1, -4, -1, -2, 2, 1, -1, 12, 1, -2, 1,
  888. -12, 4, 2, 1, -1, -4, 4, -4, 2, -2, -1, 1, 7, -1, -1, -7,
  889. -1, -3, 1, 3, 1, 5, 2, 1, -1, -5, 13, -2, -1, 2, -2, -13,
  890. 1, -1, 5, 6, 5, -5, 1, 1, -6, 1, -1, -1, -5, -1, 14, 2,
  891. -2, 1, -14, -1, 8, 1, -1, -8, 1, 5, 1, 5, -5, 1, -1, 1,
  892. -5, -1, 15, 1, -1, -1, -1, 3, -15, -3, 6, 1, 16, -1, 6, -6,
  893. -6, 1, -1, 1, -16, 1, 7, -1, 1, -1, -6, -3, 6, -7, 3, -1}
  894. },{
  895. /* MapTab3 */
  896. 0, /* eob_sym */
  897. 35, /* esc_sym */
  898. /* run table */
  899. {0, 1, 1, 2, 2, 3, 3, 4, 4, 1, 1, 5, 5, 6, 6, 7,
  900. 7, 8, 8, 9, 9, 2, 2, 10, 10, 1, 1, 11, 11, 12, 12, 3,
  901. 3, 13, 13, 0, 14, 14, 16, 15, 16, 15, 4, 4, 17, 1, 17, 1,
  902. 5, 5, 18, 18, 2, 2, 6, 6, 8, 19, 7, 8, 7, 19, 20, 20,
  903. 21, 21, 22, 24, 22, 24, 23, 23, 1, 1, 25, 25, 3, 3, 26, 26,
  904. 9, 9, 27, 27, 28, 28, 33, 29, 4, 33, 29, 1, 4, 1, 32, 32,
  905. 2, 2, 31, 10, 30, 10, 30, 31, 34, 34, 5, 5, 36, 36, 35, 41,
  906. 35, 11, 41, 11, 37, 1, 8, 8, 37, 6, 1, 6, 40, 7, 7, 40,
  907. 12, 38, 12, 39, 39, 38, 49, 13, 49, 13, 3, 42, 3, 42, 16, 16,
  908. 43, 43, 14, 14, 1, 1, 44, 15, 44, 15, 2, 2, 57, 48, 50, 48,
  909. 57, 50, 4, 45, 45, 4, 46, 47, 47, 46, 1, 51, 1, 17, 17, 51,
  910. 8, 9, 9, 5, 58, 8, 58, 5, 52, 52, 55, 56, 53, 56, 55, 59,
  911. 59, 53, 54, 1, 6, 54, 7, 7, 6, 1, 2, 3, 2, 3, 64, 60,
  912. 60, 10, 10, 64, 61, 62, 61, 63, 1, 63, 62, 1, 18, 24, 18, 4,
  913. 25, 4, 8, 21, 21, 1, 24, 22, 25, 22, 8, 11, 19, 11, 23, 1,
  914. 20, 23, 19, 20, 5, 12, 5, 1, 16, 2, 12, 13, 2, 13, 1, 16},
  915. /* value table */
  916. { 0, 1, -1, 1, -1, 1, -1, 1, -1, 2, -2, 1, -1, 1, -1, 1,
  917. -1, 1, -1, 1, -1, 2, -2, 1, -1, 3, -3, 1, -1, 1, -1, 2,
  918. -2, 1, -1, 0, 1, -1, 1, 1, -1, -1, 2, -2, 1, 4, -1, -4,
  919. 2, -2, 1, -1, -3, 3, 2, -2, 2, 1, 2, -2, -2, -1, 1, -1,
  920. 1, -1, 1, 1, -1, -1, 1, -1, 5, -5, 1, -1, 3, -3, 1, -1,
  921. 2, -2, 1, -1, 1, -1, 1, 1, 3, -1, -1, 6, -3, -6, -1, 1,
  922. 4, -4, 1, 2, 1, -2, -1, -1, 1, -1, 3, -3, 1, -1, 1, 1,
  923. -1, 2, -1, -2, 1, 7, -3, 3, -1, 3, -7, -3, 1, -3, 3, -1,
  924. 2, 1, -2, 1, -1, -1, 1, 2, -1, -2, -4, -1, 4, 1, 2, -2,
  925. 1, -1, -2, 2, 8, -8, -1, 2, 1, -2, -5, 5, 1, -1, -1, 1,
  926. -1, 1, 4, -1, 1, -4, -1, -1, 1, 1, 9, 1, -9, 2, -2, -1,
  927. -4, 3, -3, -4, -1, 4, 1, 4, 1, -1, 1, -1, 1, 1, -1, 1,
  928. -1, -1, -1, 10, 4, 1, 4, -4, -4, -10, 6, 5, -6, -5, 1, -1,
  929. 1, 3, -3, -1, 1, -1, -1, -1, 11, 1, 1, -11, -2, -2, 2, 5,
  930. -2, -5, -5, 2, -2, 12, 2, -2, 2, 2, 5, -3, -2, 3, -2, -12,
  931. -2, 2, 2, 2, -5, 3, 5, 13, -3, 7, -3, -3, -7, 3, -13, 3}
  932. },{
  933. /* MapTab4 */
  934. 0, /* eob_sym */
  935. 34, /* esc_sym */
  936. /* run table */
  937. {0, 1, 1, 1, 2, 2, 1, 3, 3, 1, 1, 1, 4, 4, 1, 5,
  938. 2, 1, 5, 2, 1, 1, 6, 6, 1, 1, 1, 1, 1, 7, 3, 1,
  939. 2, 3, 0, 1, 2, 7, 1, 1, 1, 8, 1, 1, 8, 1, 1, 1,
  940. 9, 1, 9, 1, 2, 1, 1, 2, 1, 1, 10, 4, 1, 10, 1, 4,
  941. 1, 1, 1, 1, 1, 3, 1, 1, 1, 3, 2, 1, 5, 1, 1, 1,
  942. 2, 5, 1, 11, 1, 11, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
  943. 2, 1, 6, 1, 6, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 12,
  944. 3, 1, 12, 1, 1, 1, 2, 1, 1, 3, 1, 1, 1, 1, 1, 1,
  945. 4, 1, 1, 1, 2, 1, 1, 4, 1, 1, 1, 1, 1, 1, 2, 1,
  946. 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 2, 1, 1, 5,
  947. 1, 1, 1, 1, 1, 7, 1, 7, 1, 1, 2, 3, 1, 1, 1, 1,
  948. 5, 1, 1, 1, 1, 1, 1, 2, 13, 1, 1, 1, 1, 1, 1, 1,
  949. 1, 1, 1, 1, 1, 1, 1, 1, 13, 2, 1, 1, 4, 1, 1, 1,
  950. 3, 1, 6, 1, 1, 1, 14, 1, 1, 1, 1, 1, 14, 6, 1, 1,
  951. 1, 1, 15, 2, 4, 1, 2, 3, 15, 1, 1, 1, 8, 1, 1, 8,
  952. 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1},
  953. /* value table */
  954. { 0, 1, -1, 2, 1, -1, -2, 1, -1, 3, -3, 4, 1, -1, -4, 1,
  955. 2, 5, -1, -2, -5, 6, 1, -1, -6, 7, -7, 8, -8, 1, 2, 9,
  956. 3, -2, 0, -9, -3, -1, 10, -10, 11, 1, -11, 12, -1, -12, 13, -13,
  957. 1, 14, -1, -14, 4, 15, -15, -4, 16, -16, 1, 2, 17, -1, -17, -2,
  958. 18, -18, 19, -19, 20, 3, -20, 21, -21, -3, 5, 22, 2, -22, -23, 23,
  959. -5, -2, 24, 1, -24, -1, 25, -25, 26, -26, -27, 27, 28, 29, -28, -29,
  960. 6, 30, 2, -31, -2, -30, 31, -6, -32, 32, 33, -33, 34, -35, -34, 1,
  961. 4, -36, -1, 35, 37, 36, 7, -37, 38, -4, -38, 39, 41, 40, -40, -39,
  962. 3, 42, -43, -41, -7, -42, 43, -3, 44, -44, 45, -45, 46, 47, 8, -47,
  963. -48, -46, 50, -50, 48, 49, 51, -49, 52, -52, 5, -51, -8, -53, 53, 3,
  964. -56, 56, 55, 54, -54, 2, 60, -2, -55, 58, 9, -5, 59, 57, -57, -63,
  965. -3, -58, -60, -61, 61, -59, -62, -9, 1, 64, 62, 69, -64, 63, 65, -67,
  966. -68, 66, -65, 68, -66, -69, 67, -70, -1, 10, 71, -71, 4, 73, 72, 70,
  967. 6, -76, -3, 74, -78, -74, 1, 78, 80, -72, -75, 76, -1, 3, -73, 79,
  968. 75, 77, 1, 11, -4, -79, -10, -6, -1, -77, -83, -80, 2, 81, -84, -2,
  969. 83, -81, 82, -82, 84, -87, -86, 85, -11, -85, 86, -89, 87, -88, 88, 89}
  970. },{
  971. /* MapTab5 */
  972. 2, /* eob_sym */
  973. 33, /* esc_sym */
  974. /* run table */
  975. {1, 1, 0, 2, 1, 2, 1, 3, 3, 1, 1, 4, 4, 2, 2, 1,
  976. 1, 5, 5, 6, 1, 6, 1, 7, 7, 3, 3, 2, 8, 2, 8, 1,
  977. 1, 0, 9, 9, 1, 1, 10, 4, 10, 4, 11, 11, 2, 1, 2, 1,
  978. 12, 12, 3, 3, 1, 1, 13, 5, 5, 13, 14, 1, 1, 14, 2, 2,
  979. 6, 6, 15, 1, 1, 15, 16, 4, 7, 16, 4, 7, 1, 1, 3, 3,
  980. 8, 8, 2, 2, 1, 1, 17, 17, 1, 1, 18, 18, 5, 5, 2, 2,
  981. 1, 1, 9, 19, 9, 19, 20, 3, 3, 20, 1, 10, 21, 1, 10, 4,
  982. 4, 21, 22, 6, 6, 22, 1, 1, 23, 24, 2, 2, 23, 24, 11, 1,
  983. 1, 11, 7, 25, 7, 1, 1, 25, 8, 8, 3, 26, 3, 1, 12, 2,
  984. 2, 26, 1, 12, 5, 5, 27, 4, 1, 4, 1, 27, 28, 1, 28, 13,
  985. 1, 13, 2, 29, 2, 1, 32, 6, 1, 30, 14, 29, 14, 6, 3, 31,
  986. 3, 1, 30, 1, 32, 31, 33, 9, 33, 1, 1, 7, 9, 7, 2, 2,
  987. 1, 1, 4, 36, 34, 4, 5, 10, 10, 5, 34, 1, 1, 35, 8, 8,
  988. 36, 3, 35, 1, 15, 3, 2, 1, 16, 15, 16, 2, 37, 1, 37, 1,
  989. 1, 1, 6, 6, 38, 1, 38, 11, 1, 39, 39, 40, 11, 2, 41, 4,
  990. 40, 1, 2, 4, 1, 1, 1, 41, 3, 1, 3, 1, 5, 7, 5, 7},
  991. /* value table */
  992. { 1, -1, 0, 1, 2, -1, -2, 1, -1, 3, -3, 1, -1, 2, -2, 4,
  993. -4, 1, -1, 1, 5, -1, -5, 1, -1, 2, -2, 3, 1, -3, -1, 6,
  994. -6, 0, 1, -1, 7, -7, 1, 2, -1, -2, 1, -1, 4, 8, -4, -8,
  995. 1, -1, 3, -3, 9, -9, 1, 2, -2, -1, 1, 10, -10, -1, 5, -5,
  996. 2, -2, 1, 11, -11, -1, 1, 3, 2, -1, -3, -2, 12, -12, 4, -4,
  997. 2, -2, -6, 6, 13, -13, 1, -1, 14, -14, 1, -1, 3, -3, 7, -7,
  998. 15, -15, 2, 1, -2, -1, 1, 5, -5, -1, -16, 2, 1, 16, -2, 4,
  999. -4, -1, 1, 3, -3, -1, 17, -17, 1, 1, -8, 8, -1, -1, 2, 18,
  1000. -18, -2, 3, 1, -3, 19, -19, -1, 3, -3, 6, 1, -6, 20, 2, 9,
  1001. -9, -1, -20, -2, 4, -4, 1, -5, 21, 5, -21, -1, 1, -22, -1, 2,
  1002. 22, -2, 10, 1, -10, 23, 1, 4, -23, 1, 2, -1, -2, -4, -7, 1,
  1003. 7, -24, -1, 24, -1, -1, 1, 3, -1, -25, 25, 4, -3, -4, 11, -11,
  1004. 26, -26, 6, 1, 1, -6, -5, -3, 3, 5, -1, -27, 27, 1, 4, -4,
  1005. -1, -8, -1, 28, 2, 8, -12, -28, -2, -2, 2, 12, -1, 29, 1, -29,
  1006. 30, -30, 5, -5, 1, -31, -1, 3, 31, -1, 1, 1, -3, -13, 1, -7,
  1007. -1, -32, 13, 7, 32, 33, -33, -1, -9, -34, 9, 34, -6, 5, 6, -5}
  1008. },{
  1009. /* MapTab6 */
  1010. 2, /* eob_sym */
  1011. 13, /* esc_sym */
  1012. /* run table */
  1013. {1, 1, 0, 1, 1, 2, 2, 1, 1, 3, 3, 1, 1, 0, 2, 2,
  1014. 4, 1, 4, 1, 1, 1, 5, 5, 1, 1, 6, 6, 2, 2, 1, 1,
  1015. 3, 3, 7, 7, 1, 1, 8, 8, 1, 1, 2, 2, 1, 9, 1, 9,
  1016. 4, 4, 10, 1, 1, 10, 1, 1, 11, 11, 3, 3, 1, 2, 1, 2,
  1017. 1, 1, 12, 12, 5, 5, 1, 1, 13, 1, 1, 13, 2, 2, 1, 1,
  1018. 6, 6, 1, 1, 4, 14, 4, 14, 3, 1, 3, 1, 1, 1, 15, 7,
  1019. 15, 2, 2, 7, 1, 1, 1, 8, 1, 8, 16, 16, 1, 1, 1, 1,
  1020. 2, 1, 1, 2, 1, 1, 3, 5, 5, 3, 4, 1, 1, 4, 1, 1,
  1021. 17, 17, 9, 1, 1, 9, 2, 2, 1, 1, 10, 10, 1, 6, 1, 1,
  1022. 6, 18, 1, 1, 18, 1, 1, 1, 2, 2, 3, 1, 3, 1, 1, 1,
  1023. 4, 1, 19, 1, 19, 7, 1, 1, 20, 1, 4, 20, 1, 7, 11, 2,
  1024. 1, 11, 21, 2, 8, 5, 1, 8, 1, 5, 21, 1, 1, 1, 22, 1,
  1025. 1, 22, 1, 1, 3, 3, 1, 23, 2, 12, 24, 1, 1, 2, 1, 1,
  1026. 12, 23, 1, 1, 24, 1, 1, 1, 4, 1, 1, 1, 2, 1, 6, 6,
  1027. 4, 2, 1, 1, 1, 1, 1, 1, 1, 14, 13, 3, 1, 25, 9, 25,
  1028. 14, 1, 9, 3, 13, 1, 1, 1, 1, 1, 10, 1, 1, 2, 10, 2},
  1029. /* value table */
  1030. {-20, -1, 0, 2, -2, 1, -1, 3, -3, 1, -1, 4, -4, 0, 2, -2,
  1031. 1, 5, -1, -5, 6, -6, 1, -1, 7, -7, 1, -1, 3, -3, 8, -8,
  1032. 2, -2, 1, -1, 9, -9, 1, -1, 10, -10, 4, -4, 11, 1, -11, -1,
  1033. 2, -2, 1, 12, -12, -1, 13, -13, 1, -1, 3, -3, 14, 5, -14, -5,
  1034. -15, 15, -1, 1, 2, -2, 16, -16, 1, 17, -17, -1, 6, -6, 18, -18,
  1035. 2, -2, -19, 19, -3, 1, 3, -1, 4, 20, -4, 1, -21, 21, 1, 2,
  1036. -1, -7, 7, -2, 22, -22, 23, 2, -23, -2, 1, -1, -24, 24, -25, 25,
  1037. -8, -26, 26, 8, -27, 27, 5, 3, -3, -5, -4, 28, -28, 4, 29, -29,
  1038. 1, -1, -2, -30, 30, 2, 9, -9, -31, 31, 2, -2, -32, 3, 32, -33,
  1039. -3, 1, 33, -34, -1, 34, -35, 35, -10, 10, -6, 36, 6, -36, 37, -37,
  1040. -5, 38, 1, -38, -1, 3, 39, -39, -1, 40, 5, 1, -40, -3, 2, -11,
  1041. -41, -2, 1, 11, -3, -4, 41, 3, 42, 4, -1, -43, -42, 43, 1, -44,
  1042. 45, -1, 44, -45, -7, 7, -46, 1, -12, 2, 1, -47, 46, 12, 47, 48,
  1043. -2, -1, -48, 49, -1, -50, -49, 50, -6, -51, 51, 52, -13, 53, -4, 4,
  1044. 6, 13, -53, -52, -54, 55, 54, -55, -56, -2, 2, -8, 56, 1, -3, -1,
  1045. 2, 58, 3, 8, -2, 57, -58, -60, -59, -57, -3, 60, 59, -14, 3, 14}
  1046. },{
  1047. /* MapTab7 */
  1048. 2, /* eob_sym */
  1049. 38, /* esc_sym */
  1050. /* run table */
  1051. {1, 1, 0, 2, 2, 1, 1, 3, 3, 4, 4, 5, 5, 1, 1, 6,
  1052. 6, 2, 2, 7, 7, 8, 8, 1, 1, 3, 3, 9, 9, 10, 10, 1,
  1053. 1, 2, 2, 4, 4, 11, 0, 11, 12, 12, 13, 13, 1, 1, 5, 5,
  1054. 14, 14, 15, 16, 15, 16, 3, 3, 1, 6, 1, 6, 2, 2, 7, 7,
  1055. 8, 8, 17, 17, 1, 1, 4, 4, 18, 18, 2, 2, 1, 19, 1, 20,
  1056. 19, 20, 21, 21, 3, 3, 22, 22, 5, 5, 24, 1, 1, 23, 9, 23,
  1057. 24, 9, 2, 2, 10, 1, 1, 10, 6, 6, 25, 4, 4, 25, 7, 7,
  1058. 26, 8, 1, 8, 3, 1, 26, 3, 11, 11, 27, 27, 2, 28, 1, 2,
  1059. 28, 1, 12, 12, 5, 5, 29, 13, 13, 29, 32, 1, 1, 33, 31, 30,
  1060. 32, 4, 30, 33, 4, 31, 3, 14, 1, 1, 3, 34, 34, 2, 2, 14,
  1061. 6, 6, 35, 36, 35, 36, 1, 15, 1, 16, 16, 15, 7, 9, 7, 9,
  1062. 37, 8, 8, 37, 1, 1, 39, 2, 38, 39, 2, 40, 5, 38, 40, 5,
  1063. 3, 3, 4, 4, 10, 10, 1, 1, 1, 1, 41, 2, 41, 2, 6, 6,
  1064. 1, 1, 11, 42, 11, 43, 3, 42, 3, 17, 4, 43, 1, 17, 7, 1,
  1065. 8, 44, 4, 7, 44, 5, 8, 2, 5, 1, 2, 48, 45, 1, 12, 45,
  1066. 12, 48, 13, 13, 1, 9, 9, 46, 1, 46, 47, 47, 49, 18, 18, 49},
  1067. /* value table */
  1068. { 1, -1, 0, 1, -1, 2, -2, 1, -1, 1, -1, 1, -1, 3, -3, 1,
  1069. -1, -2, 2, 1, -1, 1, -1, 4, -4, -2, 2, 1, -1, 1, -1, 5,
  1070. -5, -3, 3, 2, -2, 1, 0, -1, 1, -1, 1, -1, 6, -6, 2, -2,
  1071. 1, -1, 1, 1, -1, -1, -3, 3, 7, 2, -7, -2, -4, 4, 2, -2,
  1072. 2, -2, 1, -1, 8, -8, 3, -3, 1, -1, -5, 5, 9, 1, -9, 1,
  1073. -1, -1, 1, -1, -4, 4, 1, -1, 3, -3, 1, -10, 10, 1, 2, -1,
  1074. -1, -2, 6, -6, 2, 11, -11, -2, 3, -3, 1, -4, 4, -1, 3, -3,
  1075. 1, 3, 12, -3, -5, -12, -1, 5, 2, -2, 1, -1, -7, 1, 13, 7,
  1076. -1, -13, 2, -2, 4, -4, 1, 2, -2, -1, 1, 14, -14, 1, 1, 1,
  1077. -1, -5, -1, -1, 5, -1, -6, 2, -15, 15, 6, 1, -1, -8, 8, -2,
  1078. -4, 4, 1, 1, -1, -1, 16, 2, -16, -2, 2, -2, 4, 3, -4, -3,
  1079. -1, -4, 4, 1, -17, 17, -1, -9, 1, 1, 9, 1, -5, -1, -1, 5,
  1080. -7, 7, 6, -6, 3, -3, 18, -18, 19, -19, 1, -10, -1, 10, -5, 5,
  1081. 20, -20, -3, 1, 3, 1, 8, -1, -8, 2, 7, -1, -21, -2, 5, 21,
  1082. 5, -1, -7, -5, 1, -6, -5, -11, 6, 22, 11, 1, 1, -22, -3, -1,
  1083. 3, -1, 3, -3, -23, 4, -4, 1, 23, -1, 1, -1, 1, -2, 2, -1}
  1084. },{
  1085. /* MapTab8 */
  1086. 4, /* eob_sym */
  1087. 11, /* esc_sym */
  1088. /* run table */
  1089. {1, 1, 1, 1, 0, 2, 2, 1, 1, 3, 3, 0, 1, 1, 2, 2,
  1090. 4, 4, 1, 1, 5, 5, 1, 1, 2, 2, 3, 3, 6, 6, 1, 1,
  1091. 7, 7, 8, 1, 8, 2, 2, 1, 4, 4, 1, 3, 1, 3, 9, 9,
  1092. 2, 2, 1, 5, 1, 5, 10, 10, 1, 1, 11, 11, 3, 6, 3, 4,
  1093. 4, 6, 2, 2, 1, 12, 1, 12, 7, 13, 7, 13, 1, 1, 8, 8,
  1094. 2, 2, 14, 14, 16, 15, 16, 5, 5, 1, 3, 15, 1, 3, 4, 4,
  1095. 1, 1, 17, 17, 2, 2, 6, 6, 1, 18, 1, 18, 22, 21, 22, 21,
  1096. 25, 24, 25, 19, 9, 20, 9, 23, 19, 24, 20, 3, 23, 7, 3, 1,
  1097. 1, 7, 28, 26, 29, 5, 28, 26, 5, 8, 29, 4, 8, 27, 2, 2,
  1098. 4, 27, 1, 1, 10, 36, 10, 33, 33, 36, 30, 1, 32, 32, 1, 30,
  1099. 6, 31, 31, 35, 3, 6, 11, 11, 3, 2, 35, 2, 34, 1, 34, 1,
  1100. 37, 37, 12, 7, 12, 5, 41, 5, 4, 7, 1, 8, 13, 4, 1, 41,
  1101. 13, 38, 8, 38, 9, 1, 40, 40, 9, 1, 39, 2, 2, 49, 39, 42,
  1102. 3, 3, 14, 16, 49, 14, 16, 42, 43, 43, 6, 6, 15, 1, 1, 15,
  1103. 44, 44, 1, 1, 50, 48, 4, 5, 4, 7, 5, 2, 10, 10, 48, 7,
  1104. 50, 45, 2, 1, 45, 8, 8, 1, 46, 46, 3, 47, 47, 3, 1, 1},
  1105. /* value table */
  1106. { 1, -1, 2, -2, 0, 1, -1, 3, -3, 1, -1, 0, 4, -4, 2, -2,
  1107. 1, -1, 5, -5, 1, -1, 6, -6, 3, -3, 2, -2, 1, -1, 7, -7,
  1108. 1, -1, 1, 8, -1, 4, -4, -8, 2, -2, 9, 3, -9, -3, 1, -1,
  1109. 5, -5, 10, 2, -10, -2, 1, -1, 11, -11, 1, -1, -4, 2, 4, 3,
  1110. -3, -2, 6, -6, 12, 1, -12, -1, 2, 1, -2, -1, 13, -13, 2, -2,
  1111. 7, -7, 1, -1, 1, 1, -1, 3, -3, 14, 5, -1, -14, -5, 4, -4,
  1112. 15, -15, 1, -1, 8, -8, -3, 3, 16, 1, -16, -1, 1, 1, -1, -1,
  1113. 1, 1, -1, 1, 2, 1, -2, 1, -1, -1, -1, 6, -1, 3, -6, 17,
  1114. -17, -3, 1, 1, 1, 4, -1, -1, -4, 3, -1, 5, -3, -1, -9, 9,
  1115. -5, 1, 18, -18, 2, 1, -2, 1, -1, -1, 1, 19, -1, 1, -19, -1,
  1116. 4, 1, -1, 1, 7, -4, -2, 2, -7, 10, -1, -10, 1, 20, -1, -20,
  1117. 1, -1, 2, 4, -2, 5, 1, -5, 6, -4, 21, 4, 2, -6, -21, -1,
  1118. -2, 1, -4, -1, -3, 22, -1, 1, 3, -22, -1, 11, -11, 1, 1, 1,
  1119. 8, -8, 2, 2, -1, -2, -2, -1, 1, -1, -5, 5, 2, 23, -23, -2,
  1120. 1, -1, 24, -24, -1, -1, 7, 6, -7, 5, -6, 12, -3, 3, 1, -5,
  1121. 1, 1, -12, 25, -1, -5, 5, -25, -1, 1, 9, 1, -1, -9, 26, -26}
  1122. }
  1123. };