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