projects / nihav.git / blob
? search:
summary | shortlog | log | commit | commitdiff | tree
history | raw | HEAD
add VivoActive (imperfect) format family support
[nihav.git] / nihav-vivo / src / codecs / g723_1.rs
1 use nihav_core::formats::SND_S16_FORMAT;
2 use nihav_core::frame::*;
3 use nihav_core::codecs::*;
4 use nihav_core::io::bitreader::*;
5 use std::str::FromStr;
6
7 const SAMPLES: usize = 240;
8 const SUBFRAMES: usize = 4;
9 const SUBFRAME_LEN: usize = SAMPLES / SUBFRAMES;
10 const LPC_ORDER: usize = 10;
11 const MIN_PITCH_LAG: usize = 18;
12 const MAX_PITCH: usize = MIN_PITCH_LAG + 127;
13 const MAX_ERASED_FRAMES:u8 = 3;
14 const MAX_PULSES: usize = 6;
15
16 #[derive(Clone,Copy,PartialEq)]
17 enum G7231FrameType {
18 Active,
19 SID,
20 Untransmitted,
21 }
22
23 #[derive(Clone,Copy,Default)]
24 struct Subframe {
25 pitch_lag: usize,
26 acb_lag: usize,
27 acb_gain: usize,
28 dirac_train: bool,
29 pulse_sign: u8,
30 grid_index: usize,
31 amp_index: usize,
32 pulse_pos: usize,
33
34 ppf_index: usize,
35 ppf_opt_gain: i32,
36 ppf_sc_gain: i32,
37 }
38
39 impl Subframe {
40 fn read_acb_params(&mut self, br: &mut BitReader, is_6300: bool) -> DecoderResult<()> {
41 let mut val = br.read(12)? as usize;
42 let acb_len;
43 if is_6300 && (self.pitch_lag < SUBFRAME_LEN - 2) {
44 self.dirac_train = (val >> 11) != 0;
45 val &= 0x7FF;
46 acb_len = 85;
47 } else {
48 self.dirac_train = false;
49 acb_len = 170;
50 }
51 self.acb_gain = val / 24;
52 self.amp_index = val % 24;
53 validate!(self.acb_gain < acb_len);
54
55 Ok(())
56 }
57 fn gen_fcb_excitation(&self, dst: &mut [i16], is_6300: bool, max_pos: usize, npulses: usize) {
58 for el in dst.iter_mut().take(SUBFRAME_LEN) {
59 *el = 0;
60 }
61 let fcb_gain = FIXED_CB_GAIN[self.amp_index];
62 if is_6300 {
63 if self.pulse_pos >= max_pos {
64 return;
65 }
66 let mut pulse = MAX_PULSES - npulses;
67 let mut ppos = self.pulse_pos;
68 for i in 0..SUBFRAME_LEN/2 {
69 let ret = ppos.checked_sub(PULSE_POSITIONS[pulse][i] as usize);
70 if let Some(new_ppos) = ret {
71 ppos = new_ppos;
72 continue;
73 }
74 pulse += 1;
75 let sign = (self.pulse_sign & (1 << (MAX_PULSES - pulse))) != 0;
76 dst[i * 2 + self.grid_index] = if !sign { fcb_gain } else { -fcb_gain };
77 if pulse == MAX_PULSES {
78 break;
79 }
80 }
81 if self.dirac_train {
82 let mut orig = [0; SUBFRAME_LEN];
83 orig.copy_from_slice(&dst[..SUBFRAME_LEN]);
84 for i in (self.pitch_lag..SUBFRAME_LEN).step_by(self.pitch_lag) {
85 for j in 0..SUBFRAME_LEN - i {
86 dst[i + j] += orig[j];
87 }
88 }
89 }
90 } else {
91 let mut cb_pos = self.pulse_pos;
92 let mut cb_sign = self.pulse_sign;
93 for i in (0..8).step_by(2) {
94 let off = (cb_pos & 7) * 8 + i + self.grid_index;
95 dst[off] = if (cb_sign & 1) != 0 { fcb_gain } else { -fcb_gain };
96 cb_pos >>= 3;
97 cb_sign >>= 1;
98 }
99
100 let lag = (PITCH_CONTRIBUTION[self.acb_gain * 2] as isize + self.pitch_lag as isize + self.acb_lag as isize - 1) as usize;
101 let beta = PITCH_CONTRIBUTION[self.acb_gain * 2 + 1];
102 if lag < SUBFRAME_LEN - 2 {
103 for i in lag..SUBFRAME_LEN {
104 dst[i] += ((beta * i32::from(dst[i - lag])) >> 15) as i16;
105 }
106 }
107 }
108 }
109 fn gen_acb_excitation(&mut self, acb_vector: &mut [i16; SUBFRAME_LEN], excitation: &[i16], is_6300: bool) {
110 let lag = self.pitch_lag + self.acb_lag - 1;
111
112 let mut residual = [0; SUBFRAME_LEN + 5 - 1];
113 residual[0] = excitation[MAX_PITCH - 2 - lag];
114 residual[1] = excitation[MAX_PITCH - 1 - lag];
115 for (i, dst) in residual.iter_mut().skip(2).enumerate() {
116 *dst = excitation[MAX_PITCH - lag + i % lag];
117 }
118
119 let codebook = if is_6300 && self.pitch_lag < SUBFRAME_LEN - 2 {
120 &ACB_GAIN85[self.acb_gain]
121 } else {
122 &ACB_GAIN170[self.acb_gain]
123 };
124 for i in 0..SUBFRAME_LEN {
125 let sum = dot_product(&residual[i..], codebook, 5);
126 acb_vector[i] = (sum.saturating_add(sum).saturating_add(1 << 15) >> 16) as i16;
127 }
128 }
129 fn compute_ppf_coeffs(&mut self, src: &[i16], offset: usize, is_6300: bool) {
130 self.ppf_index = offset;
131 self.ppf_opt_gain = 0;
132 self.ppf_sc_gain = 0x7FFF;
133
134 let (fwd_lag, fwd_energy) = autocorr_max(src, offset, SUBFRAME_LEN, self.pitch_lag, true);
135 let (bwd_lag, bwd_energy) = autocorr_max(src, offset, SUBFRAME_LEN, self.pitch_lag, false);
136 if fwd_lag == 0 && bwd_lag == 0 {
137 return;
138 }
139
140 let tgt_energy = dot_product(&src[offset..], &src[offset..], SUBFRAME_LEN);
141 let fwd_res_energy = if fwd_lag != 0 {
142 dot_product(&src[offset + fwd_lag..], &src[offset + fwd_lag..], SUBFRAME_LEN)
143 } else { 0 };
144 let bwd_res_energy = if bwd_lag != 0 {
145 dot_product(&src[offset - bwd_lag..], &src[offset - bwd_lag..], SUBFRAME_LEN)
146 } else { 0 };
147
148 let max_energy = tgt_energy.max(fwd_energy).max(fwd_res_energy).max(bwd_energy).max(bwd_res_energy);
149 let scale = norm_bits(max_energy, 31);
150 let tgt_energy = tgt_energy << scale >> 16;
151 let fwd_energy = fwd_energy << scale >> 16;
152 let fwd_res_energy = fwd_res_energy << scale >> 16;
153 let bwd_energy = bwd_energy << scale >> 16;
154 let bwd_res_energy = bwd_res_energy << scale >> 16;
155
156 let use_fwd = if fwd_lag != 0 && bwd_lag == 0 {
157 true
158 } else if fwd_lag == 0 {
159 false
160 } else {
161 let tmp1 = bwd_res_energy * ((fwd_energy * fwd_energy + (1 << 14)) >> 15);
162 let tmp2 = fwd_res_energy * ((bwd_energy * bwd_energy + (1 << 14)) >> 15);
163 tmp1 >= tmp2
164 };
165 if use_fwd {
166 self.compute_ppf_gains(offset + fwd_lag, is_6300, tgt_energy, fwd_energy, fwd_res_energy);
167 } else {
168 self.compute_ppf_gains(offset - bwd_lag, is_6300, tgt_energy, bwd_energy, bwd_res_energy);
169 }
170 }
171 fn compute_ppf_gains(&mut self, offset: usize, is_6300: bool, tgt_energy: i32, corr_energy: i32, res_energy: i32) {
172 self.ppf_index = offset;
173
174 let tmp1 = tgt_energy * res_energy >> 1;
175 let tmp2 = corr_energy * corr_energy * 2;
176 if tmp1 >= tmp2 {
177 return;
178 }
179 let gain_weight = if is_6300 { 0x1800 } else { 0x2000 };
180 self.ppf_opt_gain = if corr_energy >= res_energy {
181 gain_weight
182 } else {
183 (corr_energy << 15) / res_energy * gain_weight >> 15
184 };
185
186 let tmp1 = (tgt_energy << 15) + (corr_energy * self.ppf_opt_gain * 2);
187 let tmp2 = (self.ppf_opt_gain * self.ppf_opt_gain >> 15) * res_energy;
188
189 let residual = tmp1.saturating_add(tmp2).saturating_add(1 << 15) >> 16;
190
191 if tgt_energy >= residual * 2 {
192 self.ppf_sc_gain = square_root_i32(0x7FFF0000);
193 } else {
194 let val = (tgt_energy << 14) / residual;
195 self.ppf_sc_gain = square_root_i32(val << 16);
196 }
197
198 self.ppf_opt_gain = i32::from(clip16(self.ppf_opt_gain * self.ppf_sc_gain >> 15));
199 }
200 }
201
202 fn square_root_i32(val: i32) -> i32 {
203 let mut res = 0;
204 let mut exp = 1 << 14;
205
206 for _ in 0..14 {
207 let res_exp = res + exp;
208 if val >= res_exp * res_exp * 2 {
209 res += exp;
210 }
211 exp >>= 1;
212 }
213
214 res
215 }
216
217 struct PRNG {
218 seed: i32,
219 }
220
221 const CNG_RND_SEED: i32 = 12345;
222 impl PRNG {
223 fn new() -> Self { Self { seed: 0 } }
224 fn new_cng() -> Self { Self { seed: CNG_RND_SEED } }
225 fn reset(&mut self) { self.seed = CNG_RND_SEED; }
226 fn next(&mut self) -> i32 {
227 self.seed = self.seed.wrapping_mul(521).wrapping_add(259);
228 self.seed
229 }
230 fn next_range(&mut self, range: usize) -> usize {
231 let val = (self.next() & 0x7FFF) as usize;
232 val * range >> 15
233 }
234 }
235
236 struct G7231Decoder {
237 chmap: NAChannelMap,
238 ainfo: NAAudioInfo,
239 info: NACodecInfoRef,
240
241 is_6300: bool,
242 prev_lsp: [i16; LPC_ORDER],
243 sid_lsp: [i16; LPC_ORDER],
244 prev_ftype: G7231FrameType,
245 cur_ftype: G7231FrameType,
246 lsp_index: [usize; 3],
247 subframe: [Subframe; SUBFRAMES],
248 excitation: [i16; MAX_PITCH + SAMPLES + 4],
249 prev_excitation: [i16; MAX_PITCH],
250 lpc: [[i16; LPC_ORDER]; SUBFRAMES],
251
252 filt_mem: [i16; LPC_ORDER],
253
254 interp_index: usize,
255 interp_gain: i16,
256 sid_gain: i32,
257 cur_gain: i32,
258
259 fir_mem: [i16; LPC_ORDER],
260 iir_mem: [i32; LPC_ORDER],
261 reflection_coef: i32,
262 pf_gain: i32,
263
264 cng_rnd: PRNG,
265 rnd: PRNG,
266 erased_frames: u8,
267
268 synth_buf: [i16; SAMPLES + LPC_ORDER + MAX_PITCH + 4],
269 }
270
271 macro_rules! weighted_sum {
272 ($dst: expr, $src1: expr, $weight1: expr, $src2: expr, $weight2: expr, $shift: expr) => {
273 let bias = 1 << ($shift - 1);
274 for (dst, (src1, src2)) in $dst.iter_mut().zip($src1.iter().zip($src2.iter())) {
275 let val = (i32::from(*src1) * $weight1 + i32::from(*src2) * $weight2 + bias) >> $shift;
276 *dst = clip16(val);
277 }
278 }
279 }
280
281 impl G7231Decoder {
282 fn new() -> Self {
283 let prev_lsp = DC_LSP;
284 let sid_lsp = DC_LSP;
285 G7231Decoder {
286 chmap: NAChannelMap::from_str("C").unwrap(),
287 ainfo: NAAudioInfo::new(8000, 1, SND_S16_FORMAT, SAMPLES),
288 info: NACodecInfo::new_dummy(),
289
290 is_6300: true,
291 prev_lsp, sid_lsp,
292 prev_ftype: G7231FrameType::SID,
293 cur_ftype: G7231FrameType::Untransmitted,
294 lsp_index: [0; 3],
295 subframe: [Subframe::default(); SUBFRAMES],
296 excitation: [0; MAX_PITCH + SAMPLES + 4],
297 prev_excitation: [0; MAX_PITCH],
298 lpc: [[0; LPC_ORDER]; SUBFRAMES],
299
300 filt_mem: [0; LPC_ORDER],
301
302 interp_index: 0,
303 interp_gain: 0,
304 sid_gain: 0,
305 cur_gain: 0,
306
307 fir_mem: [0; LPC_ORDER],
308 iir_mem: [0; LPC_ORDER],
309 reflection_coef: 0,
310 pf_gain: 1 << 12,
311
312 cng_rnd: PRNG::new_cng(),
313 rnd: PRNG::new(),
314 erased_frames: 0,
315
316 synth_buf: [0; SAMPLES + LPC_ORDER + MAX_PITCH + 4],
317 }
318 }
319 fn unpack_frame(&mut self, br: &mut BitReader) -> DecoderResult<()> {
320 const FRAME_SIZE: [isize; 4] = [ 24, 20, 4, 1 ];
321 let ftype = br.read(2)? as usize;
322 validate!(br.left() + 2 >= FRAME_SIZE[ftype] * 8);
323 if ftype == 3 {
324 self.cur_ftype = G7231FrameType::Untransmitted;
325 return Ok(());
326 }
327
328 for lsp_index in self.lsp_index.iter_mut().rev() {
329 *lsp_index = br.read(8)? as usize;
330 }
331
332 if ftype == 2 {
333 self.cur_ftype = G7231FrameType::Untransmitted;
334 self.subframe[0].amp_index = br.read(8)? as usize;
335 return Ok(());
336 }
337
338 self.cur_ftype = G7231FrameType::Active;
339 self.is_6300 = ftype == 0;
340
341 for i in (0..4).step_by(2) {
342 self.subframe[i].pitch_lag = (br.read(7)? as usize) + MIN_PITCH_LAG;
343 validate!(self.subframe[i].pitch_lag < MIN_PITCH_LAG + 124);
344 self.subframe[i + 1].pitch_lag = self.subframe[i].pitch_lag;
345 self.subframe[i].acb_lag = 1;
346 self.subframe[i + 1].acb_lag = br.read(2)? as usize;
347 }
348
349 for subframe in self.subframe.iter_mut() {
350 subframe.read_acb_params(br, self.is_6300)?;
351 }
352 for subframe in self.subframe.iter_mut() {
353 subframe.grid_index = br.read(1)? as usize;
354 }
355 if self.is_6300 {
356 br.read(1)?;
357 let mut ppos = br.read(13)? as usize;
358 self.subframe[0].pulse_pos = ppos / 810;
359 ppos = ppos % 810;
360 self.subframe[1].pulse_pos = ppos / 90;
361 ppos = ppos % 90;
362 self.subframe[2].pulse_pos = ppos / 9;
363 self.subframe[3].pulse_pos = ppos % 9;
364
365 for (i, subframe) in self.subframe.iter_mut().enumerate() {
366 let bits = if (i & 1) == 0 { 16 } else { 14 };
367 let val = br.read(bits)? as usize;
368 subframe.pulse_pos = (subframe.pulse_pos << bits) | val;
369 }
370 for (i, subframe) in self.subframe.iter_mut().enumerate() {
371 let bits = if (i & 1) == 0 { 6 } else { 5 };
372 subframe.pulse_sign = br.read(bits)? as u8;
373 }
374 } else {
375 for subframe in self.subframe.iter_mut() {
376 subframe.pulse_pos = br.read(12)? as usize;
377 }
378 for subframe in self.subframe.iter_mut() {
379 subframe.pulse_sign = br.read(4)? as u8;
380 }
381 }
382
383 Ok(())
384 }
385 fn synth_frame_active(&mut self, dst: &mut [i16], bad_frame: bool) {
386 const FCB_MAX_POS: [usize; 4] = [ 593775, 142506, 593775, 142506 ];
387 if !bad_frame {
388 self.erased_frames = 0;
389 } else {
390 self.erased_frames = (self.erased_frames + 1).min(MAX_ERASED_FRAMES);
391 }
392 if bad_frame {
393 self.lsp_index = [0; 3];
394 }
395 let mut cur_lsp = [0; LPC_ORDER];
396 Self::inverse_quant(&self.prev_lsp, &mut cur_lsp, &self.lsp_index, bad_frame);
397 Self::interpolate_lsp(&mut self.lpc, &cur_lsp, &self.prev_lsp);
398 self.prev_lsp.copy_from_slice(&cur_lsp);
399 (&mut self.excitation[..MAX_PITCH]).copy_from_slice(&self.prev_excitation);
400 if self.erased_frames == 0 {
401 let mut acb_vector = [0; SUBFRAME_LEN];
402 self.interp_gain = FIXED_CB_GAIN[(self.subframe[2].amp_index + self.subframe[3].amp_index) >> 1];
403 let mut exc_start = MAX_PITCH;
404 for (i, subframe) in self.subframe.iter_mut().enumerate() {
405 subframe.gen_fcb_excitation(&mut self.excitation[exc_start..], self.is_6300, FCB_MAX_POS[i], if (i & 1) == 0 { 6 } else { 5 });
406 subframe.gen_acb_excitation(&mut acb_vector, &self.excitation[SUBFRAME_LEN * i..], self.is_6300);
407 for i in 0..SUBFRAME_LEN {
408 let val = self.excitation[exc_start + i];
409 self.excitation[exc_start + i] = val.saturating_add(val).saturating_add(acb_vector[i]);
410 }
411 exc_start += SUBFRAME_LEN;
412 }
413 self.compute_interp_index();
414 let mut offset = MAX_PITCH;
415 for subframe in self.subframe.iter_mut() {
416 subframe.compute_ppf_coeffs(&self.synth_buf[LPC_ORDER..], offset, self.is_6300);
417 offset += SUBFRAME_LEN;
418 }
419 for i in 0..SUBFRAMES {
420 let src1 = &self.excitation[MAX_PITCH + i * SUBFRAME_LEN..][..SUBFRAME_LEN];
421 let src2 = &self.excitation[self.subframe[i].ppf_index..][..SUBFRAME_LEN];
422 let dst = &mut self.synth_buf[LPC_ORDER + i * SUBFRAME_LEN..][..SUBFRAME_LEN];
423 weighted_sum!(dst, src1, self.subframe[i].ppf_sc_gain, src2, self.subframe[i].ppf_opt_gain, 15);
424 }
425 self.prev_excitation.copy_from_slice(&self.excitation[SAMPLES..][..MAX_PITCH]);
426 } else {
427 self.interp_gain = (self.interp_gain * 3 + 2) >> 2;
428 if self.erased_frames == MAX_ERASED_FRAMES {
429 for el in self.excitation.iter_mut() {
430 *el = 0;
431 }
432 for el in self.prev_excitation.iter_mut() {
433 *el = 0;
434 }
435 for el in dst.iter_mut() {
436 *el = 0;
437 }
438 } else {
439 if self.interp_index != 0 {
440 for i in 0..self.interp_index {
441 let sample = i32::from(self.excitation[MAX_PITCH + i]) * 3 >> 2;
442 self.synth_buf[LPC_ORDER + i] = sample as i16;
443 }
444 for i in self.interp_index..SAMPLES {
445 self.synth_buf[LPC_ORDER + i] = self.synth_buf[LPC_ORDER + i - self.interp_index];
446 }
447 } else {
448 for i in 0..SAMPLES {
449 let sample = self.rnd.next().wrapping_mul(i32::from(self.interp_gain)) >> 15;
450 self.synth_buf[LPC_ORDER + i] = sample as i16;
451 }
452 }
453 self.prev_excitation.copy_from_slice(&self.synth_buf[LPC_ORDER + SAMPLES - MAX_PITCH..][..MAX_PITCH]);
454 }
455 }
456 self.cng_rnd.reset();
457 }
458 fn amp_index_to_sid_gain(val: usize) -> i32 {
459 if val < 0x10 {
460 (val as i32) << 6
461 } else if val < 0x20 {
462 ((val as i32) - 8) << 7
463 } else {
464 ((val as i32) - 20) << 8
465 }
466 }
467 fn estimate_sid_gain(cur_gain: i32, sid_gain: i32) -> i32 {
468 const CNG_BSEG: [i32; 3] = [ 2048, 18432, 231233 ];
469 let shift = 16 - cur_gain * 2;
470 let t = if shift > 0 { sid_gain << shift } else { sid_gain >> -shift };
471 let x = t * 273 >> 16;
472 if x >= CNG_BSEG[2] {
473 return 63;
474 }
475 let (shift, seg) = if x >= CNG_BSEG[1] { (4, 3) } else { (3, if x >= CNG_BSEG[0] { 1 } else { 0 }) };
476 let seg2 = seg.min(3);
477
478 let mut val = 1 << shift;
479 let mut val_add = val >> 1;
480 for _ in 0..shift {
481 let t = seg * 32 + (val << seg2);
482 let t2 = t * t;
483 if x >= t2 {
484 val += val_add;
485 } else {
486 val -= val_add;
487 }
488 val_add >>= 1;
489 }
490
491 let t = seg * 32 + (val << seg2);
492 let y = t * t - x;
493 if y <= 0 {
494 let t = seg * 32 + ((val + 1) << seg2);
495 let t2 = t * t - x;
496 let val = (seg2 - 1) * 16 + val;
497 if t2 >= y {
498 val + 1
499 } else {
500 val
501 }
502 } else {
503 let t = seg * 32 + ((val - 1) << seg2);
504 let t2 = t * t - x;
505 let val = (seg2 - 1) * 16 + val;
506 if t2 >= y {
507 val - 1
508 } else {
509 val
510 }
511 }
512 }
513 fn synth_frame_other(&mut self) {
514 if self.cur_ftype == G7231FrameType::SID {
515 self.sid_gain = Self::amp_index_to_sid_gain(self.subframe[0].amp_index);
516 Self::inverse_quant(&self.prev_lsp, &mut self.sid_lsp, &self.lsp_index, false);
517 } else if self.prev_ftype == G7231FrameType::Active {
518 self.sid_gain = Self::estimate_sid_gain(self.cur_gain, self.sid_gain);
519 }
520 if self.prev_ftype == G7231FrameType::Active {
521 self.cur_gain = self.sid_gain;
522 } else {
523 self.cur_gain = (self.cur_gain * 7 + self.sid_gain) >> 3;
524 }
525 self.generate_noise();
526 Self::interpolate_lsp(&mut self.lpc, &self.sid_lsp, &self.prev_lsp);
527 self.prev_lsp = self.sid_lsp;
528 }
529 fn generate_noise(&mut self) {
530 const ADAPTIVE_LAG: [usize; SUBFRAMES] = [1, 0, 1, 3];
531 self.subframe[0].pitch_lag = self.cng_rnd.next_range(21) + 123;
532 self.subframe[1].pitch_lag = self.subframe[0].pitch_lag;
533 self.subframe[2].pitch_lag = self.cng_rnd.next_range(19) + 123;
534 self.subframe[3].pitch_lag = self.subframe[2].pitch_lag;
535
536 for (i, subframe) in self.subframe.iter_mut().enumerate() {
537 subframe.acb_gain = self.cng_rnd.next_range(50) + 1;
538 subframe.acb_lag = ADAPTIVE_LAG[i];
539 }
540
541 let mut off = [0; 4];
542 let mut signs = [[0; 11]; 2];
543 for i in (0..4).step_by(2) {
544 let t = self.cng_rnd.next_range(1 << 13);
545 off[i] = t & 1;
546 off[i + 1] = ((t >> 1) & 1) + SUBFRAME_LEN;
547 for j in 0..11 {
548 signs[i/2][j] = ((((t >> (j + 2)) & 1) * 2 - 1) << 14) as i32;
549 }
550 }
551
552 let mut pos = [0; 11 * SUBFRAMES];
553 let mut pidx = 0;
554 let mut tmp = [0; SUBFRAME_LEN / 2];
555 for i in 0..SUBFRAMES {
556 let npulses = if (i & 1) == 0 { 6 } else { 5 };
557 for j in 0..npulses {
558 let idx = self.cng_rnd.next_range(SUBFRAME_LEN / 2 - j);
559 pos[pidx] = tmp[idx] * 2 + off[i];
560 pidx += 1;
561 tmp[idx] = tmp[SUBFRAME_LEN / 2 - 1 - j];
562 }
563 }
564
565 (&mut self.synth_buf[LPC_ORDER..][..MAX_PITCH]).copy_from_slice(&self.prev_excitation);
566 let mut acb_vec = [0; SUBFRAME_LEN];
567 let mut tmp = [0; SUBFRAME_LEN * 2];
568 for i in (0..SUBFRAMES).step_by(2) {
569 let buf = &mut self.synth_buf[LPC_ORDER + SUBFRAME_LEN * i..];
570 self.subframe[i].gen_acb_excitation(&mut acb_vec, buf, self.is_6300);
571 (&mut buf[..SUBFRAME_LEN]).copy_from_slice(&acb_vec);
572 self.subframe[i + 1].gen_acb_excitation(&mut acb_vec, &buf[SUBFRAME_LEN..], self.is_6300);
573 (&mut buf[SUBFRAME_LEN..][..SUBFRAME_LEN]).copy_from_slice(&acb_vec);
574
575 let mut max = 0;
576 for j in 0..SUBFRAME_LEN*2 {
577 max |= i32::from(buf[j]).abs();
578 }
579 let shift = if max == 0 { 0 } else {
580 (-10 + (32 - max.min(0x7FFF).leading_zeros()) as i32).max(-2)
581 };
582 let mut sum = 0;
583 if shift < 0 {
584 for j in 0..SUBFRAME_LEN*2 {
585 let val = buf[j] << -shift;
586 sum += i64::from(val) * i64::from(val);
587 tmp[j] = val;
588 }
589 } else {
590 for j in 0..SUBFRAME_LEN*2 {
591 let val = buf[j] >> shift;
592 sum += i64::from(val) * i64::from(val);
593 tmp[j] = val;
594 }
595 }
596
597 let mut b0 = 0;
598 for j in 0..11 {
599 b0 += i32::from(tmp[pos[i / 2 * 11 + j]]) * signs[i / 2][j];
600 }
601 b0 = ((i64::from(b0) * 2 * 2979 + (1 << 29)) >> 30) as i32;
602
603 let mut c = self.cur_gain * (self.cur_gain * (SUBFRAME_LEN as i32) >> 5);
604 if shift * 2 + 3 >= 0 {
605 c >>= shift * 2 + 3;
606 } else {
607 c <<= -(shift * 2 + 3);
608 }
609 c = (i64::from(clip32(sum * 2) - c) * 2979 >> 15) as i32;
610
611 let delta = b0 * b0 * 2 - c;
612 let x = if delta <= 0 {
613 -b0
614 } else {
615 let d0 = i32::from(square_root_i32(delta));
616 let x0 = d0 - b0;
617 let t = d0 + b0;
618 if t.abs() < x0.abs() {
619 -t
620 } else {
621 x0
622 }
623 };
624 let shift = shift + 1;
625 let x = (if shift >= 0 { x << shift } else { x >> -shift }).min(10000).max(-10000);
626 for j in 0..11 {
627 let val = x * signs[i / 2][j] >> 15;
628 buf[pos[i / 2 * 11 + j]] = buf[pos[i / 2 * 11 + j]].saturating_add(val as i16);
629 }
630
631 for j in 0..SUBFRAME_LEN*2 {
632 buf[MAX_PITCH + j] = buf[j];
633 }
634 }
635 self.prev_excitation.copy_from_slice(&self.synth_buf[LPC_ORDER + SAMPLES..][..MAX_PITCH]);
636 }
637 fn compute_interp_index(&mut self) {
638 let pitch_lag = self.subframe[3].pitch_lag;
639
640 self.cur_gain = Self::scale_vector(&mut self.synth_buf[LPC_ORDER..][..SAMPLES + MAX_PITCH], &self.excitation[..SAMPLES + MAX_PITCH]);
641
642 let (pos, cre) = autocorr_max(&self.synth_buf[LPC_ORDER..], MAX_PITCH + SUBFRAME_LEN * 2, SUBFRAME_LEN * 2, pitch_lag, false);
643 let corr_energy = cre.saturating_add(1 << 15) >> 16;
644
645 let tgt_energy = dot_product(&self.synth_buf[LPC_ORDER + MAX_PITCH + SUBFRAME_LEN * 2..], &self.synth_buf[LPC_ORDER + MAX_PITCH + SUBFRAME_LEN * 2..], SUBFRAME_LEN * 2);
646 self.sid_gain = tgt_energy.saturating_add(1 << 15) >> 16;
647
648 if corr_energy <= 0 {
649 self.interp_index = 0;
650 return;
651 }
652
653 let best_energy = dot_product(&self.synth_buf[LPC_ORDER + MAX_PITCH + SUBFRAME_LEN * 2 - pos..], &self.synth_buf[LPC_ORDER + MAX_PITCH + SUBFRAME_LEN * 2 - pos..], SUBFRAME_LEN * 2);
654 let best_energy = best_energy.saturating_add(1 << 15) >> 16;
655
656 let tmp = best_energy * self.sid_gain >> 3;
657
658 if tmp < corr_energy * corr_energy {
659 self.interp_index = pos;
660 } else {
661 self.interp_index = 0;
662 }
663 }
664 fn scale_vector(dst: &mut [i16], src: &[i16]) -> i32 {
665 let mut max = 0;
666 for el in src.iter() {
667 max = max.max(i32::from(*el).abs());
668 }
669 max = max.min(0x7FFF);
670 let shift = norm_bits(max, 15);
671 for (dst, src) in dst.iter_mut().zip(src.iter()) {
672 *dst = *src << shift >> 3;
673 }
674 i32::from(shift) - 3
675 }
676 fn inverse_quant(prev_lsp: &[i16; LPC_ORDER], lsp: &mut [i16; LPC_ORDER], lsp_index: &[usize; 3], bad_frame: bool) {
677 let (min_dist, pred) = if !bad_frame { (0x100, 12288) } else { (0x200, 23552) };
678 lsp[0] = LSP_CODEBOOK0[lsp_index[0]][0];
679 lsp[1] = LSP_CODEBOOK0[lsp_index[0]][1];
680 lsp[2] = LSP_CODEBOOK0[lsp_index[0]][2];
681 lsp[3] = LSP_CODEBOOK1[lsp_index[1]][0];
682 lsp[4] = LSP_CODEBOOK1[lsp_index[1]][1];
683 lsp[5] = LSP_CODEBOOK1[lsp_index[1]][2];
684 lsp[6] = LSP_CODEBOOK2[lsp_index[2]][0];
685 lsp[7] = LSP_CODEBOOK2[lsp_index[2]][1];
686 lsp[8] = LSP_CODEBOOK2[lsp_index[2]][2];
687 lsp[9] = LSP_CODEBOOK2[lsp_index[2]][3];
688
689 for i in 0..LPC_ORDER {
690 let diff = ((i32::from(prev_lsp[i]) - i32::from(DC_LSP[i])) * pred + (1 << 14)) >> 15;
691 lsp[i] = (i32::from(lsp[i]) + i32::from(DC_LSP[i]) + diff) as i16;
692 }
693
694 let mut stable = false;
695
696 for _ in 0..LPC_ORDER {
697 lsp[0] = lsp[0].max(0x100);
698 lsp[LPC_ORDER - 1] = lsp[LPC_ORDER - 1].min(0x7E00);
699
700 for i in 1..LPC_ORDER {
701 let mut val = min_dist + lsp[i - 1] - lsp[i];
702 if val > 0 {
703 val >>= 1;
704 lsp[i - 1] -= val;
705 lsp[i] += val;
706 }
707 }
708 stable = true;
709 for i in 1..LPC_ORDER {
710 let val = lsp[i - 1] + min_dist - lsp[i] - 4;
711 if val > 0 {
712 stable = false;
713 break;
714 }
715 }
716 if stable {
717 break;
718 }
719 }
720 if !stable {
721 lsp.copy_from_slice(prev_lsp);
722 }
723 }
724 fn interpolate_lsp(lpc: &mut [[i16; LPC_ORDER]; SUBFRAMES], lsp: &[i16; LPC_ORDER], prev_lsp: &[i16; LPC_ORDER]) {
725 weighted_sum!(lpc[0], lsp, 0x1000, prev_lsp, 0x3000, 14);
726 weighted_sum!(lpc[1], lsp, 0x2000, prev_lsp, 0x2000, 14);
727 weighted_sum!(lpc[2], lsp, 0x3000, prev_lsp, 0x1000, 14);
728 lpc[3].copy_from_slice(lsp);
729 for clpc in lpc.iter_mut() {
730 Self::lsp2lpc(clpc);
731 }
732 }
733 fn lsp2lpc(lpc: &mut [i16; LPC_ORDER]) {
734 let mut tmp1 = [0; LPC_ORDER/2 + 1];
735 let mut tmp2 = [0; LPC_ORDER/2 + 1];
736
737 for lpc in lpc.iter_mut() {
738 let index = ((*lpc >> 7) & 0x1FF) as usize;
739 let offset = (i32::from(*lpc & 0x7F) << 8) + 0x80;
740 let val1 = i32::from(COS_TAB[index]);
741 let val2 = (i32::from(COS_TAB[index + 1]) - val1) * offset * 2;
742 *lpc = -(((val1 << 17) + val2 * 2).saturating_add(1 << 15) >> 16) as i16;
743 }
744
745 tmp1[0] = 1 << 28;
746 tmp1[1] = (i32::from(lpc[0]) << 14) + (i32::from(lpc[2]) << 14);
747 tmp1[2] = i32::from(lpc[0]) * i32::from(lpc[2]) + (2 << 28);
748
749 tmp2[0] = 1 << 28;
750 tmp2[1] = (i32::from(lpc[1]) << 14) + (i32::from(lpc[3]) << 14);
751 tmp2[2] = i32::from(lpc[1]) * i32::from(lpc[3]) + (2 << 28);
752
753 for i in 2..LPC_ORDER/2 {
754 tmp1[i + 1] = tmp1[i - 1] + mul16(tmp1[i], lpc[i * 2]);
755 tmp2[i + 1] = tmp2[i - 1] + mul16(tmp2[i], lpc[i * 2 + 1]);
756
757 for j in (2..=i).rev() {
758 tmp1[j] = mul16(tmp1[j - 1], lpc[i * 2]) + (tmp1[j] >> 1) + (tmp1[j - 2] >> 1);
759 tmp2[j] = mul16(tmp2[j - 1], lpc[i * 2 + 1]) + (tmp2[j] >> 1) + (tmp2[j - 2] >> 1);
760 }
761
762 tmp1[0] >>= 1;
763 tmp2[0] >>= 1;
764 tmp1[1] = (tmp1[1] + (i32::from(lpc[2 * i]) << 16 >> i)) >> 1;
765 tmp2[1] = (tmp2[1] + (i32::from(lpc[2 * i + 1]) << 16 >> i)) >> 1;
766 }
767
768 for i in 0..LPC_ORDER/2 {
769 let c0 = i64::from(tmp1[i + 1]) + i64::from(tmp1[i]);
770 let c1 = i64::from(tmp2[i + 1]) - i64::from(tmp2[i]);
771 lpc[i] = (clip32((c0 + c1) * 8 + (1 << 15)) >> 16) as i16;
772 lpc[LPC_ORDER - i - 1] = (clip32((c0 - c1) * 8 + (1 << 15)) >> 16) as i16;
773 }
774 }
775 fn do_lpc(buf: &mut [i16], offset: usize, lpc: &[i16; LPC_ORDER]) {
776 for i in 0..SUBFRAME_LEN {
777 let mut sum = 0i32;
778 for j in 0..LPC_ORDER {
779 sum = sum.wrapping_sub(i32::from(buf[offset + i - j - 1]) * i32::from(lpc[j]));
780 }
781 let pred = (sum + (1 << 12)) >> 12;
782 buf[offset + i] = clip16((i32::from(buf[offset + i]) + pred) >> 1);
783 }
784 }
785 fn formant_postfilter(&mut self, dst: &mut [i16]) {
786 (&mut self.synth_buf[..LPC_ORDER]).copy_from_slice(&self.fir_mem);
787 let mut filter_data = [0; LPC_ORDER + SAMPLES];
788 (&mut filter_data[..LPC_ORDER]).copy_from_slice(&self.iir_mem);
789
790 let mut filter_coef = [[0; LPC_ORDER]; 2];
791 for i in 0..SUBFRAMES {
792 for j in 0..LPC_ORDER {
793 filter_coef[0][j] = (i32::from(-self.lpc[i][j]) * i32::from(POSTFILTER_COEFFS[0][j]) + (1 << 14)) >> 15;
794 filter_coef[1][j] = (i32::from(-self.lpc[i][j]) * i32::from(POSTFILTER_COEFFS[1][j]) + (1 << 14)) >> 15;
795 }
796 Self::iir_filter(&filter_coef, &self.synth_buf, &mut filter_data, LPC_ORDER + i * SUBFRAME_LEN);
797 }
798 self.fir_mem.copy_from_slice(&self.synth_buf[SAMPLES..][..LPC_ORDER]);
799 self.iir_mem.copy_from_slice(&filter_data[SAMPLES..]);
800 let mut offset = 0;
801 for _ in 0..SUBFRAMES {
802 let scale = Self::scale_vector(&mut dst[offset..][..SUBFRAME_LEN], &self.synth_buf[offset + LPC_ORDER..][..SUBFRAME_LEN]);
803 let ac1 = dot_product(&dst[offset..], &dst[offset + 1..], SUBFRAME_LEN - 1);
804 let ac0 = dot_product(&dst[offset..], &dst[offset..], SUBFRAME_LEN);
805 let tmp = if (ac0 >> 16) != 0 { (ac1 >> 2) / (ac0 >> 16) } else { 0 };
806 self.reflection_coef = (3 * self.reflection_coef + tmp + 2) >> 2;
807 let gain = (-self.reflection_coef >> 1) & !3;
808
809 for i in 0..SUBFRAME_LEN {
810 let val = (filter_data[offset + LPC_ORDER + i - 1] >> 16) * gain;
811 dst[offset + i] = (filter_data[offset + LPC_ORDER + i].saturating_add(val).saturating_add(val) >> 16) as i16;
812 }
813
814 let shift = 2 * scale + 4;
815 let energy = if shift < 0 { clip32(i64::from(ac0) << -shift) } else { ac0 >> shift };
816 Self::gain_scale(&mut self.pf_gain, &mut dst[offset..], energy);
817
818 offset += SUBFRAME_LEN;
819 }
820 }
821 fn iir_filter(coef: &[[i32; LPC_ORDER]; 2], src: &[i16], dst: &mut [i32], offset: usize) {
822 for i in 0..SUBFRAME_LEN {
823 let mut sum = 0;
824 for j in 0..LPC_ORDER {
825 sum -= i64::from(coef[0][j]) * i64::from(src[offset - 1 + i - j]) -
826 i64::from(coef[1][j]) * i64::from(dst[offset - 1 + i - j] >> 16);
827 }
828 dst[offset + i] = clip32((i64::from(src[offset + i]) << 16) + (sum << 3) + (1 << 15));
829 }
830 }
831 fn gain_scale(pf_gain: &mut i32, buf: &mut [i16], energy: i32) {
832 let mut den = 0i32;
833 for i in 0..SUBFRAME_LEN {
834 let val = i32::from(buf[i] >> 2);
835 den = den.saturating_add(val * val).saturating_add(val * val);
836 }
837 let mut num = energy;
838 let gain;
839 if (num != 0) && (den != 0) {
840 let bits1 = norm_bits(num, 31);
841 let bits2 = norm_bits(den, 31);
842 num = num << bits1 >> 1;
843 den <<= bits2;
844 let shift = (5 + bits1 - bits2).max(0);
845
846 gain = square_root_i32(((num >> 1) / (den >> 16)) << 16 >> shift);
847 } else {
848 gain = 1 << 12;
849 }
850
851 for i in 0..SUBFRAME_LEN {
852 *pf_gain = (15 * *pf_gain + gain + (1 << 3)) >> 4;
853 buf[i] = clip16((i32::from(buf[i]) * (*pf_gain + (*pf_gain >> 4)) + (1 << 10)) >> 11);
854 }
855 }
856 }
857
858 fn clip16(a: i32) -> i16 {
859 a.min(i32::from(std::i16::MAX)).max(i32::from(std::i16::MIN)) as i16
860 }
861
862 fn clip32(a: i64) -> i32 {
863 a.min(i64::from(std::i32::MAX)).max(i64::from(std::i32::MIN)) as i32
864 }
865
866 fn mul16(a: i32, b: i16) -> i32 {
867 let b = i32::from(b);
868 (a >> 16) * b * 2 + (((a & 0xFFFF) * b) >> 15)
869 }
870
871 fn norm_bits(val: i32, target: u8) -> u8 {
872 if val == 0 {
873 target
874 } else {
875 target - (32 - val.leading_zeros()) as u8
876 }
877 }
878
879 fn autocorr_max(src: &[i16], offset: usize, length: usize, pitch_lag: usize, forward: bool) -> (usize, i32) {
880 let pitch_lag = pitch_lag.min(MAX_PITCH - 3);
881 let mut max_energy = 0;
882 let mut lag = 0;
883
884 if forward {
885 let end = (pitch_lag + 3).min(SAMPLES + MAX_PITCH - offset - length);
886 for i in pitch_lag-3..=end {
887 let energy = dot_product(&src[offset..], &src[offset + i..], length);
888 if max_energy < energy {
889 max_energy = energy;
890 lag = i;
891 }
892 }
893 } else {
894 for i in pitch_lag-3..=pitch_lag+3 {
895 let energy = dot_product(&src[offset..], &src[offset - i..], length);
896 if max_energy < energy {
897 max_energy = energy;
898 lag = i;
899 }
900 }
901 }
902 (lag, max_energy)
903 }
904
905 fn dot_product(src1: &[i16], src2: &[i16], length: usize) -> i32 {
906 let mut sum = 0;
907 for (a, b) in src1.iter().zip(src2.iter()).take(length) {
908 sum += i64::from(*a) * i64::from(*b)
909 }
910 clip32(sum << 1)
911 }
912
913 impl NADecoder for G7231Decoder {
914 fn init(&mut self, _supp: &mut NADecoderSupport, info: NACodecInfoRef) -> DecoderResult<()> {
915 if let NACodecTypeInfo::Audio(_ainfo) = info.get_properties() {
916 self.info = info.replace_info(NACodecTypeInfo::Audio(self.ainfo));
917 Ok(())
918 } else {
919 Err(DecoderError::InvalidData)
920 }
921 }
922 fn decode(&mut self, _supp: &mut NADecoderSupport, pkt: &NAPacket) -> DecoderResult<NAFrameRef> {
923 let info = pkt.get_stream().get_info();
924 validate!(info.get_properties().is_audio());
925 let src = pkt.get_buffer();
926
927 let mut bad_frame = false;
928
929 let mut br = BitReader::new(src.as_slice(), BitReaderMode::LE);
930 if let Err(_) = self.unpack_frame(&mut br) {
931 bad_frame = true;
932 self.cur_ftype = if self.prev_ftype == G7231FrameType::Active {
933 G7231FrameType::Active
934 } else {
935 G7231FrameType::Untransmitted
936 };
937 }
938
939 let abuf = alloc_audio_buffer(self.ainfo, SAMPLES, self.chmap.clone())?;
940 let mut adata = abuf.get_abuf_i16().unwrap();
941 let mut asamples = adata.get_data_mut().unwrap();
942 if self.cur_ftype == G7231FrameType::Active {
943 self.synth_frame_active(&mut asamples, bad_frame);
944 } else {
945 self.synth_frame_other();
946 }
947 self.prev_ftype = self.cur_ftype;
948 (&mut self.synth_buf[..LPC_ORDER]).copy_from_slice(&self.filt_mem);
949 for i in 0..SUBFRAMES {
950 Self::do_lpc(&mut self.synth_buf, LPC_ORDER + i * SUBFRAME_LEN, &self.lpc[i]);
951 }
952 self.filt_mem.copy_from_slice(&self.synth_buf[SAMPLES..][..LPC_ORDER]);
953 self.formant_postfilter(asamples);
954
955 let mut frm = NAFrame::new_from_pkt(pkt, self.info.clone(), abuf);
956 frm.set_keyframe(true);
957 Ok(frm.into_ref())
958 }
959 fn flush(&mut self) {
960 }
961 }
962
963 pub fn get_decoder() -> Box<dyn NADecoder + Send> {
964 Box::new(G7231Decoder::new())
965 }
966
967 #[cfg(test)]
968 mod test {
969 use nihav_core::codecs::RegisteredDecoders;
970 use nihav_core::demuxers::RegisteredDemuxers;
971 use nihav_codec_support::test::dec_video::*;
972 use crate::vivo_register_all_codecs;
973 use crate::vivo_register_all_demuxers;
974 #[test]
975 fn test_g723_1() {
976 let mut dmx_reg = RegisteredDemuxers::new();
977 vivo_register_all_demuxers(&mut dmx_reg);
978 let mut dec_reg = RegisteredDecoders::new();
979 vivo_register_all_codecs(&mut dec_reg);
980
981 let file = "assets/Misc/adalogo.viv";
982 //let file = "assets/Misc/gr_al.viv";
983 //test_decode_audio("vivo", file, Some(1500), None/*Some("g7231")*/, &dmx_reg, &dec_reg);
984 test_decoding("vivo", "g723.1", file, None, &dmx_reg, &dec_reg,
985 ExpectedTestResult::MD5([0x94e60e2f, 0x578493e7, 0xe7ab1f1c, 0xae155977]));
986 }
987 }
988
989 const DC_LSP: [i16; LPC_ORDER] = [
990 0x0c3b, 0x1271, 0x1e0a, 0x2a36, 0x3630,
991 0x406f, 0x4d28, 0x56f4, 0x638c, 0x6c46
992 ];
993
994 const LSP_CODEBOOK0: [[i16; 3]; 256] = [
995 [ 0, 0, 0], [ -270, -1372, -1032], [ -541, -1650, -1382],
996 [ -723, -2011, -2213], [ -941, -1122, -1942], [ -780, -1145, -2454],
997 [ -884, -1309, -1373], [-1051, -1523, -1766], [-1083, -1622, -2300],
998 [ -777, -1377, -2147], [ -935, -1467, -2763], [ -802, -1327, -3471],
999 [ -935, -1959, -3999], [ -240, -89, 222], [ -661, -257, -160],
1000 [ -994, -466, -419], [ -188, -164, -278], [ -342, -512, -415],
1001 [ -607, -511, -797], [ 16, 19, -716], [ 374, 425, -972],
1002 [ -346, 245, -282], [ -265, 506, -754], [ -620, -147, 1955],
1003 [ -742, -860, 2597], [ -150, -352, 2704], [ 305, 880, 1954],
1004 [ 123, 731, 2766], [ -348, 765, 3327], [ 618, 221, 3258],
1005 [ -178, -47, 4219], [ 393, 1304, 3842], [ 698, 1702, 4801],
1006 [ 63, -584, 1229], [ -215, -732, 1704], [ 172, -335, 1909],
1007 [ -2, 216, 1797], [ 353, 127, 2205], [-1208, 188, 11],
1008 [ -513, -75, -683], [ -973, 222, -646], [ -616, -843, -388],
1009 [ -950, -1113, -359], [-1431, -623, -705], [-1398, -1063, -178],
1010 [ -45, -461, 35], [ -9, -657, -216], [ 127, -1078, 95],
1011 [ -950, -1156, 584], [-1480, -1494, 449], [ -120, -705, 516],
1012 [ -368, -961, 727], [ -378, -526, 973], [ -793, -614, 676],
1013 [ -801, -755, 1287], [-1476, -340, 1636], [ -505, -1254, 1543],
1014 [-1243, -1622, 1532], [ -776, -1477, -655], [-1151, -1296, -823],
1015 [-1153, -1672, -1124], [-1291, -2003, -1702], [ -622, -1283, 57],
1016 [ -471, -1611, 509], [-1060, -1570, -139], [ -873, -2156, -536],
1017 [-1716, -2021, -364], [-2150, -3218, -1291], [-1248, -1945, -2904],
1018 [-1215, -2633, -2855], [ 167, -244, 84], [ 349, -412, -217],
1019 [ -40, -352, 632], [ 227, -529, 405], [ 68, -383, -443],
1020 [ 167, -558, -706], [ -275, -854, -14], [ -351, -1089, -449],
1021 [ 341, -72, -289], [ 603, -106, -474], [ 322, -219, -649],
1022 [ 179, -317, -998], [ 450, -291, -996], [ 555, 195, -525],
1023 [ 784, 272, -831], [ -148, -384, -849], [ 82, -536, -1357],
1024 [ 238, -172, -1354], [ 422, -268, -1841], [ 297, -737, -2079],
1025 [ -111, -801, -598], [ 1, -668, -984], [ -131, -818, -1299],
1026 [ -329, -521, -1310], [ -151, -778, -1834], [ -93, -352, -1746],
1027 [ -568, -640, -1821], [ -509, -941, -2183], [ 464, -815, -1250],
1028 [ 79, -1133, -1597], [ -184, -1353, -2123], [ -196, -410, -2427],
1029 [ -192, -833, -2810], [ -259, -1382, -3045], [ -217, 4, -1166],
1030 [ -800, -325, -1219], [ -363, -830, -898], [ -661, -1134, -960],
1031 [ -386, -980, -1501], [ -627, -1159, -1722], [ -903, -829, -855],
1032 [ -685, -829, -1313], [-1065, -959, -1405], [ 441, 25, -847],
1033 [ 655, -27, -1181], [ 1159, -110, -705], [ 856, 253, -1671],
1034 [ 415, 404, -1], [ 322, 903, -398], [ 670, 499, -292],
1035 [ 803, 591, -610], [ 1144, 591, -814], [ 717, 183, 393],
1036 [ 857, 381, 106], [ 609, 62, -27], [ 792, 198, -325],
1037 [ 735, 805, 88], [ 1142, 812, 78], [ 1028, 366, -292],
1038 [ 1309, 743, -237], [ 1615, 589, -79], [ 1010, 639, -243],
1039 [ 999, 964, -311], [ 1500, 1137, -615], [ 988, 357, 646],
1040 [ 1227, 667, 683], [ 1164, 1565, 894], [ 1392, 2015, 477],
1041 [ 1138, 533, 250], [ 1437, 896, 391], [ 1765, 1118, 99],
1042 [ 1112, 1090, 802], [ 1596, 846, 1134], [ 937, 1161, 279],
1043 [ 1719, 1254, 683], [ 1338, 1086, 35], [ 1419, 1324, 428],
1044 [ 1428, 1524, 40], [ 2108, 1594, 89], [ 1015, 544, 1222],
1045 [ 1121, 925, 1263], [ 1030, 1318, 1485], [ 1295, 789, 1817],
1046 [ 1323, 1272, 1909], [ 1724, 1237, 1803], [ 1797, 1689, 858],
1047 [ 2149, 1367, 1301], [ 2302, 1867, 761], [ 2863, 2351, 1053],
1048 [ 52, 163, -76], [ 230, 309, -492], [ -71, 619, 39],
1049 [ -218, 856, 499], [ -654, 736, -207], [ -535, 1259, 155],
1050 [ -480, 1476, 643], [ 262, 1081, 102], [ 309, 1592, -182],
1051 [ 627, 1629, 534], [ 337, 643, 456], [ 758, 670, 713],
1052 [ 202, 1126, 658], [ 612, 1131, 666], [ 686, 1223, 1136],
1053 [ -131, 377, 525], [ 42, 708, 907], [ 87, 1488, 1035],
1054 [ 432, 2117, 904], [ 137, 981, 1332], [ -447, 1014, 1136],
1055 [ -839, 1793, 1246], [ -559, 297, 198], [ -850, 685, 446],
1056 [-1273, 632, 826], [ -401, -544, 173], [ -753, -793, 144],
1057 [ -436, -9, 772], [ -115, -243, 1310], [ -670, -269, 374],
1058 [-1027, -13, 639], [ -887, -81, 1137], [-1277, -455, 158],
1059 [-1411, -720, 736], [ 172, 88, 403], [ 386, 255, 756],
1060 [ -500, 522, 910], [ -958, 659, 1388], [ -395, 301, 1344],
1061 [ -356, 768, 1813], [ -613, 841, 2419], [ 445, -122, 252],
1062 [ 629, -87, 723], [ 283, -253, 870], [ 456, -116, 1381],
1063 [ 757, 180, 1059], [ 532, 408, 1509], [ 947, 288, 1806],
1064 [ 1325, 994, 2524], [ 892, 1219, 3023], [ 1397, 1596, 3406],
1065 [ 1143, 1552, 2546], [ 1850, 1433, 2710], [ -10, 134, 1002],
1066 [ 154, 499, 1323], [ 508, 792, 1117], [ 509, 1340, 1616],
1067 [ 762, 862, 1608], [ 787, 740, 2320], [ 794, 1727, 1283],
1068 [ 465, 2108, 1660], [ -120, 1451, 1613], [ -386, 2016, 2169],
1069 [ 891, 1225, 2050], [ 456, 1480, 2185], [ 1493, 1283, 1209],
1070 [ 1397, 1636, 1518], [ 1776, 1738, 1552], [ 1572, 1698, 2141],
1071 [ 1389, 2126, 1271], [ 1959, 2413, 1119], [ 1365, 2892, 1505],
1072 [ 2206, 1971, 1623], [ 2076, 1950, 2280], [ 1717, 2291, 1867],
1073 [ 2366, 2515, 1953], [ 2865, 2838, 2522], [ 2535, 3465, 2011],
1074 [ 3381, 4127, 2638], [ 836, 2667, 2289], [ 1761, 2773, 2337],
1075 [ 1415, 3325, 2911], [ 2354, 3138, 3126], [ 2659, 4192, 4010],
1076 [ 1048, 1786, 1818], [ 1242, 2111, 2240], [ 1512, 2079, 2780],
1077 [ 1573, 2491, 3138], [ 2230, 2377, 2782], [ 416, 1773, 2704],
1078 [ 725, 2336, 3297], [ 1252, 2373, 3978], [ 2094, 2268, 3568],
1079 [ 2011, 2712, 4528], [ 1341, 3507, 3876], [ 1216, 3919, 4922],
1080 [ 1693, 4793, 6012]
1081 ];
1082
1083 const LSP_CODEBOOK1: [[i16; 3]; 256] = [
1084 [ 0, 0, 0], [-2114, -1302, 76], [-2652, -1278, -1368],
1085 [-2847, -828, -349], [-3812, -2190, -349], [-3946, -364, -449],
1086 [-2725, -4492, -3607], [-3495, -4764, -1744], [ -51, -756, 84],
1087 [ -153, -1191, 504], [ 108, -1418, 1167], [ -835, -896, 390],
1088 [ -569, -1702, 87], [-1151, -1818, 933], [-1826, -2547, 411],
1089 [-1842, -1818, 1451], [-2438, -1611, 781], [-2747, -2477, 1311],
1090 [ -940, 1252, 477], [-1629, 1688, 602], [-1202, 617, 280],
1091 [-1737, 393, 580], [-1528, 1077, 1199], [-2165, -161, 1408],
1092 [-2504, -1087, 2371], [-3458, -175, 1395], [-1397, -98, -843],
1093 [-2252, -177, -1149], [-1489, -726, -1283], [-1558, -265, -1744],
1094 [-1867, -821, -1897], [-2062, -1516, -2340], [-2595, -1142, -2861],
1095 [ 170, 46, -819], [ -193, -204, -1151], [ 326, -196, -1532],
1096 [ 780, 329, -816], [ 201, 369, -1243], [ 650, -209, -1060],
1097 [ 1144, -15, -1216], [ 1203, -259, -1867], [ -890, -564, -1430],
1098 [ -638, -852, -1921], [ 177, -739, -1358], [ -261, -526, -1666],
1099 [ 206, -407, -2255], [ 338, -526, -822], [ 421, -1095, -1009],
1100 [ 765, -607, -1408], [ 825, -1295, -2004], [ 357, -905, -1815],
1101 [ -58, -1248, -1588], [ -596, -1436, -2046], [ -73, -1159, -2116],
1102 [ -115, -1382, -2581], [ -160, -1723, -1952], [ -6, -2196, -2954],
1103 [ -649, -1705, -2603], [ -617, -1453, -3282], [ -949, -2019, -3102],
1104 [ -812, 1544, 1937], [-1854, 574, 2000], [-1463, 1140, 2649],
1105 [-2683, 1748, 1452], [-2486, 2241, 2523], [ 783, 1910, 1435],
1106 [ 581, 2682, 1376], [ 236, 2197, 1885], [ -453, 2943, 2057],
1107 [ -682, 2178, 2565], [-1342, 3201, 3328], [ -288, -184, 262],
1108 [ 121, -149, -183], [ 758, -412, 206], [ 1038, -204, 853],
1109 [ 1577, -457, 700], [ 937, -640, -567], [ 1508, -528, -1024],
1110 [ -225, -527, -427], [ -564, -1095, -332], [ -742, -353, -186],
1111 [-1288, -459, 84], [-1853, -484, -274], [-1554, -731, 825],
1112 [-2425, -234, 382], [-1722, 293, -271], [-2515, 425, -564],
1113 [-2599, 818, 464], [ -358, 118, -375], [ -613, 198, -874],
1114 [ -690, 683, -324], [-1352, 1155, -168], [-1093, 129, -324],
1115 [-1184, 611, -858], [ 433, 386, -372], [ -120, 486, -634],
1116 [ 234, 851, -631], [ 602, 128, 46], [ 1099, 410, 159],
1117 [ 715, -145, -424], [ 1198, -85, -593], [ 1390, 367, -358],
1118 [ 1683, 362, -964], [ 1711, 622, 45], [ 2033, 833, -383],
1119 [ 2890, 549, -506], [ 7, 401, 52], [ 72, 811, 415],
1120 [ 566, 668, 41], [ 467, 1218, 130], [ 68, 957, -187],
1121 [ -25, 1649, -103], [ -661, 260, 214], [ -925, -94, 612],
1122 [ -321, -422, 965], [ -788, -672, 1783], [ 400, -673, 779],
1123 [ 741, -595, 1635], [ -161, 307, 657], [ -382, 836, 871],
1124 [ -814, 400, 1223], [ 364, 606, 1247], [ 57, 75, 1571],
1125 [ 151, 471, 2287], [ -81, 1021, 1502], [ 227, 1470, 1097],
1126 [ 658, 1275, 1653], [ 664, 1478, 2377], [ 263, -127, 444],
1127 [ 264, 89, 969], [ 794, 171, 576], [ 821, 186, 1226],
1128 [ 404, 462, 517], [ 339, 918, 794], [ 1280, 1423, 196],
1129 [ 1453, 2019, 365], [ 1615, 1481, 672], [ 2394, 1708, 508],
1130 [ 806, 1238, 573], [ 713, 1158, 1078], [ 1285, 1436, 1232],
1131 [ 1790, 1188, 1141], [ 765, 643, 864], [ 1032, 797, 1279],
1132 [ 900, 563, 1827], [ 1514, 673, 2312], [ 1544, 1129, 3240],
1133 [ 1469, 1050, 1594], [ 1945, 1318, 1988], [ 2397, 2026, 2060],
1134 [ 3538, 2057, 2620], [ 1249, -118, 74], [ 1727, 194, 421],
1135 [ 2078, -50, -463], [ 970, 688, -432], [ 1149, 952, -110],
1136 [ 1254, 1275, -651], [ 1386, 929, 401], [ 1960, 1167, 232],
1137 [ 407, -752, -243], [ 859, -1118, 172], [ -227, -860, -992],
1138 [ -796, -1175, -1380], [ 8, -1282, -388], [ 353, -1781, -1037],
1139 [ -732, -397, -807], [ -853, -28, -1342], [-1229, -1207, -1959],
1140 [-1015, -1125, -2543], [-1452, -1791, -2725], [-1891, -2416, -3269],
1141 [ -918, -1629, -783], [ -580, -2155, -698], [-1097, -2364, -96],
1142 [-1387, -1513, 7], [-1588, -2076, -664], [-1473, -2740, -784],
1143 [-2378, -3149, -56], [-2856, -2092, -169], [-3391, -3708, 316],
1144 [-1176, -890, -614], [-1944, -1061, -800], [ -299, -1517, -1000],
1145 [ -640, -1850, -1526], [-1454, -1536, -1233], [-1890, -1955, -1756],
1146 [-1086, -1921, -2122], [ -750, -2325, -2260], [-1325, -2413, -2673],
1147 [-1114, -2542, -3459], [-1341, -2901, -3963], [-1160, -2226, -1393],
1148 [-1001, -2772, -1573], [-1594, -2641, -1978], [-1534, -3046, -2624],
1149 [-2224, -2196, -675], [-2807, -3054, -1102], [-2008, -2840, -1186],
1150 [-1980, -3332, -1695], [-1715, -3562, -505], [-2527, -4000, -1887],
1151 [-2333, -2734, -2296], [-3440, -2401, -3211], [-2008, -3528, -3337],
1152 [-2247, -3291, -4510], [ -475, 949, 155], [ -149, 1365, 545],
1153 [ -757, 1644, 1083], [ -217, 2053, 1353], [-1433, 2301, 1462],
1154 [ 495, 1661, 529], [ 10, 2037, 740], [ 2082, 1898, 978],
1155 [ 2831, 2294, 911], [ 842, 793, 420], [ 1223, 1023, 863],
1156 [ 1237, 451, 780], [ 1744, 708, 822], [ 1533, 284, 1384],
1157 [ 2135, 609, 1538], [ 2305, 626, 540], [ 2368, 1187, 955],
1158 [ 2586, 1255, -7], [ 3116, 1131, 726], [ 3431, 1730, 428],
1159 [ 2734, 1648, 1307], [ 2988, 1231, 2010], [ 3523, 2024, 1488],
1160 [ 1034, 1657, 871], [ 1206, 2163, 1036], [ 1807, 2372, 1233],
1161 [ 1808, 1769, 1493], [ 1573, 2332, 1779], [ 1216, 1609, 1866],
1162 [ 1480, 1898, 2513], [ 465, 2708, 2776], [ 771, 3638, 3338],
1163 [ 1869, 2599, 2623], [ 2825, 2745, 2468], [ 2638, 2439, 1585],
1164 [ 2094, 2970, 1308], [ 2022, 3057, 1999], [ 3428, 2912, 1816],
1165 [ 4536, 2974, 2129], [ 1046, 2563, 2086], [ 1363, 3562, 2318],
1166 [ 2511, 1891, 2984], [ 1866, 2306, 3986], [ 3272, 2924, 3682],
1167 [ 3146, 3564, 2272], [ 3592, 3968, 2822], [ 2431, 3369, 3069],
1168 [ 1931, 4709, 3090], [ 2629, 4220, 3986], [ 4639, 4056, 3664],
1169 [ 4035, 5334, 4912]
1170 ];
1171
1172 const LSP_CODEBOOK2: [[i16; 4]; 256] = [
1173 [ 0, 0, 0, 0], [ 601, 512, -542, 334],
1174 [ 428, 1087, -484, -132], [ 652, 622, -391, -572],
1175 [ 378, 799, 141, -860], [ 1040, 409, 112, -554],
1176 [ 1123, 670, -75, -847], [ 1421, 494, -315, -1095],
1177 [ 787, 1001, 114, -460], [ 988, 1672, 216, -681],
1178 [ 1007, 1241, -132, -1247], [ 1073, 399, 186, -5],
1179 [ 1262, 193, -694, -129], [ 325, 196, 51, -641],
1180 [ 861, -59, 350, -458], [ 1261, 567, 586, -346],
1181 [ 1532, 885, 210, -517], [ 2027, 937, 113, -792],
1182 [ 1383, 1064, 334, 38], [ 1964, 1468, 459, 133],
1183 [ 2062, 1186, -98, -121], [ 2577, 1445, 506, -373],
1184 [ 2310, 1682, -2, -960], [ 2876, 1939, 765, 138],
1185 [ 3581, 2360, 649, -414], [ 219, 176, -398, -309],
1186 [ 434, -78, -435, -880], [ -344, 301, 265, -552],
1187 [ -915, 470, 657, -380], [ 419, -432, -163, -453],
1188 [ 351, -953, 8, -562], [ 789, -43, 20, -958],
1189 [ 302, -594, -352, -1159], [ 1040, 108, -668, -924],
1190 [ 1333, 210, -1217, -1663], [ 483, 589, -350, -1140],
1191 [ 1003, 824, -802, -1184], [ 745, 58, -589, -1443],
1192 [ 346, 247, -915, -1683], [ 270, 796, -720, -2043],
1193 [ 1208, 722, -222, -193], [ 1486, 1180, -412, -672],
1194 [ 1722, 179, -69, -521], [ 2047, 860, -666, -1410],
1195 [ -146, 222, -281, -805], [ -189, 90, -114, -1307],
1196 [ -152, 1086, -241, -764], [ -439, 733, -601, -1302],
1197 [ -833, -167, -351, -601], [ -856, -422, -411, -1059],
1198 [ -747, -355, -582, -1644], [ -837, 210, -916, -1144],
1199 [-1800, 32, -878, -1687], [ -48, -23, -1146, 52],
1200 [ -350, -409, -1656, -364], [ 265, -728, -858, -577],
1201 [ 458, -247, -1141, -997], [ 691, -407, -1988, -1161],
1202 [ -66, -104, -705, -1249], [ -431, -93, -1191, -1844],
1203 [ 203, -732, -1000, -1693], [ 10, -832, -1846, -1819],
1204 [ 493, -128, -1436, -1768], [ 488, -311, -1730, -2540],
1205 [ -653, -532, -1150, -1172], [-1086, -289, -1706, -1533],
1206 [ -699, -1205, -1216, -1766], [-1032, -1481, -2074, -1523],
1207 [ -721, -1220, -2277, -2600], [ 12, -539, -1484, -1131],
1208 [ -40, -911, -2106, -441], [ -471, -484, -2267, -1549],
1209 [ -141, -988, -3006, -1721], [-1545, -2102, -583, 342],
1210 [-1383, -2772, -386, -13], [-2118, -2589, -1205, 72],
1211 [-2147, -3231, -965, 390], [-2949, -3300, -621, 637],
1212 [-3907, -4138, -865, 803], [-1287, -845, -375, -548],
1213 [-1416, -1169, -487, -1277], [-1400, -1690, -1027, -418],
1214 [-2018, -1909, -1188, -1260], [-1418, -2222, -2029, -128],
1215 [-2067, -2998, -2693, -310], [ -950, -1028, -1538, 185],
1216 [-1616, -915, -2205, -549], [ 19, -821, -1145, 352],
1217 [ 184, -1175, -1356, -627], [ -547, -1088, -1661, -911],
1218 [ -216, -1502, -2197, -948], [ -795, -1306, -2374, -451],
1219 [ -924, -1889, -2796, -680], [ -600, -1614, -3609, -885],
1220 [-2392, -2528, 319, 303], [-2908, -2095, -310, 573],
1221 [-3460, -2141, 49, -113], [-2231, -448, 675, -146],
1222 [-2805, -532, 1231, 479], [-2684, -486, -200, 611],
1223 [-3525, -971, -198, 704], [-3707, 173, 349, 254],
1224 [-4734, -1447, -34, 880], [ 777, -512, 114, -10],
1225 [ 1250, -66, 442, -5], [ 604, 613, 452, -352],
1226 [ 1224, 777, 675, -1014], [-1372, -79, -1208, -238],
1227 [-2389, -17, -1157, -818], [-1504, -673, -1133, -1060],
1228 [-1984, -799, -2005, -1973], [-2037, -798, -1068, -105],
1229 [-3190, -899, -1817, -194], [ -156, -886, 394, -318],
1230 [ -258, -1283, 551, 202], [ -536, -1729, 910, 331],
1231 [ -847, -1109, 795, -163], [-1171, -1128, 715, 519],
1232 [-1080, -1319, 1685, 668], [-1000, -1921, 96, 211],
1233 [-1487, -2148, 831, 174], [-1139, -374, 414, -4],
1234 [-1517, -1383, 396, -352], [-1012, 439, -59, -967],
1235 [-1812, 706, -440, -1030], [-1971, -329, -34, -827],
1236 [-2472, -1588, -151, -606], [-2161, 374, -281, 76],
1237 [-3012, 231, -15, -690], [ 1104, 566, 721, 209],
1238 [ 1685, 564, 383, 98], [ 1898, 750, 792, -97],
1239 [ 556, -64, 561, -93], [ 876, 162, 913, -22],
1240 [ 961, 675, 1296, 140], [ 756, -396, 851, 544],
1241 [ 360, -303, 1341, 396], [ 878, -22, 1464, 863],
1242 [ -309, -273, 642, -129], [ -686, -82, 842, 454],
1243 [ -5, -47, 1069, 998], [ -94, 967, 1277, 298],
1244 [ -489, 385, 1473, 746], [ -369, -717, 1333, 242],
1245 [ 281, -993, 1726, 924], [ 464, 601, 1575, 1376],
1246 [ -250, 206, 2339, 1175], [ -438, 377, -597, -285],
1247 [-1020, 787, -790, -287], [ -458, -410, 215, 295],
1248 [ -589, -860, -121, 797], [-1175, 122, -437, 466],
1249 [-1480, -121, 367, 924], [ 234, 323, 770, -555],
1250 [ 145, 30, 996, 26], [ 66, 849, 93, -145],
1251 [ -117, 1261, 474, -399], [-1495, 1051, 218, -506],
1252 [-1390, 694, 994, 88], [ 616, 7, 78, 304],
1253 [ 1060, 52, -62, 835], [ 833, 454, 649, 1359],
1254 [ -770, 464, 47, 93], [ -574, 1199, -39, 379],
1255 [ 114, -98, 488, 485], [ 727, 244, 606, 696],
1256 [ -76, 455, 671, 546], [ -565, -13, 145, 819],
1257 [ -376, 569, 448, 1128], [ 218, 122, 265, 1167],
1258 [ 230, 738, 932, 1003], [ 138, 477, 36, 450],
1259 [ 404, 787, -73, 1000], [ 497, 1259, 387, 1231],
1260 [ 17, 207, 195, -79], [ 562, 358, 53, -158],
1261 [ 493, 387, 478, 189], [ 678, 831, 640, 558],
1262 [ -197, 523, 613, 57], [ 429, 894, 769, 111],
1263 [ 67, 1174, 568, 511], [ 1242, 824, 251, 840],
1264 [ 1419, 1074, 864, 481], [ 924, 1474, 669, 724],
1265 [ 1539, 1879, 654, 1590], [ 445, 337, 1111, 541],
1266 [ 472, 1421, 1264, 1094], [ 794, 735, 1103, 668],
1267 [ 1055, 863, 1192, 1020], [ 778, 1105, 806, 1798],
1268 [ 1052, 1527, 1587, 2151], [ 881, 1552, 1265, 391],
1269 [ 726, 872, 1812, 601], [ 1469, 280, 1008, 616],
1270 [ 1403, 577, 1803, 1244], [ 1650, 1314, 1148, 1072],
1271 [ 1297, 1669, 1911, 1026], [ 2093, 1044, 2115, 1189],
1272 [ 1644, 1961, 2587, 1512], [ 25, -315, -9, -106],
1273 [ 290, -339, 428, -444], [ -68, -783, 735, 772],
1274 [ 245, -555, 468, 47], [ 334, -895, 814, 146],
1275 [ 235, 368, -964, -959], [ -203, 315, -1566, -1217],
1276 [ 801, 17, -276, -354], [ 894, -495, -789, -635],
1277 [ 716, 291, -1189, -357], [ 560, -260, -733, -2],
1278 [ 679, -508, -1429, 211], [ -51, -62, -428, 557],
1279 [ 322, -638, -211, 614], [ -878, -1057, -84, -71],
1280 [ -388, -1415, -167, -318], [ -754, -1574, 214, -539],
1281 [-1419, -2004, -92, -787], [ -47, -856, -347, -255],
1282 [ 23, -1211, -173, 320], [ -658, -487, -893, 353],
1283 [ -783, -1587, -584, 507], [-1420, -859, -378, 441],
1284 [-2095, -1491, -137, 439], [ -321, -1450, -1288, -12],
1285 [ -359, -2113, -553, -8], [ -831, -1918, -1561, 32],
1286 [-1014, -2487, -1359, -939], [ -475, -311, -169, -236],
1287 [ -907, -426, 276, -611], [ -96, -400, 50, -710],
1288 [ -426, -1022, -10, -985], [ -197, -258, -744, -575],
1289 [ -611, -930, -771, -394], [ -267, -776, -612, -939],
1290 [ -256, -1346, -802, -1122], [ -796, -1570, -825, -754],
1291 [ 712, 876, 141, 227], [ 981, 1509, 85, 124],
1292 [ 1462, 1228, 979, -39], [ 1734, 999, 1481, 440],
1293 [ 2293, 1116, 769, 440], [ 2504, 1480, 1241, 356],
1294 [ 2474, 1909, 1558, 810], [ 917, 1134, 607, -134],
1295 [ 509, 1809, 781, -123], [ 1712, 1506, 559, -423],
1296 [ 2037, 2317, 726, -155], [ 3031, 2676, 1203, 331],
1297 [ 3664, 3274, 1768, 531], [ 1610, 1839, 867, 183],
1298 [ 1774, 1972, 1538, 97], [ 1822, 2158, 1282, 659],
1299 [ 2222, 2758, 1818, 900], [ 3251, 2124, 1723, 996],
1300 [ 3633, 2336, 2408, 1453], [ 2923, 3517, 2567, 1318],
1301 ];
1302
1303 const COS_TAB: [i16; 513] = [
1304 16384, 16383, 16379, 16373, 16364, 16353, 16340, 16324,
1305 16305, 16284, 16261, 16235, 16207, 16176, 16143, 16107,
1306 16069, 16029, 15986, 15941, 15893, 15843, 15791, 15736,
1307 15679, 15619, 15557, 15493, 15426, 15357, 15286, 15213,
1308 15137, 15059, 14978, 14896, 14811, 14724, 14635, 14543,
1309 14449, 14354, 14256, 14155, 14053, 13949, 13842, 13733,
1310 13623, 13510, 13395, 13279, 13160, 13039, 12916, 12792,
1311 12665, 12537, 12406, 12274, 12140, 12004, 11866, 11727,
1312 11585, 11442, 11297, 11151, 11003, 10853, 10702, 10549,
1313 10394, 10238, 10080, 9921, 9760, 9598, 9434, 9269,
1314 9102, 8935, 8765, 8595, 8423, 8250, 8076, 7900,
1315 7723, 7545, 7366, 7186, 7005, 6823, 6639, 6455,
1316 6270, 6084, 5897, 5708, 5520, 5330, 5139, 4948,
1317 4756, 4563, 4370, 4176, 3981, 3786, 3590, 3393,
1318 3196, 2999, 2801, 2603, 2404, 2205, 2006, 1806,
1319 1606, 1406, 1205, 1005, 804, 603, 402, 201,
1320 0, -201, -402, -603, -804, -1005, -1205, -1406,
1321 -1606, -1806, -2006, -2205, -2404, -2603, -2801, -2999,
1322 -3196, -3393, -3590, -3786, -3981, -4176, -4370, -4563,
1323 -4756, -4948, -5139, -5330, -5520, -5708, -5897, -6084,
1324 -6270, -6455, -6639, -6823, -7005, -7186, -7366, -7545,
1325 -7723, -7900, -8076, -8250, -8423, -8595, -8765, -8935,
1326 -9102, -9269, -9434, -9598, -9760, -9921, -10080, -10238,
1327 -10394, -10549, -10702, -10853, -11003, -11151, -11297, -11442,
1328 -11585, -11727, -11866, -12004, -12140, -12274, -12406, -12537,
1329 -12665, -12792, -12916, -13039, -13160, -13279, -13395, -13510,
1330 -13623, -13733, -13842, -13949, -14053, -14155, -14256, -14354,
1331 -14449, -14543, -14635, -14724, -14811, -14896, -14978, -15059,
1332 -15137, -15213, -15286, -15357, -15426, -15493, -15557, -15619,
1333 -15679, -15736, -15791, -15843, -15893, -15941, -15986, -16029,
1334 -16069, -16107, -16143, -16176, -16207, -16235, -16261, -16284,
1335 -16305, -16324, -16340, -16353, -16364, -16373, -16379, -16383,
1336 -16384, -16383, -16379, -16373, -16364, -16353, -16340, -16324,
1337 -16305, -16284, -16261, -16235, -16207, -16176, -16143, -16107,
1338 -16069, -16029, -15986, -15941, -15893, -15843, -15791, -15736,
1339 -15679, -15619, -15557, -15493, -15426, -15357, -15286, -15213,
1340 -15137, -15059, -14978, -14896, -14811, -14724, -14635, -14543,
1341 -14449, -14354, -14256, -14155, -14053, -13949, -13842, -13733,
1342 -13623, -13510, -13395, -13279, -13160, -13039, -12916, -12792,
1343 -12665, -12537, -12406, -12274, -12140, -12004, -11866, -11727,
1344 -11585, -11442, -11297, -11151, -11003, -10853, -10702, -10549,
1345 -10394, -10238, -10080, -9921, -9760, -9598, -9434, -9269,
1346 -9102, -8935, -8765, -8595, -8423, -8250, -8076, -7900,
1347 -7723, -7545, -7366, -7186, -7005, -6823, -6639, -6455,
1348 -6270, -6084, -5897, -5708, -5520, -5330, -5139, -4948,
1349 -4756, -4563, -4370, -4176, -3981, -3786, -3590, -3393,
1350 -3196, -2999, -2801, -2603, -2404, -2205, -2006, -1806,
1351 -1606, -1406, -1205, -1005, -804, -603, -402, -201,
1352 0, 201, 402, 603, 804, 1005, 1205, 1406,
1353 1606, 1806, 2006, 2205, 2404, 2603, 2801, 2999,
1354 3196, 3393, 3590, 3786, 3981, 4176, 4370, 4563,
1355 4756, 4948, 5139, 5330, 5520, 5708, 5897, 6084,
1356 6270, 6455, 6639, 6823, 7005, 7186, 7366, 7545,
1357 7723, 7900, 8076, 8250, 8423, 8595, 8765, 8935,
1358 9102, 9269, 9434, 9598, 9760, 9921, 10080, 10238,
1359 10394, 10549, 10702, 10853, 11003, 11151, 11297, 11442,
1360 11585, 11727, 11866, 12004, 12140, 12274, 12406, 12537,
1361 12665, 12792, 12916, 13039, 13160, 13279, 13395, 13510,
1362 13623, 13733, 13842, 13949, 14053, 14155, 14256, 14354,
1363 14449, 14543, 14635, 14724, 14811, 14896, 14978, 15059,
1364 15137, 15213, 15286, 15357, 15426, 15493, 15557, 15619,
1365 15679, 15736, 15791, 15843, 15893, 15941, 15986, 16029,
1366 16069, 16107, 16143, 16176, 16207, 16235, 16261, 16284,
1367 16305, 16324, 16340, 16353, 16364, 16373, 16379, 16383,
1368 16384
1369 ];
1370
1371 const FIXED_CB_GAIN: [i16; 24] = [
1372 1, 2, 3, 4, 6, 9, 13, 18,
1373 26, 38, 55, 80, 115, 166, 240, 348,
1374 502, 726, 1050, 1517, 2193, 3170, 4582, 6623
1375 ];
1376
1377 const ACB_GAIN170: [[i16; 20]; 170] = [
1378 [
1379 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1380 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
1381 ], [
1382 776, 212, 715, 670, 809, -36, -2, -31, -27, -39,
1383 -10, -33, -9, -31, -8, -29, -38, -10, -35, -33
1384 ], [
1385 1296, 1316, -168, -320, -815, -102, -105, -1, -6, -40,
1386 -104, 13, 13, 25, 25, -3, 64, 65, -8, -15
1387 ], [
1388 -589, 680, 2478, 308, -596, -21, -28, -375, -5, -21,
1389 24, 89, -102, 11, -12, -46, -21, 24, 90, 11
1390 ], [
1391 -735, -487, -5, 2948, 468, -33, -14, 0, -530, -13,
1392 -21, 0, 0, 132, 87, 0, 21, 13, 0, -84
1393 ], [
1394 1042, 1730, 1068, 333, 626, -66, -182, -69, -6, -23,
1395 -110, -67, -112, -21, -35, -21, -39, -66, -40, -12
1396 ], [
1397 486, -769, 4074, 2825, -1107, -14, -36, -1013, -487, -74,
1398 22, -120, 191, -83, 132, -702, 32, -52, 275, 191
1399 ], [
1400 1521, -767, -124, 4320, 1026, -141, -35, 0, -1139, -64,
1401 71, 11, -5, -401, 202, 32, -95, 48, 7, -270
1402 ], [
1403 2425, 1267, 3439, -91, -1166, -359, -98, -722, 0, -83,
1404 -187, -509, -266, 13, 7, 19, 172, 90, 244, -6
1405 ], [
1406 -1251, 975, 173, 4039, 2005, -95, -58, -1, -996, -245,
1407 74, 13, -10, 308, -240, -42, 153, -119, -21, -494
1408 ], [
1409 1820, 632, 1322, 2062, 1031, -202, -24, -106, -259, -64,
1410 -70, -146, -51, -229, -79, -166, -114, -39, -83, -129
1411 ], [
1412 -447, 4904, 244, -315, -2038, -12, -1467, -3, -6, -253,
1413 134, 6, -73, -8, 94, 4, -55, 610, 30, -39
1414 ], [
1415 -208, -1102, 463, -448, 5653, -2, -74, -13, -12, -1950,
1416 -14, 5, 31, -5, -30, 12, 71, 380, -159, 154
1417 ], [
1418 4739, 2600, -1864, 856, -1554, -1371, -412, -212, -44, -147,
1419 -752, 539, 295, -247, -135, 97, 449, 246, -176, 81
1420 ], [
1421 1894, 3533, 35, -26, 2145, -219, -762, 0, 0, -280,
1422 -408, -4, -7, 3, 5, 0, -248, -462, -4, 3
1423 ], [
1424 -2699, 1841, 4072, 2443, 1582, -444, -207, -1012, -364, -152,
1425 303, 670, -457, 402, -274, -607, 260, -177, -393, -236
1426 ], [
1427 -844, 3358, 6106, -1059, -537, -43, -688, -2275, -68, -17,
1428 173, 314, -1251, -54, 217, 395, -27, 110, 200, -34
1429 ], [
1430 1251, 1016, 3020, 2210, 1445, -95, -63, -556, -298, -127,
1431 -77, -230, -187, -168, -137, -407, -110, -89, -266, -194
1432 ], [
1433 2099, 2277, 4038, 3533, -2870, -269, -316, -995, -762, -503,
1434 -291, -517, -561, -452, -491, -871, 367, 399, 707, 619
1435 ], [
1436 400, -1114, 8516, 2422, -1117, -9, -75, -4426, -358, -76,
1437 27, -208, 579, -59, 164, -1259, 27, -75, 580, 165
1438 ], [
1439 -4398, -2011, 3912, -2407, 2258, -1180, -247, -934, -353, -311,
1440 -540, 1050, 480, -646, -295, 575, 606, 277, -539, 331
1441 ], [
1442 1767, -1447, 4240, 6160, -757, -190, -127, -1097, -2316, -35,
1443 156, -457, 374, -664, 544, -1594, 81, -66, 195, 284
1444 ], [
1445 1594, -1463, 1035, 6938, 1920, -155, -130, -65, -2938, -225,
1446 142, -100, 92, -675, 619, -438, -186, 171, -121, -813
1447 ], [
1448 -562, 4716, 4085, -591, 2421, -19, -1357, -1018, -21, -357,
1449 162, 140, -1175, -20, 170, 147, 83, -696, -603, 87
1450 ], [
1451 1552, 8778, -935, 354, -1424, -147, -4703, -53, -7, -123,
1452 -831, 88, 501, -33, -189, 20, 134, 763, -81, 30
1453 ], [
1454 4831, -4431, 41, -1479, -2976, -1424, -1198, 0, -133, -540,
1455 1306, -12, 11, 436, -400, 3, 877, -804, 7, -268
1456 ], [
1457 2090, 1192, 1006, 1645, 4853, -266, -86, -61, -165, -1437,
1458 -152, -128, -73, -210, -119, -101, -619, -353, -298, -487
1459 ], [
1460 2386, 5712, 1426, -94, 1350, -347, -1991, -124, 0, -111,
1461 -832, -207, -497, 13, 32, 8, -196, -470, -117, 7
1462 ], [
1463 -1349, 1091, 1659, 8891, 313, -111, -72, -168, -4825, -5,
1464 89, 136, -110, 732, -592, -900, 25, -20, -31, -170
1465 ], [
1466 9980, 916, -381, -808, 88, -6080, -51, -8, -39, 0,
1467 -558, 232, 21, 492, 45, -18, -53, -4, 2, 4
1468 ], [
1469 2338, -1031, -248, 3928, 6484, -333, -64, -3, -942, -2566,
1470 147, 35, -15, -560, 247, 59, -925, 408, 98, -1555
1471 ], [
1472 6166, -1240, -337, 3672, -1277, -2320, -93, -6, -823, -99,
1473 466, 126, -25, -1382, 278, 75, 480, -96, -26, 286
1474 ], [
1475 4377, -132, -2588, 1701, 4865, -1169, -1, -409, -176, -1444,
1476 35, 691, -20, -454, 13, 268, -1299, 39, 768, -505
1477 ], [
1478 2594, 3295, 3944, 1481, 682, -410, -662, -949, -133, -28,
1479 -521, -624, -793, -234, -297, -356, -108, -137, -164, -61
1480 ], [
1481 4151, 624, 815, 4485, 2229, -1052, -23, -40, -1228, -303,
1482 -158, -206, -31, -1136, -170, -223, -565, -84, -111, -610
1483 ], [
1484 -3575, -361, 4924, 2791, 4698, -780, -7, -1480, -475, -1347,
1485 -78, 1074, 108, 609, 61, -839, 1025, 103, -1412, -800
1486 ], [
1487 -2518, 3791, 8623, 315, 2465, -387, -877, -4538, -6, -370,
1488 582, 1325, -1995, 48, -73, -166, 378, -570, -1297, -47
1489 ], [
1490 -691, 2989, 9957, -421, -1142, -29, -545, -6051, -10, -79,
1491 126, 420, -1817, -17, 76, 256, -48, 208, 694, -29
1492 ], [
1493 -1918, 104, -3190, -3410, -4440, -224, 0, -621, -709, -1203,
1494 12, -373, 20, -399, 21, -664, -519, 28, -864, -924
1495 ], [
1496 -3359, -1668, 1854, 6939, 1430, -688, -169, -209, -2939, -124,
1497 -341, 380, 188, 1422, 706, -785, 293, 145, -161, -606
1498 ], [
1499 42, 9706, 3164, -952, 907, 0, -5750, -611, -55, -50,
1500 -25, -8, -1874, 2, 564, 183, -2, -537, -175, 52
1501 ], [
1502 1607, 785, 2862, 4327, 3307, -157, -37, -500, -1143, -667,
1503 -77, -280, -137, -424, -207, -756, -324, -158, -577, -873
1504 ], [
1505 6801, 3416, 2227, 1682, -3217, -2823, -712, -302, -172, -631,
1506 -1418, -924, -464, -698, -350, -228, 1335, 670, 437, 330
1507 ], [
1508 3459, 3898, 364, 7841, -2640, -730, -927, -8, -3753, -425,
1509 -823, -76, -86, -1655, -1865, -174, 557, 628, 58, 1263
1510 ], [
1511 -5902, -3458, -2465, -1886, 4334, -2126, -730, -371, -217, -1146,
1512 -1245, -888, -520, -679, -398, -283, 1561, 915, 652, 499
1513 ], [
1514 -3710, 1133, 7849, 3443, -215, -840, -78, -3760, -723, -2,
1515 256, 1777, -543, 779, -238, -1649, -48, 14, 103, 45
1516 ], [
1517 4132, 2828, 2, -4212, -4116, -1042, -488, 0, -1083, -1034,
1518 -713, 0, 0, 1062, 727, 0, 1038, 710, 0, -1058
1519 ], [
1520 5875, 8496, -1796, 1376, -1786, -2107, -4406, -197, -115, -194,
1521 -3047, 644, 931, -493, -713, 150, 640, 926, -195, 150
1522 ], [
1523 3143, 3483, 3546, -793, 4489, -603, -740, -767, -38, -1230,
1524 -668, -680, -754, 152, 168, 171, -861, -954, -971, 217
1525 ], [
1526 2845, 7965, 3695, -5432, 3978, -494, -3873, -833, -1801, -966,
1527 -1383, -641, -1796, 943, 2641, 1225, -691, -1934, -897, 1319
1528 ], [
1529 1538, 150, 7139, 2049, 3097, -144, -1, -3110, -256, -585,
1530 -14, -670, -65, -192, -18, -892, -290, -28, -1349, -387
1531 ], [
1532 618, 7520, 4729, -238, -3373, -23, -3452, -1365, -3, -694,
1533 -283, -178, -2170, 8, 109, 68, 127, 1548, 973, -49
1534 ], [
1535 2965, -3013, 7912, 7076, -1997, -536, -554, -3821, -3056, -243,
1536 545, -1431, 1455, -1280, 1301, -3417, 361, -367, 964, 862
1537 ], [
1538 2443, -929, -1113, 9677, 4138, -364, -52, -75, -5716, -1045,
1539 138, 166, -63, -1443, 549, 657, -617, 234, 281, -2444
1540 ], [
1541 1966, 3309, 10085, -3399, 2105, -236, -668, -6207, -705, -270,
1542 -397, -1210, -2037, 408, 686, 2092, -252, -425, -1295, 436
1543 ], [
1544 -112, -1368, 8868, 4822, 2048, 0, -114, -4800, -1419, -256,
1545 -9, 61, 740, 33, 402, -2610, 14, 171, -1108, -602
1546 ], [
1547 -2597, 438, -1839, 6229, 7266, -411, -11, -206, -2368, -3223,
1548 69, -291, 49, 987, -166, 699, 1152, -194, 816, -2763
1549 ], [
1550 3454, 553, 9127, 4946, -5596, -728, -18, -5084, -1493, -1911,
1551 -116, -1924, -308, -1042, -166, -2755, 1179, 188, 3117, 1689
1552 ], [
1553 -532, -663, 12262, 2495, -1004, -17, -26, -9177, -380, -61,
1554 -21, 398, 496, 81, 101, -1867, -32, -40, 751, 152
1555 ], [
1556 -2100, 1317, -1509, 11425, 2997, -269, -105, -139, -7967, -548,
1557 168, -193, 121, 1464, -918, 1052, 384, -240, 276, -2090
1558 ], [
1559 1193, -2697, 11259, 5373, -763, -86, -444, -7737, -1762, -35,
1560 196, -819, 1853, -391, 884, -3692, 55, -125, 525, 250
1561 ], [
1562 2405, -471, 11079, 203, 782, -353, -13, -7491, -2, -37,
1563 69, -1626, 318, -29, 5, -137, -114, 22, -529, -9
1564 ], [
1565 -1871, 5685, 11290, -2662, 1353, -213, -1972, -7780, -432, -111,
1566 649, 1289, -3917, -304, 923, 1834, 154, -469, -932, 220
1567 ], [
1568 -3768, 5927, -3093, 5041, 5212, -866, -2144, -584, -1551, -1658,
1569 1363, -711, 1119, 1159, -1824, 951, 1198, -1885, 984, -1603
1570 ], [
1571 -2546, 9502, 5969, -2440, 1928, -395, -5511, -2175, -363, -226,
1572 1477, 927, -3462, -379, 1415, 889, 299, -1118, -702, 287
1573 ], [
1574 -4963, 3568, 4592, 5508, 3451, -1503, -777, -1287, -1851, -727,
1575 1080, 1391, -1000, 1668, -1199, -1543, 1045, -751, -967, -1160
1576 ], [
1577 1745, -2586, 3983, 10899, -1551, -186, -408, -968, -7250, -146,
1578 275, -424, 628, -1161, 1720, -2649, 165, -244, 377, 1032
1579 ], [
1580 867, -456, -727, 3369, 11822, -45, -12, -32, -692, -8531,
1581 24, 38, -20, -178, 93, 149, -625, 329, 525, -2431
1582 ], [
1583 7535, 2422, 1926, 1405, 1599, -3466, -358, -226, -120, -156,
1584 -1114, -886, -284, -646, -207, -165, -735, -236, -188, -137
1585 ], [
1586 1041, -735, -142, 13209, 1515, -66, -33, -1, -10649, -140,
1587 46, 9, -6, -839, 593, 114, -96, 68, 13, -1222
1588 ], [
1589 7950, 6745, -1444, -1008, 2721, -3857, -2777, -127, -62, -452,
1590 -3273, 700, 594, 489, 415, -88, -1320, -1120, 239, 167
1591 ], [
1592 -4754, -1379, 4522, -578, -5733, -1379, -116, -1248, -20, -2006,
1593 -400, 1312, 380, -167, -48, 159, -1663, -482, 1582, -202
1594 ], [
1595 3220, 5978, 5923, 2430, -2689, -633, -2181, -2141, -360, -441,
1596 -1175, -1164, -2161, -477, -886, -878, 528, 981, 972, 398
1597 ], [
1598 377, 1312, 13978, -1470, 677, -8, -105, -11925, -132, -28,
1599 -30, -321, -1119, 33, 117, 1254, -15, -54, -577, 60
1600 ], [
1601 -3435, 6770, 314, -885, 5686, -720, -2797, -6, -47, -1973,
1602 1419, 65, -129, -185, 366, 16, 1192, -2349, -109, 307
1603 ], [
1604 3171, 8774, -2260, 2679, 3069, -613, -4699, -312, -438, -575,
1605 -1698, 437, 1210, -518, -1435, 369, -594, -1643, 423, -501
1606 ], [
1607 5557, 1509, 5407, -125, -7386, -1884, -139, -1784, 0, -3330,
1608 -511, -1834, -498, 42, 11, 41, 2505, 680, 2438, -56
1609 ], [
1610 -2838, 2595, 13228, 271, 1793, -491, -411, -10680, -4, -196,
1611 449, 2291, -2095, 47, -42, -219, 310, -284, -1447, -29
1612 ], [
1613 664, -278, 14966, 951, -711, -26, -4, -13672, -55, -30,
1614 11, -606, 253, -38, 16, -869, 28, -12, 650, 41
1615 ], [
1616 808, 1770, 8658, 5863, -1486, -39, -191, -4576, -2098, -134,
1617 -87, -427, -935, -289, -633, -3098, 73, 160, 785, 531
1618 ], [
1619 3063, 1539, 2000, -542, 9576, -572, -144, -244, -17, -5597,
1620 -287, -374, -188, 101, 51, 66, -1790, -900, -1169, 317
1621 ], [
1622 514, 14083, -323, 896, -891, -16, -12106, -6, -49, -48,
1623 -442, 10, 277, -28, -770, 17, 27, 766, -17, 48
1624 ], [
1625 892, 158, 5237, 11057, -1603, -48, -1, -1674, -7462, -156,
1626 -8, -285, -50, -602, -106, -3534, 87, 15, 512, 1082
1627 ], [
1628 -1612, 2564, -4296, 12526, 5710, -158, -401, -1126, -9576, -1990,
1629 252, -422, 672, 1232, -1960, 3284, 561, -893, 1497, -4365
1630 ], [
1631 4889, -6878, 612, 6109, 4753, -1459, -2887, -22, -2277, -1379,
1632 2052, -182, 257, -1823, 2564, -228, -1418, 1995, -177, -1772
1633 ], [
1634 3053, -506, 2403, 9625, 1322, -569, -15, -352, -5655, -106,
1635 94, -448, 74, -1794, 297, -1412, -246, 40, -194, -777
1636 ], [
1637 -754, 12904, 4480, -2113, 1471, -34, -10163, -1225, -272, -132,
1638 594, 206, -3529, -97, 1664, 577, 67, -1159, -402, 189
1639 ], [
1640 4255, 1476, 5055, 2393, 2912, -1105, -132, -1559, -349, -517,
1641 -383, -1313, -455, -621, -215, -738, -756, -262, -898, -425
1642 ], [
1643 -1371, 535, 1417, 14604, -997, -114, -17, -122, -13017, -60,
1644 44, 118, -46, 1222, -477, -1263, -83, 32, 86, 888
1645 ], [
1646 5368, -1744, 4083, -1236, 3753, -1758, -185, -1017, -93, -860,
1647 571, -1338, 434, 405, -131, 308, -1229, 399, -935, 283
1648 ], [
1649 1588, -3097, 14415, 3699, -1171, -154, -585, -12683, -835, -83,
1650 300, -1397, 2725, -358, 699, -3255, 113, -221, 1030, 264
1651 ], [
1652 212, 7989, 9471, -3344, 2009, -2, -3895, -5475, -682, -246,
1653 -103, -123, -4618, 43, 1630, 1933, -26, -979, -1161, 410
1654 ], [
1655 856, 2294, -627, 6930, 6929, -44, -321, -24, -2931, -2930,
1656 -119, 32, 87, -362, -970, 265, -362, -970, 265, -2931
1657 ], [
1658 2357, -4187, 7162, 7683, 3371, -339, -1070, -3131, -3603, -693,
1659 602, -1030, 1830, -1105, 1963, -3359, -485, 861, -1474, -1581
1660 ], [
1661 350, 4585, 14053, -3819, 1218, -7, -1283, -12054, -890, -90,
1662 -97, -300, -3933, 81, 1068, 3275, -26, -341, -1045, 284
1663 ], [
1664 -3248, 3531, 475, 2137, 11711, -644, -761, -13, -278, -8372,
1665 700, 94, -102, 423, -460, -62, 2322, -2524, -340, -1528
1666 ], [
1667 -3017, 3852, 1725, 8440, 5257, -555, -905, -181, -4348, -1686,
1668 709, 317, -405, 1554, -1984, -889, 968, -1236, -553, -2708
1669 ], [
1670 -909, 3196, 15512, -2528, 1066, -50, -623, -14686, -390, -69,
1671 177, 861, -3026, -140, 493, 2393, 59, -208, -1009, 164
1672 ], [
1673 959, -3370, 9617, 9545, -1761, -56, -693, -5645, -5561, -189,
1674 197, -563, 1978, -558, 1963, -5603, 103, -362, 1034, 1026
1675 ], [
1676 7575, 11796, -4845, 3252, -1703, -3502, -8493, -1433, -645, -177,
1677 -5454, 2240, 3488, -1503, -2341, 961, 787, 1226, -503, 338
1678 ], [
1679 6409, 1722, 1764, -4191, 6015, -2507, -181, -189, -1072, -2208,
1680 -673, -690, -185, 1639, 440, 451, -2353, -632, -647, 1538
1681 ], [
1682 -2420, 12161, 5038, 1286, -2098, -357, -9027, -1549, -100, -268,
1683 1796, 744, -3740, 190, -954, -395, -310, 1557, 645, 164
1684 ], [
1685 -2232, -1341, 7246, 9470, -1977, -304, -109, -3204, -5474, -238,
1686 -182, 987, 593, 1290, 775, -4188, -269, -161, 874, 1143
1687 ], [
1688 1030, 7034, 4231, 1551, 3077, -64, -3019, -1093, -146, -577,
1689 -442, -266, -1816, -97, -666, -400, -193, -1321, -794, -291
1690 ], [
1691 5121, 11835, -477, -1749, 2298, -1601, -8549, -13, -186, -322,
1692 -3699, 149, 344, 546, 1264, -50, -718, -1660, 66, 245
1693 ], [
1694 -3328, 3827, 5921, 9976, -1045, -676, -894, -2140, -6075, -66,
1695 777, 1203, -1383, 2027, -2330, -3605, -212, 244, 377, 636
1696 ], [
1697 3813, 5718, -4666, -3412, 5674, -887, -1995, -1329, -710, -1965,
1698 -1331, 1086, 1628, 794, 1191, -972, -1320, -1980, 1616, 1181
1699 ], [
1700 1348, -3672, 13154, 6938, -1690, -110, -823, -10561, -2938, -174,
1701 302, -1082, 2948, -570, 1555, -5570, 139, -379, 1357, 716
1702 ], [
1703 2151, -3586, 6949, 12131, -1224, -282, -785, -2947, -8982, -91,
1704 470, -912, 1521, -1592, 2655, -5145, 160, -268, 519, 906
1705 ], [
1706 -2889, 9647, 10276, -2728, 995, -509, -5680, -6445, -454, -60,
1707 1701, 1812, -6051, -481, 1606, 1711, 175, -586, -624, 165
1708 ], [
1709 6177, 2184, 555, 1985, 6589, -2329, -291, -18, -240, -2650,
1710 -823, -209, -74, -748, -264, -67, -2484, -878, -223, -798
1711 ], [
1712 -492, 391, 17166, -681, 240, -14, -9, -17987, -28, -3,
1713 11, 515, -410, -20, 16, 713, 7, -5, -252, 10
1714 ], [
1715 12628, 5448, -2630, 3011, -2695, -9733, -1811, -422, -553, -443,
1716 -4199, 2027, 874, -2321, -1001, 483, 2077, 896, -432, 495
1717 ], [
1718 -3628, -534, 3447, 7002, 6751, -803, -17, -725, -2992, -2782,
1719 -118, 763, 112, 1550, 228, -1473, 1495, 220, -1420, -2885
1720 ], [
1721 -5239, 5901, 8107, 3650, 4846, -1675, -2125, -4012, -813, -1433,
1722 1887, 2592, -2920, 1167, -1315, -1806, 1550, -1745, -2398, -1080
1723 ], [
1724 6157, 6678, 4099, -1074, 2348, -2314, -2722, -1025, -70, -336,
1725 -2509, -1540, -1670, 403, 437, 268, -882, -957, -587, 153
1726 ], [
1727 1079, 16099, 242, -881, 1690, -71, -15820, -3, -47, -174,
1728 -1060, -16, -238, 58, 865, 13, -111, -1661, -25, 90
1729 ], [
1730 -278, 227, -1039, 1636, 16945, -4, -3, -65, -163, -17526,
1731 3, -17, 14, 27, -22, 103, 287, -234, 1074, -1693
1732 ], [
1733 15778, -1454, 574, -603, -107, -15195, -129, -20, -22, 0,
1734 1400, -553, 51, 581, -53, 21, 103, -9, 3, -3
1735 ], [
1736 2406, -836, 13224, 7993, -4266, -353, -42, -10673, -3899, -1111,
1737 122, -1942, 674, -1174, 407, -6451, 626, -217, 3443, 2081
1738 ], [
1739 3184, 14368, -3336, 2255, -1801, -619, -12600, -679, -310, -198,
1740 -2793, 648, 2926, -438, -1977, 459, 350, 1580, -366, 247
1741 ], [
1742 -1698, 17076, 2504, -539, -646, -176, -17798, -382, -17, -25,
1743 1770, 259, -2610, -55, 561, 82, -67, 673, 98, -21
1744 ], [
1745 2375, -797, -2696, 14483, 5383, -344, -38, -443, -12803, -1769,
1746 115, 391, -131, -2100, 705, 2384, -780, 262, 886, -4759
1747 ], [
1748 -2691, 2554, -4520, 9573, 10655, -442, -398, -1247, -5594, -6930,
1749 419, -742, 704, 1572, -1492, 2641, 1750, -1661, 2939, -6226
1750 ], [
1751 -4332, -4399, -1657, 4880, 7375, -1145, -1181, -167, -1453, -3319,
1752 -1163, -438, -444, 1290, 1310, 493, 1950, 1980, 745, -2196
1753 ], [
1754 -3498, 7405, 9955, 2693, -2971, -746, -3347, -6049, -442, -538,
1755 1581, 2125, -4499, 575, -1217, -1636, -634, 1342, 1805, 488
1756 ], [
1757 6717, -3792, 7739, 2798, 3489, -2754, -877, -3655, -477, -743,
1758 1554, -3173, 1791, -1147, 647, -1321, -1430, 807, -1648, -595
1759 ], [
1760 5263, 9770, 3463, 1069, -3971, -1690, -5826, -732, -69, -962,
1761 -3138, -1112, -2065, -343, -637, -226, 1275, 2368, 839, 259
1762 ], [
1763 1243, -2634, 16772, 1871, 332, -94, -423, -17169, -213, -6,
1764 199, -1273, 2696, -142, 300, -1915, -25, 53, -339, -37
1765 ], [
1766 2691, 2836, 3105, 5711, 4817, -442, -491, -588, -1991, -1416,
1767 -465, -510, -537, -938, -988, -1082, -791, -834, -913, -1679
1768 ], [
1769 4366, 2944, 7210, 3627, 1161, -1163, -529, -3172, -803, -82,
1770 -784, -1921, -1295, -966, -651, -1596, -309, -208, -511, -257
1771 ], [
1772 13888, 3951, -671, -2305, 3354, -11773, -953, -27, -324, -686,
1773 -3349, 569, 161, 1954, 556, -94, -2843, -809, 137, 472
1774 ], [
1775 7053, 5847, 2929, 8378, -4794, -3036, -2086, -523, -4284, -1403,
1776 -2517, -1261, -1045, -3607, -2990, -1498, 2064, 1711, 857, 2451
1777 ], [
1778 -2191, 12838, 9182, -3915, 1617, -293, -10059, -5146, -935, -159,
1779 1717, 1228, -7195, -523, 3068, 2194, 216, -1267, -906, 386
1780 ], [
1781 -4881, 13114, 5767, -435, 4155, -1454, -10498, -2030, -11, -1054,
1782 3907, 1718, -4616, -129, 348, 153, 1238, -3326, -1462, 110
1783 ], [
1784 7843, -1250, 210, 7106, -5203, -3754, -95, -2, -3082, -1652,
1785 598, -100, 16, -3402, 542, -91, 2491, -397, 66, 2257
1786 ], [
1787 -2463, 8168, 14551, -3908, 1828, -370, -4072, -12923, -932, -204,
1788 1228, 2188, -7254, -587, 1948, 3471, 274, -911, -1623, 436
1789 ], [
1790 -1579, 347, -272, -2735, 16031, -152, -7, -4, -456, -15686,
1791 33, -26, 5, -263, 58, -45, 1545, -340, 266, 2676
1792 ], [
1793 -6327, 1328, 5093, -5079, 7617, -2443, -107, -1583, -1574, -3541,
1794 513, 1967, -413, -1961, 411, 1578, 2941, -617, -2367, 2361
1795 ], [
1796 3286, -4509, 11306, 11025, -2623, -659, -1241, -7802, -7419, -420,
1797 904, -2267, 3112, -2211, 3034, -7608, 526, -722, 1810, 1765
1798 ], [
1799 5567, 17853, -3754, 1166, -519, -1892, -19455, -860, -83, -16,
1800 -6067, 1275, 4090, -396, -1271, 267, 176, 566, -119, 37
1801 ], [
1802 -2136, -424, 15292, 5108, -1648, -278, -10, -14273, -1593, -165,
1803 -55, 1993, 396, 666, 132, -4768, -214, -42, 1538, 514
1804 ], [
1805 2267, -3297, 2549, 16563, -791, -313, -663, -396, -16745, -38,
1806 456, -352, 513, -2291, 3333, -2576, 109, -159, 123, 799
1807 ], [
1808 3655, 1899, -3364, 6279, 12510, -815, -220, -690, -2406, -9552,
1809 -423, 750, 390, -1400, -728, 1289, -2791, -1450, 2568, -4794
1810 ], [
1811 8052, 2285, -6193, 5138, 6003, -3957, -318, -2341, -1611, -2199,
1812 -1123, 3044, 864, -2525, -716, 1942, -2950, -837, 2269, -1882
1813 ], [
1814 -386, -2291, 7679, 15387, -2723, -9, -320, -3599, -14452, -452,
1815 -54, 181, 1074, 362, 2152, -7212, -64, -380, 1276, 2557
1816 ], [
1817 2777, -1173, 3984, 13079, 2508, -470, -84, -969, -10440, -384,
1818 198, -675, 285, -2217, 936, -3180, -425, 179, -610, -2002
1819 ], [
1820 -1879, 1771, -2684, 16705, 1833, -215, -191, -439, -17032, -205,
1821 203, -308, 290, 1916, -1805, 2736, 210, -198, 300, -1869
1822 ], [
1823 1052, 4495, 15519, 1467, -4032, -67, -1233, -14700, -131, -992,
1824 -288, -997, -4257, -94, -402, -1389, 259, 1106, 3819, 361
1825 ], [
1826 3010, 2544, 6969, 7559, 1996, -553, -395, -2964, -3487, -243,
1827 -467, -1280, -1082, -1388, -1174, -3215, -366, -310, -849, -921
1828 ], [
1829 -5209, -1867, 8713, 10351, 1549, -1656, -212, -4634, -6540, -146,
1830 -593, 2770, 993, 3291, 1180, -5505, 492, 176, -824, -979
1831 ], [
1832 -4314, 8513, 913, 7547, -2723, -1135, -4423, -50, -3476, -452,
1833 2241, 240, -474, 1987, -3921, -420, -717, 1415, 151, 1254
1834 ], [
1835 12929, -1219, 2448, 1757, 6303, -10204, -90, -365, -188, -2425,
1836 962, -1932, 182, -1386, 130, -262, -4974, 469, -941, -676
1837 ], [
1838 6465, 4132, 3167, 3160, 5697, -2551, -1042, -612, -609, -1981,
1839 -1630, -1249, -798, -1247, -797, -611, -2248, -1437, -1101, -1099
1840 ], [
1841 -3636, 4859, 18914, -1335, 810, -807, -1441, -21836, -108, -40,
1842 1078, 4198, -5609, -296, 396, 1541, 179, -240, -936, 66
1843 ], [
1844 8844, 7864, 654, -4063, -5680, -4774, -3774, -26, -1007, -1969,
1845 -4245, -353, -314, 2193, 1950, 162, 3066, 2726, 226, -1408
1846 ], [
1847 1859, 2634, 9228, 996, 9464, -211, -423, -5197, -60, -5467,
1848 -299, -1047, -1483, -113, -160, -561, -1074, -1521, -5330, -575
1849 ], [
1850 2949, 12260, 10290, -497, -3943, -530, -9174, -6463, -15, -949,
1851 -2206, -1852, -7700, 89, 372, 312, 709, 2950, 2476, -119
1852 ], [
1853 -2903, 1552, 14867, 9970, -496, -514, -147, -13491, -6068, -15,
1854 275, 2634, -1408, 1766, -944, -9047, -87, 47, 450, 302
1855 ], [
1856 3243, 8234, 7586, 3373, 2151, -642, -4138, -3512, -694, -282,
1857 -1630, -1501, -3812, -667, -1695, -1561, -425, -1081, -996, -442
1858 ], [
1859 -9631, 60, 3501, 5359, 10150, -5662, 0, -748, -1752, -6288,
1860 35, 2058, -12, 3150, -19, -1145, 5967, -37, -2169, -3320
1861 ], [
1862 -6874, -2553, -5446, -2195, -7841, -2884, -397, -1810, -294, -3753,
1863 -1071, -2285, -848, -921, -342, -729, -3290, -1221, -2606, -1050
1864 ], [
1865 -3413, -1141, 4630, 13612, 7897, -711, -79, -1308, -11310, -3806,
1866 -237, 964, 322, 2836, 948, -3847, 1645, 550, -2231, -6561
1867 ], [
1868 4410, -5678, 8006, -3992, 3811, -1187, -1968, -3912, -973, -886,
1869 1528, -2155, 2775, 1074, -1383, 1951, -1025, 1321, -1862, 928
1870 ], [
1871 5659, 11535, 2203, -452, 7169, -1954, -8121, -296, -12, -3137,
1872 -3984, -761, -1551, 156, 318, 60, -2476, -5048, -964, 197
1873 ], [
1874 2914, -2914, 3485, -3965, 13675, -518, -518, -741, -959, -11414,
1875 518, -620, 620, 705, -705, 843, -2433, 2432, -2909, 3310
1876 ], [
1877 7843, 1907, 1022, 8882, 7972, -3755, -222, -63, -4815, -3879,
1878 -913, -489, -119, -4252, -1034, -554, -3816, -928, -497, -4322
1879 ], [
1880 13807, 9531, 1436, 1612, 1779, -11636, -5544, -125, -158, -193,
1881 -8032, -1210, -835, -1358, -938, -141, -1499, -1035, -156, -175
1882 ], [
1883 13620, -5337, 5450, -2263, 1723, -11322, -1738, -1813, -312, -181,
1884 4436, -4531, 1775, 1881, -737, 752, -1432, 561, -573, 238
1885 ], [
1886 5297, 8374, 8872, 7694, 6538, -1712, -4280, -4804, -3613, -2609,
1887 -2707, -2868, -4534, -2487, -3932, -4166, -2113, -3341, -3540, -3070
1888 ]
1889 ];
1890 const ACB_GAIN85: [[i16; 20]; 85] = [
1891 [
1892 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1893 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
1894 ], [
1895 800, 1496, 167, -256, -338, -39, -136, -1, -4, -6,
1896 -73, -8, -15, 12, 23, 2, 16, 30, 3, -5
1897 ], [
1898 -462, -686, 493, 2575, 311, -13, -28, -14, -404, -5,
1899 -19, 13, 20, 72, 107, -77, 8, 13, -9, -48
1900 ], [
1901 1483, 144, 784, 928, 1243, -134, -1, -37, -52, -94,
1902 -13, -71, -6, -84, -8, -44, -112, -10, -59, -70
1903 ], [
1904 -77, 275, 3522, 1056, -1254, 0, -4, -757, -68, -95,
1905 1, 16, -59, 4, -17, -227, -5, 21, 269, 80
1906 ], [
1907 -125, -40, -264, 381, 5027, 0, 0, -4, -8, -1542,
1908 0, -2, 0, 2, 0, 6, 38, 12, 81, -117
1909 ], [
1910 138, 332, 2215, 2574, 1339, -1, -6, -299, -404, -109,
1911 -2, -18, -44, -21, -52, -348, -11, -27, -181, -210
1912 ], [
1913 3685, 2883, -887, 866, -1639, -828, -507, -48, -45, -164,
1914 -648, 199, 156, -194, -152, 46, 368, 288, -88, 86
1915 ], [
1916 1396, 2146, 2235, 345, 942, -118, -281, -305, -7, -54,
1917 -182, -190, -292, -29, -45, -47, -80, -123, -128, -19
1918 ], [
1919 13, 4475, 3549, -804, -655, 0, -1222, -768, -39, -26,
1920 -3, -2, -969, 0, 219, 174, 0, 179, 141, -32
1921 ], [
1922 -724, 254, 242, 6049, 2462, -32, -3, -3, -2233, -370,
1923 11, 10, -3, 267, -94, -89, 108, -38, -36, -909
1924 ], [
1925 626, -1713, 6121, 4561, -1061, -23, -179, -2287, -1270, -68,
1926 65, -233, 640, -174, 477, -1704, 40, -111, 396, 295
1927 ], [
1928 -350, 1391, 7985, 511, -405, -7, -118, -3892, -15, -10,
1929 29, 170, -678, 10, -43, -249, -8, 34, 197, 12
1930 ], [
1931 3144, -529, 608, 2530, 3878, -603, -17, -22, -390, -918,
1932 101, -116, 19, -485, 81, -93, -744, 125, -144, -599
1933 ], [
1934 2589, -689, 3045, 5603, -404, -409, -29, -566, -1916, -10,
1935 108, -481, 128, -885, 235, -1041, 63, -17, 75, 138
1936 ], [
1937 3107, 513, 1374, -3594, -4922, -589, -16, -115, -788, -1478,
1938 -97, -260, -43, 681, 112, 301, 933, 154, 413, -1079
1939 ], [
1940 2468, 6010, 1107, -390, 1961, -372, -2204, -74, -9, -234,
1941 -905, -166, -406, 58, 143, 26, -295, -719, -132, 46
1942 ], [
1943 4773, 2766, 2368, 4862, -4044, -1390, -467, -342, -1443, -998,
1944 -806, -690, -399, -1416, -821, -702, 1178, 682, 584, 1200
1945 ], [
1946 1665, -1879, 1443, 1701, 8562, -169, -215, -127, -176, -4475,
1947 190, -146, 165, -172, 195, -149, -870, 982, -754, -889
1948 ], [
1949 2716, 9011, -1007, 755, -1785, -450, -4956, -61, -34, -194,
1950 -1493, 167, 554, -125, -415, 46, 296, 982, -109, 82
1951 ], [
1952 -2727, 7548, 1285, 938, 3420, -453, -3478, -100, -53, -714,
1953 1256, 213, -592, 156, -432, -73, 569, -1576, -268, -196
1954 ], [
1955 3677, 882, 4050, 1202, 2323, -825, -47, -1001, -88, -329,
1956 -198, -909, -218, -269, -64, -297, -521, -125, -574, -170
1957 ], [
1958 2046, -753, 122, 10102, 603, -255, -34, 0, -6229, -22,
1959 94, -15, 5, -1261, 464, -75, -75, 27, -4, -372
1960 ], [
1961 449, -1815, 10690, 3870, -527, -12, -201, -6976, -914, -16,
1962 49, -293, 1184, -106, 428, -2525, 14, -58, 344, 124
1963 ], [
1964 -941, 2352, 5049, 3650, 2637, -54, -337, -1556, -813, -424,
1965 135, 290, -725, 209, -524, -1125, 151, -378, -812, -587
1966 ], [
1967 -1879, 796, 3117, 9569, -404, -215, -38, -593, -5589, -9,
1968 91, 357, -151, 1097, -464, -1821, -46, 19, 76, 236
1969 ], [
1970 -1715, 2043, -2096, 9946, 4001, -179, -254, -268, -6038, -977,
1971 213, -219, 261, 1041, -1240, 1272, 418, -498, 511, -2429
1972 ], [
1973 -5772, -618, -3921, 284, -3155, -2033, -23, -938, -4, -607,
1974 -218, -1381, -148, 100, 10, 68, -1111, -119, -755, 54
1975 ], [
1976 382, 4748, 8003, -2064, 2198, -8, -1376, -3909, -260, -294,
1977 -110, -186, -2319, 48, 598, 1008, -51, -637, -1073, 277
1978 ], [
1979 -867, 3015, 11926, -1675, 947, -45, -555, -8681, -171, -54,
1980 159, 631, -2195, -88, 308, 1219, 50, -174, -690, 96
1981 ], [
1982 -4933, -432, 6757, 3771, 1352, -1485, -11, -2786, -867, -111,
1983 -130, 2034, 178, 1135, 99, -1555, 407, 35, -557, -311
1984 ], [
1985 152, 9726, 4231, -1928, 1490, -1, -5774, -1092, -226, -135,
1986 -90, -39, -2511, 17, 1144, 498, -13, -884, -384, 175
1987 ], [
1988 2512, 193, 9033, 5361, -3148, -385, -2, -4980, -1754, -605,
1989 -29, -1385, -106, -822, -63, -2956, 482, 37, 1735, 1030
1990 ], [
1991 8464, 2844, 12, 549, 2132, -4373, -493, 0, -18, -277,
1992 -1469, -6, -2, -284, -95, 0, -1101, -370, -1, -71
1993 ], [
1994 2141, -2602, 7166, 9046, -1350, -279, -413, -3134, -4994, -111,
1995 340, -936, 1138, -1182, 1436, -3957, 176, -214, 590, 745
1996 ], [
1997 -244, 278, 13307, 1227, -161, -3, -4, -10808, -91, -1,
1998 4, 198, -226, 18, -20, -997, -2, 2, 131, 12
1999 ], [
2000 -1947, 8217, 6269, 917, -2559, -231, -4121, -2399, -51, -399,
2001 976, 745, -3144, 108, -460, -350, -304, 1283, 979, 143
2002 ], [
2003 -1810, 2061, -2781, 6056, 10058, -200, -259, -472, -2238, -6174,
2004 227, -307, 349, 669, -761, 1028, 1111, -1265, 1707, -3717
2005 ], [
2006 7827, 9161, -3409, 2473, -1510, -3739, -5122, -709, -373, -139,
2007 -4376, 1628, 1906, -1181, -1382, 514, 721, 844, -314, 228
2008 ], [
2009 -1430, 8313, 9541, -2955, 1626, -124, -4218, -5556, -533, -161,
2010 725, 832, -4841, -257, 1499, 1721, 142, -825, -947, 293
2011 ], [
2012 2819, -4247, 5391, 8673, 2756, -485, -1101, -1774, -4591, -463,
2013 730, -927, 1397, -1492, 2248, -2854, -474, 714, -907, -1459
2014 ], [
2015 141, 14552, 690, 257, -112, -1, -12926, -29, -4, 0,
2016 -125, -5, -613, -2, -228, -10, 0, 99, 4, 1
2017 ], [
2018 11938, -1859, 1806, -962, -884, -8699, -211, -199, -56, -47,
2019 1355, -1316, 205, 701, -109, 106, 644, -100, 97, -51
2020 ], [
2021 3728, 1982, 2264, 4584, 3131, -848, -239, -312, -1282, -598,
2022 -451, -515, -273, -1043, -554, -633, -712, -378, -432, -876
2023 ], [
2024 -1181, 766, 720, 14303, -216, -85, -35, -31, -12486, -2,
2025 55, 51, -33, 1031, -668, -628, -15, 10, 9, 189
2026 ], [
2027 -4385, 4826, 10112, 1569, 3388, -1173, -1421, -6242, -150, -700,
2028 1291, 2706, -2979, 420, -462, -969, 906, -998, -2091, -324
2029 ], [
2030 -448, 1932, 15591, -1842, 657, -12, -227, -14837, -207, -26,
2031 52, 427, -1838, -50, 217, 1753, 18, -77, -626, 74
2032 ], [
2033 -4141, 1844, 3962, 5517, 6220, -1046, -207, -958, -1858, -2361,
2034 466, 1001, -446, 1394, -621, -1334, 1572, -700, -1504, -2094
2035 ], [
2036 729, -2299, 14755, 3657, -952, -32, -322, -13288, -816, -55,
2037 102, -656, 2071, -162, 513, -3294, 42, -133, 857, 212
2038 ], [
2039 -1385, 5801, 13339, -3137, 1344, -117, -2054, -10861, -600, -110,
2040 490, 1127, -4723, -265, 1111, 2554, 113, -476, -1094, 257
2041 ], [
2042 4710, 9661, 1073, -2467, 3274, -1354, -5697, -70, -371, -654,
2043 -2777, -308, -633, 709, 1455, 161, -941, -1930, -214, 493
2044 ], [
2045 1843, -3624, 12422, 6898, -1559, -207, -802, -9419, -2904, -148,
2046 407, -1397, 2748, -775, 1526, -5230, 175, -344, 1182, 656
2047 ], [
2048 1433, 2394, 2507, 1380, 8780, -125, -349, -383, -116, -4705,
2049 -209, -219, -366, -120, -201, -211, -768, -1283, -1343, -740
2050 ], [
2051 -1712, 12915, 5883, -2197, 991, -179, -10181, -2112, -294, -60,
2052 1350, 615, -4638, -229, 1732, 789, 103, -781, -356, 133
2053 ], [
2054 15072, 2158, -1245, 910, -496, -13865, -284, -94, -50, -15,
2055 -1986, 1145, 164, -837, -119, 69, 456, 65, -37, 27
2056 ], [
2057 4655, 7319, 4916, 586, -3381, -1322, -3270, -1475, -20, -697,
2058 -2079, -1396, -2196, -166, -261, -175, 960, 1510, 1014, 120
2059 ], [
2060 1191, -2140, 5120, 13498, -1418, -86, -279, -1600, -11121, -122,
2061 155, -372, 669, -981, 1763, -4218, 103, -185, 443, 1168
2062 ], [
2063 -1530, -817, 8191, 9632, -1452, -143, -40, -4095, -5663, -128,
2064 -76, 765, 408, 900, 480, -4815, -135, -72, 726, 854
2065 ], [
2066 -3236, 607, 1696, -2106, 11485, -639, -22, -175, -270, -8051,
2067 119, 335, -62, -416, 78, 218, 2268, -425, -1189, 1476
2068 ], [
2069 3203, -1903, -837, 9679, 7057, -626, -221, -42, -5718, -3039,
2070 372, 163, -97, -1892, 1124, 494, -1380, 819, 360, -4169
2071 ], [
2072 213, -655, 17015, 620, -384, -2, -26, -17671, -23, -9,
2073 8, -221, 681, -8, 24, -644, 5, -15, 399, 14
2074 ], [
2075 5088, 35, -3339, 3726, 8488, -1580, 0, -680, -847, -4397,
2076 -10, 1037, 7, -1157, -8, 759, -2636, -18, 1730, -1930
2077 ], [
2078 -988, 1454, -2688, 15039, 2682, -59, -129, -441, -13805, -439,
2079 87, -162, 238, 907, -1335, 2467, 161, -238, 440, -2462
2080 ], [
2081 -4865, -2842, -53, 5495, 6523, -1445, -493, 0, -1843, -2597,
2082 -844, -16, -9, 1632, 953, 18, 1937, 1131, 21, -2188
2083 ], [
2084 3076, 15069, -2914, 1810, -971, -577, -13860, -518, -200, -57,
2085 -2829, 547, 2680, -339, -1665, 322, 182, 893, -172, 107
2086 ], [
2087 1311, 5355, 11054, 2299, -3654, -105, -1750, -7458, -322, -814,
2088 -428, -885, -3613, -184, -751, -1551, 292, 1194, 2465, 512
2089 ], [
2090 4035, 5619, 4618, 1815, 1912, -994, -1927, -1301, -201, -223,
2091 -1384, -1137, -1583, -447, -622, -511, -471, -656, -539, -211
2092 ], [
2093 -2131, 2754, -4501, 12879, 7432, -277, -463, -1236, -10124, -3371,
2094 358, -585, 756, 1675, -2165, 3538, 967, -1249, 2042, -5842
2095 ], [
2096 5618, -515, 3219, -4149, 4857, -1926, -16, -632, -1050, -1440,
2097 176, -1104, 101, 1422, -130, 815, -1666, 152, -954, 1230
2098 ], [
2099 1838, -1709, 1139, 16867, 716, -206, -178, -79, -17366, -31,
2100 191, -127, 118, -1892, 1759, -1173, -80, 74, -49, -737
2101 ], [
2102 1978, -3845, 10050, 11854, -2492, -238, -902, -6164, -8576, -379,
2103 464, -1213, 2358, -1431, 2782, -7271, 301, -585, 1529, 1803
2104 ], [
2105 -2600, 11246, 11289, -3647, 1463, -412, -7720, -7778, -812, -130,
2106 1784, 1791, -7749, -578, 2504, 2513, 232, -1004, -1008, 325
2107 ], [
2108 3442, 907, 2725, 8970, 3638, -723, -50, -453, -4911, -808,
2109 -190, -572, -150, -1884, -496, -1492, -764, -201, -605, -1992
2110 ], [
2111 -126, 17498, 3481, -2003, 1090, 0, -18689, -739, -244, -72,
2112 135, 26, -3717, -15, 2139, 425, 8, -1165, -231, 133
2113 ], [
2114 -1814, 1048, -2164, 4070, 16272, -200, -67, -285, -1011, -16160,
2115 116, -239, 138, 450, -260, 537, 1801, -1041, 2149, -4042
2116 ], [
2117 9354, 12580, -1883, 962, -617, -5341, -9660, -216, -56, -23,
2118 -7183, 1075, 1446, -549, -738, 110, 352, 474, -71, 36
2119 ], [
2120 1708, 4199, 7387, 6335, 1003, -178, -1076, -3330, -2449, -61,
2121 -437, -770, -1893, -660, -1623, -2856, -104, -257, -452, -388
2122 ], [
2123 -2624, 5623, 17310, -2353, 592, -420, -1930, -18288, -338, -21,
2124 900, 2772, -5941, -376, 807, 2486, 94, -203, -625, 85
2125 ], [
2126 1211, -850, 1193, -1926, 15992, -89, -44, -86, -226, -15609,
2127 62, -88, 61, 142, -100, 140, -1182, 830, -1165, 1880
2128 ], [
2129 3983, -2054, 11506, -19, 3622, -968, -257, -8080, 0, -801,
2130 499, -2797, 1442, 4, -2, 13, -880, 454, -2544, 4
2131 ], [
2132 -786, -1354, 16092, 7246, -1665, -37, -111, -15805, -3205, -169,
2133 -65, 772, 1330, 348, 599, -7117, -80, -137, 1636, 736
2134 ], [
2135 -4316, -511, 6674, 11665, 4633, -1137, -15, -2719, -8305, -1310,
2136 -134, 1758, 208, 3073, 364, -4752, 1220, 144, -1887, -3299
2137 ], [
2138 7912, 4557, 1937, 1885, 7037, -3821, -1267, -229, -216, -3022,
2139 -2200, -935, -538, -910, -524, -222, -3398, -1957, -832, -809
2140 ], [
2141 3434, 2967, 5867, 8196, 8766, -720, -537, -2101, -4100, -4690,
2142 -622, -1230, -1062, -1718, -1484, -2935, -1837, -1588, -3139, -4385
2143 ], [
2144 5881, 9176, 8119, 3934, 3355, -2111, -5139, -4023, -944, -687,
2145 -3294, -2914, -4547, -1412, -2203, -1949, -1204, -1879, -1662, -805
2146 ]
2147 ];
2148
2149 const PULSE_POSITIONS: [[u32; SUBFRAME_LEN/2]; MAX_PULSES] = [
2150 [
2151 118755, 98280, 80730, 65780, 53130, 42504, 33649, 26334, 20349, 15504,
2152 11628, 8568, 6188, 4368, 3003, 2002, 1287, 792, 462, 252,
2153 126, 56, 21, 6, 1, 0, 0, 0, 0, 0
2154 ], [
2155 23751, 20475, 17550, 14950, 12650, 10626, 8855, 7315, 5985, 4845,
2156 3876, 3060, 2380, 1820, 1365, 1001, 715, 495, 330, 210,
2157 126, 70, 35, 15, 5, 1, 0, 0, 0, 0
2158 ], [
2159 3654, 3276, 2925, 2600, 2300, 2024, 1771, 1540, 1330, 1140,
2160 969, 816, 680, 560, 455, 364, 286, 220, 165, 120,
2161 84, 56, 35, 20, 10, 4, 1, 0, 0, 0
2162 ], [
2163 406, 378, 351, 325, 300, 276, 253, 231, 210, 190,
2164 171, 153, 136, 120, 105, 91, 78, 66, 55, 45,
2165 36, 28, 21, 15, 10, 6, 3, 1, 0, 0
2166 ], [
2167 29, 28, 27, 26, 25, 24, 23, 22, 21, 20,
2168 19, 18, 17, 16, 15, 14, 13, 12, 11, 10,
2169 9, 8, 7, 6, 5, 4, 3, 2, 1, 0
2170 ], [
2171 1; SUBFRAME_LEN/2
2172 ]
2173 ];
2174
2175 const PITCH_CONTRIBUTION: [i32; 340] = [
2176 60, 0, 0, 2489, 60, 0, 0, 5217, 1, 6171, 0, 3953, 0, 10364, 1, 9357,
2177 -1, 8843, 1, 9396, 0, 5794, -1, 10816, 2, 11606, -2, 12072, 0, 8616, 1, 12170,
2178 0, 14440, 0, 7787, -1, 13721, 0, 18205, 0, 14471, 0, 15807, 1, 15275, 0, 13480,
2179 -1, 18375, -1, 0, 1, 11194, -1, 13010, 1, 18836, -2, 20354, 1, 16233, -1, 0,
2180 60, 0, 0, 12130, 0, 13385, 1, 17834, 1, 20875, 0, 21996, 1, 0, 1, 18277,
2181 -1, 21321, 1, 13738, -1, 19094, -1, 20387, -1, 0, 0, 21008, 60, 0, -2, 22807,
2182 0, 15900, 1, 0, 0, 17989, -1, 22259, 1, 24395, 1, 23138, 0, 23948, 1, 22997,
2183 2, 22604, -1, 25942, 0, 26246, 1, 25321, 0, 26423, 0, 24061, 0, 27247, 60, 0,
2184 -1, 25572, 1, 23918, 1, 25930, 2, 26408, -1, 19049, 1, 27357, -1, 24538, 60, 0,
2185 -1, 25093, 0, 28549, 1, 0, 0, 22793, -1, 25659, 0, 29377, 0, 30276, 0, 26198,
2186 1, 22521, -1, 28919, 0, 27384, 1, 30162, -1, 0, 0, 24237, -1, 30062, 0, 21763,
2187 1, 30917, 60, 0, 0, 31284, 0, 29433, 1, 26821, 1, 28655, 0, 31327, 2, 30799,
2188 1, 31389, 0, 32322, 1, 31760, -2, 31830, 0, 26936, -1, 31180, 1, 30875, 0, 27873,
2189 -1, 30429, 1, 31050, 0, 0, 0, 31912, 1, 31611, 0, 31565, 0, 25557, 0, 31357,
2190 60, 0, 1, 29536, 1, 28985, -1, 26984, -1, 31587, 2, 30836, -2, 31133, 0, 30243,
2191 -1, 30742, -1, 32090, 60, 0, 2, 30902, 60, 0, 0, 30027, 0, 29042, 60, 0,
2192 0, 31756, 0, 24553, 0, 25636, -2, 30501, 60, 0, -1, 29617, 0, 30649, 60, 0,
2193 0, 29274, 2, 30415, 0, 27480, 0, 31213, -1, 28147, 0, 30600, 1, 31652, 2, 29068,
2194 60, 0, 1, 28571, 1, 28730, 1, 31422, 0, 28257, 0, 24797, 60, 0, 0, 0,
2195 60, 0, 0, 22105, 0, 27852, 60, 0, 60, 0, -1, 24214, 0, 24642, 0, 23305,
2196 60, 0, 60, 0, 1, 22883, 0, 21601, 60, 0, 2, 25650, 60, 0, -2, 31253,
2197 -2, 25144, 0, 17998
2198 ];
2199
2200 const POSTFILTER_COEFFS: [[i16; LPC_ORDER]; 2] = [
2201 [ 21299, 13844, 8999, 5849, 3802, 2471, 1606, 1044, 679, 441 ],
2202 [ 24576, 18432, 13824, 10368, 7776, 5832, 4374, 3281, 2460, 1845 ]
2203 ];
NihAV is a research multimedia framework