186d7f5d3SJohn Marino ///////////////////////////////////////////////////////////////////////////////
286d7f5d3SJohn Marino //
386d7f5d3SJohn Marino /// \file price.h
486d7f5d3SJohn Marino /// \brief Probability price calculation
586d7f5d3SJohn Marino //
686d7f5d3SJohn Marino // Author: Igor Pavlov
786d7f5d3SJohn Marino //
886d7f5d3SJohn Marino // This file has been put into the public domain.
986d7f5d3SJohn Marino // You can do whatever you want with this file.
1086d7f5d3SJohn Marino //
1186d7f5d3SJohn Marino ///////////////////////////////////////////////////////////////////////////////
1286d7f5d3SJohn Marino
1386d7f5d3SJohn Marino #ifndef LZMA_PRICE_H
1486d7f5d3SJohn Marino #define LZMA_PRICE_H
1586d7f5d3SJohn Marino
1686d7f5d3SJohn Marino
1786d7f5d3SJohn Marino #define RC_MOVE_REDUCING_BITS 4
1886d7f5d3SJohn Marino #define RC_BIT_PRICE_SHIFT_BITS 4
1986d7f5d3SJohn Marino #define RC_PRICE_TABLE_SIZE (RC_BIT_MODEL_TOTAL >> RC_MOVE_REDUCING_BITS)
2086d7f5d3SJohn Marino
2186d7f5d3SJohn Marino #define RC_INFINITY_PRICE (UINT32_C(1) << 30)
2286d7f5d3SJohn Marino
2386d7f5d3SJohn Marino
2486d7f5d3SJohn Marino /// Lookup table for the inline functions defined in this file.
2586d7f5d3SJohn Marino extern const uint8_t lzma_rc_prices[RC_PRICE_TABLE_SIZE];
2686d7f5d3SJohn Marino
2786d7f5d3SJohn Marino
2886d7f5d3SJohn Marino static inline uint32_t
rc_bit_price(const probability prob,const uint32_t bit)2986d7f5d3SJohn Marino rc_bit_price(const probability prob, const uint32_t bit)
3086d7f5d3SJohn Marino {
3186d7f5d3SJohn Marino return lzma_rc_prices[(prob ^ ((UINT32_C(0) - bit)
3286d7f5d3SJohn Marino & (RC_BIT_MODEL_TOTAL - 1))) >> RC_MOVE_REDUCING_BITS];
3386d7f5d3SJohn Marino }
3486d7f5d3SJohn Marino
3586d7f5d3SJohn Marino
3686d7f5d3SJohn Marino static inline uint32_t
rc_bit_0_price(const probability prob)3786d7f5d3SJohn Marino rc_bit_0_price(const probability prob)
3886d7f5d3SJohn Marino {
3986d7f5d3SJohn Marino return lzma_rc_prices[prob >> RC_MOVE_REDUCING_BITS];
4086d7f5d3SJohn Marino }
4186d7f5d3SJohn Marino
4286d7f5d3SJohn Marino
4386d7f5d3SJohn Marino static inline uint32_t
rc_bit_1_price(const probability prob)4486d7f5d3SJohn Marino rc_bit_1_price(const probability prob)
4586d7f5d3SJohn Marino {
4686d7f5d3SJohn Marino return lzma_rc_prices[(prob ^ (RC_BIT_MODEL_TOTAL - 1))
4786d7f5d3SJohn Marino >> RC_MOVE_REDUCING_BITS];
4886d7f5d3SJohn Marino }
4986d7f5d3SJohn Marino
5086d7f5d3SJohn Marino
5186d7f5d3SJohn Marino static inline uint32_t
rc_bittree_price(const probability * const probs,const uint32_t bit_levels,uint32_t symbol)5286d7f5d3SJohn Marino rc_bittree_price(const probability *const probs,
5386d7f5d3SJohn Marino const uint32_t bit_levels, uint32_t symbol)
5486d7f5d3SJohn Marino {
5586d7f5d3SJohn Marino uint32_t price = 0;
5686d7f5d3SJohn Marino symbol += UINT32_C(1) << bit_levels;
5786d7f5d3SJohn Marino
5886d7f5d3SJohn Marino do {
5986d7f5d3SJohn Marino const uint32_t bit = symbol & 1;
6086d7f5d3SJohn Marino symbol >>= 1;
6186d7f5d3SJohn Marino price += rc_bit_price(probs[symbol], bit);
6286d7f5d3SJohn Marino } while (symbol != 1);
6386d7f5d3SJohn Marino
6486d7f5d3SJohn Marino return price;
6586d7f5d3SJohn Marino }
6686d7f5d3SJohn Marino
6786d7f5d3SJohn Marino
6886d7f5d3SJohn Marino static inline uint32_t
rc_bittree_reverse_price(const probability * const probs,uint32_t bit_levels,uint32_t symbol)6986d7f5d3SJohn Marino rc_bittree_reverse_price(const probability *const probs,
7086d7f5d3SJohn Marino uint32_t bit_levels, uint32_t symbol)
7186d7f5d3SJohn Marino {
7286d7f5d3SJohn Marino uint32_t price = 0;
7386d7f5d3SJohn Marino uint32_t model_index = 1;
7486d7f5d3SJohn Marino
7586d7f5d3SJohn Marino do {
7686d7f5d3SJohn Marino const uint32_t bit = symbol & 1;
7786d7f5d3SJohn Marino symbol >>= 1;
7886d7f5d3SJohn Marino price += rc_bit_price(probs[model_index], bit);
7986d7f5d3SJohn Marino model_index = (model_index << 1) + bit;
8086d7f5d3SJohn Marino } while (--bit_levels != 0);
8186d7f5d3SJohn Marino
8286d7f5d3SJohn Marino return price;
8386d7f5d3SJohn Marino }
8486d7f5d3SJohn Marino
8586d7f5d3SJohn Marino
8686d7f5d3SJohn Marino static inline uint32_t
rc_direct_price(const uint32_t bits)8786d7f5d3SJohn Marino rc_direct_price(const uint32_t bits)
8886d7f5d3SJohn Marino {
8986d7f5d3SJohn Marino return bits << RC_BIT_PRICE_SHIFT_BITS;
9086d7f5d3SJohn Marino }
9186d7f5d3SJohn Marino
9286d7f5d3SJohn Marino #endif
93