xref: /dflybsd-src/contrib/xz/src/liblzma/rangecoder/range_decoder.h (revision 2940b44dd4f550115e00e03e0f6149cebf62ded6)
1*2940b44dSPeter Avalos ///////////////////////////////////////////////////////////////////////////////
2*2940b44dSPeter Avalos //
3*2940b44dSPeter Avalos /// \file       range_decoder.h
4*2940b44dSPeter Avalos /// \brief      Range Decoder
5*2940b44dSPeter Avalos ///
6*2940b44dSPeter Avalos //  Authors:    Igor Pavlov
7*2940b44dSPeter Avalos //              Lasse Collin
8*2940b44dSPeter Avalos //
9*2940b44dSPeter Avalos //  This file has been put into the public domain.
10*2940b44dSPeter Avalos //  You can do whatever you want with this file.
11*2940b44dSPeter Avalos //
12*2940b44dSPeter Avalos ///////////////////////////////////////////////////////////////////////////////
13*2940b44dSPeter Avalos 
14*2940b44dSPeter Avalos #ifndef LZMA_RANGE_DECODER_H
15*2940b44dSPeter Avalos #define LZMA_RANGE_DECODER_H
16*2940b44dSPeter Avalos 
17*2940b44dSPeter Avalos #include "range_common.h"
18*2940b44dSPeter Avalos 
19*2940b44dSPeter Avalos 
20*2940b44dSPeter Avalos typedef struct {
21*2940b44dSPeter Avalos 	uint32_t range;
22*2940b44dSPeter Avalos 	uint32_t code;
23*2940b44dSPeter Avalos 	uint32_t init_bytes_left;
24*2940b44dSPeter Avalos } lzma_range_decoder;
25*2940b44dSPeter Avalos 
26*2940b44dSPeter Avalos 
27*2940b44dSPeter Avalos /// Reads the first five bytes to initialize the range decoder.
28*2940b44dSPeter Avalos static inline bool
29*2940b44dSPeter Avalos rc_read_init(lzma_range_decoder *rc, const uint8_t *restrict in,
30*2940b44dSPeter Avalos 		size_t *restrict in_pos, size_t in_size)
31*2940b44dSPeter Avalos {
32*2940b44dSPeter Avalos 	while (rc->init_bytes_left > 0) {
33*2940b44dSPeter Avalos 		if (*in_pos == in_size)
34*2940b44dSPeter Avalos 			return false;
35*2940b44dSPeter Avalos 
36*2940b44dSPeter Avalos 		rc->code = (rc->code << 8) | in[*in_pos];
37*2940b44dSPeter Avalos 		++*in_pos;
38*2940b44dSPeter Avalos 		--rc->init_bytes_left;
39*2940b44dSPeter Avalos 	}
40*2940b44dSPeter Avalos 
41*2940b44dSPeter Avalos 	return true;
42*2940b44dSPeter Avalos }
43*2940b44dSPeter Avalos 
44*2940b44dSPeter Avalos 
45*2940b44dSPeter Avalos /// Makes local copies of range decoder and *in_pos variables. Doing this
46*2940b44dSPeter Avalos /// improves speed significantly. The range decoder macros expect also
47*2940b44dSPeter Avalos /// variables `in' and `in_size' to be defined.
48*2940b44dSPeter Avalos #define rc_to_local(range_decoder, in_pos) \
49*2940b44dSPeter Avalos 	lzma_range_decoder rc = range_decoder; \
50*2940b44dSPeter Avalos 	size_t rc_in_pos = (in_pos); \
51*2940b44dSPeter Avalos 	uint32_t rc_bound
52*2940b44dSPeter Avalos 
53*2940b44dSPeter Avalos 
54*2940b44dSPeter Avalos /// Stores the local copes back to the range decoder structure.
55*2940b44dSPeter Avalos #define rc_from_local(range_decoder, in_pos) \
56*2940b44dSPeter Avalos do { \
57*2940b44dSPeter Avalos 	range_decoder = rc; \
58*2940b44dSPeter Avalos 	in_pos = rc_in_pos; \
59*2940b44dSPeter Avalos } while (0)
60*2940b44dSPeter Avalos 
61*2940b44dSPeter Avalos 
62*2940b44dSPeter Avalos /// Resets the range decoder structure.
63*2940b44dSPeter Avalos #define rc_reset(range_decoder) \
64*2940b44dSPeter Avalos do { \
65*2940b44dSPeter Avalos 	(range_decoder).range = UINT32_MAX; \
66*2940b44dSPeter Avalos 	(range_decoder).code = 0; \
67*2940b44dSPeter Avalos 	(range_decoder).init_bytes_left = 5; \
68*2940b44dSPeter Avalos } while (0)
69*2940b44dSPeter Avalos 
70*2940b44dSPeter Avalos 
71*2940b44dSPeter Avalos /// When decoding has been properly finished, rc.code is always zero unless
72*2940b44dSPeter Avalos /// the input stream is corrupt. So checking this can catch some corrupt
73*2940b44dSPeter Avalos /// files especially if they don't have any other integrity check.
74*2940b44dSPeter Avalos #define rc_is_finished(range_decoder) \
75*2940b44dSPeter Avalos 	((range_decoder).code == 0)
76*2940b44dSPeter Avalos 
77*2940b44dSPeter Avalos 
78*2940b44dSPeter Avalos /// Read the next input byte if needed. If more input is needed but there is
79*2940b44dSPeter Avalos /// no more input available, "goto out" is used to jump out of the main
80*2940b44dSPeter Avalos /// decoder loop.
81*2940b44dSPeter Avalos #define rc_normalize(seq) \
82*2940b44dSPeter Avalos do { \
83*2940b44dSPeter Avalos 	if (rc.range < RC_TOP_VALUE) { \
84*2940b44dSPeter Avalos 		if (unlikely(rc_in_pos == in_size)) { \
85*2940b44dSPeter Avalos 			coder->sequence = seq; \
86*2940b44dSPeter Avalos 			goto out; \
87*2940b44dSPeter Avalos 		} \
88*2940b44dSPeter Avalos 		rc.range <<= RC_SHIFT_BITS; \
89*2940b44dSPeter Avalos 		rc.code = (rc.code << RC_SHIFT_BITS) | in[rc_in_pos++]; \
90*2940b44dSPeter Avalos 	} \
91*2940b44dSPeter Avalos } while (0)
92*2940b44dSPeter Avalos 
93*2940b44dSPeter Avalos 
94*2940b44dSPeter Avalos /// Start decoding a bit. This must be used together with rc_update_0()
95*2940b44dSPeter Avalos /// and rc_update_1():
96*2940b44dSPeter Avalos ///
97*2940b44dSPeter Avalos ///     rc_if_0(prob, seq) {
98*2940b44dSPeter Avalos ///         rc_update_0(prob);
99*2940b44dSPeter Avalos ///         // Do something
100*2940b44dSPeter Avalos ///     } else {
101*2940b44dSPeter Avalos ///         rc_update_1(prob);
102*2940b44dSPeter Avalos ///         // Do something else
103*2940b44dSPeter Avalos ///     }
104*2940b44dSPeter Avalos ///
105*2940b44dSPeter Avalos #define rc_if_0(prob, seq) \
106*2940b44dSPeter Avalos 	rc_normalize(seq); \
107*2940b44dSPeter Avalos 	rc_bound = (rc.range >> RC_BIT_MODEL_TOTAL_BITS) * (prob); \
108*2940b44dSPeter Avalos 	if (rc.code < rc_bound)
109*2940b44dSPeter Avalos 
110*2940b44dSPeter Avalos 
111*2940b44dSPeter Avalos /// Update the range decoder state and the used probability variable to
112*2940b44dSPeter Avalos /// match a decoded bit of 0.
113*2940b44dSPeter Avalos #define rc_update_0(prob) \
114*2940b44dSPeter Avalos do { \
115*2940b44dSPeter Avalos 	rc.range = rc_bound; \
116*2940b44dSPeter Avalos 	prob += (RC_BIT_MODEL_TOTAL - (prob)) >> RC_MOVE_BITS; \
117*2940b44dSPeter Avalos } while (0)
118*2940b44dSPeter Avalos 
119*2940b44dSPeter Avalos 
120*2940b44dSPeter Avalos /// Update the range decoder state and the used probability variable to
121*2940b44dSPeter Avalos /// match a decoded bit of 1.
122*2940b44dSPeter Avalos #define rc_update_1(prob) \
123*2940b44dSPeter Avalos do { \
124*2940b44dSPeter Avalos 	rc.range -= rc_bound; \
125*2940b44dSPeter Avalos 	rc.code -= rc_bound; \
126*2940b44dSPeter Avalos 	prob -= (prob) >> RC_MOVE_BITS; \
127*2940b44dSPeter Avalos } while (0)
128*2940b44dSPeter Avalos 
129*2940b44dSPeter Avalos 
130*2940b44dSPeter Avalos /// Decodes one bit and runs action0 or action1 depending on the decoded bit.
131*2940b44dSPeter Avalos /// This macro is used as the last step in bittree reverse decoders since
132*2940b44dSPeter Avalos /// those don't use "symbol" for anything else than indexing the probability
133*2940b44dSPeter Avalos /// arrays.
134*2940b44dSPeter Avalos #define rc_bit_last(prob, action0, action1, seq) \
135*2940b44dSPeter Avalos do { \
136*2940b44dSPeter Avalos 	rc_if_0(prob, seq) { \
137*2940b44dSPeter Avalos 		rc_update_0(prob); \
138*2940b44dSPeter Avalos 		action0; \
139*2940b44dSPeter Avalos 	} else { \
140*2940b44dSPeter Avalos 		rc_update_1(prob); \
141*2940b44dSPeter Avalos 		action1; \
142*2940b44dSPeter Avalos 	} \
143*2940b44dSPeter Avalos } while (0)
144*2940b44dSPeter Avalos 
145*2940b44dSPeter Avalos 
146*2940b44dSPeter Avalos /// Decodes one bit, updates "symbol", and runs action0 or action1 depending
147*2940b44dSPeter Avalos /// on the decoded bit.
148*2940b44dSPeter Avalos #define rc_bit(prob, action0, action1, seq) \
149*2940b44dSPeter Avalos 	rc_bit_last(prob, \
150*2940b44dSPeter Avalos 		symbol <<= 1; action0, \
151*2940b44dSPeter Avalos 		symbol = (symbol << 1) + 1; action1, \
152*2940b44dSPeter Avalos 		seq);
153*2940b44dSPeter Avalos 
154*2940b44dSPeter Avalos 
155*2940b44dSPeter Avalos /// Like rc_bit() but add "case seq:" as a prefix. This makes the unrolled
156*2940b44dSPeter Avalos /// loops more readable because the code isn't littered with "case"
157*2940b44dSPeter Avalos /// statements. On the other hand this also makes it less readable, since
158*2940b44dSPeter Avalos /// spotting the places where the decoder loop may be restarted is less
159*2940b44dSPeter Avalos /// obvious.
160*2940b44dSPeter Avalos #define rc_bit_case(prob, action0, action1, seq) \
161*2940b44dSPeter Avalos 	case seq: rc_bit(prob, action0, action1, seq)
162*2940b44dSPeter Avalos 
163*2940b44dSPeter Avalos 
164*2940b44dSPeter Avalos /// Decode a bit without using a probability.
165*2940b44dSPeter Avalos #define rc_direct(dest, seq) \
166*2940b44dSPeter Avalos do { \
167*2940b44dSPeter Avalos 	rc_normalize(seq); \
168*2940b44dSPeter Avalos 	rc.range >>= 1; \
169*2940b44dSPeter Avalos 	rc.code -= rc.range; \
170*2940b44dSPeter Avalos 	rc_bound = UINT32_C(0) - (rc.code >> 31); \
171*2940b44dSPeter Avalos 	rc.code += rc.range & rc_bound; \
172*2940b44dSPeter Avalos 	dest = (dest << 1) + (rc_bound + 1); \
173*2940b44dSPeter Avalos } while (0)
174*2940b44dSPeter Avalos 
175*2940b44dSPeter Avalos 
176*2940b44dSPeter Avalos // NOTE: No macros are provided for bittree decoding. It seems to be simpler
177*2940b44dSPeter Avalos // to just write them open in the code.
178*2940b44dSPeter Avalos 
179*2940b44dSPeter Avalos #endif
180