FFmpeg coverage


Directory: ../../../ffmpeg/
File: src/libavcodec/aacenc.c
Date: 2026-07-22 00:55:30
Exec Total Coverage
Lines: 764 859 88.9%
Functions: 28 30 93.3%
Branches: 558 724 77.1%

Line Branch Exec Source
1 /*
2 * AAC encoder
3 * Copyright (C) 2008 Konstantin Shishkov
4 *
5 * This file is part of FFmpeg.
6 *
7 * FFmpeg is free software; you can redistribute it and/or
8 * modify it under the terms of the GNU Lesser General Public
9 * License as published by the Free Software Foundation; either
10 * version 2.1 of the License, or (at your option) any later version.
11 *
12 * FFmpeg is distributed in the hope that it will be useful,
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15 * Lesser General Public License for more details.
16 *
17 * You should have received a copy of the GNU Lesser General Public
18 * License along with FFmpeg; if not, write to the Free Software
19 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
20 */
21
22 /**
23 * @file
24 * AAC encoder
25 */
26
27 /***********************************
28 * TODOs:
29 * add sane pulse detection
30 ***********************************/
31 #include <float.h>
32
33 #include "libavutil/channel_layout.h"
34 #include "libavutil/libm.h"
35 #include "libavutil/float_dsp.h"
36 #include "libavutil/mem.h"
37 #include "libavutil/opt.h"
38 #include "avcodec.h"
39 #include "codec_internal.h"
40 #include "encode.h"
41 #include "put_bits.h"
42 #include "mpeg4audio.h"
43 #include "sinewin.h"
44 #include "profiles.h"
45 #include "version.h"
46
47 #include "aac.h"
48 #include "aactab.h"
49 #include "aacenc.h"
50 #include "aacenctab.h"
51 #include "aacenc_utils.h"
52
53 #include "psymodel.h"
54
55 /**
56 * List of PCE (Program Configuration Element) for the channel layouts listed
57 * in channel_layout.h
58 *
59 * For those wishing in the future to add other layouts:
60 *
61 * - num_ele: number of elements in each group of front, side, back, lfe channels
62 * (an element is of type SCE (single channel), CPE (channel pair) for
63 * the first 3 groups; and is LFE for LFE group).
64 *
65 * - pairing: 0 for an SCE element or 1 for a CPE; does not apply to LFE group
66 *
67 * - index: there are three independent indices for SCE, CPE and LFE;
68 * they are incremented irrespective of the group to which the element belongs;
69 * they are not reset when going from one group to another
70 *
71 * Example: for 7.0 channel layout,
72 * .pairing = { { 1, 0 }, { 1 }, { 1 }, }, (3 CPE and 1 SCE in front group)
73 * .index = { { 0, 0 }, { 1 }, { 2 }, },
74 * (index is 0 for the single SCE but goes from 0 to 2 for the CPEs)
75 *
76 * The index order impacts the channel ordering. But is otherwise arbitrary
77 * (the sequence could have been 2, 0, 1 instead of 0, 1, 2).
78 *
79 * Spec allows for discontinuous indices, e.g. if one has a total of two SCE,
80 * SCE.0 SCE.15 is OK per spec; BUT it won't be decoded by our AAC decoder
81 * which at this time requires that indices fully cover some range starting
82 * from 0 (SCE.1 SCE.0 is OK but not SCE.0 SCE.15).
83 *
84 * - config_map: total number of elements and their types. Beware, the way the
85 * types are ordered impacts the final channel ordering.
86 *
87 * - reorder_map: reorders the channels.
88 *
89 */
90 static const AACPCEInfo aac_pce_configs[] = {
91 {
92 .layout = AV_CHANNEL_LAYOUT_MONO,
93 .num_ele = { 1, 0, 0, 0 },
94 .pairing = { { 0 }, },
95 .index = { { 0 }, },
96 .config_map = { 1, TYPE_SCE, },
97 .reorder_map = { 0 },
98 },
99 {
100 .layout = AV_CHANNEL_LAYOUT_STEREO,
101 .num_ele = { 1, 0, 0, 0 },
102 .pairing = { { 1 }, },
103 .index = { { 0 }, },
104 .config_map = { 1, TYPE_CPE, },
105 .reorder_map = { 0, 1 },
106 },
107 {
108 .layout = AV_CHANNEL_LAYOUT_2POINT1,
109 .num_ele = { 1, 0, 0, 1 },
110 .pairing = { { 1 }, },
111 .index = { { 0 },{ 0 },{ 0 },{ 0 } },
112 .config_map = { 2, TYPE_CPE, TYPE_LFE },
113 .reorder_map = { 0, 1, 2 },
114 },
115 {
116 .layout = AV_CHANNEL_LAYOUT_2_1,
117 .num_ele = { 1, 0, 1, 0 },
118 .pairing = { { 1 },{ 0 },{ 0 } },
119 .index = { { 0 },{ 0 },{ 0 }, },
120 .config_map = { 2, TYPE_CPE, TYPE_SCE },
121 .reorder_map = { 0, 1, 2 },
122 },
123 {
124 .layout = AV_CHANNEL_LAYOUT_SURROUND,
125 .num_ele = { 2, 0, 0, 0 },
126 .pairing = { { 0, 1 }, },
127 .index = { { 0, 0 }, },
128 .config_map = { 2, TYPE_SCE, TYPE_CPE },
129 .reorder_map = { 2, 0, 1 },
130 },
131 {
132 .layout = AV_CHANNEL_LAYOUT_3POINT1,
133 .num_ele = { 2, 0, 0, 1 },
134 .pairing = { { 0, 1 }, },
135 .index = { { 0, 0 }, { 0 }, { 0 }, { 0 }, },
136 .config_map = { 3, TYPE_SCE, TYPE_CPE, TYPE_LFE },
137 .reorder_map = { 2, 0, 1, 3 },
138 },
139 {
140 .layout = AV_CHANNEL_LAYOUT_4POINT0,
141 .num_ele = { 2, 0, 1, 0 },
142 .pairing = { { 0, 1 }, { 0 }, { 0 }, },
143 .index = { { 0, 0 }, { 0 }, { 1 } },
144 .config_map = { 3, TYPE_SCE, TYPE_CPE, TYPE_SCE },
145 .reorder_map = { 2, 0, 1, 3 },
146 },
147 {
148 .layout = AV_CHANNEL_LAYOUT_4POINT1,
149 .num_ele = { 2, 0, 1, 1 },
150 .pairing = { { 0, 1 }, { 0 }, { 0 }, },
151 .index = { { 0, 0 }, { 0 }, { 1 }, { 0 } },
152 .config_map = { 4, TYPE_SCE, TYPE_CPE, TYPE_SCE, TYPE_LFE },
153 .reorder_map = { 2, 0, 1, 4, 3 },
154 },
155 {
156 .layout = AV_CHANNEL_LAYOUT_2_2,
157 .num_ele = { 1, 0, 1, 0 },
158 .pairing = { { 1 }, { 0 }, { 1 }, },
159 .index = { { 0 }, { 0 }, { 1 } },
160 .config_map = { 2, TYPE_CPE, TYPE_CPE },
161 .reorder_map = { 0, 1, 2, 3 },
162 },
163 {
164 .layout = AV_CHANNEL_LAYOUT_QUAD,
165 .num_ele = { 1, 0, 1, 0 },
166 .pairing = { { 1 }, { 0 }, { 1 }, },
167 .index = { { 0 }, { 0 }, { 1 } },
168 .config_map = { 2, TYPE_CPE, TYPE_CPE },
169 .reorder_map = { 0, 1, 2, 3 },
170 },
171 {
172 .layout = AV_CHANNEL_LAYOUT_5POINT0,
173 .num_ele = { 2, 0, 1, 0 },
174 .pairing = { { 0, 1 }, { 0 }, { 1 } },
175 .index = { { 0, 0 }, { 0 }, { 1 } },
176 .config_map = { 3, TYPE_SCE, TYPE_CPE, TYPE_CPE },
177 .reorder_map = { 2, 0, 1, 3, 4 },
178 },
179 {
180 .layout = AV_CHANNEL_LAYOUT_5POINT1,
181 .num_ele = { 2, 0, 1, 1 },
182 .pairing = { { 0, 1 }, { 0 }, { 1 }, },
183 .index = { { 0, 0 }, { 0 }, { 1 }, { 0 } },
184 .config_map = { 4, TYPE_SCE, TYPE_CPE, TYPE_CPE, TYPE_LFE },
185 .reorder_map = { 2, 0, 1, 4, 5, 3 },
186 },
187 {
188 .layout = AV_CHANNEL_LAYOUT_5POINT0_BACK,
189 .num_ele = { 2, 0, 1, 0 },
190 .pairing = { { 0, 1 }, { 0 }, { 1 } },
191 .index = { { 0, 0 }, { 0 }, { 1 } },
192 .config_map = { 3, TYPE_SCE, TYPE_CPE, TYPE_CPE },
193 .reorder_map = { 2, 0, 1, 3, 4 },
194 },
195 {
196 .layout = AV_CHANNEL_LAYOUT_5POINT1_BACK,
197 .num_ele = { 2, 0, 1, 1 },
198 .pairing = { { 0, 1 }, { 0 }, { 1 }, },
199 .index = { { 0, 0 }, { 0 }, { 1 }, { 0 } },
200 .config_map = { 4, TYPE_SCE, TYPE_CPE, TYPE_CPE, TYPE_LFE },
201 .reorder_map = { 2, 0, 1, 4, 5, 3 },
202 },
203 {
204 .layout = AV_CHANNEL_LAYOUT_6POINT0,
205 .num_ele = { 2, 0, 2, 0 },
206 .pairing = { { 0, 1 }, { 0 }, { 1, 0 } },
207 .index = { { 0, 0 }, { 0 }, { 1, 1 } },
208 .config_map = { 4, TYPE_SCE, TYPE_CPE, TYPE_CPE, TYPE_SCE },
209 .reorder_map = { 2, 0, 1, 4, 5, 3 },
210 },
211 {
212 .layout = AV_CHANNEL_LAYOUT_6POINT0_FRONT,
213 .num_ele = { 2, 0, 1, 0 },
214 .pairing = { { 1, 1 }, { 0 }, { 1 } },
215 .index = { { 0, 1 }, { 0 }, { 2 }, },
216 .config_map = { 3, TYPE_CPE, TYPE_CPE, TYPE_CPE, },
217 .reorder_map = { 2, 3, 0, 1, 4, 5 },
218 },
219 {
220 .layout = AV_CHANNEL_LAYOUT_HEXAGONAL,
221 .num_ele = { 2, 0, 2, 0 },
222 .pairing = { { 0, 1 }, { 0 }, { 1, 0 } },
223 .index = { { 0, 0 }, { 0 }, { 1, 1 } },
224 .config_map = { 4, TYPE_SCE, TYPE_CPE, TYPE_CPE, TYPE_SCE },
225 .reorder_map = { 2, 0, 1, 3, 4, 5 },
226 },
227 {
228 .layout = AV_CHANNEL_LAYOUT_6POINT1,
229 .num_ele = { 2, 0, 2, 1 },
230 .pairing = { { 0, 1 }, { 0 }, { 1, 0 }, },
231 .index = { { 0, 0 }, { 0 }, { 1, 1 }, { 0 } },
232 .config_map = { 5, TYPE_SCE, TYPE_CPE, TYPE_CPE, TYPE_SCE, TYPE_LFE },
233 .reorder_map = { 2, 0, 1, 5, 6, 4, 3 },
234 },
235 {
236 .layout = AV_CHANNEL_LAYOUT_6POINT1_BACK,
237 .num_ele = { 2, 0, 2, 1 },
238 .pairing = { { 0, 1 },{ 0 },{ 1, 0 }, },
239 .index = { { 0, 0 },{ 0 },{ 1, 1 },{ 0 } },
240 .config_map = { 5, TYPE_SCE, TYPE_CPE, TYPE_CPE, TYPE_SCE, TYPE_LFE },
241 .reorder_map = { 2, 0, 1, 4, 5, 6, 3 },
242 },
243 {
244 .layout = AV_CHANNEL_LAYOUT_6POINT1_FRONT,
245 .num_ele = { 2, 0, 1, 1 },
246 .pairing = { { 1, 1 }, { 0 }, { 1 }, },
247 .index = { { 0, 1 }, { 0 }, { 2 }, { 0 }, },
248 .config_map = { 4, TYPE_CPE, TYPE_CPE, TYPE_CPE, TYPE_LFE, },
249 .reorder_map = { 3, 4, 0, 1, 5, 6, 2 },
250 },
251 {
252 .layout = AV_CHANNEL_LAYOUT_7POINT0,
253 .num_ele = { 2, 0, 2, 0 },
254 .pairing = { { 0, 1 }, { 0 }, { 1, 1 }, },
255 .index = { { 0, 0 }, { 0 }, { 2, 1 }, },
256 .config_map = { 4, TYPE_SCE, TYPE_CPE, TYPE_CPE, TYPE_CPE },
257 .reorder_map = { 2, 0, 1, 3, 4, 5, 6 },
258 },
259 {
260 .layout = AV_CHANNEL_LAYOUT_7POINT0_FRONT,
261 .num_ele = { 3, 0, 1, 0 },
262 .pairing = { { 0, 1, 1 }, { 0 }, { 1 }, },
263 .index = { { 0, 0, 1 }, { 0 }, { 2 }, },
264 .config_map = { 4, TYPE_SCE, TYPE_CPE, TYPE_CPE, TYPE_CPE },
265 .reorder_map = { 2, 3, 4, 0, 1, 5, 6 },
266 },
267 {
268 .layout = AV_CHANNEL_LAYOUT_7POINT1,
269 .num_ele = { 2, 0, 2, 1 },
270 .pairing = { { 0, 1 }, { 0 }, { 1, 1 }, },
271 .index = { { 0, 0 }, { 0 }, { 2, 1 }, { 0 } },
272 .config_map = { 5, TYPE_SCE, TYPE_CPE, TYPE_CPE, TYPE_CPE, TYPE_LFE },
273 .reorder_map = { 2, 0, 1, 4, 5, 6, 7, 3 },
274 },
275 {
276 .layout = AV_CHANNEL_LAYOUT_7POINT1_WIDE,
277 .num_ele = { 3, 0, 1, 1 },
278 .pairing = { { 0, 1, 1 }, { 0 }, { 1 }, },
279 .index = { { 0, 0, 1 }, { 0 }, { 2 }, { 0 }, },
280 .config_map = { 5, TYPE_SCE, TYPE_CPE, TYPE_CPE, TYPE_CPE, TYPE_LFE },
281 .reorder_map = { 2, 4, 5, 0, 1, 6, 7, 3 },
282 },
283 {
284 .layout = AV_CHANNEL_LAYOUT_7POINT1_WIDE_BACK,
285 .num_ele = { 3, 0, 1, 1 },
286 .pairing = { { 0, 1, 1 }, { 0 }, { 1 } },
287 .index = { { 0, 0, 1 }, { 0 }, { 2 }, { 0 } },
288 .config_map = { 5, TYPE_SCE, TYPE_CPE, TYPE_CPE, TYPE_CPE, TYPE_LFE },
289 .reorder_map = { 2, 6, 7, 0, 1, 4, 5, 3 },
290 },
291 {
292 .layout = AV_CHANNEL_LAYOUT_OCTAGONAL,
293 .num_ele = { 2, 0, 3, 0 },
294 .pairing = { { 0, 1 }, { 0 }, { 1, 1, 0 }, },
295 .index = { { 0, 0 }, { 0 }, { 1, 2, 1 }, },
296 .config_map = { 5, TYPE_SCE, TYPE_CPE, TYPE_CPE, TYPE_CPE, TYPE_SCE },
297 .reorder_map = { 2, 0, 1, 6, 7, 3, 4, 5 },
298 },
299 };
300
301 static void put_pce(PutBitContext *pb, AVCodecContext *avctx)
302 {
303 int i, j;
304 AACEncContext *s = avctx->priv_data;
305 AACPCEInfo *pce = &s->pce;
306 const int bitexact = avctx->flags & AV_CODEC_FLAG_BITEXACT;
307 const char *aux_data = bitexact ? "Lavc" : LIBAVCODEC_IDENT;
308
309 put_bits(pb, 4, 0);
310
311 put_bits(pb, 2, avctx->profile);
312 put_bits(pb, 4, s->samplerate_index);
313
314 put_bits(pb, 4, pce->num_ele[0]); /* Front */
315 put_bits(pb, 4, pce->num_ele[1]); /* Side */
316 put_bits(pb, 4, pce->num_ele[2]); /* Back */
317 put_bits(pb, 2, pce->num_ele[3]); /* LFE */
318 put_bits(pb, 3, 0); /* Assoc data */
319 put_bits(pb, 4, 0); /* CCs */
320
321 put_bits(pb, 1, 0); /* Stereo mixdown */
322 put_bits(pb, 1, 0); /* Mono mixdown */
323 put_bits(pb, 1, 0); /* Something else */
324
325 for (i = 0; i < 4; i++) {
326 for (j = 0; j < pce->num_ele[i]; j++) {
327 if (i < 3)
328 put_bits(pb, 1, pce->pairing[i][j]);
329 put_bits(pb, 4, pce->index[i][j]);
330 }
331 }
332
333 align_put_bits(pb);
334 put_bits(pb, 8, strlen(aux_data));
335 ff_put_string(pb, aux_data, 0);
336 }
337
338 /**
339 * Make AAC audio config object.
340 * @see 1.6.2.1 "Syntax - AudioSpecificConfig"
341 */
342 16 static int put_audio_specific_config(AVCodecContext *avctx, int chcfg)
343 {
344 PutBitContext pb;
345 16 AACEncContext *s = avctx->priv_data;
346 16 const int max_size = 32;
347
348 16 avctx->extradata = av_mallocz(max_size);
349
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 16 times.
16 if (!avctx->extradata)
350 return AVERROR(ENOMEM);
351
352 16 init_put_bits(&pb, avctx->extradata, max_size);
353 16 put_bits(&pb, 5, s->profile+1); //profile
354 16 put_bits(&pb, 4, s->samplerate_index); //sample rate index
355 16 put_bits(&pb, 4, chcfg);
356 //GASpecificConfig
357 16 put_bits(&pb, 1, 0); //frame length - 1024 samples
358 16 put_bits(&pb, 1, 0); //does not depend on core coder
359 16 put_bits(&pb, 1, 0); //is not extension
360
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 16 times.
16 if (s->needs_pce)
361 put_pce(&pb, avctx);
362
363 //Explicitly Mark SBR absent
364 16 put_bits(&pb, 11, 0x2b7); //sync extension
365 16 put_bits(&pb, 5, AOT_SBR);
366 16 put_bits(&pb, 1, 0);
367 16 flush_put_bits(&pb);
368 16 avctx->extradata_size = put_bytes_output(&pb);
369
370 16 return 0;
371 }
372
373 10012 void ff_quantize_band_cost_cache_init(struct AACEncContext *s)
374 {
375 10012 ++s->quantize_band_cost_cache_generation;
376
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 10012 times.
10012 if (s->quantize_band_cost_cache_generation == 0) {
377 memset(s->quantize_band_cost_cache, 0, sizeof(s->quantize_band_cost_cache));
378 s->quantize_band_cost_cache_generation = 1;
379 }
380 10012 }
381
382 #define WINDOW_FUNC(type) \
383 static void apply_ ##type ##_window(AVFloatDSPContext *fdsp, \
384 SingleChannelElement *sce, \
385 const float *audio)
386
387 5915 WINDOW_FUNC(only_long)
388 {
389
2/2
✓ Branch 0 taken 5867 times.
✓ Branch 1 taken 48 times.
5915 const float *lwindow = sce->ics.use_kb_window[0] ? ff_aac_kbd_long_1024 : ff_sine_1024;
390
2/2
✓ Branch 0 taken 5852 times.
✓ Branch 1 taken 63 times.
5915 const float *pwindow = sce->ics.use_kb_window[1] ? ff_aac_kbd_long_1024 : ff_sine_1024;
391 5915 float *out = sce->ret_buf;
392
393 5915 fdsp->vector_fmul (out, audio, lwindow, 1024);
394 5915 fdsp->vector_fmul_reverse(out + 1024, audio + 1024, pwindow, 1024);
395 5915 }
396
397 194 WINDOW_FUNC(long_start)
398 {
399
2/2
✓ Branch 0 taken 175 times.
✓ Branch 1 taken 19 times.
194 const float *lwindow = sce->ics.use_kb_window[1] ? ff_aac_kbd_long_1024 : ff_sine_1024;
400
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 194 times.
194 const float *swindow = sce->ics.use_kb_window[0] ? ff_aac_kbd_short_128 : ff_sine_128;
401 194 float *out = sce->ret_buf;
402
403 194 fdsp->vector_fmul(out, audio, lwindow, 1024);
404 194 memcpy(out + 1024, audio + 1024, sizeof(out[0]) * 448);
405 194 fdsp->vector_fmul_reverse(out + 1024 + 448, audio + 1024 + 448, swindow, 128);
406 194 memset(out + 1024 + 576, 0, sizeof(out[0]) * 448);
407 194 }
408
409 194 WINDOW_FUNC(long_stop)
410 {
411
1/2
✓ Branch 0 taken 194 times.
✗ Branch 1 not taken.
194 const float *lwindow = sce->ics.use_kb_window[0] ? ff_aac_kbd_long_1024 : ff_sine_1024;
412
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 194 times.
194 const float *swindow = sce->ics.use_kb_window[1] ? ff_aac_kbd_short_128 : ff_sine_128;
413 194 float *out = sce->ret_buf;
414
415 194 memset(out, 0, sizeof(out[0]) * 448);
416 194 fdsp->vector_fmul(out + 448, audio + 448, swindow, 128);
417 194 memcpy(out + 576, audio + 576, sizeof(out[0]) * 448);
418 194 fdsp->vector_fmul_reverse(out + 1024, audio + 1024, lwindow, 1024);
419 194 }
420
421 229 WINDOW_FUNC(eight_short)
422 {
423
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 229 times.
229 const float *swindow = sce->ics.use_kb_window[0] ? ff_aac_kbd_short_128 : ff_sine_128;
424
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 229 times.
229 const float *pwindow = sce->ics.use_kb_window[1] ? ff_aac_kbd_short_128 : ff_sine_128;
425 229 const float *in = audio + 448;
426 229 float *out = sce->ret_buf;
427 int w;
428
429
2/2
✓ Branch 0 taken 1832 times.
✓ Branch 1 taken 229 times.
2061 for (w = 0; w < 8; w++) {
430
2/2
✓ Branch 0 taken 1603 times.
✓ Branch 1 taken 229 times.
1832 fdsp->vector_fmul (out, in, w ? pwindow : swindow, 128);
431 1832 out += 128;
432 1832 in += 128;
433 1832 fdsp->vector_fmul_reverse(out, in, swindow, 128);
434 1832 out += 128;
435 }
436 229 }
437
438 static void (*const apply_window[4])(AVFloatDSPContext *fdsp,
439 SingleChannelElement *sce,
440 const float *audio) = {
441 [ONLY_LONG_SEQUENCE] = apply_only_long_window,
442 [LONG_START_SEQUENCE] = apply_long_start_window,
443 [EIGHT_SHORT_SEQUENCE] = apply_eight_short_window,
444 [LONG_STOP_SEQUENCE] = apply_long_stop_window
445 };
446
447 6532 static void apply_window_and_mdct(AACEncContext *s, SingleChannelElement *sce,
448 float *audio)
449 {
450 int i;
451 6532 float *output = sce->ret_buf;
452
453 6532 apply_window[sce->ics.window_sequence[0]](s->fdsp, sce, audio);
454
455
2/2
✓ Branch 0 taken 6303 times.
✓ Branch 1 taken 229 times.
6532 if (sce->ics.window_sequence[0] != EIGHT_SHORT_SEQUENCE)
456 6303 s->mdct1024_fn(s->mdct1024, sce->coeffs, output, sizeof(float));
457 else
458
2/2
✓ Branch 0 taken 1832 times.
✓ Branch 1 taken 229 times.
2061 for (i = 0; i < 1024; i += 128)
459 1832 s->mdct128_fn(s->mdct128, &sce->coeffs[i], output + i*2, sizeof(float));
460 6532 memcpy(audio, audio + 1024, sizeof(audio[0]) * 1024);
461 6532 memcpy(sce->pcoeffs, sce->coeffs, sizeof(sce->pcoeffs));
462 6532 }
463
464 /**
465 * Encode ics_info element.
466 * @see Table 4.6 (syntax of ics_info)
467 */
468 5070 static void put_ics_info(AACEncContext *s, IndividualChannelStream *info)
469 {
470 int w;
471
472 5070 put_bits(&s->pb, 1, 0); // ics_reserved bit
473 5070 put_bits(&s->pb, 2, info->window_sequence[0]);
474 5070 put_bits(&s->pb, 1, info->use_kb_window[0]);
475
2/2
✓ Branch 0 taken 4948 times.
✓ Branch 1 taken 122 times.
5070 if (info->window_sequence[0] != EIGHT_SHORT_SEQUENCE) {
476 4948 put_bits(&s->pb, 6, info->max_sfb);
477 4948 put_bits(&s->pb, 1, 0); /* No predictor present */
478 } else {
479 122 put_bits(&s->pb, 4, info->max_sfb);
480
2/2
✓ Branch 0 taken 854 times.
✓ Branch 1 taken 122 times.
976 for (w = 1; w < 8; w++)
481 854 put_bits(&s->pb, 1, !info->group_len[w]);
482 }
483 5070 }
484
485 /**
486 * Encode MS data.
487 * @see 4.6.8.1 "Joint Coding - M/S Stereo"
488 */
489 4388 static void encode_ms_info(PutBitContext *pb, ChannelElement *cpe)
490 {
491 int i, w;
492
493 4388 put_bits(pb, 2, cpe->ms_mode);
494
2/2
✓ Branch 0 taken 366 times.
✓ Branch 1 taken 4022 times.
4388 if (cpe->ms_mode == 1)
495
2/2
✓ Branch 0 taken 404 times.
✓ Branch 1 taken 366 times.
770 for (w = 0; w < cpe->ch[0].ics.num_windows; w += cpe->ch[0].ics.group_len[w])
496
2/2
✓ Branch 0 taken 17996 times.
✓ Branch 1 taken 404 times.
18400 for (i = 0; i < cpe->ch[0].ics.max_sfb; i++)
497 17996 put_bits(pb, 1, cpe->ms_mask[w*16 + i]);
498 4388 }
499
500 /**
501 * Produce integer coefficients from scalefactors provided by the model.
502 */
503 5070 static void adjust_frame_information(ChannelElement *cpe, int chans)
504 {
505 int i, w, w2, g, ch;
506 int maxsfb, cmaxsfb;
507
508
2/2
✓ Branch 0 taken 9458 times.
✓ Branch 1 taken 5070 times.
14528 for (ch = 0; ch < chans; ch++) {
509 9458 IndividualChannelStream *ics = &cpe->ch[ch].ics;
510 9458 maxsfb = 0;
511 9458 cpe->ch[ch].pulse.num_pulse = 0;
512
2/2
✓ Branch 0 taken 10146 times.
✓ Branch 1 taken 9458 times.
19604 for (w = 0; w < ics->num_windows; w += ics->group_len[w]) {
513
4/4
✓ Branch 0 taken 18309 times.
✓ Branch 1 taken 103 times.
✓ Branch 2 taken 8266 times.
✓ Branch 3 taken 10043 times.
18412 for (cmaxsfb = ics->num_swb; cmaxsfb > 0 && cpe->ch[ch].zeroes[w*16+cmaxsfb-1]; cmaxsfb--)
514 ;
515 10146 maxsfb = FFMAX(maxsfb, cmaxsfb);
516 }
517 9458 ics->max_sfb = maxsfb;
518
519 //adjust zero bands for window groups
520
2/2
✓ Branch 0 taken 10146 times.
✓ Branch 1 taken 9458 times.
19604 for (w = 0; w < ics->num_windows; w += ics->group_len[w]) {
521
2/2
✓ Branch 0 taken 452675 times.
✓ Branch 1 taken 10146 times.
462821 for (g = 0; g < ics->max_sfb; g++) {
522 452675 i = 1;
523
2/2
✓ Branch 0 taken 452765 times.
✓ Branch 1 taken 1960 times.
454725 for (w2 = w; w2 < w + ics->group_len[w]; w2++) {
524
2/2
✓ Branch 0 taken 450715 times.
✓ Branch 1 taken 2050 times.
452765 if (!cpe->ch[ch].zeroes[w2*16 + g]) {
525 450715 i = 0;
526 450715 break;
527 }
528 }
529 452675 cpe->ch[ch].zeroes[w*16 + g] = i;
530 }
531 }
532 }
533
534
3/4
✓ Branch 0 taken 4388 times.
✓ Branch 1 taken 682 times.
✓ Branch 2 taken 4388 times.
✗ Branch 3 not taken.
5070 if (chans > 1 && cpe->common_window) {
535 4388 IndividualChannelStream *ics0 = &cpe->ch[0].ics;
536 4388 IndividualChannelStream *ics1 = &cpe->ch[1].ics;
537 4388 int msc = 0;
538 4388 ics0->max_sfb = FFMAX(ics0->max_sfb, ics1->max_sfb);
539 4388 ics1->max_sfb = ics0->max_sfb;
540
2/2
✓ Branch 0 taken 5221 times.
✓ Branch 1 taken 4388 times.
9609 for (w = 0; w < ics0->num_windows*16; w += 16)
541
2/2
✓ Branch 0 taken 220551 times.
✓ Branch 1 taken 5221 times.
225772 for (i = 0; i < ics0->max_sfb; i++)
542
2/2
✓ Branch 0 taken 28479 times.
✓ Branch 1 taken 192072 times.
220551 if (cpe->ms_mask[w+i])
543 28479 msc++;
544
3/4
✓ Branch 0 taken 852 times.
✓ Branch 1 taken 3536 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 852 times.
4388 if (msc == 0 || ics0->max_sfb == 0)
545 3536 cpe->ms_mode = 0;
546 else
547
2/2
✓ Branch 0 taken 366 times.
✓ Branch 1 taken 486 times.
852 cpe->ms_mode = msc < ics0->max_sfb * ics0->num_windows ? 1 : 2;
548 }
549 5070 }
550
551 740 static void apply_intensity_stereo(ChannelElement *cpe)
552 {
553 int w, w2, g, i;
554 740 IndividualChannelStream *ics = &cpe->ch[0].ics;
555
2/2
✓ Branch 0 taken 164 times.
✓ Branch 1 taken 576 times.
740 if (!cpe->common_window)
556 164 return;
557
2/2
✓ Branch 0 taken 633 times.
✓ Branch 1 taken 576 times.
1209 for (w = 0; w < ics->num_windows; w += ics->group_len[w]) {
558
2/2
✓ Branch 0 taken 716 times.
✓ Branch 1 taken 633 times.
1349 for (w2 = 0; w2 < ics->group_len[w]; w2++) {
559 716 int start = (w+w2) * 128;
560
2/2
✓ Branch 0 taken 29484 times.
✓ Branch 1 taken 716 times.
30200 for (g = 0; g < ics->num_swb; g++) {
561 29484 int p = -1 + 2 * (cpe->ch[1].band_type[w*16+g] - 14);
562 29484 float scale = cpe->ch[0].is_ener[w*16+g];
563
2/2
✓ Branch 0 taken 20331 times.
✓ Branch 1 taken 9153 times.
29484 if (!cpe->is_mask[w*16 + g]) {
564 20331 start += ics->swb_sizes[g];
565 20331 continue;
566 }
567
2/2
✓ Branch 0 taken 2365 times.
✓ Branch 1 taken 6788 times.
9153 if (cpe->ms_mask[w*16 + g])
568 2365 p *= -1;
569
2/2
✓ Branch 0 taken 297772 times.
✓ Branch 1 taken 9153 times.
306925 for (i = 0; i < ics->swb_sizes[g]; i++) {
570 297772 float sum = (cpe->ch[0].coeffs[start+i] + p*cpe->ch[1].coeffs[start+i])*scale;
571 297772 cpe->ch[0].coeffs[start+i] = sum;
572 297772 cpe->ch[1].coeffs[start+i] = 0.0f;
573 }
574 9153 start += ics->swb_sizes[g];
575 }
576 }
577 }
578 }
579
580 /* I/S acceptance level for the image-error EMA at full rate pressure */
581 #define NMR_IS_IMG_GATE 8000.0f
582
583 /* Frequency in Hz for the lower limit of intensity stereo */
584 #define NMR_IS_LOW_LIMIT 6100
585
586 /* M/S adoption: es < 0.5*em, content-driven and rate-free */
587 #define NMR_MS_EQUIV 0.5f
588 #define NMR_MS_MASK 0.0f
589
590 /* Pair decouple threshold on the joint-tool candidacy fraction EMA: pairs
591 * whose joint tools are mostly dead (diffuse decorrelated content) window
592 * per-channel and skip M/S; recouple above 1.3x. */
593 #define NMR_DECORR_LO 0.20f
594
595 /* Stereo-decision hysteresis: leaving a joint mode costs a margin. */
596 #define NMR_STICKY 2.0f
597
598 /* Decision statistics are EMA-smoothed across frames. */
599 #define NMR_SDEC_EMA 0.75f
600
601 /* PNS-stereo gate: substitute only clearly-decorrelated (wide) bands. */
602 #define NMR_PNS_STEREO_DECORR 0.6f
603
604 /* Recode one band's window group as mid+side in place. */
605 4958 static void nmr_apply_ms_band(AACEncContext *s, ChannelElement *cpe,
606 int w, int g, int start, int len, int gl)
607 {
608 4958 SingleChannelElement *sce0 = &cpe->ch[0];
609 4958 SingleChannelElement *sce1 = &cpe->ch[1];
610 4958 cpe->ms_mask[w*16+g] = 1;
611
2/2
✓ Branch 0 taken 4986 times.
✓ Branch 1 taken 4958 times.
9944 for (int w2 = 0; w2 < gl; w2++) {
612 4986 FFPsyBand *b0 = &s->psy.ch[s->cur_channel+0].psy_bands[(w+w2)*16+g];
613 4986 FFPsyBand *b1 = &s->psy.ch[s->cur_channel+1].psy_bands[(w+w2)*16+g];
614 4986 float *L = sce0->coeffs + start + (w+w2)*128;
615 4986 float *R = sce1->coeffs + start + (w+w2)*128;
616 4986 float em = 0.0f, es = 0.0f;
617
2/2
✓ Branch 0 taken 89848 times.
✓ Branch 1 taken 4986 times.
94834 for (int i = 0; i < len; i++) {
618 89848 float m = (L[i] + R[i]) * 0.5f;
619 89848 R[i] = m - R[i]; L[i] = m;
620 89848 em += L[i]*L[i]; es += R[i]*R[i];
621 }
622
2/2
✓ Branch 0 taken 190 times.
✓ Branch 1 taken 4796 times.
4986 b0->threshold = FFMIN(b0->threshold, b1->threshold) * 0.5f;
623 4986 b1->threshold = b0->threshold;
624 4986 b0->energy = em; b1->energy = es;
625 }
626 4958 }
627
628 /* I/S perceptual test: reconstruction image error vs the pair's masks. */
629 2423 static int nmr_is_image_masked(AACEncContext *s, ChannelElement *cpe,
630 int w, int g, int start, int len, int gl,
631 float ener0, float ener1, float dot,
632 float minthr0, float minthr1, float *ratio_out,
633 float *scale_out, float *sr_out, int *p_out)
634 {
635
2/2
✓ Branch 0 taken 2363 times.
✓ Branch 1 taken 60 times.
2423 int p = dot >= 0.0f ? 1 : -1;
636 2423 float ener01 = ener0 + ener1 + 2*p*dot; /* energy of L + p*R */
637 2423 *ratio_out = FLT_MAX;
638
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2423 times.
2423 if (ener01 <= FLT_MIN)
639 return 0;
640 2423 float scale = sqrtf(ener0 / ener01); /* carrier = (L + p*R)*scale */
641 2423 float sr_ = sqrtf(ener1 / ener0); /* decoder: R = p*sr_*carrier */
642 2423 float img0 = 0.0f, img1 = 0.0f;
643
2/2
✓ Branch 0 taken 2451 times.
✓ Branch 1 taken 2423 times.
4874 for (int w2 = 0; w2 < gl; w2++) {
644 2451 const float *L = cpe->ch[0].coeffs + start + (w+w2)*128;
645 2451 const float *R = cpe->ch[1].coeffs + start + (w+w2)*128;
646
2/2
✓ Branch 0 taken 85236 times.
✓ Branch 1 taken 2451 times.
87687 for (int i = 0; i < len; i++) {
647 85236 float c = (L[i] + p*R[i]) * scale;
648 85236 float dl = L[i] - c, dr = R[i] - p*sr_*c;
649 85236 img0 += dl*dl; img1 += dr*dr;
650 }
651 }
652
6/10
✓ Branch 0 taken 2423 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2423 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 285 times.
✓ Branch 5 taken 2138 times.
✓ Branch 6 taken 285 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2138 times.
✗ Branch 9 not taken.
2423 *ratio_out = FFMAX(img0 / FFMAX(minthr0 * gl, FLT_MIN),
653 img1 / FFMAX(minthr1 * gl, FLT_MIN));
654 2423 *scale_out = scale; *sr_out = sr_; *p_out = p;
655 2423 return 1;
656 }
657
658 /* Recode one band's window group as intensity stereo in place: replace L with the
659 * carrier, zero R, signal the phase via the side channel's band type, and fold the
660 * pair's masking into the surviving (carrier) channel. */
661 798 static void nmr_apply_is_band(AACEncContext *s, ChannelElement *cpe,
662 int w, int g, int start, int len, int gl,
663 float scale, float sr_, int p,
664 float ener0, float ener1)
665 {
666 798 cpe->is_mask[w*16+g] = 1;
667 798 cpe->ch[0].is_ener[w*16+g] = scale;
668 798 cpe->ch[1].is_ener[w*16+g] = ener0 / ener1;
669
2/2
✓ Branch 0 taken 748 times.
✓ Branch 1 taken 50 times.
798 cpe->ch[1].band_type[w*16+g] = p > 0 ? INTENSITY_BT : INTENSITY_BT2;
670
2/2
✓ Branch 0 taken 826 times.
✓ Branch 1 taken 798 times.
1624 for (int w2 = 0; w2 < gl; w2++) {
671 826 FFPsyBand *b0 = &s->psy.ch[s->cur_channel+0].psy_bands[(w+w2)*16+g];
672 826 FFPsyBand *b1 = &s->psy.ch[s->cur_channel+1].psy_bands[(w+w2)*16+g];
673 826 float *L = cpe->ch[0].coeffs + start + (w+w2)*128;
674 826 float *R = cpe->ch[1].coeffs + start + (w+w2)*128;
675 826 float ec = 0.0f;
676
2/2
✓ Branch 0 taken 27796 times.
✓ Branch 1 taken 826 times.
28622 for (int i = 0; i < len; i++) {
677 27796 L[i] = (L[i] + p*R[i]) * scale;
678 27796 R[i] = 0.0f;
679 27796 ec += L[i]*L[i];
680 }
681
4/6
✓ Branch 0 taken 826 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 136 times.
✓ Branch 3 taken 690 times.
✓ Branch 4 taken 136 times.
✗ Branch 5 not taken.
826 b0->threshold = FFMIN(b0->threshold, b1->threshold / FFMAX(sr_*sr_, 1e-9f));
682 826 b0->energy = ec; b1->energy = 0.0f;
683 }
684 798 }
685
686 /*
687 * Per-band stereo-mode decision (L/R vs M/S vs intensity) for the NMR coder,
688 * made before quantization from the psychoacoustic model alone, so the
689 * quantizer search allocates natively on the spectra that are actually coded.
690 */
691 137 static void nmr_decide_stereo(AACEncContext *s, ChannelElement *cpe)
692 {
693 137 SingleChannelElement *sce0 = &cpe->ch[0];
694 137 SingleChannelElement *sce1 = &cpe->ch[1];
695 137 IndividualChannelStream *ics = &sce0->ics;
696 137 const AVCodecContext *avctx = s->psy.avctx;
697 137 const float freq_mult = avctx->sample_rate / (1024.0f / ics->num_windows) / 2.0f;
698 137 int is_count = 0;
699
700
1/2
✓ Branch 0 taken 137 times.
✗ Branch 1 not taken.
137 if (s->nmr) {
701 137 int pi = (s->cur_channel >> 1) & 7;
702 137 pi = pi * 2 + (ics->num_windows == 8); /* per-grid state bank */
703
2/2
✓ Branch 0 taken 8 times.
✓ Branch 1 taken 129 times.
137 if (!s->nmr->sinit[pi]) {
704 /* one-time init; per-grid banks persist across window switches
705 * (wiping them churned stereo modes audibly) */
706 8 memset(s->nmr->smode[pi], 0, sizeof(s->nmr->smode[pi]));
707
2/2
✓ Branch 0 taken 1024 times.
✓ Branch 1 taken 8 times.
1032 for (int b = 0; b < 128; b++) {
708 1024 s->nmr->sema_em[pi][b] = 0.0f;
709 1024 s->nmr->sema_img[pi][b] = -1.0f;
710 }
711 8 s->nmr->sinit[pi] = 1;
712 }
713 }
714
715 /* Per-band stereo decision (L/R vs M/S vs I/S), made pre-quantization from
716 * the psy model so the trellis allocates on the coded spectra. */
717
718 /* I/S engages under SUSTAINED strain only: rate pressure gated by the
719 * lambda floor (pressure spikes at a comfortable operating point must
720 * not admit it). Unengaged candidates fall back to M/S. */
721 411 float is_ramp = s->nmr ? s->nmr->press *
722
1/2
✓ Branch 0 taken 137 times.
✗ Branch 1 not taken.
137 av_clipf((s->nmr->lam_floor - 40.0f) / (120.0f - 40.0f), 0.0f, 1.0f) : 0.0f;
723
3/4
✓ Branch 0 taken 137 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 51 times.
✓ Branch 3 taken 86 times.
137 const int allow_is = s->options.intensity_stereo && is_ramp > 0.0f;
724
725 137 const int pidx = (s->cur_channel >> 1) & 15;
726 137 const int decoupled = s->psy.pair_decoupled[pidx];
727 137 int njoint = 0, nbands = 0; /* joint-tool candidacy census, decouple feed */
728
729
2/2
✓ Branch 0 taken 140 times.
✓ Branch 1 taken 137 times.
277 for (int w = 0; w < ics->num_windows; w += ics->group_len[w]) {
730 140 int start = 0;
731
2/2
✓ Branch 0 taken 6720 times.
✓ Branch 1 taken 140 times.
6860 for (int g = 0; g < ics->num_swb; start += ics->swb_sizes[g++]) {
732 6720 int len = ics->swb_sizes[g], gl = ics->group_len[w];
733 6720 float ener0 = 0.0f, ener1 = 0.0f, dot = 0.0f, es_tot = 0.0f, em_tot = 0.0f;
734 6720 float minthr0 = FLT_MAX, minthr1 = FLT_MAX;
735
736 6720 cpe->is_mask[w*16+g] = 0;
737 6720 cpe->ms_mask[w*16+g] = 0;
738
739
2/2
✓ Branch 0 taken 6776 times.
✓ Branch 1 taken 6720 times.
13496 for (int w2 = 0; w2 < gl; w2++) {
740 6776 FFPsyBand *b0 = &s->psy.ch[s->cur_channel+0].psy_bands[(w+w2)*16+g];
741 6776 FFPsyBand *b1 = &s->psy.ch[s->cur_channel+1].psy_bands[(w+w2)*16+g];
742 6776 const float *L = sce0->coeffs + start + (w+w2)*128;
743 6776 const float *R = sce1->coeffs + start + (w+w2)*128;
744 6776 float el = 0.0f, er = 0.0f, em = 0.0f, es = 0.0f, d = 0.0f;
745
2/2
✓ Branch 0 taken 140288 times.
✓ Branch 1 taken 6776 times.
147064 for (int i = 0; i < len; i++) {
746 140288 float m = (L[i] + R[i]) * 0.5f;
747 140288 float sv = m - R[i];
748 140288 el += L[i]*L[i]; er += R[i]*R[i];
749 140288 em += m*m; es += sv*sv; d += L[i]*R[i];
750 }
751 6776 ener0 += el; ener1 += er; dot += d; es_tot += es; em_tot += em;
752
2/2
✓ Branch 0 taken 6745 times.
✓ Branch 1 taken 31 times.
6776 minthr0 = FFMIN(minthr0, b0->threshold);
753
2/2
✓ Branch 0 taken 6745 times.
✓ Branch 1 taken 31 times.
6776 minthr1 = FFMIN(minthr1, b1->threshold);
754 }
755
2/2
✓ Branch 0 taken 274 times.
✓ Branch 1 taken 6446 times.
6720 float thr_g = FFMIN(minthr0, minthr1) * gl; /* group masking budget */
756
757 /* PNS-stereo reservation: keep clearly-wide noise bands for PNS. */
758 6720 const int sidx = w*16+g;
759 {
760 6720 float es_w = es_tot, em_w = em_tot;
761
1/2
✓ Branch 0 taken 6720 times.
✗ Branch 1 not taken.
6720 if (s->nmr) {
762 6720 int pi_ = ((s->cur_channel >> 1) & 7) * 2 + (cpe->ch[0].ics.num_windows == 8);
763 6720 float pe = s->nmr->sema_es[pi_][sidx];
764 6720 float pm = s->nmr->sema_em[pi_][sidx];
765
2/2
✓ Branch 0 taken 5821 times.
✓ Branch 1 taken 899 times.
6720 if (pm > 0.0f) {
766 5821 es_w = NMR_SDEC_EMA * pe + (1.0f - NMR_SDEC_EMA) * es_tot;
767 5821 em_w = NMR_SDEC_EMA * pm + (1.0f - NMR_SDEC_EMA) * em_tot;
768 }
769 }
770
3/4
✓ Branch 0 taken 3014 times.
✓ Branch 1 taken 3706 times.
✓ Branch 2 taken 3014 times.
✗ Branch 3 not taken.
6720 if (cpe->ch[0].can_pns[w*16+g] && cpe->ch[1].can_pns[w*16+g] &&
771
2/2
✓ Branch 0 taken 415 times.
✓ Branch 1 taken 2599 times.
3014 es_w > NMR_PNS_STEREO_DECORR * em_w)
772 415 continue;
773 }
774 6305 cpe->ch[0].can_pns[w*16+g] = cpe->ch[1].can_pns[w*16+g] = 0;
775
776 6305 int pi = ((s->cur_channel >> 1) & 7) * 2 + (cpe->ch[0].ics.num_windows == 8);
777
1/2
✓ Branch 0 taken 6305 times.
✗ Branch 1 not taken.
6305 uint8_t *pmode = s->nmr ? s->nmr->smode[pi] : NULL;
778
1/2
✓ Branch 0 taken 6305 times.
✗ Branch 1 not taken.
6305 int prev = pmode ? pmode[sidx] : 0;
779
2/2
✓ Branch 0 taken 4756 times.
✓ Branch 1 taken 1549 times.
6305 float eqgate = NMR_MS_EQUIV * (prev == 1 ? 1.5f : 1.0f); /* stay-until es>0.75em */
780 /* I/S = lossy economy: image-error budget scales with pressure */
781
2/2
✓ Branch 0 taken 690 times.
✓ Branch 1 taken 5615 times.
6305 float imgate = NMR_IS_IMG_GATE * is_ramp * (prev == 2 ? NMR_STICKY : 1.0f);
782 6305 float es_d = es_tot, em_d = em_tot;
783
1/2
✓ Branch 0 taken 6305 times.
✗ Branch 1 not taken.
6305 if (s->nmr) {
784 6305 float *ees = &s->nmr->sema_es[pi][sidx];
785 6305 float *eem = &s->nmr->sema_em[pi][sidx];
786
2/2
✓ Branch 0 taken 624 times.
✓ Branch 1 taken 5681 times.
6305 if (*eem <= 0.0f) { *ees = es_tot; *eem = em_tot; }
787 else {
788 5681 *ees = NMR_SDEC_EMA * *ees + (1.0f - NMR_SDEC_EMA) * es_tot;
789 5681 *eem = NMR_SDEC_EMA * *eem + (1.0f - NMR_SDEC_EMA) * em_tot;
790 }
791 6305 es_d = *ees; em_d = *eem;
792 }
793
1/2
✓ Branch 0 taken 6305 times.
✗ Branch 1 not taken.
12610 int ms_would = s->options.mid_side &&
794
1/2
✓ Branch 0 taken 6305 times.
✗ Branch 1 not taken.
6305 (s->options.mid_side == 1 ||
795
2/2
✓ Branch 0 taken 639 times.
✓ Branch 1 taken 5666 times.
6305 es_d < eqgate * em_d ||
796
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 639 times.
639 es_tot < NMR_MS_MASK * thr_g);
797
3/4
✓ Branch 0 taken 5666 times.
✓ Branch 1 taken 639 times.
✓ Branch 2 taken 5666 times.
✗ Branch 3 not taken.
6305 int ms_ok = ms_would && !decoupled;
798 float scale, sr_, imgratio; int p;
799 /* I/S competes with M/S above the frequency limit (candidacy must
800 * not be gated on !ms_ok - that leaves only unrenderable bands) */
801
2/2
✓ Branch 0 taken 2423 times.
✓ Branch 1 taken 126 times.
2549 int is_cand = start * freq_mult > NMR_IS_LOW_LIMIT &&
802
4/6
✓ Branch 0 taken 2549 times.
✓ Branch 1 taken 3756 times.
✓ Branch 2 taken 2423 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 2423 times.
✗ Branch 5 not taken.
11277 ener0 > FLT_MIN && ener1 > FLT_MIN &&
803 2423 nmr_is_image_masked(s, cpe, w, g, start, len, gl,
804 ener0, ener1, dot, minthr0, minthr1,
805 &imgratio, &scale, &sr_, &p);
806 6305 int is_ok = is_cand;
807
3/4
✓ Branch 0 taken 6305 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2549 times.
✓ Branch 3 taken 3756 times.
6305 if (s->nmr && start * freq_mult > NMR_IS_LOW_LIMIT) {
808 /* smoothed image-error; updated only while candidate (fail-value
809 * feeding jammed it permanently high) */
810 2549 float *eim = &s->nmr->sema_img[pi][sidx];
811
2/2
✓ Branch 0 taken 2423 times.
✓ Branch 1 taken 126 times.
2549 if (is_cand) {
812 /* seed from first measurement; freeze when not candidate */
813
2/2
✓ Branch 0 taken 139 times.
✓ Branch 1 taken 2284 times.
2423 if (*eim < 0.0f) *eim = imgratio;
814
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2284 times.
2284 else *eim = NMR_SDEC_EMA * *eim + (1.0f - NMR_SDEC_EMA) * FFMIN(imgratio, 100.0f * NMR_IS_IMG_GATE);
815 }
816
5/6
✓ Branch 0 taken 2423 times.
✓ Branch 1 taken 126 times.
✓ Branch 2 taken 2423 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 798 times.
✓ Branch 5 taken 1625 times.
2549 is_ok = is_cand && *eim >= 0.0f && *eim < imgate;
817 }
818
819
4/4
✓ Branch 0 taken 639 times.
✓ Branch 1 taken 5666 times.
✓ Branch 2 taken 90 times.
✓ Branch 3 taken 549 times.
6305 njoint += ms_would || is_ok; nbands++;
820
1/2
✓ Branch 0 taken 6305 times.
✗ Branch 1 not taken.
6305 if (pmode) {
821
6/8
✓ Branch 0 taken 798 times.
✓ Branch 1 taken 5507 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 798 times.
✓ Branch 4 taken 549 times.
✓ Branch 5 taken 4958 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 549 times.
6305 int m_ = (is_ok && allow_is) ? 2 : ms_ok ? 1 :
822 (is_ok && s->options.mid_side) ? 1 : 0;
823 6305 pmode[sidx] = m_;
824 6305 s->nmr->smode_band[(s->cur_channel >> 1) & 7][w*16+g] = m_;
825 }
826
3/4
✓ Branch 0 taken 798 times.
✓ Branch 1 taken 5507 times.
✓ Branch 2 taken 798 times.
✗ Branch 3 not taken.
6305 if (is_ok && allow_is) {
827 798 nmr_apply_is_band(s, cpe, w, g, start, len, gl,
828 scale, sr_, p, ener0, ener1);
829 798 is_count++;
830
3/6
✓ Branch 0 taken 549 times.
✓ Branch 1 taken 4958 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 549 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
5507 } else if (ms_ok || (is_ok && s->options.mid_side)) {
831 4958 nmr_apply_ms_band(s, cpe, w, g, start, len, gl);
832 }
833 /* else: keep full L/R stereo */
834 }
835 }
836 137 cpe->is_mode = !!is_count;
837
838
1/2
✓ Branch 0 taken 137 times.
✗ Branch 1 not taken.
137 if (nbands > 0) {
839 /* Pair joint-tool value, read next frame by the psy pair-synced window
840 * decision and the M/S candidacy above. Measured as CANDIDACY (not
841 * adoption) so decoupling cannot starve its own signal and self-lock. */
842 137 float r = (float)njoint / nbands;
843 137 float *pj = &s->psy.pair_joint[pidx];
844
2/2
✓ Branch 0 taken 125 times.
✓ Branch 1 taken 12 times.
137 *pj = *pj > 0.0f ? 0.95f * *pj + 0.05f * r : r;
845 137 s->psy.pair_decoupled[pidx] = *pj <
846
2/2
✓ Branch 0 taken 5 times.
✓ Branch 1 taken 132 times.
137 (s->psy.pair_decoupled[pidx] ? 1.3f * NMR_DECORR_LO : NMR_DECORR_LO);
847 }
848 137 }
849
850 411 static void apply_mid_side_stereo(ChannelElement *cpe)
851 {
852 int w, w2, g, i;
853 411 IndividualChannelStream *ics = &cpe->ch[0].ics;
854
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 411 times.
411 if (!cpe->common_window)
855 return;
856
2/2
✓ Branch 0 taken 462 times.
✓ Branch 1 taken 411 times.
873 for (w = 0; w < ics->num_windows; w += ics->group_len[w]) {
857
2/2
✓ Branch 0 taken 537 times.
✓ Branch 1 taken 462 times.
999 for (w2 = 0; w2 < ics->group_len[w]; w2++) {
858 537 int start = (w+w2) * 128;
859
2/2
✓ Branch 0 taken 21273 times.
✓ Branch 1 taken 537 times.
21810 for (g = 0; g < ics->num_swb; g++) {
860 /* ms_mask can be used for other purposes in PNS and I/S,
861 * so must not apply M/S if any band uses either, even if
862 * ms_mask is set.
863 */
864
2/4
✓ Branch 0 taken 21273 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 21273 times.
✗ Branch 3 not taken.
21273 if (!cpe->ms_mask[w*16 + g] || cpe->is_mask[w*16 + g]
865
1/2
✓ Branch 0 taken 21273 times.
✗ Branch 1 not taken.
21273 || cpe->ch[0].band_type[w*16 + g] >= NOISE_BT
866
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 21273 times.
21273 || cpe->ch[1].band_type[w*16 + g] >= NOISE_BT) {
867 start += ics->swb_sizes[g];
868 continue;
869 }
870
2/2
✓ Branch 0 taken 420864 times.
✓ Branch 1 taken 21273 times.
442137 for (i = 0; i < ics->swb_sizes[g]; i++) {
871 420864 float L = (cpe->ch[0].coeffs[start+i] + cpe->ch[1].coeffs[start+i]) * 0.5f;
872 420864 float R = L - cpe->ch[1].coeffs[start+i];
873 420864 cpe->ch[0].coeffs[start+i] = L;
874 420864 cpe->ch[1].coeffs[start+i] = R;
875 }
876 21273 start += ics->swb_sizes[g];
877 }
878 }
879 }
880 }
881
882 /**
883 * Encode scalefactor band coding type.
884 */
885 9458 static void encode_band_info(AACEncContext *s, SingleChannelElement *sce)
886 {
887 int w;
888
889
1/2
✓ Branch 0 taken 9458 times.
✗ Branch 1 not taken.
9458 if (s->coder->set_special_band_scalefactors)
890 9458 s->coder->set_special_band_scalefactors(s, sce);
891
892
2/2
✓ Branch 0 taken 10146 times.
✓ Branch 1 taken 9458 times.
19604 for (w = 0; w < sce->ics.num_windows; w += sce->ics.group_len[w])
893 10146 s->coder->encode_window_bands_info(s, sce, w, sce->ics.group_len[w], s->lambda);
894 9458 }
895
896 /**
897 * Encode scalefactors.
898 */
899 9458 static void encode_scale_factors(AVCodecContext *avctx, AACEncContext *s,
900 SingleChannelElement *sce)
901 {
902 9458 int diff, off_sf = sce->sf_idx[0], off_pns = sce->sf_idx[0] - NOISE_OFFSET;
903 9458 int off_is = 0, noise_flag = 1;
904 int i, w;
905
906
2/2
✓ Branch 0 taken 10146 times.
✓ Branch 1 taken 9458 times.
19604 for (w = 0; w < sce->ics.num_windows; w += sce->ics.group_len[w]) {
907
2/2
✓ Branch 0 taken 456327 times.
✓ Branch 1 taken 10146 times.
466473 for (i = 0; i < sce->ics.max_sfb; i++) {
908
2/2
✓ Branch 0 taken 444210 times.
✓ Branch 1 taken 12117 times.
456327 if (!sce->zeroes[w*16 + i]) {
909
2/2
✓ Branch 0 taken 5957 times.
✓ Branch 1 taken 438253 times.
444210 if (sce->band_type[w*16 + i] == NOISE_BT) {
910 5957 diff = sce->sf_idx[w*16 + i] - off_pns;
911 5957 off_pns = sce->sf_idx[w*16 + i];
912
2/2
✓ Branch 0 taken 1025 times.
✓ Branch 1 taken 4932 times.
5957 if (noise_flag-- > 0) {
913 1025 put_bits(&s->pb, NOISE_PRE_BITS, diff + NOISE_PRE);
914 1025 continue;
915 }
916
2/2
✓ Branch 0 taken 430934 times.
✓ Branch 1 taken 7319 times.
438253 } else if (sce->band_type[w*16 + i] == INTENSITY_BT ||
917
2/2
✓ Branch 0 taken 2129 times.
✓ Branch 1 taken 428805 times.
430934 sce->band_type[w*16 + i] == INTENSITY_BT2) {
918 9448 diff = sce->sf_idx[w*16 + i] - off_is;
919 9448 off_is = sce->sf_idx[w*16 + i];
920 } else {
921 428805 diff = sce->sf_idx[w*16 + i] - off_sf;
922 428805 off_sf = sce->sf_idx[w*16 + i];
923 }
924 443185 diff += SCALE_DIFF_ZERO;
925
2/4
✓ Branch 0 taken 443185 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 443185 times.
443185 av_assert0(diff >= 0 && diff <= 120);
926 443185 put_bits(&s->pb, ff_aac_scalefactor_bits[diff], ff_aac_scalefactor_code[diff]);
927 }
928 }
929 }
930 9458 }
931
932 /**
933 * Encode pulse data.
934 */
935 9458 static void encode_pulses(AACEncContext *s, Pulse *pulse)
936 {
937 int i;
938
939 9458 put_bits(&s->pb, 1, !!pulse->num_pulse);
940
1/2
✓ Branch 0 taken 9458 times.
✗ Branch 1 not taken.
9458 if (!pulse->num_pulse)
941 9458 return;
942
943 put_bits(&s->pb, 2, pulse->num_pulse - 1);
944 put_bits(&s->pb, 6, pulse->start);
945 for (i = 0; i < pulse->num_pulse; i++) {
946 put_bits(&s->pb, 5, pulse->pos[i]);
947 put_bits(&s->pb, 4, pulse->amp[i]);
948 }
949 }
950
951 /**
952 * Encode spectral coefficients processed by psychoacoustic model.
953 */
954 9458 static void encode_spectral_coeffs(AACEncContext *s, SingleChannelElement *sce)
955 {
956 int start, i, w, w2;
957
958
2/2
✓ Branch 0 taken 10146 times.
✓ Branch 1 taken 9458 times.
19604 for (w = 0; w < sce->ics.num_windows; w += sce->ics.group_len[w]) {
959 10146 start = 0;
960
2/2
✓ Branch 0 taken 456327 times.
✓ Branch 1 taken 10146 times.
466473 for (i = 0; i < sce->ics.max_sfb; i++) {
961
2/2
✓ Branch 0 taken 12117 times.
✓ Branch 1 taken 444210 times.
456327 if (sce->zeroes[w*16 + i]) {
962 12117 start += sce->ics.swb_sizes[i];
963 12117 continue;
964 }
965
2/2
✓ Branch 0 taken 456649 times.
✓ Branch 1 taken 444210 times.
900859 for (w2 = w; w2 < w + sce->ics.group_len[w]; w2++) {
966 456649 s->coder->quantize_and_encode_band(s, &s->pb,
967 456649 &sce->coeffs[start + w2*128],
968 456649 NULL, sce->ics.swb_sizes[i],
969 456649 sce->sf_idx[w*16 + i],
970 456649 sce->band_type[w*16 + i],
971 s->lambda,
972 456649 sce->ics.window_clipping[w]);
973 }
974 444210 start += sce->ics.swb_sizes[i];
975 }
976 }
977 9458 }
978
979 /**
980 * Downscale spectral coefficients for near-clipping windows to avoid artifacts
981 */
982 6532 static void avoid_clipping(AACEncContext *s, SingleChannelElement *sce)
983 {
984 int start, i, j, w;
985
986
2/2
✓ Branch 0 taken 157 times.
✓ Branch 1 taken 6375 times.
6532 if (sce->ics.clip_avoidance_factor < 1.0f) {
987
2/2
✓ Branch 0 taken 220 times.
✓ Branch 1 taken 157 times.
377 for (w = 0; w < sce->ics.num_windows; w++) {
988 220 start = 0;
989
2/2
✓ Branch 0 taken 8115 times.
✓ Branch 1 taken 220 times.
8335 for (i = 0; i < sce->ics.max_sfb; i++) {
990 8115 float *swb_coeffs = &sce->coeffs[start + w*128];
991
2/2
✓ Branch 0 taken 156312 times.
✓ Branch 1 taken 8115 times.
164427 for (j = 0; j < sce->ics.swb_sizes[i]; j++)
992 156312 swb_coeffs[j] *= sce->ics.clip_avoidance_factor;
993 8115 start += sce->ics.swb_sizes[i];
994 }
995 }
996 }
997 6532 }
998
999 /**
1000 * Encode one channel of audio data.
1001 */
1002 9458 static int encode_individual_channel(AVCodecContext *avctx, AACEncContext *s,
1003 SingleChannelElement *sce,
1004 int common_window)
1005 {
1006 9458 put_bits(&s->pb, 8, sce->sf_idx[0]);
1007
2/2
✓ Branch 0 taken 682 times.
✓ Branch 1 taken 8776 times.
9458 if (!common_window)
1008 682 put_ics_info(s, &sce->ics);
1009 9458 encode_band_info(s, sce);
1010 9458 encode_scale_factors(avctx, s, sce);
1011 9458 encode_pulses(s, &sce->pulse);
1012 9458 put_bits(&s->pb, 1, !!sce->tns.present);
1013
1/2
✓ Branch 0 taken 9458 times.
✗ Branch 1 not taken.
9458 if (s->coder->encode_tns_info)
1014 9458 s->coder->encode_tns_info(s, sce);
1015 9458 put_bits(&s->pb, 1, 0); //ssr
1016 9458 encode_spectral_coeffs(s, sce);
1017 9458 return 0;
1018 }
1019
1020 /**
1021 * Write some auxiliary information about the created AAC file.
1022 */
1023 static void put_bitstream_info(AACEncContext *s, const char *name)
1024 {
1025 int i, namelen, padbits;
1026
1027 namelen = strlen(name) + 2;
1028 put_bits(&s->pb, 3, TYPE_FIL);
1029 put_bits(&s->pb, 4, FFMIN(namelen, 15));
1030 if (namelen >= 15)
1031 put_bits(&s->pb, 8, namelen - 14);
1032 put_bits(&s->pb, 4, 0); //extension type - filler
1033 padbits = -put_bits_count(&s->pb) & 7;
1034 align_put_bits(&s->pb);
1035 for (i = 0; i < namelen - 2; i++)
1036 put_bits(&s->pb, 8, name[i]);
1037 put_bits(&s->pb, 12 - padbits, 0);
1038 }
1039
1040 /*
1041 * Copy input samples.
1042 * Channels are reordered from libavcodec's default order to AAC order.
1043 */
1044 3445 static void copy_input_samples(AACEncContext *s, const AVFrame *frame)
1045 {
1046 int ch;
1047
2/2
✓ Branch 0 taken 3413 times.
✓ Branch 1 taken 32 times.
3445 int end = 2048 + (frame ? frame->nb_samples : 0);
1048 3445 const uint8_t *channel_map = s->reorder_map;
1049
1050 /* copy and remap input samples */
1051
2/2
✓ Branch 0 taken 6567 times.
✓ Branch 1 taken 3445 times.
10012 for (ch = 0; ch < s->channels; ch++) {
1052 /* copy last 1024 samples of previous frame to the start of the current frame */
1053 6567 memcpy(&s->planar_samples[ch][1024], &s->planar_samples[ch][2048], 1024 * sizeof(s->planar_samples[0][0]));
1054
1055 /* copy new samples and zero any remaining samples */
1056
2/2
✓ Branch 0 taken 6497 times.
✓ Branch 1 taken 70 times.
6567 if (frame) {
1057 6497 memcpy(&s->planar_samples[ch][2048],
1058 6497 frame->extended_data[channel_map[ch]],
1059 6497 frame->nb_samples * sizeof(s->planar_samples[0][0]));
1060 }
1061 6567 memset(&s->planar_samples[ch][end], 0,
1062 6567 (3072 - end) * sizeof(s->planar_samples[0][0]));
1063 }
1064 3445 }
1065
1066 3461 static int aac_encode_frame(AVCodecContext *avctx, AVPacket *avpkt,
1067 const AVFrame *frame, int *got_packet_ptr)
1068 {
1069 3461 AACEncContext *s = avctx->priv_data;
1070 3461 float **samples = s->planar_samples, *samples2, *la, *overlap;
1071 ChannelElement *cpe;
1072 SingleChannelElement *sce;
1073 IndividualChannelStream *ics;
1074 int i, its, ch, w, chans, tag, start_ch, ret, frame_bits;
1075 int target_bits, rate_bits, too_many_bits, too_few_bits;
1076 3461 int ms_mode = 0, is_mode = 0, tns_mode = 0, pred_mode = 0;
1077 int chan_el_counter[4];
1078 FFPsyWindowInfo windows[AAC_MAX_CHANNELS];
1079
1080 /* add current frame to queue */
1081
2/2
✓ Branch 0 taken 3413 times.
✓ Branch 1 taken 48 times.
3461 if (frame) {
1082
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 3413 times.
3413 if ((ret = ff_af_queue_add(&s->afq, frame)) < 0)
1083 return ret;
1084 } else {
1085
3/6
✓ Branch 0 taken 32 times.
✓ Branch 1 taken 16 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 32 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
48 if (!s->afq.remaining_samples || (!s->afq.frame_alloc && !s->afq.frame_count))
1086 16 return 0;
1087 }
1088
1089 3445 copy_input_samples(s, frame);
1090
1091
2/2
✓ Branch 0 taken 16 times.
✓ Branch 1 taken 3429 times.
3445 if (!avctx->frame_num)
1092 16 return 0;
1093
1094 3429 start_ch = 0;
1095
2/2
✓ Branch 0 taken 3573 times.
✓ Branch 1 taken 3429 times.
7002 for (i = 0; i < s->chan_map[0]; i++) {
1096 3573 FFPsyWindowInfo* wi = windows + start_ch;
1097 3573 tag = s->chan_map[i+1];
1098
2/2
✓ Branch 0 taken 2959 times.
✓ Branch 1 taken 614 times.
3573 chans = tag == TYPE_CPE ? 2 : 1;
1099 3573 cpe = &s->cpe[i];
1100 {
1101 3573 int wi_paired = 0;
1102 /* Synced pair windows: decide both channels of a CPE together so
1103 * their block switching never diverges (see psy window_pair). */
1104
6/8
✓ Branch 0 taken 2959 times.
✓ Branch 1 taken 614 times.
✓ Branch 2 taken 2959 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 2959 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2927 times.
✓ Branch 7 taken 32 times.
3573 if (chans == 2 && tag != TYPE_LFE && s->psy.model->window_pair && frame) {
1105 2927 const float *ov0 = &samples[start_ch][0], *ov1 = &samples[start_ch + 1][0];
1106 2927 s->psy.model->window_pair(&s->psy,
1107 ov0 + 1024, ov0 + 1024 + 448 + 64,
1108 ov1 + 1024, ov1 + 1024 + 448 + 64,
1109 start_ch, start_ch + 1,
1110 2927 cpe->ch[0].ics.window_sequence[0],
1111 2927 cpe->ch[1].ics.window_sequence[0],
1112 wi);
1113 2927 wi_paired = 1;
1114 }
1115
2/2
✓ Branch 0 taken 6532 times.
✓ Branch 1 taken 3573 times.
10105 for (ch = 0; ch < chans; ch++) {
1116 int k;
1117 float clip_avoidance_factor;
1118 6532 sce = &cpe->ch[ch];
1119 6532 ics = &sce->ics;
1120 6532 s->cur_channel = start_ch + ch;
1121 6532 overlap = &samples[s->cur_channel][0];
1122 6532 samples2 = overlap + 1024;
1123 6532 la = samples2 + (448+64);
1124
2/2
✓ Branch 0 taken 70 times.
✓ Branch 1 taken 6462 times.
6532 if (!frame)
1125 70 la = NULL;
1126
2/2
✓ Branch 0 taken 48 times.
✓ Branch 1 taken 6484 times.
6532 if (tag == TYPE_LFE) {
1127 48 wi[ch].window_type[0] = wi[ch].window_type[1] = ONLY_LONG_SEQUENCE;
1128 48 wi[ch].window_shape = 0;
1129 48 wi[ch].num_windows = 1;
1130 48 wi[ch].grouping[0] = 1;
1131 48 wi[ch].clipping[0] = 0;
1132
1133 /* Only the lowest 12 coefficients are used in a LFE channel.
1134 * The expression below results in only the bottom 8 coefficients
1135 * being used for 11.025kHz to 16kHz sample rates.
1136 */
1137
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 48 times.
48 ics->num_swb = s->samplerate_index >= 8 ? 1 : 3;
1138
2/2
✓ Branch 0 taken 630 times.
✓ Branch 1 taken 5854 times.
6484 } else if (!wi_paired) {
1139 630 wi[ch] = s->psy.model->window(&s->psy, samples2, la, s->cur_channel,
1140 630 ics->window_sequence[0]);
1141 }
1142 6532 ics->window_sequence[1] = ics->window_sequence[0];
1143 6532 ics->window_sequence[0] = wi[ch].window_type[0];
1144 6532 ics->use_kb_window[1] = ics->use_kb_window[0];
1145 6532 ics->use_kb_window[0] = wi[ch].window_shape;
1146 6532 ics->num_windows = wi[ch].num_windows;
1147
2/2
✓ Branch 0 taken 229 times.
✓ Branch 1 taken 6303 times.
6532 ics->swb_sizes = s->psy.bands [ics->num_windows == 8];
1148
4/4
✓ Branch 0 taken 48 times.
✓ Branch 1 taken 6484 times.
✓ Branch 2 taken 229 times.
✓ Branch 3 taken 6255 times.
6532 ics->num_swb = tag == TYPE_LFE ? ics->num_swb : s->psy.num_bands[ics->num_windows == 8];
1149 6532 ics->max_sfb = FFMIN(ics->max_sfb, ics->num_swb);
1150 13064 ics->swb_offset = wi[ch].window_type[0] == EIGHT_SHORT_SEQUENCE ?
1151
2/2
✓ Branch 0 taken 229 times.
✓ Branch 1 taken 6303 times.
6532 ff_swb_offset_128 [s->samplerate_index]:
1152 6303 ff_swb_offset_1024[s->samplerate_index];
1153 13064 ics->tns_max_bands = wi[ch].window_type[0] == EIGHT_SHORT_SEQUENCE ?
1154
2/2
✓ Branch 0 taken 229 times.
✓ Branch 1 taken 6303 times.
6532 ff_tns_max_bands_128 [s->samplerate_index]:
1155 6303 ff_tns_max_bands_1024[s->samplerate_index];
1156
1157
2/2
✓ Branch 0 taken 8135 times.
✓ Branch 1 taken 6532 times.
14667 for (w = 0; w < ics->num_windows; w++)
1158 8135 ics->group_len[w] = wi[ch].grouping[w];
1159
1160 /* Calculate input sample maximums and evaluate clipping risk */
1161 6532 clip_avoidance_factor = 0.0f;
1162
2/2
✓ Branch 0 taken 8135 times.
✓ Branch 1 taken 6532 times.
14667 for (w = 0; w < ics->num_windows; w++) {
1163 8135 const float *wbuf = overlap + w * 128;
1164 8135 const int wlen = 2048 / ics->num_windows;
1165 8135 float max = 0;
1166 int j;
1167 /* mdct input is 2 * output */
1168
2/2
✓ Branch 0 taken 13377536 times.
✓ Branch 1 taken 8135 times.
13385671 for (j = 0; j < wlen; j++)
1169
2/2
✓ Branch 0 taken 13061957 times.
✓ Branch 1 taken 315579 times.
13377536 max = FFMAX(max, fabsf(wbuf[j]));
1170 8135 wi[ch].clipping[w] = max;
1171 }
1172
2/2
✓ Branch 0 taken 8135 times.
✓ Branch 1 taken 6532 times.
14667 for (w = 0; w < ics->num_windows; w++) {
1173
2/2
✓ Branch 0 taken 175 times.
✓ Branch 1 taken 7960 times.
8135 if (wi[ch].clipping[w] > CLIP_AVOIDANCE_FACTOR) {
1174 175 ics->window_clipping[w] = 1;
1175
2/2
✓ Branch 0 taken 5 times.
✓ Branch 1 taken 170 times.
175 clip_avoidance_factor = FFMAX(clip_avoidance_factor, wi[ch].clipping[w]);
1176 } else {
1177 7960 ics->window_clipping[w] = 0;
1178 }
1179 }
1180
2/2
✓ Branch 0 taken 157 times.
✓ Branch 1 taken 6375 times.
6532 if (clip_avoidance_factor > CLIP_AVOIDANCE_FACTOR) {
1181 157 ics->clip_avoidance_factor = CLIP_AVOIDANCE_FACTOR / clip_avoidance_factor;
1182 } else {
1183 6375 ics->clip_avoidance_factor = 1.0f;
1184 }
1185
1186 6532 apply_window_and_mdct(s, sce, overlap);
1187
1188
2/2
✓ Branch 0 taken 6688768 times.
✓ Branch 1 taken 6532 times.
6695300 for (k = 0; k < 1024; k++) {
1189
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 6688768 times.
6688768 if (!(fabs(cpe->ch[ch].coeffs[k]) < 1E16)) { // Ensure headroom for energy calculation
1190 av_log(avctx, AV_LOG_ERROR, "Input contains (near) NaN/+-Inf\n");
1191 return AVERROR(EINVAL);
1192 }
1193 }
1194 6532 avoid_clipping(s, sce);
1195 }
1196 }
1197 3573 start_ch += chans;
1198 }
1199
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 3429 times.
3429 if ((ret = ff_alloc_packet(avctx, avpkt, 8192 * s->channels)) < 0)
1200 return ret;
1201 3429 frame_bits = its = 0;
1202 do {
1203 4824 init_put_bits(&s->pb, avpkt->data, avpkt->size);
1204
1205
3/4
✓ Branch 0 taken 35 times.
✓ Branch 1 taken 4789 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 35 times.
4824 if ((avctx->frame_num & 0xFF)==1 && !(avctx->flags & AV_CODEC_FLAG_BITEXACT))
1206 put_bitstream_info(s, LIBAVCODEC_IDENT);
1207 4824 start_ch = 0;
1208 4824 target_bits = 0;
1209 4824 memset(chan_el_counter, 0, sizeof(chan_el_counter));
1210
2/2
✓ Branch 0 taken 5070 times.
✓ Branch 1 taken 4824 times.
9894 for (i = 0; i < s->chan_map[0]; i++) {
1211 5070 FFPsyWindowInfo* wi = windows + start_ch;
1212 const float *coeffs[2];
1213 5070 tag = s->chan_map[i+1];
1214
2/2
✓ Branch 0 taken 4388 times.
✓ Branch 1 taken 682 times.
5070 chans = tag == TYPE_CPE ? 2 : 1;
1215 5070 cpe = &s->cpe[i];
1216 5070 cpe->common_window = 0;
1217 5070 memset(cpe->is_mask, 0, sizeof(cpe->is_mask));
1218 5070 memset(cpe->ms_mask, 0, sizeof(cpe->ms_mask));
1219 5070 put_bits(&s->pb, 3, tag);
1220 5070 put_bits(&s->pb, 4, chan_el_counter[tag]++);
1221
2/2
✓ Branch 0 taken 9458 times.
✓ Branch 1 taken 5070 times.
14528 for (ch = 0; ch < chans; ch++) {
1222 9458 sce = &cpe->ch[ch];
1223 9458 coeffs[ch] = sce->coeffs;
1224 9458 memset(&sce->tns, 0, sizeof(TemporalNoiseShaping));
1225
2/2
✓ Branch 0 taken 1210624 times.
✓ Branch 1 taken 9458 times.
1220082 for (w = 0; w < 128; w++)
1226
2/2
✓ Branch 0 taken 15242 times.
✓ Branch 1 taken 1195382 times.
1210624 if (sce->band_type[w] > RESERVED_BT)
1227 15242 sce->band_type[w] = 0;
1228 }
1229 5070 s->psy.bitres.alloc = -1;
1230 5070 s->psy.bitres.bits = s->last_frame_pb_count / s->channels;
1231 5070 s->psy.model->analyze(&s->psy, start_ch, coeffs, wi);
1232
1/2
✓ Branch 0 taken 5070 times.
✗ Branch 1 not taken.
5070 if (s->psy.bitres.alloc > 0) {
1233 /* Lambda unused here on purpose, we need to take psy's unscaled allocation */
1234 10140 target_bits += s->psy.bitres.alloc
1235
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 5070 times.
5070 * (s->lambda / (avctx->global_quality ? avctx->global_quality : 120));
1236 5070 s->psy.bitres.alloc /= chans;
1237 }
1238 5070 s->cur_type = tag;
1239
2/2
✓ Branch 0 taken 4388 times.
✓ Branch 1 taken 682 times.
5070 if (chans > 1
1240
1/2
✓ Branch 0 taken 4388 times.
✗ Branch 1 not taken.
4388 && wi[0].window_type[0] == wi[1].window_type[0]
1241
1/2
✓ Branch 0 taken 4388 times.
✗ Branch 1 not taken.
4388 && wi[0].window_shape == wi[1].window_shape) {
1242
1243 4388 cpe->common_window = 1;
1244
2/2
✓ Branch 0 taken 5221 times.
✓ Branch 1 taken 4388 times.
9609 for (w = 0; w < wi[0].num_windows; w++) {
1245
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 5221 times.
5221 if (wi[0].grouping[w] != wi[1].grouping[w]) {
1246 cpe->common_window = 0;
1247 break;
1248 }
1249 }
1250 }
1251
1252
3/4
✓ Branch 0 taken 1394 times.
✓ Branch 1 taken 3676 times.
✓ Branch 2 taken 1394 times.
✗ Branch 3 not taken.
6464 const int use_tns = s->options.tns && s->coder->search_for_tns &&
1253
1/2
✓ Branch 0 taken 1394 times.
✗ Branch 1 not taken.
1394 s->coder->apply_tns_filt;
1254
1255 /* The NMR coder rate-controls itself and never re-quantizes, so TNS must run
1256 * before the quantizer */
1257 5070 const int tns_first = s->options.coder == AAC_CODER_NMR;
1258
3/4
✓ Branch 0 taken 655 times.
✓ Branch 1 taken 4415 times.
✓ Branch 2 taken 655 times.
✗ Branch 3 not taken.
5070 if (tns_first && use_tns) {
1259
2/2
✓ Branch 0 taken 792 times.
✓ Branch 1 taken 655 times.
1447 for (ch = 0; ch < chans; ch++) {
1260 792 sce = &cpe->ch[ch];
1261 792 s->cur_channel = start_ch + ch;
1262 /* mono: mark_pns before TNS so the region cap sees PNS bands. Stereo
1263 * PNS is marked in its own block (below) after the stereo decision. */
1264
4/6
✓ Branch 0 taken 518 times.
✓ Branch 1 taken 274 times.
✓ Branch 2 taken 518 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 518 times.
✗ Branch 5 not taken.
792 if (chans == 1 && s->options.pns && s->coder->mark_pns)
1265 518 s->coder->mark_pns(s, avctx, sce);
1266 792 s->coder->search_for_tns(s, sce);
1267 792 s->coder->apply_tns_filt(s, sce);
1268
2/2
✓ Branch 0 taken 24 times.
✓ Branch 1 taken 768 times.
792 if (sce->tns.present)
1269 24 tns_mode = 1;
1270 }
1271 }
1272
1273 /* NMR stereo PNS (imaging-safe). Mark each channel's noise-like bands on the
1274 * original L/R psy, then keep PNS only where BOTH channels are noise-like. */
1275
5/6
✓ Branch 0 taken 4388 times.
✓ Branch 1 taken 682 times.
✓ Branch 2 taken 4388 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 137 times.
✓ Branch 5 taken 4251 times.
5070 if (chans == 2 && cpe->common_window && tns_first &&
1276
2/4
✓ Branch 0 taken 137 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 137 times.
✗ Branch 3 not taken.
137 s->options.pns && s->coder->mark_pns) {
1277 137 s->cur_channel = start_ch; s->coder->mark_pns(s, avctx, &cpe->ch[0]);
1278 137 s->cur_channel = start_ch + 1; s->coder->mark_pns(s, avctx, &cpe->ch[1]);
1279
2/2
✓ Branch 0 taken 17536 times.
✓ Branch 1 taken 137 times.
17673 for (int b = 0; b < 128; b++)
1280
4/4
✓ Branch 0 taken 3029 times.
✓ Branch 1 taken 14507 times.
✓ Branch 2 taken 15 times.
✓ Branch 3 taken 3014 times.
17536 if (!cpe->ch[0].can_pns[b] || !cpe->ch[1].can_pns[b])
1281 14522 cpe->ch[0].can_pns[b] = cpe->ch[1].can_pns[b] = 0;
1282 }
1283
1284 /* The NMR coder decides I/S and M/S BEFORE quantization, from the psy model,
1285 * and the trellis then allocates natively on the coeffs actually coded. */
1286
5/6
✓ Branch 0 taken 4388 times.
✓ Branch 1 taken 682 times.
✓ Branch 2 taken 4388 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 137 times.
✓ Branch 5 taken 4251 times.
5070 if (chans == 2 && cpe->common_window && s->options.coder == AAC_CODER_NMR &&
1287
1/4
✗ Branch 0 not taken.
✓ Branch 1 taken 137 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
137 (s->options.mid_side || s->options.intensity_stereo)) {
1288 137 s->cur_channel = start_ch;
1289 137 nmr_decide_stereo(s, cpe);
1290 }
1291 /* NMR pools the CPE bit budget: both channels of a pair are solved
1292 * jointly under one shared lambda (see aaccoder_nmr.h). */
1293
3/4
✓ Branch 0 taken 655 times.
✓ Branch 1 taken 4415 times.
✓ Branch 2 taken 655 times.
✗ Branch 3 not taken.
5070 if (s->options.coder == AAC_CODER_NMR && s->nmr)
1294 655 s->nmr->pair = (chans == 2);
1295
2/2
✓ Branch 0 taken 9458 times.
✓ Branch 1 taken 5070 times.
14528 for (ch = 0; ch < chans; ch++) {
1296 9458 s->cur_channel = start_ch + ch;
1297 /* NMR PNS is mono-only */
1298
5/6
✓ Branch 0 taken 2106 times.
✓ Branch 1 taken 7352 times.
✓ Branch 2 taken 2106 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 1314 times.
✓ Branch 5 taken 792 times.
9458 if (s->options.pns && s->coder->mark_pns && !tns_first)
1299 1314 s->coder->mark_pns(s, avctx, &cpe->ch[ch]);
1300 9458 s->coder->search_for_quantizers(avctx, s, &cpe->ch[ch], s->lambda);
1301 }
1302
2/2
✓ Branch 0 taken 9458 times.
✓ Branch 1 taken 5070 times.
14528 for (ch = 0; ch < chans; ch++) { /* TNS (non-NMR) and PNS */
1303 9458 sce = &cpe->ch[ch];
1304 9458 s->cur_channel = start_ch + ch;
1305
4/4
✓ Branch 0 taken 8666 times.
✓ Branch 1 taken 792 times.
✓ Branch 2 taken 1314 times.
✓ Branch 3 taken 7352 times.
9458 if (!tns_first && use_tns) {
1306 1314 s->coder->search_for_tns(s, sce);
1307 1314 s->coder->apply_tns_filt(s, sce);
1308
2/2
✓ Branch 0 taken 54 times.
✓ Branch 1 taken 1260 times.
1314 if (sce->tns.present)
1309 54 tns_mode = 1;
1310 }
1311
4/4
✓ Branch 0 taken 2106 times.
✓ Branch 1 taken 7352 times.
✓ Branch 2 taken 1314 times.
✓ Branch 3 taken 792 times.
9458 if (s->options.pns && s->coder->search_for_pns)
1312 1314 s->coder->search_for_pns(s, avctx, sce);
1313 }
1314 5070 s->cur_channel = start_ch;
1315
2/2
✓ Branch 0 taken 1395 times.
✓ Branch 1 taken 3675 times.
5070 if (s->options.intensity_stereo) { /* Intensity Stereo */
1316
2/2
✓ Branch 0 taken 740 times.
✓ Branch 1 taken 655 times.
1395 if (s->options.coder != AAC_CODER_NMR) { /* NMR: decided pre-search */
1317
1/2
✓ Branch 0 taken 740 times.
✗ Branch 1 not taken.
740 if (s->coder->search_for_is)
1318 740 s->coder->search_for_is(s, avctx, cpe);
1319 740 apply_intensity_stereo(cpe);
1320 }
1321
2/2
✓ Branch 0 taken 467 times.
✓ Branch 1 taken 928 times.
1395 if (cpe->is_mode) is_mode = 1;
1322 }
1323
4/4
✓ Branch 0 taken 1066 times.
✓ Branch 1 taken 4004 times.
✓ Branch 2 taken 411 times.
✓ Branch 3 taken 655 times.
5070 if (s->options.mid_side && s->options.coder != AAC_CODER_NMR) { /* Mid/Side stereo */
1324
1/4
✗ Branch 0 not taken.
✓ Branch 1 taken 411 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
411 if (s->options.mid_side == -1 && s->coder->search_for_ms)
1325 s->coder->search_for_ms(s, cpe);
1326
1/2
✓ Branch 0 taken 411 times.
✗ Branch 1 not taken.
411 else if (cpe->common_window)
1327 411 memset(cpe->ms_mask, 1, sizeof(cpe->ms_mask));
1328 411 apply_mid_side_stereo(cpe);
1329 }
1330 5070 adjust_frame_information(cpe, chans);
1331
2/2
✓ Branch 0 taken 4388 times.
✓ Branch 1 taken 682 times.
5070 if (chans == 2) {
1332 4388 put_bits(&s->pb, 1, cpe->common_window);
1333
1/2
✓ Branch 0 taken 4388 times.
✗ Branch 1 not taken.
4388 if (cpe->common_window) {
1334 4388 put_ics_info(s, &cpe->ch[0].ics);
1335 4388 encode_ms_info(&s->pb, cpe);
1336
2/2
✓ Branch 0 taken 852 times.
✓ Branch 1 taken 3536 times.
4388 if (cpe->ms_mode) ms_mode = 1;
1337 }
1338 }
1339
2/2
✓ Branch 0 taken 9458 times.
✓ Branch 1 taken 5070 times.
14528 for (ch = 0; ch < chans; ch++) {
1340 9458 s->cur_channel = start_ch + ch;
1341 9458 encode_individual_channel(avctx, s, &cpe->ch[ch], cpe->common_window);
1342 }
1343 5070 start_ch += chans;
1344 }
1345
1346
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 4824 times.
4824 if (avctx->flags & AV_CODEC_FLAG_QSCALE) {
1347 /* When using a constant Q-scale, don't mess with lambda */
1348 break;
1349 }
1350
1351 4824 frame_bits = put_bits_count(&s->pb);
1352
1353 /* The NMR coder rate-controls itself (global-lambda reservoir servo):
1354 * per-frame bits intentionally float around the nominal rate, so skip
1355 * the lambda rate loop and only intervene on a hard overflow. */
1356
3/4
✓ Branch 0 taken 655 times.
✓ Branch 1 taken 4169 times.
✓ Branch 2 taken 655 times.
✗ Branch 3 not taken.
4824 if (s->options.coder == AAC_CODER_NMR && avctx->bit_rate_tolerance != 0 &&
1357
1/2
✓ Branch 0 taken 655 times.
✗ Branch 1 not taken.
655 frame_bits < 6144 * s->channels - 3)
1358 655 break;
1359
1360 /* rate control stuff
1361 * allow between the nominal bitrate, and what psy's bit reservoir says to target
1362 * but drift towards the nominal bitrate always
1363 */
1364 4169 rate_bits = avctx->bit_rate * 1024 / avctx->sample_rate;
1365
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 4169 times.
4169 rate_bits = FFMIN(rate_bits, 6144 * s->channels - 3);
1366 4169 too_many_bits = FFMAX(target_bits, rate_bits);
1367
2/2
✓ Branch 0 taken 2025 times.
✓ Branch 1 taken 2144 times.
4169 too_many_bits = FFMIN(too_many_bits, 6144 * s->channels - 3);
1368 4169 too_few_bits = FFMIN(FFMAX(rate_bits - rate_bits/4, target_bits), too_many_bits);
1369
1370 /* When strict bit-rate control is demanded */
1371
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 4169 times.
4169 if (avctx->bit_rate_tolerance == 0) {
1372 if (rate_bits < frame_bits) {
1373 float ratio = ((float)rate_bits) / frame_bits;
1374 s->lambda *= FFMIN(0.9f, ratio);
1375 continue;
1376 }
1377 /* reset lambda when solution is found */
1378 s->lambda = avctx->global_quality > 0 ? avctx->global_quality : 120;
1379 break;
1380 }
1381
1382 /* When using ABR, be strict (but only for increasing) */
1383 4169 too_few_bits = too_few_bits - too_few_bits/8;
1384 4169 too_many_bits = too_many_bits + too_many_bits/2;
1385
1386
2/2
✓ Branch 0 taken 1395 times.
✓ Branch 1 taken 2774 times.
4169 if ( its == 0 /* for steady-state Q-scale tracking */
1387
5/6
✓ Branch 0 taken 1117 times.
✓ Branch 1 taken 278 times.
✓ Branch 2 taken 1 times.
✓ Branch 3 taken 1116 times.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
1395 || (its < 5 && (frame_bits < too_few_bits || frame_bits > too_many_bits))
1388
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 279 times.
279 || frame_bits >= 6144 * s->channels - 3 )
1389 {
1390 3890 float ratio = ((float)rate_bits) / frame_bits;
1391
1392
3/4
✓ Branch 0 taken 2165 times.
✓ Branch 1 taken 1725 times.
✓ Branch 2 taken 2165 times.
✗ Branch 3 not taken.
3890 if (frame_bits >= too_few_bits && frame_bits <= too_many_bits) {
1393 /*
1394 * This path is for steady-state Q-scale tracking
1395 * When frame bits fall within the stable range, we still need to adjust
1396 * lambda to maintain it like so in a stable fashion (large jumps in lambda
1397 * create artifacts and should be avoided), but slowly
1398 */
1399 2165 ratio = sqrtf(sqrtf(ratio));
1400 2165 ratio = av_clipf(ratio, 0.9f, 1.1f);
1401 } else {
1402 /* Not so fast though */
1403 1725 ratio = sqrtf(ratio);
1404 }
1405 3890 s->lambda = av_clipf(s->lambda * ratio, FLT_EPSILON, 65536.f);
1406
1407 /* Keep iterating if we must reduce and lambda is in the sky */
1408
4/4
✓ Branch 0 taken 3862 times.
✓ Branch 1 taken 28 times.
✓ Branch 2 taken 1367 times.
✓ Branch 3 taken 2495 times.
3890 if (ratio > 0.9f && ratio < 1.1f) {
1409 break;
1410 } else {
1411
6/8
✓ Branch 0 taken 1394 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 1394 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 1382 times.
✓ Branch 5 taken 12 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 1382 times.
1395 if (is_mode || ms_mode || tns_mode || pred_mode) {
1412
2/2
✓ Branch 0 taken 49 times.
✓ Branch 1 taken 13 times.
62 for (i = 0; i < s->chan_map[0]; i++) {
1413 // Must restore coeffs
1414
2/2
✓ Branch 0 taken 1 times.
✓ Branch 1 taken 48 times.
49 chans = tag == TYPE_CPE ? 2 : 1;
1415 49 cpe = &s->cpe[i];
1416
2/2
✓ Branch 0 taken 50 times.
✓ Branch 1 taken 49 times.
99 for (ch = 0; ch < chans; ch++)
1417 50 memcpy(cpe->ch[ch].coeffs, cpe->ch[ch].pcoeffs, sizeof(cpe->ch[ch].coeffs));
1418 }
1419 }
1420 1395 its++;
1421 }
1422 } else {
1423 break;
1424 }
1425 } while (1);
1426
1427 /* tool-usage stats over the final per-band decisions of this frame */
1428
2/2
✓ Branch 0 taken 3573 times.
✓ Branch 1 taken 3429 times.
7002 for (i = 0; i < s->chan_map[0]; i++) {
1429
2/2
✓ Branch 0 taken 2959 times.
✓ Branch 1 taken 614 times.
3573 int etag = s->chan_map[i + 1], echans = etag == TYPE_CPE ? 2 : 1;
1430 3573 ChannelElement *ce = &s->cpe[i];
1431 3573 IndividualChannelStream *ics = &ce->ch[0].ics;
1432
2/2
✓ Branch 0 taken 6532 times.
✓ Branch 1 taken 3573 times.
10105 for (ch = 0; ch < echans; ch++) { /* per-channel frame stats */
1433 6532 int is_short = ce->ch[ch].ics.window_sequence[0] == EIGHT_SHORT_SEQUENCE;
1434 6532 s->stat_chans++;
1435
2/2
✓ Branch 0 taken 229 times.
✓ Branch 1 taken 6303 times.
6532 if (is_short)
1436 229 s->stat_short++;
1437
2/2
✓ Branch 0 taken 64 times.
✓ Branch 1 taken 6468 times.
6532 if (ce->ch[ch].tns.present) {
1438
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 64 times.
64 if (is_short) s->stat_tns_short++;
1439 64 else s->stat_tns_long++;
1440 }
1441 }
1442
2/2
✓ Branch 0 taken 3903 times.
✓ Branch 1 taken 3573 times.
7476 for (w = 0; w < ics->num_windows; w += ics->group_len[w]) {
1443
2/2
✓ Branch 0 taken 173429 times.
✓ Branch 1 taken 3903 times.
177332 for (int g = 0; g < ics->num_swb; g++) {
1444 173429 int idx = w*16 + g, coded = 0;
1445
2/2
✓ Branch 0 taken 318973 times.
✓ Branch 1 taken 173429 times.
492402 for (ch = 0; ch < echans; ch++) {
1446 318973 SingleChannelElement *sce = &ce->ch[ch];
1447
4/4
✓ Branch 0 taken 12735 times.
✓ Branch 1 taken 306238 times.
✓ Branch 2 taken 10853 times.
✓ Branch 3 taken 1882 times.
318973 if (sce->zeroes[idx] && sce->band_type[idx] == 0)
1448 10853 continue;
1449 308120 s->stat_ch_bands++;
1450
2/2
✓ Branch 0 taken 5263 times.
✓ Branch 1 taken 302857 times.
308120 if (sce->band_type[idx] == NOISE_BT)
1451 5263 s->stat_pns++;
1452 308120 coded = 1;
1453 }
1454
4/4
✓ Branch 0 taken 145544 times.
✓ Branch 1 taken 27885 times.
✓ Branch 2 taken 142966 times.
✓ Branch 3 taken 2578 times.
173429 if (etag == TYPE_CPE && coded) {
1455 142966 s->stat_cpe_bands++;
1456
2/2
✓ Branch 0 taken 26224 times.
✓ Branch 1 taken 116742 times.
142966 if (ce->ms_mask[idx]) s->stat_ms++;
1457
2/2
✓ Branch 0 taken 9424 times.
✓ Branch 1 taken 133542 times.
142966 if (ce->is_mask[idx]) s->stat_is++;
1458 }
1459 }
1460 }
1461 }
1462
1463 3429 put_bits(&s->pb, 3, TYPE_END);
1464 3429 flush_put_bits(&s->pb);
1465
1466 3429 s->last_frame_pb_count = put_bits_count(&s->pb);
1467
1468 /* NMR rate accounting: how many bits the frame really took beyond what the
1469 * trellis counted; feeds the next frame's budget correction */
1470
2/2
✓ Branch 0 taken 655 times.
✓ Branch 1 taken 2774 times.
3429 if (s->nmr) {
1471 655 int counted = 0;
1472
2/2
✓ Branch 0 taken 792 times.
✓ Branch 1 taken 655 times.
1447 for (i = 0; i < s->channels; i++)
1473 792 counted += s->nmr->counted[i];
1474
2/2
✓ Branch 0 taken 649 times.
✓ Branch 1 taken 6 times.
655 if (counted > 0) {
1475 649 float side = (float)s->last_frame_pb_count - counted;
1476
2/2
✓ Branch 0 taken 642 times.
✓ Branch 1 taken 7 times.
649 if (s->nmr->side_inited) {
1477 642 s->nmr->side_ema += 0.125f * (side - s->nmr->side_ema);
1478 } else {
1479 7 s->nmr->side_ema = side;
1480 7 s->nmr->side_inited = 1;
1481 }
1482 }
1483 }
1484 3429 avpkt->size = put_bytes_output(&s->pb);
1485
1486
4/4
✓ Branch 0 taken 655 times.
✓ Branch 1 taken 2774 times.
✓ Branch 2 taken 649 times.
✓ Branch 3 taken 6 times.
3429 s->lambda_sum += (s->nmr && s->nmr->lam_rc > 0.0f) ? s->nmr->lam_rc : s->lambda;
1487 3429 s->lambda_count++;
1488
1489 3429 ret = ff_af_queue_remove(&s->afq, avctx->frame_size, avpkt);
1490
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 3429 times.
3429 if (ret < 0)
1491 return ret;
1492
1493 3429 avpkt->flags |= AV_PKT_FLAG_KEY;
1494
1495 3429 *got_packet_ptr = 1;
1496 3429 return 0;
1497 }
1498
1499 16 static av_cold int aac_encode_end(AVCodecContext *avctx)
1500 {
1501 16 AACEncContext *s = avctx->priv_data;
1502
1503 27 av_log(avctx, AV_LOG_INFO,
1504 "Qavg: %.3f Tr: %.1f%% TNS(L): %.1f%% TNS(S): %.1f%% M/S: %.1f%% I/S: %.1f%% PNS: %.1f%%\n",
1505
1/2
✓ Branch 0 taken 16 times.
✗ Branch 1 not taken.
16 s->lambda_count ? s->lambda_sum / s->lambda_count : NAN,
1506
1/2
✓ Branch 0 taken 16 times.
✗ Branch 1 not taken.
16 s->stat_chans ? 100.0 * s->stat_short / s->stat_chans : 0.0,
1507
1/2
✓ Branch 0 taken 16 times.
✗ Branch 1 not taken.
16 s->stat_chans - s->stat_short ? 100.0 * s->stat_tns_long / (s->stat_chans - s->stat_short) : 0.0,
1508
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 6 times.
16 s->stat_short ? 100.0 * s->stat_tns_short / s->stat_short : 0.0,
1509
2/2
✓ Branch 0 taken 14 times.
✓ Branch 1 taken 2 times.
16 s->stat_cpe_bands ? 100.0 * s->stat_ms / s->stat_cpe_bands : 0.0,
1510
2/2
✓ Branch 0 taken 14 times.
✓ Branch 1 taken 2 times.
16 s->stat_cpe_bands ? 100.0 * s->stat_is / s->stat_cpe_bands : 0.0,
1511
2/2
✓ Branch 0 taken 15 times.
✓ Branch 1 taken 1 times.
16 s->stat_ch_bands ? 100.0 * s->stat_pns / s->stat_ch_bands : 0.0);
1512
1513 16 av_tx_uninit(&s->mdct1024);
1514 16 av_tx_uninit(&s->mdct128);
1515 16 ff_psy_end(&s->psy);
1516 16 ff_lpc_end(&s->lpc);
1517 16 av_freep(&s->buffer.samples);
1518 16 av_freep(&s->cpe);
1519 16 av_freep(&s->fdsp);
1520 16 av_freep(&s->nmr);
1521 16 ff_af_queue_close(&s->afq);
1522 16 return 0;
1523 }
1524
1525 16 static av_cold int dsp_init(AVCodecContext *avctx, AACEncContext *s)
1526 {
1527 16 int ret = 0;
1528 16 float scale = 32768.0f;
1529
1530 16 s->fdsp = avpriv_float_dsp_alloc(avctx->flags & AV_CODEC_FLAG_BITEXACT);
1531
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 16 times.
16 if (!s->fdsp)
1532 return AVERROR(ENOMEM);
1533
1534
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 16 times.
16 if ((ret = av_tx_init(&s->mdct1024, &s->mdct1024_fn, AV_TX_FLOAT_MDCT, 0,
1535 1024, &scale, 0)) < 0)
1536 return ret;
1537
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 16 times.
16 if ((ret = av_tx_init(&s->mdct128, &s->mdct128_fn, AV_TX_FLOAT_MDCT, 0,
1538 128, &scale, 0)) < 0)
1539 return ret;
1540
1541 16 return 0;
1542 }
1543
1544 16 static av_cold int alloc_buffers(AVCodecContext *avctx, AACEncContext *s)
1545 {
1546 int ch;
1547
1/2
✓ Branch 1 taken 16 times.
✗ Branch 2 not taken.
16 if (!FF_ALLOCZ_TYPED_ARRAY(s->buffer.samples, s->channels * 3 * 1024) ||
1548
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 16 times.
16 !FF_ALLOCZ_TYPED_ARRAY(s->cpe, s->chan_map[0]))
1549 return AVERROR(ENOMEM);
1550
1551
2/2
✓ Branch 0 taken 35 times.
✓ Branch 1 taken 16 times.
51 for(ch = 0; ch < s->channels; ch++)
1552 35 s->planar_samples[ch] = s->buffer.samples + 3 * 1024 * ch;
1553
1554
2/2
✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
16 if (s->options.coder == AAC_CODER_NMR) {
1555 8 s->nmr = av_mallocz(sizeof(*s->nmr));
1556
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 8 times.
8 if (!s->nmr)
1557 return AVERROR(ENOMEM);
1558 }
1559
1560 16 return 0;
1561 }
1562
1563 16 static av_cold int aac_encode_init(AVCodecContext *avctx)
1564 {
1565 16 AACEncContext *s = avctx->priv_data;
1566 16 int i, ret = 0;
1567 int chcfg;
1568 const uint8_t *sizes[2];
1569 uint8_t grouping[AAC_MAX_CHANNELS];
1570 int lengths[2];
1571
1572 /* Constants */
1573 16 s->last_frame_pb_count = 0;
1574 16 avctx->frame_size = 1024;
1575 16 avctx->initial_padding = 1024;
1576
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 16 times.
16 s->lambda = avctx->global_quality > 0 ? avctx->global_quality : 120;
1577
1578 /* Channel map and unspecified bitrate guessing */
1579 16 s->channels = avctx->ch_layout.nb_channels;
1580
1581 16 s->needs_pce = 1;
1582
1/2
✓ Branch 0 taken 35 times.
✗ Branch 1 not taken.
35 for (chcfg = 1; chcfg < FF_ARRAY_ELEMS(aac_normal_chan_layouts); chcfg++) {
1583
2/2
✓ Branch 1 taken 16 times.
✓ Branch 2 taken 19 times.
35 if (!av_channel_layout_compare(&avctx->ch_layout, &aac_normal_chan_layouts[chcfg])) {
1584 16 s->needs_pce = s->options.pce;
1585 16 break;
1586 }
1587 }
1588
1589
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 16 times.
16 if (s->needs_pce) {
1590 char buf[64];
1591 for (i = 0; i < FF_ARRAY_ELEMS(aac_pce_configs); i++)
1592 if (!av_channel_layout_compare(&avctx->ch_layout, &aac_pce_configs[i].layout))
1593 break;
1594 av_channel_layout_describe(&avctx->ch_layout, buf, sizeof(buf));
1595 if (i == FF_ARRAY_ELEMS(aac_pce_configs)) {
1596 av_log(avctx, AV_LOG_ERROR, "Unsupported channel layout \"%s\"\n", buf);
1597 return AVERROR(EINVAL);
1598 }
1599 av_log(avctx, AV_LOG_INFO, "Using a PCE to encode channel layout \"%s\"\n", buf);
1600 s->pce = aac_pce_configs[i];
1601 s->reorder_map = s->pce.reorder_map;
1602 s->chan_map = s->pce.config_map;
1603 chcfg = 0;
1604 } else {
1605 16 s->reorder_map = aac_chan_maps[chcfg - 1];
1606 16 s->chan_map = aac_chan_configs[chcfg - 1];
1607 }
1608
1609
2/2
✓ Branch 0 taken 9 times.
✓ Branch 1 taken 7 times.
16 if (!avctx->bit_rate) {
1610
2/2
✓ Branch 0 taken 12 times.
✓ Branch 1 taken 9 times.
21 for (i = 1; i <= s->chan_map[0]; i++) {
1611
2/2
✓ Branch 0 taken 3 times.
✓ Branch 1 taken 9 times.
15 avctx->bit_rate += s->chan_map[i] == TYPE_CPE ? 128000 : /* Pair */
1612
2/2
✓ Branch 0 taken 1 times.
✓ Branch 1 taken 2 times.
3 s->chan_map[i] == TYPE_LFE ? 16000 : /* LFE */
1613 69000 ; /* SCE */
1614 }
1615 }
1616
1617 /* Samplerate */
1618 16 for (int i = 0;; i++) {
1619 63 av_assert1(i < 13);
1620
2/2
✓ Branch 0 taken 16 times.
✓ Branch 1 taken 63 times.
79 if (avctx->sample_rate == ff_mpeg4audio_sample_rates[i]) {
1621 16 s->samplerate_index = i;
1622 16 break;
1623 }
1624 }
1625
1626 /* Bitrate limiting */
1627
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 16 times.
16 WARN_IF(1024.0 * avctx->bit_rate / avctx->sample_rate > 6144 * s->channels,
1628 "Too many bits %f > %d per frame requested, clamping to max\n",
1629 1024.0 * avctx->bit_rate / avctx->sample_rate,
1630 6144 * s->channels);
1631
1/2
✓ Branch 0 taken 16 times.
✗ Branch 1 not taken.
16 avctx->bit_rate = (int64_t)FFMIN(6144 * s->channels / 1024.0 * avctx->sample_rate,
1632 avctx->bit_rate);
1633
1634 /* Profile and option setting */
1635
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 16 times.
16 avctx->profile = avctx->profile == AV_PROFILE_UNKNOWN ? AV_PROFILE_AAC_LOW :
1636 avctx->profile;
1637
1/2
✓ Branch 0 taken 16 times.
✗ Branch 1 not taken.
16 for (i = 0; i < FF_ARRAY_ELEMS(aacenc_profiles); i++)
1638
1/2
✓ Branch 0 taken 16 times.
✗ Branch 1 not taken.
16 if (avctx->profile == aacenc_profiles[i])
1639 16 break;
1640
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 16 times.
16 ERROR_IF(i == FF_ARRAY_ELEMS(aacenc_profiles), "Profile not supported!\n");
1641
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 16 times.
16 if (avctx->profile == AV_PROFILE_MPEG2_AAC_LOW) {
1642 avctx->profile = AV_PROFILE_AAC_LOW;
1643 WARN_IF(s->options.pns,
1644 "PNS unavailable in the \"mpeg2_aac_low\" profile, turning off\n");
1645 s->options.pns = 0;
1646 }
1647 16 s->profile = avctx->profile;
1648
1649 /* Coder limitations */
1650 16 s->coder = &ff_aac_coders[s->options.coder];
1651
1652 /* M/S introduces horrible artifacts with multichannel files, this is temporary */
1653
2/2
✓ Branch 0 taken 1 times.
✓ Branch 1 taken 15 times.
16 if (s->channels > 3)
1654 1 s->options.mid_side = 0;
1655
1656 /* Coding bandwidth, fixed at init time */
1657
2/2
✓ Branch 0 taken 5 times.
✓ Branch 1 taken 11 times.
16 if (avctx->cutoff > 0) {
1658 5 s->bandwidth = avctx->cutoff;
1659 } else {
1660
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 11 times.
11 int frame_br = (avctx->flags & AV_CODEC_FLAG_QSCALE) ?
1661 (avctx->bit_rate / 2.0f * (s->lambda / 120.f) * 1.5f) :
1662 11 (avctx->bit_rate / avctx->ch_layout.nb_channels);
1663
1664
3/4
✓ Branch 0 taken 8 times.
✓ Branch 1 taken 3 times.
✓ Branch 2 taken 8 times.
✗ Branch 3 not taken.
11 if (s->options.coder == AAC_CODER_NMR && frame_br >= 24000) {
1665 static const int rates[] = { 24000, 32000, 48000, 64000, 96000, 192000 };
1666 static const int bws[] = { 14000, 14000, 18500, 20000, 21000, 22000 };
1667 8 int bw_i = 0;
1668
3/4
✓ Branch 0 taken 25 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 17 times.
✓ Branch 3 taken 8 times.
25 for (; bw_i < FF_ARRAY_ELEMS(rates) - 2 && frame_br > rates[bw_i + 1]; bw_i++);
1669 8 s->bandwidth = bws[bw_i] + (int)((int64_t)(bws[bw_i + 1] - bws[bw_i]) *
1670 8 (frame_br - rates[bw_i]) / (rates[bw_i + 1] - rates[bw_i]));
1671 8 s->bandwidth = FFMIN3(s->bandwidth, 22000, avctx->sample_rate / 2);
1672 } else {
1673
3/4
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
3 if (s->options.pns || s->options.intensity_stereo)
1674 1 frame_br *= 1.15f;
1675
5/10
✓ Branch 0 taken 3 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 3 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 3 times.
✓ Branch 6 taken 3 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 3 times.
3 s->bandwidth = FFMAX(3000, AAC_CUTOFF_FROM_BITRATE(frame_br, 1,
1676 avctx->sample_rate));
1677 }
1678
1679 11 s->bandwidth = FFMIN(FFMAX(s->bandwidth, 8000), avctx->sample_rate / 2);
1680 }
1681
1682
2/4
✓ Branch 0 taken 16 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 16 times.
✗ Branch 3 not taken.
16 if (!(avctx->flags & AV_CODEC_FLAG_QSCALE) && avctx->bit_rate > 0) {
1683 16 int bpc = avctx->bit_rate / avctx->ch_layout.nb_channels;
1684
1/4
✗ Branch 0 not taken.
✓ Branch 1 taken 16 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
16 if (bpc <= 32000 && avctx->sample_rate > 32000)
1685 av_log(avctx, AV_LOG_INFO,
1686 "%d kb/s per channel at %d Hz: consider resampling the "
1687 "input to 32000 Hz or lower for better quality.\n",
1688 bpc / 1000, avctx->sample_rate);
1689 }
1690
1691 // Initialize static tables
1692 16 ff_aac_float_common_init();
1693
1694
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 16 times.
16 if ((ret = dsp_init(avctx, s)) < 0)
1695 return ret;
1696
1697
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 16 times.
16 if ((ret = alloc_buffers(avctx, s)) < 0)
1698 return ret;
1699
1700
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 16 times.
16 if ((ret = put_audio_specific_config(avctx, chcfg)))
1701 return ret;
1702
1703 16 sizes[0] = ff_aac_swb_size_1024[s->samplerate_index];
1704 16 sizes[1] = ff_aac_swb_size_128[s->samplerate_index];
1705 16 lengths[0] = ff_aac_num_swb_1024[s->samplerate_index];
1706 16 lengths[1] = ff_aac_num_swb_128[s->samplerate_index];
1707
2/2
✓ Branch 0 taken 19 times.
✓ Branch 1 taken 16 times.
35 for (i = 0; i < s->chan_map[0]; i++)
1708 19 grouping[i] = s->chan_map[i + 1] == TYPE_CPE;
1709
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 16 times.
16 if ((ret = ff_psy_init(&s->psy, avctx, 2, sizes, lengths,
1710 16 s->chan_map[0], grouping, s->bandwidth)) < 0)
1711 return ret;
1712 16 ff_lpc_init(&s->lpc, 2*avctx->frame_size, TNS_MAX_ORDER, FF_LPC_TYPE_LEVINSON);
1713 16 s->random_state = 0x1f2e3d4c;
1714
1715 16 ff_aacenc_dsp_init(&s->aacdsp);
1716
1717 16 ff_af_queue_init(avctx, &s->afq);
1718
1719 16 return 0;
1720 }
1721
1722 #define AACENC_FLAGS AV_OPT_FLAG_ENCODING_PARAM | AV_OPT_FLAG_AUDIO_PARAM
1723 static const AVOption aacenc_options[] = {
1724 {"aac_coder", "Coding algorithm", offsetof(AACEncContext, options.coder), AV_OPT_TYPE_INT, {.i64 = AAC_CODER_NMR}, 0, AAC_CODER_NB-1, AACENC_FLAGS, .unit = "coder"},
1725 {"twoloop", "Two loop searching method", 0, AV_OPT_TYPE_CONST, {.i64 = AAC_CODER_TWOLOOP}, INT_MIN, INT_MAX, AACENC_FLAGS, .unit = "coder"},
1726 {"fast", "Fast search", 0, AV_OPT_TYPE_CONST, {.i64 = AAC_CODER_FAST}, INT_MIN, INT_MAX, AACENC_FLAGS, .unit = "coder"},
1727 {"nmr", "Noise-to-mask ratio scalefactor trellis", 0, AV_OPT_TYPE_CONST, {.i64 = AAC_CODER_NMR}, INT_MIN, INT_MAX, AACENC_FLAGS, .unit = "coder"},
1728 {"aac_ms", "Force M/S stereo coding", offsetof(AACEncContext, options.mid_side), AV_OPT_TYPE_BOOL, {.i64 = -1}, -1, 1, AACENC_FLAGS},
1729 {"aac_is", "Intensity stereo coding", offsetof(AACEncContext, options.intensity_stereo), AV_OPT_TYPE_BOOL, {.i64 = 1}, -1, 1, AACENC_FLAGS},
1730 {"aac_pns", "Perceptual noise substitution", offsetof(AACEncContext, options.pns), AV_OPT_TYPE_BOOL, {.i64 = 1}, -1, 1, AACENC_FLAGS},
1731 {"aac_tns", "Temporal noise shaping", offsetof(AACEncContext, options.tns), AV_OPT_TYPE_BOOL, {.i64 = 1}, -1, 1, AACENC_FLAGS},
1732 {"aac_pce", "Forces the use of PCEs", offsetof(AACEncContext, options.pce), AV_OPT_TYPE_BOOL, {.i64 = 0}, -1, 1, AACENC_FLAGS},
1733 {"aac_nmr_speed", "NMR coder speed level: 0 = slowest/best, higher trades quality for speed", offsetof(AACEncContext, options.nmr_speed), AV_OPT_TYPE_INT, {.i64 = 0}, 0, 4, AACENC_FLAGS},
1734 FF_AAC_PROFILE_OPTS
1735 {NULL}
1736 };
1737
1738 static const AVClass aacenc_class = {
1739 .class_name = "AAC encoder",
1740 .item_name = av_default_item_name,
1741 .option = aacenc_options,
1742 .version = LIBAVUTIL_VERSION_INT,
1743 };
1744
1745 static const FFCodecDefault aac_encode_defaults[] = {
1746 { "b", "0" },
1747 { NULL }
1748 };
1749
1750 const FFCodec ff_aac_encoder = {
1751 .p.name = "aac",
1752 CODEC_LONG_NAME("AAC (Advanced Audio Coding)"),
1753 .p.type = AVMEDIA_TYPE_AUDIO,
1754 .p.id = AV_CODEC_ID_AAC,
1755 .p.capabilities = AV_CODEC_CAP_DR1 | AV_CODEC_CAP_DELAY |
1756 AV_CODEC_CAP_SMALL_LAST_FRAME,
1757 .priv_data_size = sizeof(AACEncContext),
1758 .init = aac_encode_init,
1759 FF_CODEC_ENCODE_CB(aac_encode_frame),
1760 .close = aac_encode_end,
1761 .defaults = aac_encode_defaults,
1762 CODEC_SAMPLERATES_ARRAY(ff_mpeg4audio_sample_rates),
1763 .caps_internal = FF_CODEC_CAP_INIT_CLEANUP,
1764 CODEC_SAMPLEFMTS(AV_SAMPLE_FMT_FLTP),
1765 .p.priv_class = &aacenc_class,
1766 };
1767