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