yading@10: /* yading@10: * Copyright (c) 1990 James Ashton - Sydney University yading@10: * Copyright (c) 2012 Stefano Sabatini yading@10: * yading@10: * This file is part of FFmpeg. yading@10: * yading@10: * FFmpeg is free software; you can redistribute it and/or yading@10: * modify it under the terms of the GNU Lesser General Public yading@10: * License as published by the Free Software Foundation; either yading@10: * version 2.1 of the License, or (at your option) any later version. yading@10: * yading@10: * FFmpeg is distributed in the hope that it will be useful, yading@10: * but WITHOUT ANY WARRANTY; without even the implied warranty of yading@10: * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU yading@10: * Lesser General Public License for more details. yading@10: * yading@10: * You should have received a copy of the GNU Lesser General Public yading@10: * License along with FFmpeg; if not, write to the Free Software yading@10: * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA yading@10: */ yading@10: yading@10: /** yading@10: * @file yading@10: * X-Face common data and utilities definition. yading@10: */ yading@10: yading@10: #include "xface.h" yading@10: yading@10: void ff_big_add(BigInt *b, uint8_t a) yading@10: { yading@10: int i; yading@10: uint8_t *w; yading@10: uint16_t c; yading@10: yading@10: a &= XFACE_WORDMASK; yading@10: if (a == 0) yading@10: return; yading@10: w = b->words; yading@10: c = a; yading@10: for (i = 0; i < b->nb_words && c; i++) { yading@10: c += *w; yading@10: *w++ = c & XFACE_WORDMASK; yading@10: c >>= XFACE_BITSPERWORD; yading@10: } yading@10: if (i == b->nb_words && c) { yading@10: b->nb_words++; yading@10: *w = c & XFACE_WORDMASK; yading@10: } yading@10: } yading@10: yading@10: void ff_big_div(BigInt *b, uint8_t a, uint8_t *r) yading@10: { yading@10: int i; yading@10: uint8_t *w; yading@10: uint16_t c, d; yading@10: yading@10: a &= XFACE_WORDMASK; yading@10: if (a == 1 || b->nb_words == 0) { yading@10: *r = 0; yading@10: return; yading@10: } yading@10: yading@10: /* treat this as a == WORDCARRY and just shift everything right a WORD */ yading@10: if (a == 0) { yading@10: i = --b->nb_words; yading@10: w = b->words; yading@10: *r = *w; yading@10: while (i--) { yading@10: *w = *(w + 1); yading@10: w++; yading@10: } yading@10: *w = 0; yading@10: return; yading@10: } yading@10: i = b->nb_words; yading@10: w = b->words + i; yading@10: c = 0; yading@10: while (i--) { yading@10: c <<= XFACE_BITSPERWORD; yading@10: c += *--w; yading@10: d = c / (uint16_t)a; yading@10: c = c % (uint16_t)a; yading@10: *w = d & XFACE_WORDMASK; yading@10: } yading@10: *r = c; yading@10: if (b->words[b->nb_words - 1] == 0) yading@10: b->nb_words--; yading@10: } yading@10: yading@10: void ff_big_mul(BigInt *b, uint8_t a) yading@10: { yading@10: int i; yading@10: uint8_t *w; yading@10: uint16_t c; yading@10: yading@10: a &= XFACE_WORDMASK; yading@10: if (a == 1 || b->nb_words == 0) yading@10: return; yading@10: if (a == 0) { yading@10: /* treat this as a == WORDCARRY and just shift everything left a WORD */ yading@10: i = b->nb_words++; yading@10: w = b->words + i; yading@10: while (i--) { yading@10: *w = *(w - 1); yading@10: w--; yading@10: } yading@10: *w = 0; yading@10: return; yading@10: } yading@10: i = b->nb_words; yading@10: w = b->words; yading@10: c = 0; yading@10: while (i--) { yading@10: c += (uint16_t)*w * (uint16_t)a; yading@10: *(w++) = c & XFACE_WORDMASK; yading@10: c >>= XFACE_BITSPERWORD; yading@10: } yading@10: if (c) { yading@10: b->nb_words++; yading@10: *w = c & XFACE_WORDMASK; yading@10: } yading@10: } yading@10: yading@10: const ProbRange ff_xface_probranges_per_level[4][3] = { yading@10: // black grey white yading@10: { { 1, 255}, {251, 0}, { 4, 251} }, /* Top of tree almost always grey */ yading@10: { { 1, 255}, {200, 0}, { 55, 200} }, yading@10: { { 33, 223}, {159, 0}, { 64, 159} }, yading@10: { {131, 0}, { 0, 0}, {125, 131} }, /* Grey disallowed at bottom */ yading@10: }; yading@10: yading@10: const ProbRange ff_xface_probranges_2x2[16] = { yading@10: { 0, 0}, {38, 0}, {38, 38}, {13, 152}, yading@10: {38, 76}, {13, 165}, {13, 178}, { 6, 230}, yading@10: {38, 114}, {13, 191}, {13, 204}, { 6, 236}, yading@10: {13, 217}, { 6, 242}, { 5, 248}, { 3, 253}, yading@10: }; yading@10: yading@10: /* yading@10: * The "guess the next pixel" tables follow. Normally there are 12 yading@10: * neighbour pixels used to give 1<<12 cases as we get closer to the yading@10: * upper left corner lesser numbers of neighbours are available. yading@10: * yading@10: * Each byte in the tables represents 8 boolean values starting from yading@10: * the most significant bit. yading@10: */ yading@10: yading@10: static const uint8_t g_00[] = { yading@10: 0x00, 0x00, 0x01, 0x01, 0x00, 0x00, 0xe3, 0xdf, 0x05, 0x17, yading@10: 0x05, 0x0f, 0x00, 0x1b, 0x0f, 0xdf, 0x00, 0x04, 0x00, 0x00, yading@10: 0x0d, 0x0f, 0x03, 0x7f, 0x00, 0x00, 0x00, 0x01, 0x00, 0x1d, yading@10: 0x45, 0x2f, 0x00, 0x00, 0x00, 0x0d, 0x00, 0x0a, 0xff, 0xff, yading@10: 0x00, 0x04, 0x00, 0x05, 0x01, 0x3f, 0xcf, 0xff, 0x10, 0x01, yading@10: 0x80, 0xc9, 0x0f, 0x0f, 0xff, 0xff, 0x00, 0x00, 0x00, 0x00, yading@10: 0x1b, 0x1f, 0xff, 0xff, 0x4f, 0x54, 0x07, 0x1f, 0x57, 0x47, yading@10: 0xd7, 0x3d, 0xff, 0xff, 0x5f, 0x1f, 0x7f, 0xff, 0x7f, 0x7f, yading@10: 0x05, 0x0f, 0x01, 0x0f, 0x0f, 0x5f, 0x9b, 0xdf, 0x7f, 0xff, yading@10: 0x5f, 0x1d, 0x5f, 0xff, 0x0f, 0x1f, 0x0f, 0x5f, 0x03, 0x1f, yading@10: 0x4f, 0x5f, 0xf7, 0x7f, 0x7f, 0xff, 0x0d, 0x0f, 0xfb, 0xff, yading@10: 0xf7, 0xbf, 0x0f, 0x4f, 0xd7, 0x3f, 0x4f, 0x7f, 0xff, 0xff, yading@10: 0x67, 0xbf, 0x56, 0x25, 0x1f, 0x7f, 0x9f, 0xff, 0x00, 0x00, yading@10: 0x00, 0x05, 0x5f, 0x7f, 0x01, 0xdf, 0x14, 0x00, 0x05, 0x0f, yading@10: 0x07, 0xa2, 0x09, 0x0f, 0x00, 0x00, 0x00, 0x00, 0x0f, 0x5f, yading@10: 0x18, 0xd7, 0x94, 0x71, 0x00, 0x05, 0x1f, 0xb7, 0x0c, 0x07, yading@10: 0x0f, 0x0f, 0x00, 0x0f, 0x0f, 0x1f, 0x84, 0x8f, 0x05, 0x15, yading@10: 0x05, 0x0f, 0x4f, 0xff, 0x87, 0xdf, 0x05, 0x01, 0x10, 0x00, yading@10: 0x0f, 0x0f, 0x00, 0x08, 0x05, 0x04, 0x04, 0x01, 0x4f, 0xff, yading@10: 0x9f, 0x8f, 0x4a, 0x40, 0x5f, 0x5f, 0xff, 0xfe, 0xdf, 0xff, yading@10: 0x7f, 0xf7, 0xff, 0x7f, 0xff, 0xff, 0x7b, 0xff, 0x0f, 0xfd, yading@10: 0xd7, 0x5f, 0x4f, 0x7f, 0x7f, 0xdf, 0xff, 0xff, 0xff, 0xff, yading@10: 0xff, 0x77, 0xdf, 0x7f, 0x4f, 0xef, 0xff, 0xff, 0x77, 0xff, yading@10: 0xff, 0xff, 0x6f, 0xff, 0x0f, 0x4f, 0xff, 0xff, 0x9d, 0xff, yading@10: 0x0f, 0xef, 0xff, 0xdf, 0x6f, 0xff, 0xff, 0xff, 0x4f, 0xff, yading@10: 0xcd, 0x0f, 0x4f, 0xff, 0xff, 0xdf, 0x00, 0x00, 0x00, 0x0b, yading@10: 0x05, 0x02, 0x02, 0x0f, 0x04, 0x00, 0x00, 0x0c, 0x01, 0x06, yading@10: 0x00, 0x0f, 0x20, 0x03, 0x00, 0x00, 0x05, 0x0f, 0x40, 0x08, yading@10: 0x00, 0x00, 0x00, 0x01, 0x00, 0x01, 0x0c, 0x0f, 0x01, 0x00, yading@10: 0x80, 0x00, 0x00, 0x00, 0x80, 0x00, 0x00, 0x14, 0x01, 0x05, yading@10: 0x01, 0x15, 0xaf, 0x0f, 0x00, 0x01, 0x10, 0x00, 0x08, 0x00, yading@10: 0x46, 0x0c, 0x20, 0x00, 0x88, 0x00, 0x0f, 0x15, 0xff, 0xdf, yading@10: 0x02, 0x00, 0x00, 0x0f, 0x7f, 0x5f, 0xdb, 0xff, 0x4f, 0x3e, yading@10: 0x05, 0x0f, 0x7f, 0xf7, 0x95, 0x4f, 0x0d, 0x0f, 0x01, 0x0f, yading@10: 0x4f, 0x5f, 0x9f, 0xdf, 0x25, 0x0e, 0x0d, 0x0d, 0x4f, 0x7f, yading@10: 0x8f, 0x0f, 0x0f, 0xfa, 0x04, 0x4f, 0x4f, 0xff, 0xf7, 0x77, yading@10: 0x47, 0xed, 0x05, 0x0f, 0xff, 0xff, 0xdf, 0xff, 0x4f, 0x6f, yading@10: 0xd8, 0x5f, 0x0f, 0x7f, 0xdf, 0x5f, 0x07, 0x0f, 0x94, 0x0d, yading@10: 0x1f, 0xff, 0xff, 0xff, 0x00, 0x02, 0x00, 0x03, 0x46, 0x57, yading@10: 0x01, 0x0d, 0x01, 0x08, 0x01, 0x0f, 0x47, 0x6c, 0x0d, 0x0f, yading@10: 0x02, 0x00, 0x00, 0x00, 0x0b, 0x4f, 0x00, 0x08, 0x05, 0x00, yading@10: 0x95, 0x01, 0x0f, 0x7f, 0x0c, 0x0f, 0x01, 0x0e, 0x00, 0x00, yading@10: 0x0f, 0x41, 0x00, 0x00, 0x04, 0x24, 0x0d, 0x0f, 0x0f, 0x7f, yading@10: 0xcf, 0xdf, 0x00, 0x00, 0x00, 0x00, 0x04, 0x40, 0x00, 0x00, yading@10: 0x06, 0x26, 0xcf, 0x05, 0xcf, 0x7f, 0xdf, 0xdf, 0x00, 0x00, yading@10: 0x17, 0x5f, 0xff, 0xfd, 0xff, 0xff, 0x46, 0x09, 0x4f, 0x5f, yading@10: 0x7f, 0xfd, 0xdf, 0xff, 0x0a, 0x88, 0xa7, 0x7f, 0x7f, 0xff, yading@10: 0xff, 0xff, 0x0f, 0x04, 0xdf, 0x7f, 0x4f, 0xff, 0x9f, 0xff, yading@10: 0x0e, 0xe6, 0xdf, 0xff, 0x7f, 0xff, 0xff, 0xff, 0x0f, 0xec, yading@10: 0x8f, 0x4f, 0x7f, 0xff, 0xdf, 0xff, 0x0f, 0xcf, 0xdf, 0xff, yading@10: 0x6f, 0x7f, 0xff, 0xff, 0x03, 0x0c, 0x9d, 0x0f, 0x7f, 0xff, yading@10: 0xff, 0xff, yading@10: }; yading@10: yading@10: static const uint8_t g_01[] = { yading@10: 0x37, 0x73, 0x00, 0x19, 0x57, 0x7f, 0xf5, 0xfb, 0x70, 0x33, yading@10: 0xf0, 0xf9, 0x7f, 0xff, 0xff, 0xff, yading@10: }; yading@10: yading@10: static const uint8_t g_02[] = { yading@10: 0x50, yading@10: }; yading@10: yading@10: static const uint8_t g_10[] = { yading@10: 0x00, 0x00, 0x00, 0x00, 0x50, 0x00, 0xf3, 0x5f, 0x84, 0x04, yading@10: 0x17, 0x9f, 0x04, 0x23, 0x05, 0xff, 0x00, 0x00, 0x00, 0x02, yading@10: 0x03, 0x03, 0x33, 0xd7, 0x05, 0x03, 0x5f, 0x3f, 0x17, 0x33, yading@10: 0xff, 0xff, 0x00, 0x80, 0x02, 0x04, 0x12, 0x00, 0x11, 0x57, yading@10: 0x05, 0x25, 0x05, 0x03, 0x35, 0xbf, 0x9f, 0xff, 0x07, 0x6f, yading@10: 0x20, 0x40, 0x17, 0x06, 0xfa, 0xe8, 0x01, 0x07, 0x1f, 0x9f, yading@10: 0x1f, 0xff, 0xff, 0xff, yading@10: }; yading@10: yading@10: static const uint8_t g_20[] = { yading@10: 0x04, 0x00, 0x01, 0x01, 0x43, 0x2e, 0xff, 0x3f, yading@10: }; yading@10: yading@10: static const uint8_t g_30[] = { yading@10: 0x11, 0x11, 0x11, 0x11, 0x51, 0x11, 0x13, 0x11, 0x11, 0x11, yading@10: 0x13, 0x11, 0x11, 0x11, 0x33, 0x11, 0x13, 0x11, 0x13, 0x13, yading@10: 0x13, 0x13, 0x31, 0x31, 0x11, 0x01, 0x11, 0x11, 0x71, 0x11, yading@10: 0x11, 0x75, yading@10: }; yading@10: yading@10: static const uint8_t g_40[] = { yading@10: 0x00, 0x0f, 0x00, 0x09, 0x00, 0x0d, 0x00, 0x0d, 0x00, 0x0f, yading@10: 0x00, 0x4e, 0xe4, 0x0d, 0x10, 0x0f, 0x00, 0x0f, 0x44, 0x4f, yading@10: 0x00, 0x1e, 0x0f, 0x0f, 0xae, 0xaf, 0x45, 0x7f, 0xef, 0xff, yading@10: 0x0f, 0xff, 0x00, 0x09, 0x01, 0x11, 0x00, 0x01, 0x1c, 0xdd, yading@10: 0x00, 0x15, 0x00, 0xff, 0x00, 0x10, 0x00, 0xfd, 0x00, 0x0f, yading@10: 0x4f, 0x5f, 0x3d, 0xff, 0xff, 0xff, 0x4f, 0xff, 0x1c, 0xff, yading@10: 0xdf, 0xff, 0x8f, 0xff, 0x00, 0x0d, 0x00, 0x00, 0x00, 0x15, yading@10: 0x01, 0x07, 0x00, 0x01, 0x02, 0x1f, 0x01, 0x11, 0x05, 0x7f, yading@10: 0x00, 0x1f, 0x41, 0x57, 0x1f, 0xff, 0x05, 0x77, 0x0d, 0x5f, yading@10: 0x4d, 0xff, 0x4f, 0xff, 0x0f, 0xff, 0x00, 0x00, 0x02, 0x05, yading@10: 0x00, 0x11, 0x05, 0x7d, 0x10, 0x15, 0x2f, 0xff, 0x40, 0x50, yading@10: 0x0d, 0xfd, 0x04, 0x0f, 0x07, 0x1f, 0x07, 0x7f, 0x0f, 0xbf, yading@10: 0x0d, 0x7f, 0x0f, 0xff, 0x4d, 0x7d, 0x0f, 0xff, yading@10: }; yading@10: yading@10: static const uint8_t g_11[] = { yading@10: 0x01, 0x13, 0x03, 0x7f, yading@10: }; yading@10: yading@10: static const uint8_t g_21[] = { yading@10: 0x17, yading@10: }; yading@10: yading@10: static const uint8_t g_31[] = { yading@10: 0x55, 0x57, 0x57, 0x7f, yading@10: }; yading@10: yading@10: static const uint8_t g_41[] = { yading@10: 0x01, 0x01, 0x01, 0x1f, 0x03, 0x1f, 0x3f, 0xff, yading@10: }; yading@10: yading@10: static const uint8_t g_12[] = { yading@10: 0x40, yading@10: }; yading@10: yading@10: static const uint8_t g_22[] = { yading@10: 0x00, yading@10: }; yading@10: yading@10: static const uint8_t g_32[] = { yading@10: 0x10, yading@10: }; yading@10: yading@10: static const uint8_t g_42[] = { yading@10: 0x10, yading@10: }; yading@10: yading@10: void ff_xface_generate_face(uint8_t *dst, uint8_t * const src) yading@10: { yading@10: int h, i, j, k, l, m; yading@10: yading@10: for (j = 0; j < XFACE_HEIGHT; j++) { yading@10: for (i = 0; i < XFACE_WIDTH; i++) { yading@10: h = i + j * XFACE_WIDTH; yading@10: k = 0; yading@10: yading@10: /* yading@10: Compute k, encoding the bits *before* the current one, contained in the yading@10: image buffer. That is, given the grid: yading@10: yading@10: l i yading@10: | | yading@10: v v yading@10: +--+--+--+--+--+ yading@10: m -> | 1| 2| 3| 4| 5| yading@10: +--+--+--+--+--+ yading@10: | 6| 7| 8| 9|10| yading@10: +--+--+--+--+--+ yading@10: j -> |11|12| *| | | yading@10: +--+--+--+--+--+ yading@10: yading@10: the value k for the pixel marked as "*" will contain the bit encoding of yading@10: the values in the matrix marked from "1" to "12". In case the pixel is yading@10: near the border of the grid, the number of values contained within the yading@10: grid will be lesser than 12. yading@10: */ yading@10: yading@10: for (l = i - 2; l <= i + 2; l++) { yading@10: for (m = j - 2; m <= j; m++) { yading@10: if (l >= i && m == j) yading@10: continue; yading@10: if (l > 0 && l <= XFACE_WIDTH && m > 0) yading@10: k = 2*k + src[l + m * XFACE_WIDTH]; yading@10: } yading@10: } yading@10: yading@10: /* yading@10: Use the guess for the given position and the computed value of k. yading@10: yading@10: The following table shows the number of digits in k, depending on yading@10: the position of the pixel, and shows the corresponding guess table yading@10: to use: yading@10: yading@10: i=1 i=2 i=3 i=w-1 i=w yading@10: +----+----+----+ ... +----+----+ yading@10: j=1 | 0 | 1 | 2 | | 2 | 2 | yading@10: |g22 |g12 |g02 | |g42 |g32 | yading@10: +----+----+----+ ... +----+----+ yading@10: j=2 | 3 | 5 | 7 | | 6 | 5 | yading@10: |g21 |g11 |g01 | |g41 |g31 | yading@10: +----+----+----+ ... +----+----+ yading@10: j=3 | 5 | 9 | 12 | | 10 | 8 | yading@10: |g20 |g10 |g00 | |g40 |g30 | yading@10: +----+----+----+ ... +----+----+ yading@10: */ yading@10: yading@10: #define GEN(table) dst[h] ^= (table[k>>3]>>(7-(k&7)))&1 yading@10: yading@10: switch (i) { yading@10: case 1: yading@10: switch (j) { yading@10: case 1: GEN(g_22); break; yading@10: case 2: GEN(g_21); break; yading@10: default: GEN(g_20); break; yading@10: } yading@10: break; yading@10: case 2: yading@10: switch (j) { yading@10: case 1: GEN(g_12); break; yading@10: case 2: GEN(g_11); break; yading@10: default: GEN(g_10); break; yading@10: } yading@10: break; yading@10: case XFACE_WIDTH - 1: yading@10: switch (j) { yading@10: case 1: GEN(g_42); break; yading@10: case 2: GEN(g_41); break; yading@10: default: GEN(g_40); break; yading@10: } yading@10: break; yading@10: case XFACE_WIDTH: yading@10: switch (j) { yading@10: case 1: GEN(g_32); break; yading@10: case 2: GEN(g_31); break; yading@10: default: GEN(g_30); break; yading@10: } yading@10: break; yading@10: default: yading@10: switch (j) { yading@10: case 1: GEN(g_02); break; yading@10: case 2: GEN(g_01); break; yading@10: default: GEN(g_00); break; yading@10: } yading@10: break; yading@10: } yading@10: } yading@10: } yading@10: }