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