xref: /dflybsd-src/contrib/xz/src/liblzma/rangecoder/price.h (revision 86d7f5d305c6adaa56ff4582ece9859d73106103)
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