Line data Source code
1 : /******************************************************************************************************
2 :
3 : (C) 2022-2025 IVAS codec Public Collaboration with portions copyright Dolby International AB, Ericsson AB,
4 : Fraunhofer-Gesellschaft zur Foerderung der angewandten Forschung e.V., Huawei Technologies Co. LTD.,
5 : Koninklijke Philips N.V., Nippon Telegraph and Telephone Corporation, Nokia Technologies Oy, Orange,
6 : Panasonic Holdings Corporation, Qualcomm Technologies, Inc., VoiceAge Corporation, and other
7 : contributors to this repository. All Rights Reserved.
8 :
9 : This software is protected by copyright law and by international treaties.
10 : The IVAS codec Public Collaboration consisting of Dolby International AB, Ericsson AB,
11 : Fraunhofer-Gesellschaft zur Foerderung der angewandten Forschung e.V., Huawei Technologies Co. LTD.,
12 : Koninklijke Philips N.V., Nippon Telegraph and Telephone Corporation, Nokia Technologies Oy, Orange,
13 : Panasonic Holdings Corporation, Qualcomm Technologies, Inc., VoiceAge Corporation, and other
14 : contributors to this repository retain full ownership rights in their respective contributions in
15 : the software. This notice grants no license of any kind, including but not limited to patent
16 : license, nor is any license granted by implication, estoppel or otherwise.
17 :
18 : Contributors are required to enter into the IVAS codec Public Collaboration agreement before making
19 : contributions.
20 :
21 : This software is provided "AS IS", without any express or implied warranties. The software is in the
22 : development stage. It is intended exclusively for experts who have experience with such software and
23 : solely for the purpose of inspection. All implied warranties of non-infringement, merchantability
24 : and fitness for a particular purpose are hereby disclaimed and excluded.
25 :
26 : Any dispute, controversy or claim arising under or in relation to providing this software shall be
27 : submitted to and settled by the final, binding jurisdiction of the courts of Munich, Germany in
28 : accordance with the laws of the Federal Republic of Germany excluding its conflict of law rules and
29 : the United Nations Convention on Contracts on the International Sales of Goods.
30 :
31 : *******************************************************************************************************/
32 :
33 : #include <stdint.h>
34 : #include "options.h"
35 : #ifdef DEBUGGING
36 : #include "debug.h"
37 : #endif
38 : #include "wmc_auto.h"
39 : #include "prot.h"
40 : #include "ivas_prot.h"
41 : #include "stat_dec.h"
42 :
43 : /*-----------------------------------------------------------------------------------------*
44 : * Function ivas_ari_start_decoding_14bits_ext_1_lfe()
45 : *
46 : * Start arithemetic coding
47 : *
48 : * Similar to ari_start_decoding_14bits(), but do not read past frame boundary.
49 : *-----------------------------------------------------------------------------------------*/
50 :
51 36815531 : void ivas_ari_start_decoding_14bits_ext_1_lfe(
52 : Decoder_State *st,
53 : Tastat *s,
54 : int16_t *extra_bits_read )
55 : {
56 : int32_t val;
57 :
58 36815531 : if ( st->bits_frame >= ( st->next_bit_pos + cbitsnew ) )
59 : {
60 36815531 : val = get_next_indice( st, cbitsnew );
61 : }
62 : else
63 : {
64 : int16_t rem_bits;
65 0 : rem_bits = st->bits_frame - st->next_bit_pos;
66 0 : val = get_next_indice( st, rem_bits );
67 0 : val = ( val << ( cbitsnew - rem_bits ) );
68 0 : *extra_bits_read = *extra_bits_read + ( cbitsnew - rem_bits );
69 : }
70 :
71 36815531 : s->low = 0;
72 36815531 : s->high = ari_q4new;
73 36815531 : s->value = val;
74 :
75 36815531 : return;
76 : }
77 :
78 :
79 : /* Helper function to determine the symbol. */
80 633822108 : static uint16_t ivas_ari_get_symbol(
81 : const uint16_t *p,
82 : uint32_t range,
83 : uint32_t cum )
84 : {
85 : uint16_t symbol;
86 :
87 633822108 : symbol = 0;
88 :
89 4948487686 : while ( ( p[symbol + 1] * range ) > cum )
90 : {
91 4314665578 : symbol = symbol + 1;
92 : }
93 :
94 633822108 : return symbol;
95 : }
96 :
97 :
98 : /*---------------------------------------------------------------
99 : * Function ivas_ari_decode_14bits_bit_ext_1_lfe()
100 : *
101 : *
102 : *-------------------------------------------------------------*/
103 :
104 633822108 : uint16_t ivas_ari_decode_14bits_bit_ext_1_lfe(
105 : Decoder_State *st,
106 : Tastat *s,
107 : const uint16_t *cum_freq,
108 : int16_t *extra_bits_read )
109 : {
110 : uint16_t symbol;
111 : uint32_t low, high, range, value;
112 : uint32_t cum;
113 :
114 : /* read s->low,high,value sequentially */
115 633822108 : low = s->low;
116 633822108 : high = s->high;
117 633822108 : value = s->value;
118 :
119 633822108 : range = high - low + 1;
120 633822108 : cum = ( ( ( (uint32_t) ( value - low + 1 ) ) << stat_bitsnew ) - ( (uint32_t) 1 ) );
121 :
122 633822108 : symbol = ivas_ari_get_symbol( cum_freq, range, cum );
123 :
124 633822108 : high = low + mul_sbc_14bits( range, cum_freq[symbol] ) - 1;
125 633822108 : low += mul_sbc_14bits( range, cum_freq[symbol + 1] );
126 :
127 : for ( ;; )
128 : {
129 1899955303 : if ( high >= ari_q2new )
130 : {
131 1443320101 : if ( low >= ari_q2new )
132 : {
133 447420385 : value -= ari_q2new;
134 447420385 : low -= ari_q2new;
135 447420385 : high -= ari_q2new;
136 : }
137 : else
138 : {
139 995899716 : if ( low >= ari_q1new && high < ari_q3new )
140 : {
141 362077608 : value -= ari_q1new;
142 362077608 : low -= ari_q1new;
143 362077608 : high -= ari_q1new;
144 : }
145 : else
146 : {
147 : break;
148 : }
149 : }
150 : }
151 1266133195 : low += low;
152 1266133195 : high += high + 1;
153 :
154 1266133195 : if ( st->next_bit_pos >= st->bits_frame )
155 : {
156 0 : value = ( value << 1 );
157 0 : *extra_bits_read = *extra_bits_read + 1;
158 : }
159 : else
160 : {
161 1266133195 : value = ( value << 1 ) | get_next_indice_1( st );
162 : }
163 : }
164 :
165 633822108 : s->low = low;
166 633822108 : s->high = high;
167 633822108 : s->value = value;
168 :
169 633822108 : return symbol;
170 : }
171 :
172 :
173 : /*---------------------------------------------------------------
174 : * Function ivas_ari_done_decoding_14bits_ext_1_lfe()
175 : *
176 : *
177 : *-------------------------------------------------------------*/
178 :
179 36815531 : void ivas_ari_done_decoding_14bits_ext_1_lfe(
180 : Decoder_State *st,
181 : const int16_t extra_bits_read )
182 : {
183 36815531 : get_next_indice_tmp( st, -( cbitsnew - 2 - extra_bits_read ) );
184 :
185 36815531 : return;
186 : }
187 :
188 :
189 : /*---------------------------------------------------------------
190 : * Function ivas_ari_done_encoding_14bits()
191 : *
192 : * Arith encoding of last table entry
193 : *-------------------------------------------------------------*/
194 :
195 15883152 : void ivas_ari_done_encoding_14bits(
196 : BSTR_ENC_HANDLE hBstr,
197 : Tastat *s )
198 : {
199 : int32_t low;
200 : int32_t bits_to_follow;
201 :
202 : /* not needed, s points to s->low */
203 15883152 : low = s->low;
204 15883152 : bits_to_follow = s->bits_to_follow + 1;
205 :
206 15883152 : if ( low < ari_q1new )
207 : {
208 9948375 : push_next_indice( hBstr, 0, 1 ); /*send a zero*/
209 :
210 31663988 : for ( ; bits_to_follow > 0; bits_to_follow-- )
211 : {
212 21715613 : push_next_indice( hBstr, 1, 1 ); /*send a one*/
213 : }
214 : }
215 : else
216 : {
217 5934777 : push_next_indice( hBstr, 1, 1 ); /*send a one*/
218 17937192 : for ( ; bits_to_follow > 0; bits_to_follow-- )
219 : {
220 12002415 : push_next_indice( hBstr, 0, 1 ); /*send a zero*/
221 : }
222 : }
223 :
224 : /*It is done so no need to save values-> no counting*/
225 : /*s->low = low;
226 : s->bits_to_follow = bits_to_follow;*/
227 :
228 15883152 : return;
229 : }
230 :
231 :
232 : /*---------------------------------------------------------------
233 : * Function ivas_ari_encode_14bits_ext()
234 : *
235 : * Arith encode function for extended proba tables
236 : *-------------------------------------------------------------*/
237 :
238 279416646 : void ivas_ari_encode_14bits_ext(
239 : BSTR_ENC_HANDLE hBstr,
240 : Tastat *s,
241 : int32_t symbol,
242 : const uint16_t *cum_freq )
243 : {
244 : int32_t low, high, range;
245 : int32_t bits_to_follow;
246 :
247 : /*for all operation using bit_ptr=&ptr[bp] */
248 : /* for reading s->high,low,bits_to_follow sequentially */
249 279416646 : high = s->high;
250 279416646 : low = s->low;
251 279416646 : range = high - low + 1;
252 :
253 279416646 : high = low + mul_sbc_14bits( range, cum_freq[symbol] ) - 1;
254 279416646 : low += mul_sbc_14bits( range, cum_freq[symbol + 1] );
255 :
256 279416646 : bits_to_follow = s->bits_to_follow;
257 :
258 : for ( ;; )
259 : {
260 886690323 : if ( high < ari_q2new )
261 : {
262 218775124 : push_next_indice( hBstr, 0, 1 ); /*send a zero*/
263 300996675 : for ( ; bits_to_follow > 0; bits_to_follow-- )
264 : {
265 82221551 : push_next_indice( hBstr, 1, 1 ); /*send a one*/
266 : }
267 : }
268 : else
269 : {
270 667915199 : if ( low >= ari_q2new )
271 : {
272 214893036 : push_next_indice( hBstr, 1, 1 ); /*send a one*/
273 288349406 : for ( ; bits_to_follow > 0; bits_to_follow-- )
274 : {
275 73456370 : push_next_indice( hBstr, 0, 1 ); /*send a zero*/
276 : }
277 214893036 : low -= ari_q2new;
278 214893036 : high -= ari_q2new; /* Subtract offset to top. */
279 : }
280 : else
281 : {
282 : /* Output an opposite bit */
283 453022163 : if ( low >= ari_q1new && high < ari_q3new ) /* Output an opposite bit */
284 : {
285 : /* later if in middle half. */
286 173605517 : bits_to_follow += 1;
287 173605517 : low -= ari_q1new; /* Subtract offset to middle*/
288 173605517 : high -= ari_q1new;
289 : }
290 : else
291 : {
292 : break; /* Otherwise exit loop. */
293 : }
294 : }
295 : }
296 607273677 : low += low;
297 607273677 : high += high + 1; /* Scale up code range. */
298 : }
299 :
300 279416646 : s->low = low;
301 279416646 : s->high = high;
302 279416646 : s->bits_to_follow = bits_to_follow;
303 :
304 279416646 : return;
305 : }
|
