xref: /btstack/3rd-party/lc3-google/src/bits.h (revision 6897da5c53aac5b1f90f41b5b15d0bd43d61dfff)
19a19cd78SMatthias Ringwald /******************************************************************************
29a19cd78SMatthias Ringwald  *
3*4930cef6SMatthias Ringwald  *  Copyright 2022 Google LLC
49a19cd78SMatthias Ringwald  *
59a19cd78SMatthias Ringwald  *  Licensed under the Apache License, Version 2.0 (the "License");
69a19cd78SMatthias Ringwald  *  you may not use this file except in compliance with the License.
79a19cd78SMatthias Ringwald  *  You may obtain a copy of the License at:
89a19cd78SMatthias Ringwald  *
99a19cd78SMatthias Ringwald  *  http://www.apache.org/licenses/LICENSE-2.0
109a19cd78SMatthias Ringwald  *
119a19cd78SMatthias Ringwald  *  Unless required by applicable law or agreed to in writing, software
129a19cd78SMatthias Ringwald  *  distributed under the License is distributed on an "AS IS" BASIS,
139a19cd78SMatthias Ringwald  *  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
149a19cd78SMatthias Ringwald  *  See the License for the specific language governing permissions and
159a19cd78SMatthias Ringwald  *  limitations under the License.
169a19cd78SMatthias Ringwald  *
179a19cd78SMatthias Ringwald  ******************************************************************************/
189a19cd78SMatthias Ringwald 
199a19cd78SMatthias Ringwald /**
209a19cd78SMatthias Ringwald  * The bitstream is written by the 2 ends of the buffer :
219a19cd78SMatthias Ringwald  *
229a19cd78SMatthias Ringwald  * - Arthmetic coder put bits while increasing memory addresses
239a19cd78SMatthias Ringwald  *   in the buffer (forward)
249a19cd78SMatthias Ringwald  *
259a19cd78SMatthias Ringwald  * - Plain bits are puts starting the end of the buffer, with memeory
269a19cd78SMatthias Ringwald  *   addresses decreasing (backward)
279a19cd78SMatthias Ringwald  *
289a19cd78SMatthias Ringwald  *       .---------------------------------------------------.
299a19cd78SMatthias Ringwald  *       | > > > > > > > > > > :         : < < < < < < < < < |
309a19cd78SMatthias Ringwald  *       '---------------------------------------------------'
319a19cd78SMatthias Ringwald  *       |---------------------> - - - - - - - - - - - - - ->|
329a19cd78SMatthias Ringwald  *                              |< - - - <-------------------|
339a19cd78SMatthias Ringwald  *          Arithmetic coding                  Plain bits
349a19cd78SMatthias Ringwald  *          `lc3_put_symbol()`               `lc3_put_bits()`
359a19cd78SMatthias Ringwald  *
369a19cd78SMatthias Ringwald  * - The forward writing is protected against buffer overflow, it cannot
379a19cd78SMatthias Ringwald  *   write after the buffer, but can overwrite plain bits previously
389a19cd78SMatthias Ringwald  *   written in the buffer.
399a19cd78SMatthias Ringwald  *
409a19cd78SMatthias Ringwald  * - The backward writing is protected against overwrite of the arithmetic
419a19cd78SMatthias Ringwald  *   coder bitstream. In such way, the backward bitstream is always limited
429a19cd78SMatthias Ringwald  *   by the aritmetic coder bitstream, and can be overwritten by him.
439a19cd78SMatthias Ringwald  *
449a19cd78SMatthias Ringwald  *       .---------------------------------------------------.
459a19cd78SMatthias Ringwald  *       | > > > > > > > > > > :         : < < < < < < < < < |
469a19cd78SMatthias Ringwald  *       '---------------------------------------------------'
479a19cd78SMatthias Ringwald  *       |---------------------> - - - - - - - - - - - - - ->|
489a19cd78SMatthias Ringwald  *       |< - - - - - - - - - - -  - - - <-------------------|
499a19cd78SMatthias Ringwald  *          Arithmetic coding                  Plain bits
509a19cd78SMatthias Ringwald  *          `lc3_get_symbol()`               `lc3_get_bits()`
519a19cd78SMatthias Ringwald  *
529a19cd78SMatthias Ringwald  * - Reading is limited to read of the complementary end of the buffer.
539a19cd78SMatthias Ringwald  *
549a19cd78SMatthias Ringwald  * - The procedure `lc3_check_bits()` returns indication that read has been
559a19cd78SMatthias Ringwald  *   made crossing the other bit plane.
569a19cd78SMatthias Ringwald  *
579a19cd78SMatthias Ringwald  */
589a19cd78SMatthias Ringwald 
599a19cd78SMatthias Ringwald #ifndef __LC3_BITS_H
609a19cd78SMatthias Ringwald #define __LC3_BITS_H
619a19cd78SMatthias Ringwald 
629a19cd78SMatthias Ringwald #include "common.h"
639a19cd78SMatthias Ringwald 
649a19cd78SMatthias Ringwald 
659a19cd78SMatthias Ringwald /**
669a19cd78SMatthias Ringwald  * Bitstream mode
679a19cd78SMatthias Ringwald  */
689a19cd78SMatthias Ringwald 
699a19cd78SMatthias Ringwald enum lc3_bits_mode {
709a19cd78SMatthias Ringwald     LC3_BITS_MODE_READ,
719a19cd78SMatthias Ringwald     LC3_BITS_MODE_WRITE,
729a19cd78SMatthias Ringwald };
739a19cd78SMatthias Ringwald 
749a19cd78SMatthias Ringwald /**
759a19cd78SMatthias Ringwald  * Arithmetic coder symbol interval
769a19cd78SMatthias Ringwald  * The model split the interval in 17 symbols
779a19cd78SMatthias Ringwald  */
789a19cd78SMatthias Ringwald 
799a19cd78SMatthias Ringwald struct lc3_ac_symbol {
809a19cd78SMatthias Ringwald     uint16_t low   : 16;
819a19cd78SMatthias Ringwald     uint16_t range : 16;
829a19cd78SMatthias Ringwald };
839a19cd78SMatthias Ringwald 
849a19cd78SMatthias Ringwald struct lc3_ac_model {
859a19cd78SMatthias Ringwald     struct lc3_ac_symbol s[17];
869a19cd78SMatthias Ringwald };
879a19cd78SMatthias Ringwald 
889a19cd78SMatthias Ringwald /**
899a19cd78SMatthias Ringwald  * Bitstream context
909a19cd78SMatthias Ringwald  */
919a19cd78SMatthias Ringwald 
929a19cd78SMatthias Ringwald #define LC3_ACCU_BITS (int)(8 * sizeof(unsigned))
939a19cd78SMatthias Ringwald 
949a19cd78SMatthias Ringwald struct lc3_bits_accu {
959a19cd78SMatthias Ringwald     unsigned v;
969a19cd78SMatthias Ringwald     int n, nover;
979a19cd78SMatthias Ringwald };
989a19cd78SMatthias Ringwald 
999a19cd78SMatthias Ringwald #define LC3_AC_BITS (int)(24)
1009a19cd78SMatthias Ringwald 
1019a19cd78SMatthias Ringwald struct lc3_bits_ac {
1029a19cd78SMatthias Ringwald     unsigned low, range;
1039a19cd78SMatthias Ringwald     int cache, carry, carry_count;
1049a19cd78SMatthias Ringwald     bool error;
1059a19cd78SMatthias Ringwald };
1069a19cd78SMatthias Ringwald 
1079a19cd78SMatthias Ringwald struct lc3_bits_buffer {
1089a19cd78SMatthias Ringwald     const uint8_t *start, *end;
1099a19cd78SMatthias Ringwald     uint8_t *p_fw, *p_bw;
1109a19cd78SMatthias Ringwald };
1119a19cd78SMatthias Ringwald 
1129a19cd78SMatthias Ringwald typedef struct lc3_bits {
1139a19cd78SMatthias Ringwald     enum lc3_bits_mode mode;
1149a19cd78SMatthias Ringwald     struct lc3_bits_ac ac;
1159a19cd78SMatthias Ringwald     struct lc3_bits_accu accu;
1169a19cd78SMatthias Ringwald     struct lc3_bits_buffer buffer;
1179a19cd78SMatthias Ringwald } lc3_bits_t;
1189a19cd78SMatthias Ringwald 
1199a19cd78SMatthias Ringwald 
1209a19cd78SMatthias Ringwald /**
1219a19cd78SMatthias Ringwald  * Setup bitstream reading/writing
1229a19cd78SMatthias Ringwald  * bits            Bitstream context
1239a19cd78SMatthias Ringwald  * mode            Either READ or WRITE mode
1249a19cd78SMatthias Ringwald  * buffer, len     Output buffer and length (in bytes)
1259a19cd78SMatthias Ringwald  */
1269a19cd78SMatthias Ringwald void lc3_setup_bits(lc3_bits_t *bits,
1279a19cd78SMatthias Ringwald     enum lc3_bits_mode mode, void *buffer, int len);
1289a19cd78SMatthias Ringwald 
1299a19cd78SMatthias Ringwald /**
1309a19cd78SMatthias Ringwald  * Return number of bits left in the bitstream
1319a19cd78SMatthias Ringwald  * bits            Bitstream context
1329a19cd78SMatthias Ringwald  * return          Number of bits left
1339a19cd78SMatthias Ringwald  */
1349a19cd78SMatthias Ringwald int lc3_get_bits_left(const lc3_bits_t *bits);
1359a19cd78SMatthias Ringwald 
1369a19cd78SMatthias Ringwald /**
1379a19cd78SMatthias Ringwald  * Check if error occured on bitstream reading/writing
1389a19cd78SMatthias Ringwald  * bits            Bitstream context
1399a19cd78SMatthias Ringwald  * return          0: Ok  -1: Bitstream overflow or AC reading error
1409a19cd78SMatthias Ringwald  */
1419a19cd78SMatthias Ringwald int lc3_check_bits(const lc3_bits_t *bits);
1429a19cd78SMatthias Ringwald 
1439a19cd78SMatthias Ringwald /**
1449a19cd78SMatthias Ringwald  * Put a bit
1459a19cd78SMatthias Ringwald  * bits            Bitstream context
1469a19cd78SMatthias Ringwald  * v               Bit value, 0 or 1
1479a19cd78SMatthias Ringwald  */
1489a19cd78SMatthias Ringwald static inline void lc3_put_bit(lc3_bits_t *bits, int v);
1499a19cd78SMatthias Ringwald 
1509a19cd78SMatthias Ringwald /**
1519a19cd78SMatthias Ringwald  * Put from 1 to 32 bits
1529a19cd78SMatthias Ringwald  * bits            Bitstream context
1539a19cd78SMatthias Ringwald  * v, n            Value, in range 0 to 2^n - 1, and bits count (1 to 32)
1549a19cd78SMatthias Ringwald  */
1559a19cd78SMatthias Ringwald static inline void lc3_put_bits(lc3_bits_t *bits, unsigned v, int n);
1569a19cd78SMatthias Ringwald 
1579a19cd78SMatthias Ringwald /**
1589a19cd78SMatthias Ringwald  * Put arithmetic coder symbol
1599a19cd78SMatthias Ringwald  * bits            Bitstream context
1609a19cd78SMatthias Ringwald  * model, s        Model distribution and symbol value
1619a19cd78SMatthias Ringwald  */
1629a19cd78SMatthias Ringwald static inline void lc3_put_symbol(lc3_bits_t *bits,
1639a19cd78SMatthias Ringwald     const struct lc3_ac_model *model, unsigned s);
1649a19cd78SMatthias Ringwald 
1659a19cd78SMatthias Ringwald /**
1669a19cd78SMatthias Ringwald  * Flush and terminate bitstream writing
1679a19cd78SMatthias Ringwald  * bits            Bitstream context
1689a19cd78SMatthias Ringwald  */
1699a19cd78SMatthias Ringwald void lc3_flush_bits(lc3_bits_t *bits);
1709a19cd78SMatthias Ringwald 
1719a19cd78SMatthias Ringwald /**
1729a19cd78SMatthias Ringwald  * Get a bit
1739a19cd78SMatthias Ringwald  * bits            Bitstream context
1749a19cd78SMatthias Ringwald  */
1759a19cd78SMatthias Ringwald static inline int lc3_get_bit(lc3_bits_t *bits);
1769a19cd78SMatthias Ringwald 
1779a19cd78SMatthias Ringwald /**
1789a19cd78SMatthias Ringwald  * Get from 1 to 32 bits
1799a19cd78SMatthias Ringwald  * bits            Bitstream context
1809a19cd78SMatthias Ringwald  * n               Number of bits to read (1 to 32)
1819a19cd78SMatthias Ringwald  * return          The value read
1829a19cd78SMatthias Ringwald  */
1839a19cd78SMatthias Ringwald static inline unsigned lc3_get_bits(lc3_bits_t *bits,  int n);
1849a19cd78SMatthias Ringwald 
1859a19cd78SMatthias Ringwald /**
1869a19cd78SMatthias Ringwald  * Get arithmetic coder symbol
1879a19cd78SMatthias Ringwald  * bits            Bitstream context
1889a19cd78SMatthias Ringwald  * model           Model distribution
1899a19cd78SMatthias Ringwald  * return          The value read
1909a19cd78SMatthias Ringwald  */
1919a19cd78SMatthias Ringwald static inline unsigned lc3_get_symbol(lc3_bits_t *bits,
1929a19cd78SMatthias Ringwald     const struct lc3_ac_model *model);
1939a19cd78SMatthias Ringwald 
1949a19cd78SMatthias Ringwald 
1959a19cd78SMatthias Ringwald 
1969a19cd78SMatthias Ringwald /* ----------------------------------------------------------------------------
1979a19cd78SMatthias Ringwald  *  Inline implementations
1989a19cd78SMatthias Ringwald  * -------------------------------------------------------------------------- */
1999a19cd78SMatthias Ringwald 
2009a19cd78SMatthias Ringwald void lc3_put_bits_generic(lc3_bits_t *bits, unsigned v, int n);
2019a19cd78SMatthias Ringwald unsigned lc3_get_bits_generic(struct lc3_bits *bits, int n);
2029a19cd78SMatthias Ringwald 
2039a19cd78SMatthias Ringwald void lc3_ac_read_renorm(lc3_bits_t *bits);
2049a19cd78SMatthias Ringwald void lc3_ac_write_renorm(lc3_bits_t *bits);
2059a19cd78SMatthias Ringwald 
2069a19cd78SMatthias Ringwald 
2079a19cd78SMatthias Ringwald /**
2089a19cd78SMatthias Ringwald  * Put a bit
2099a19cd78SMatthias Ringwald  */
lc3_put_bit(lc3_bits_t * bits,int v)210*4930cef6SMatthias Ringwald LC3_HOT static inline void lc3_put_bit(lc3_bits_t *bits, int v)
2119a19cd78SMatthias Ringwald {
2129a19cd78SMatthias Ringwald     lc3_put_bits(bits, v, 1);
2139a19cd78SMatthias Ringwald }
2149a19cd78SMatthias Ringwald 
2159a19cd78SMatthias Ringwald /**
2169a19cd78SMatthias Ringwald  * Put from 1 to 32 bits
2179a19cd78SMatthias Ringwald  */
lc3_put_bits(struct lc3_bits * bits,unsigned v,int n)218*4930cef6SMatthias Ringwald LC3_HOT static inline void lc3_put_bits(
219*4930cef6SMatthias Ringwald     struct lc3_bits *bits, unsigned v, int n)
2209a19cd78SMatthias Ringwald {
2219a19cd78SMatthias Ringwald     struct lc3_bits_accu *accu = &bits->accu;
2229a19cd78SMatthias Ringwald 
2239a19cd78SMatthias Ringwald     if (accu->n + n <= LC3_ACCU_BITS) {
2249a19cd78SMatthias Ringwald         accu->v |= v << accu->n;
2259a19cd78SMatthias Ringwald         accu->n += n;
2269a19cd78SMatthias Ringwald     } else {
2279a19cd78SMatthias Ringwald         lc3_put_bits_generic(bits, v, n);
2289a19cd78SMatthias Ringwald     }
2299a19cd78SMatthias Ringwald }
2309a19cd78SMatthias Ringwald 
2319a19cd78SMatthias Ringwald /**
2329a19cd78SMatthias Ringwald  * Get a bit
2339a19cd78SMatthias Ringwald  */
lc3_get_bit(lc3_bits_t * bits)234*4930cef6SMatthias Ringwald LC3_HOT static inline int lc3_get_bit(lc3_bits_t *bits)
2359a19cd78SMatthias Ringwald {
2369a19cd78SMatthias Ringwald     return lc3_get_bits(bits, 1);
2379a19cd78SMatthias Ringwald }
2389a19cd78SMatthias Ringwald 
2399a19cd78SMatthias Ringwald /**
2409a19cd78SMatthias Ringwald  * Get from 1 to 32 bits
2419a19cd78SMatthias Ringwald  */
lc3_get_bits(struct lc3_bits * bits,int n)242*4930cef6SMatthias Ringwald LC3_HOT static inline unsigned lc3_get_bits(struct lc3_bits *bits, int n)
2439a19cd78SMatthias Ringwald {
2449a19cd78SMatthias Ringwald     struct lc3_bits_accu *accu = &bits->accu;
2459a19cd78SMatthias Ringwald 
2469a19cd78SMatthias Ringwald     if (accu->n + n <= LC3_ACCU_BITS) {
2479a19cd78SMatthias Ringwald         int v = (accu->v >> accu->n) & ((1u << n) - 1);
2489a19cd78SMatthias Ringwald         return (accu->n += n), v;
2499a19cd78SMatthias Ringwald     }
2509a19cd78SMatthias Ringwald     else {
2519a19cd78SMatthias Ringwald         return lc3_get_bits_generic(bits, n);
2529a19cd78SMatthias Ringwald     }
2539a19cd78SMatthias Ringwald }
2549a19cd78SMatthias Ringwald 
2559a19cd78SMatthias Ringwald /**
2569a19cd78SMatthias Ringwald  * Put arithmetic coder symbol
2579a19cd78SMatthias Ringwald  */
lc3_put_symbol(struct lc3_bits * bits,const struct lc3_ac_model * model,unsigned s)258*4930cef6SMatthias Ringwald LC3_HOT static inline void lc3_put_symbol(
2599a19cd78SMatthias Ringwald     struct lc3_bits *bits, const struct lc3_ac_model *model, unsigned s)
2609a19cd78SMatthias Ringwald {
2619a19cd78SMatthias Ringwald     const struct lc3_ac_symbol *symbols = model->s;
2629a19cd78SMatthias Ringwald     struct lc3_bits_ac *ac = &bits->ac;
2639a19cd78SMatthias Ringwald     unsigned range = ac->range >> 10;
2649a19cd78SMatthias Ringwald 
2659a19cd78SMatthias Ringwald     ac->low += range * symbols[s].low;
2669a19cd78SMatthias Ringwald     ac->range = range * symbols[s].range;
2679a19cd78SMatthias Ringwald 
2689a19cd78SMatthias Ringwald     ac->carry |= ac->low >> 24;
2699a19cd78SMatthias Ringwald     ac->low &= 0xffffff;
2709a19cd78SMatthias Ringwald 
2719a19cd78SMatthias Ringwald     if (ac->range < 0x10000)
2729a19cd78SMatthias Ringwald         lc3_ac_write_renorm(bits);
2739a19cd78SMatthias Ringwald }
2749a19cd78SMatthias Ringwald 
2759a19cd78SMatthias Ringwald /**
2769a19cd78SMatthias Ringwald  * Get arithmetic coder symbol
2779a19cd78SMatthias Ringwald  */
lc3_get_symbol(lc3_bits_t * bits,const struct lc3_ac_model * model)278*4930cef6SMatthias Ringwald LC3_HOT static inline unsigned lc3_get_symbol(
2799a19cd78SMatthias Ringwald     lc3_bits_t *bits, const struct lc3_ac_model *model)
2809a19cd78SMatthias Ringwald {
2819a19cd78SMatthias Ringwald     const struct lc3_ac_symbol *symbols = model->s;
2829a19cd78SMatthias Ringwald     struct lc3_bits_ac *ac = &bits->ac;
2839a19cd78SMatthias Ringwald 
2849a19cd78SMatthias Ringwald     unsigned range = (ac->range >> 10) & 0xffff;
2859a19cd78SMatthias Ringwald 
2869a19cd78SMatthias Ringwald     ac->error |= (ac->low >= (range << 10));
2879a19cd78SMatthias Ringwald     if (ac->error)
2889a19cd78SMatthias Ringwald         ac->low = 0;
2899a19cd78SMatthias Ringwald 
2909a19cd78SMatthias Ringwald     int s = 16;
2919a19cd78SMatthias Ringwald 
2929a19cd78SMatthias Ringwald     if (ac->low < range * symbols[s].low) {
2939a19cd78SMatthias Ringwald         s >>= 1;
2949a19cd78SMatthias Ringwald         s -= ac->low < range * symbols[s].low ? 4 : -4;
2959a19cd78SMatthias Ringwald         s -= ac->low < range * symbols[s].low ? 2 : -2;
2969a19cd78SMatthias Ringwald         s -= ac->low < range * symbols[s].low ? 1 : -1;
2979a19cd78SMatthias Ringwald         s -= ac->low < range * symbols[s].low;
2989a19cd78SMatthias Ringwald     }
2999a19cd78SMatthias Ringwald 
3009a19cd78SMatthias Ringwald     ac->low -= range * symbols[s].low;
3019a19cd78SMatthias Ringwald     ac->range = range * symbols[s].range;
3029a19cd78SMatthias Ringwald 
3039a19cd78SMatthias Ringwald     if (ac->range < 0x10000)
3049a19cd78SMatthias Ringwald         lc3_ac_read_renorm(bits);
3059a19cd78SMatthias Ringwald 
3069a19cd78SMatthias Ringwald     return s;
3079a19cd78SMatthias Ringwald }
3089a19cd78SMatthias Ringwald 
3099a19cd78SMatthias Ringwald #endif /* __LC3_BITS_H */
310