rangecoder.c 3.32 KB
Newer Older
Michael Niedermayer's avatar
Michael Niedermayer committed
1 2 3 4
/*
 * Range coder
 * Copyright (c) 2004 Michael Niedermayer <michaelni@gmx.at>
 *
5 6 7
 * This file is part of FFmpeg.
 *
 * FFmpeg is free software; you can redistribute it and/or
Michael Niedermayer's avatar
Michael Niedermayer committed
8 9
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
10
 * version 2.1 of the License, or (at your option) any later version.
Michael Niedermayer's avatar
Michael Niedermayer committed
11
 *
12
 * FFmpeg is distributed in the hope that it will be useful,
Michael Niedermayer's avatar
Michael Niedermayer committed
13 14 15 16 17
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public
18
 * License along with FFmpeg; if not, write to the Free Software
19
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
Michael Niedermayer's avatar
Michael Niedermayer committed
20
 */
21

Michael Niedermayer's avatar
Michael Niedermayer committed
22
/**
23
 * @file
Michael Niedermayer's avatar
Michael Niedermayer committed
24 25 26 27 28 29 30 31 32 33 34
 * Range coder.
 * based upon
 *    "Range encoding: an algorithm for removing redundancy from a digitised
 *                     message.
 *     G. N. N. Martin                  Presented in March 1979 to the Video &
 *                                      Data Recording Conference,
 *     IBM UK Scientific Center         held in Southampton July 24-27 1979."
 */

#include <string.h>

35
#include "libavutil/attributes.h"
36
#include "libavutil/avassert.h"
37 38
#include "libavutil/intreadwrite.h"

Michael Niedermayer's avatar
Michael Niedermayer committed
39 40 41
#include "avcodec.h"
#include "rangecoder.h"

42
av_cold void ff_init_range_encoder(RangeCoder *c, uint8_t *buf, int buf_size)
43 44 45 46 47 48 49 50
{
    c->bytestream_start  =
    c->bytestream        = buf;
    c->bytestream_end    = buf + buf_size;
    c->low               = 0;
    c->range             = 0xFF00;
    c->outstanding_count = 0;
    c->outstanding_byte  = -1;
Michael Niedermayer's avatar
Michael Niedermayer committed
51 52
}

53 54
av_cold void ff_init_range_decoder(RangeCoder *c, const uint8_t *buf,
                                   int buf_size)
55
{
56
    /* cast to avoid compiler warning */
57
    ff_init_range_encoder(c, (uint8_t *)buf, buf_size);
Michael Niedermayer's avatar
Michael Niedermayer committed
58

59
    c->low         = AV_RB16(c->bytestream);
60
    c->bytestream += 2;
Michael Niedermayer's avatar
Michael Niedermayer committed
61 62
}

63 64 65
void ff_build_rac_states(RangeCoder *c, int factor, int max_p)
{
    const int64_t one = 1LL << 32;
Michael Niedermayer's avatar
Michael Niedermayer committed
66
    int64_t p;
67
    int last_p8, p8, i;
Michael Niedermayer's avatar
Michael Niedermayer committed
68 69

    memset(c->zero_state, 0, sizeof(c->zero_state));
70 71 72 73 74 75 76 77 78 79 80 81 82
    memset(c->one_state, 0, sizeof(c->one_state));

    last_p8 = 0;
    p       = one / 2;
    for (i = 0; i < 128; i++) {
        p8 = (256 * p + one / 2) >> 32; // FIXME: try without the one
        if (p8 <= last_p8)
            p8 = last_p8 + 1;
        if (last_p8 && last_p8 < 256 && p8 <= max_p)
            c->one_state[last_p8] = p8;

        p      += ((one - p) * factor + one / 2) >> 32;
        last_p8 = p8;
Michael Niedermayer's avatar
Michael Niedermayer committed
83
    }
84

85 86
    for (i = 256 - max_p; i <= max_p; i++) {
        if (c->one_state[i])
Michael Niedermayer's avatar
Michael Niedermayer committed
87 88
            continue;

89 90 91 92 93 94 95 96
        p  = (i * one + 128) >> 8;
        p += ((one - p) * factor + one / 2) >> 32;
        p8 = (256 * p + one / 2) >> 32; // FIXME: try without the one
        if (p8 <= i)
            p8 = i + 1;
        if (p8 > max_p)
            p8 = max_p;
        c->one_state[i] = p8;
Michael Niedermayer's avatar
Michael Niedermayer committed
97
    }
98

99 100
    for (i = 1; i < 255; i++)
        c->zero_state[i] = 256 - c->one_state[256 - i];
Michael Niedermayer's avatar
Michael Niedermayer committed
101 102
}

103 104 105 106 107
/* Return the number of bytes written. */
int ff_rac_terminate(RangeCoder *c)
{
    c->range = 0xFF;
    c->low  += 0xFF;
Michael Niedermayer's avatar
Michael Niedermayer committed
108
    renorm_encoder(c);
109
    c->range = 0xFF;
Michael Niedermayer's avatar
Michael Niedermayer committed
110 111
    renorm_encoder(c);

112 113
    av_assert1(c->low   == 0);
    av_assert1(c->range >= 0x100);
Michael Niedermayer's avatar
Michael Niedermayer committed
114 115 116

    return c->bytestream - c->bytestream_start;
}