[FFmpeg-cvslog] r26051 - in trunk: Changelog configure doc/general.texi libavcodec/Makefile libavcodec/allcodecs.c libavcodec/amrwbdata.h libavcodec/amrwbdec.c
rbultje
subversion
Sat Dec 18 04:03:18 CET 2010
Author: rbultje
Date: Sat Dec 18 04:03:18 2010
New Revision: 26051
Log:
AMR-WB decoder, written as part of Google Summer of Code 2010 by Marcelo
Galv?o P?voa <marspeoplester gmail com>, mentored by Robert Swain <robert
dot swain gmail com>.
Added:
trunk/libavcodec/amrwbdata.h
trunk/libavcodec/amrwbdec.c
Modified:
trunk/Changelog
trunk/configure
trunk/doc/general.texi
trunk/libavcodec/Makefile
trunk/libavcodec/allcodecs.c
Modified: trunk/Changelog
==============================================================================
--- trunk/Changelog Sat Dec 18 00:42:56 2010 (r26050)
+++ trunk/Changelog Sat Dec 18 04:03:18 2010 (r26051)
@@ -65,6 +65,7 @@ version <next>:
- RTP depacketization of QCELP
- FLAC parser added
- gradfun filter added
+- AMR-WB decoder
version 0.6:
Modified: trunk/configure
==============================================================================
--- trunk/configure Sat Dec 18 00:42:56 2010 (r26050)
+++ trunk/configure Sat Dec 18 04:03:18 2010 (r26051)
@@ -1194,6 +1194,7 @@ aac_latm_decoder_select="aac_decoder aac
ac3_decoder_select="mdct ac3_parser"
alac_encoder_select="lpc"
amrnb_decoder_select="lsp"
+amrwb_decoder_select="lsp"
atrac1_decoder_select="mdct"
atrac3_decoder_select="mdct"
binkaudio_dct_decoder_select="mdct rdft dct"
Modified: trunk/doc/general.texi
==============================================================================
--- trunk/doc/general.texi Sat Dec 18 00:42:56 2010 (r26050)
+++ trunk/doc/general.texi Sat Dec 18 04:03:18 2010 (r26051)
@@ -574,8 +574,7 @@ following image formats are supported:
@item ADPCM Yamaha @tab X @tab X
@item AMR-NB @tab E @tab X
@tab encoding supported through external library libopencore-amrnb
- at item AMR-WB @tab @tab E
- @tab decoding supported through external library libopencore-amrwb
+ at item AMR-WB @tab @tab X
@item Apple lossless audio @tab X @tab X
@tab QuickTime fourcc 'alac'
@item Atrac 1 @tab @tab X
Modified: trunk/libavcodec/Makefile
==============================================================================
--- trunk/libavcodec/Makefile Sat Dec 18 00:42:56 2010 (r26050)
+++ trunk/libavcodec/Makefile Sat Dec 18 04:03:18 2010 (r26051)
@@ -62,6 +62,10 @@ OBJS-$(CONFIG_AMRNB_DECODER) +
celp_math.o acelp_filters.o \
acelp_vectors.o \
acelp_pitch_delay.o
+OBJS-$(CONFIG_AMRWB_DECODER) += amrwbdec.o celp_filters.o \
+ celp_math.o acelp_filters.o \
+ acelp_vectors.o \
+ acelp_pitch_delay.o lsp.o
OBJS-$(CONFIG_AMV_DECODER) += sp5xdec.o mjpegdec.o mjpeg.o
OBJS-$(CONFIG_ANM_DECODER) += anm.o
OBJS-$(CONFIG_ANSI_DECODER) += ansi.o cga_data.o
Modified: trunk/libavcodec/allcodecs.c
==============================================================================
--- trunk/libavcodec/allcodecs.c Sat Dec 18 00:42:56 2010 (r26050)
+++ trunk/libavcodec/allcodecs.c Sat Dec 18 04:03:18 2010 (r26051)
@@ -225,6 +225,7 @@ void avcodec_register_all(void)
REGISTER_ENCDEC (ALAC, alac);
REGISTER_DECODER (ALS, als);
REGISTER_DECODER (AMRNB, amrnb);
+ REGISTER_DECODER (AMRWB, amrwb);
REGISTER_DECODER (APE, ape);
REGISTER_DECODER (ATRAC1, atrac1);
REGISTER_DECODER (ATRAC3, atrac3);
Added: trunk/libavcodec/amrwbdata.h
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ trunk/libavcodec/amrwbdata.h Sat Dec 18 04:03:18 2010 (r26051)
@@ -0,0 +1,1889 @@
+/*
+ * AMR wideband data and definitions
+ * Copyright (c) 2010 Marcelo Galvao Povoa
+ *
+ * This file is part of FFmpeg.
+ *
+ * FFmpeg is free software; you can redistribute it and/or
+ * modify it under the terms of the GNU Lesser General Public
+ * License as published by the Free Software Foundation; either
+ * version 2.1 of the License, or (at your option) any later version.
+ *
+ * FFmpeg is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
+ * Lesser General Public License for more details.
+ *
+ * You should have received a copy of the GNU Lesser General Public
+ * License along with FFmpeg; if not, write to the Free Software
+ * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
+ */
+
+/**
+ * @file
+ * AMR wideband data and definitions
+ */
+
+#ifndef AVCODEC_AMRWBDATA_H
+#define AVCODEC_AMRWBDATA_H
+
+#include <stdint.h>
+
+#define LP_ORDER 16 ///< linear predictive coding filter order
+#define LP_ORDER_16k 20 ///< lpc filter order at 16kHz
+#define HB_FIR_SIZE 30 ///< amount of past data needed by HB filters
+#define UPS_FIR_SIZE 12 ///< upsampling filter size
+#define UPS_MEM_SIZE (2 * UPS_FIR_SIZE)
+
+#define MIN_ISF_SPACING (128.0 / 32768.0) ///< minimum isf gap
+#define PRED_FACTOR (1.0 / 3.0)
+#define MIN_ENERGY -14.0 ///< initial innnovation energy (dB)
+#define ENERGY_MEAN 30.0 ///< mean innovation energy (dB) in all modes
+#define PREEMPH_FAC 0.68 ///< factor used to de-emphasize synthesis
+
+#define AMRWB_SFR_SIZE 64 ///< samples per subframe at 12.8 kHz
+#define AMRWB_SFR_SIZE_16k 80 ///< samples per subframe at 16 kHz
+#define AMRWB_P_DELAY_MAX 231 ///< maximum pitch delay value
+#define AMRWB_P_DELAY_MIN 34
+
+/* Relative mode ordering is sensitive */
+enum Mode {
+ MODE_6k60 = 0, ///< 6.60 kbit/s
+ MODE_8k85, ///< 8.85 kbit/s
+ MODE_12k65, ///< 12.65 kbit/s
+ MODE_14k25, ///< 14.25 kbit/s
+ MODE_15k85, ///< 15.85 kbit/s
+ MODE_18k25, ///< 18.25 kbit/s
+ MODE_19k85, ///< 19.85 kbit/s
+ MODE_23k05, ///< 23.05 kbit/s
+ MODE_23k85, ///< 23.85 kbit/s
+ MODE_SID, ///< comfort noise frame
+ /* 10-13: Future use */
+ SP_LOST = 14, ///< speech lost
+ NO_DATA ///< no transmission
+};
+
+/* All decoded parameters in these structs must be 2 bytes long
+ * because of the direct indexing at the frame parsing */
+typedef struct {
+ uint16_t adap; ///< adaptive codebook index
+ uint16_t ltp; ///< ltp-filtering flag
+ uint16_t vq_gain; ///< VQ adaptive and innovative gains
+ uint16_t hb_gain; ///< high-band energy index (mode 23k85 only)
+ uint16_t pul_ih[4]; ///< MSBs part of codebook index (high modes only)
+ uint16_t pul_il[4]; ///< LSBs part of codebook index
+} AMRWBSubFrame;
+
+typedef struct {
+ uint16_t vad; ///< voice activity detection flag
+ uint16_t isp_id[7]; ///< index of ISP subvectors
+ AMRWBSubFrame subframe[4]; ///< data for subframes
+} AMRWBFrame;
+
+/** The index of a frame parameter */
+#define AMR_BIT(field) (offsetof(AMRWBFrame, field) >> 1)
+/** The index of a subframe-specific parameter */
+#define AMR_OF(frame_num, variable) AMR_BIT(subframe[frame_num].variable)
+
+//As defined in 3GPP TS 26.201 V9.0.0
+//Tables for bit parsing in Core Frame speech frames
+//The reordered bits are in order of decreasing importance and
+//may be contiguously separated in Class A, B and C bits.
+
+// Each field in AMRWBFrame is stored as:
+// * one byte for the number of bits in the field
+// * one byte for the field index
+// * then, one byte for each bit of the field (from most-significant to least)
+// of the position of that bit in the AMR frame.
+static const uint16_t order_MODE_6k60[] = {
+ 1, AMR_BIT(vad), 7,
+ 8, AMR_BIT(isp_id[0]), 24, 33, 39, 12, 6, 5, 4, 13,
+ 8, AMR_BIT(isp_id[1]), 65, 79, 64, 78, 51, 61, 71, 70,
+ 7, AMR_BIT(isp_id[2]), 52, 55, 44, 54, 53, 43, 42,
+ 7, AMR_BIT(isp_id[3]), 60, 59, 58, 57, 56, 75, 74,
+ 6, AMR_BIT(isp_id[4]), 73, 72, 86, 87, 85, 84,
+ 8, AMR_OF(0, adap), 11, 10, 9, 8, 28, 27, 49, 69,
+ 6, AMR_OF(0, pul_il[0]), 83, 91, 99, 107, 115, 123,
+ 6, AMR_OF(0, pul_il[1]), 82, 103, 111, 119, 127, 135,
+ 6, AMR_OF(0, vq_gain), 38, 23, 34, 19, 3, 15,
+ 5, AMR_OF(1, adap), 32, 41, 63, 67, 77,
+ 6, AMR_OF(1, pul_il[0]), 81, 90, 98, 106, 114, 122,
+ 6, AMR_OF(1, pul_il[1]), 80, 102, 110, 118, 126, 134,
+ 6, AMR_OF(1, vq_gain), 26, 22, 36, 18, 2, 14,
+ 5, AMR_OF(2, adap), 45, 40, 50, 48, 68,
+ 6, AMR_OF(2, pul_il[0]), 95, 89, 97, 105, 113, 121,
+ 6, AMR_OF(2, pul_il[1]), 94, 101, 109, 117, 125, 133,
+ 6, AMR_OF(2, vq_gain), 37, 21, 35, 17, 1, 31,
+ 5, AMR_OF(3, adap), 47, 46, 62, 66, 76,
+ 6, AMR_OF(3, pul_il[0]), 93, 88, 96, 104, 112, 120,
+ 6, AMR_OF(3, pul_il[1]), 92, 100, 108, 116, 124, 132,
+ 6, AMR_OF(3, vq_gain), 25, 20, 29, 16, 0, 30,
+ 0
+};
+
+static const uint16_t order_MODE_8k85[] = {
+ 1, AMR_BIT(vad), 7,
+ 8, AMR_BIT(isp_id[0]), 47, 32, 2, 6, 3, 5, 4, 60,
+ 8, AMR_BIT(isp_id[1]), 69, 50, 67, 41, 51, 49, 59, 53,
+ 6, AMR_BIT(isp_id[2]), 40, 55, 43, 54, 42, 62,
+ 7, AMR_BIT(isp_id[3]), 63, 48, 52, 61, 77, 78, 72,
+ 7, AMR_BIT(isp_id[4]), 85, 56, 86, 68, 74, 73, 81,
+ 5, AMR_BIT(isp_id[5]), 82, 95, 80, 94, 91,
+ 5, AMR_BIT(isp_id[6]), 90, 89, 88, 103, 87,
+ 8, AMR_OF(0, adap), 1, 0, 15, 35, 33, 58, 64, 84,
+ 5, AMR_OF(0, pul_il[0]), 102, 118, 134, 150, 166,
+ 5, AMR_OF(0, pul_il[1]), 101, 114, 130, 146, 162,
+ 5, AMR_OF(0, pul_il[2]), 100, 126, 142, 158, 174,
+ 5, AMR_OF(0, pul_il[3]), 99, 122, 138, 154, 170,
+ 6, AMR_OF(0, vq_gain), 11, 39, 19, 31, 27, 23,
+ 5, AMR_OF(1, adap), 46, 71, 66, 76, 93,
+ 5, AMR_OF(1, pul_il[0]), 98, 117, 133, 149, 165,
+ 5, AMR_OF(1, pul_il[1]), 97, 113, 129, 145, 161,
+ 5, AMR_OF(1, pul_il[2]), 96, 125, 141, 157, 173,
+ 5, AMR_OF(1, pul_il[3]), 111, 121, 137, 153, 169,
+ 6, AMR_OF(1, vq_gain), 10, 38, 18, 30, 26, 22,
+ 8, AMR_OF(2, adap), 14, 13, 12, 34, 45, 57, 79, 83,
+ 5, AMR_OF(2, pul_il[0]), 110, 116, 132, 148, 164,
+ 5, AMR_OF(2, pul_il[1]), 109, 112, 128, 144, 160,
+ 5, AMR_OF(2, pul_il[2]), 108, 124, 140, 156, 172,
+ 5, AMR_OF(2, pul_il[3]), 107, 120, 136, 152, 168,
+ 6, AMR_OF(2, vq_gain), 9, 37, 17, 29, 25, 21,
+ 5, AMR_OF(3, adap), 44, 70, 65, 75, 92,
+ 5, AMR_OF(3, pul_il[0]), 106, 115, 131, 147, 163,
+ 5, AMR_OF(3, pul_il[1]), 105, 127, 143, 159, 175,
+ 5, AMR_OF(3, pul_il[2]), 104, 123, 139, 155, 171,
+ 5, AMR_OF(3, pul_il[3]), 119, 135, 151, 167, 183,
+ 6, AMR_OF(3, vq_gain), 8, 36, 16, 28, 24, 20,
+ 0
+};
+
+static const uint16_t order_MODE_12k65[] = {
+ 1, AMR_BIT(vad), 7,
+ 8, AMR_BIT(isp_id[0]), 55, 40, 14, 6, 15, 5, 0, 68,
+ 8, AMR_BIT(isp_id[1]), 77, 58, 75, 49, 59, 57, 67, 61,
+ 6, AMR_BIT(isp_id[2]), 48, 63, 51, 62, 50, 70,
+ 7, AMR_BIT(isp_id[3]), 71, 56, 60, 69, 85, 86, 80,
+ 7, AMR_BIT(isp_id[4]), 93, 64, 94, 76, 82, 81, 89,
+ 5, AMR_BIT(isp_id[5]), 90, 103, 88, 102, 99,
+ 5, AMR_BIT(isp_id[6]), 98, 97, 96, 111, 95,
+ 9, AMR_OF(0, adap), 13, 12, 11, 10, 9, 41, 66, 72,
+ 92,
+ 1, AMR_OF(0, ltp), 110,
+ 9, AMR_OF(0, pul_il[0]), 106, 122, 154, 186, 218, 134, 166, 198,
+ 230,
+ 9, AMR_OF(0, pul_il[1]), 105, 130, 162, 194, 226, 142, 174, 206,
+ 238,
+ 9, AMR_OF(0, pul_il[2]), 104, 138, 170, 202, 234, 150, 182, 214,
+ 246,
+ 9, AMR_OF(0, pul_il[3]), 119, 146, 178, 210, 242, 158, 190, 222,
+ 254,
+ 7, AMR_OF(0, vq_gain), 4, 19, 45, 27, 39, 33, 31,
+ 6, AMR_OF(1, adap), 35, 54, 79, 74, 84, 101,
+ 1, AMR_OF(1, ltp), 109,
+ 9, AMR_OF(1, pul_il[0]), 118, 121, 153, 185, 217, 133, 165, 197,
+ 229,
+ 9, AMR_OF(1, pul_il[1]), 117, 129, 161, 193, 225, 141, 173, 205,
+ 237,
+ 9, AMR_OF(1, pul_il[2]), 116, 137, 169, 201, 233, 149, 181, 213,
+ 245,
+ 9, AMR_OF(1, pul_il[3]), 115, 145, 177, 209, 241, 157, 189, 221,
+ 253,
+ 7, AMR_OF(1, vq_gain), 3, 18, 44, 26, 38, 32, 30,
+ 9, AMR_OF(2, adap), 8, 23, 22, 21, 20, 52, 65, 87,
+ 91,
+ 1, AMR_OF(2, ltp), 108,
+ 9, AMR_OF(2, pul_il[0]), 114, 120, 152, 184, 216, 132, 164, 196,
+ 228,
+ 9, AMR_OF(2, pul_il[1]), 113, 128, 160, 192, 224, 140, 172, 204,
+ 236,
+ 9, AMR_OF(2, pul_il[2]), 112, 136, 168, 200, 232, 148, 180, 212,
+ 244,
+ 9, AMR_OF(2, pul_il[3]), 127, 144, 176, 208, 240, 156, 188, 220,
+ 252,
+ 7, AMR_OF(2, vq_gain), 2, 17, 43, 25, 37, 47, 29,
+ 6, AMR_OF(3, adap), 34, 53, 78, 73, 83, 100,
+ 1, AMR_OF(3, ltp), 107,
+ 9, AMR_OF(3, pul_il[0]), 126, 135, 167, 199, 231, 131, 163, 195,
+ 227,
+ 9, AMR_OF(3, pul_il[1]), 125, 143, 175, 207, 239, 139, 171, 203,
+ 235,
+ 9, AMR_OF(3, pul_il[2]), 124, 151, 183, 215, 247, 147, 179, 211,
+ 243,
+ 9, AMR_OF(3, pul_il[3]), 123, 159, 191, 223, 255, 155, 187, 219,
+ 251,
+ 7, AMR_OF(3, vq_gain), 1, 16, 42, 24, 36, 46, 28,
+ 0
+};
+
+static const uint16_t order_MODE_14k25[] = {
+ 1, AMR_BIT(vad), 7,
+ 8, AMR_BIT(isp_id[0]), 55, 40, 14, 6, 15, 5, 0, 68,
+ 8, AMR_BIT(isp_id[1]), 77, 58, 75, 49, 59, 57, 67, 61,
+ 6, AMR_BIT(isp_id[2]), 48, 63, 51, 62, 50, 70,
+ 7, AMR_BIT(isp_id[3]), 71, 56, 60, 69, 85, 86, 80,
+ 7, AMR_BIT(isp_id[4]), 93, 64, 94, 76, 82, 81, 89,
+ 5, AMR_BIT(isp_id[5]), 90, 103, 88, 102, 99,
+ 5, AMR_BIT(isp_id[6]), 98, 97, 96, 111, 95,
+ 9, AMR_OF(0, adap), 13, 12, 11, 10, 9, 41, 66, 72,
+ 92,
+ 1, AMR_OF(0, ltp), 110,
+ 13, AMR_OF(0, pul_il[0]), 114, 186, 210, 234, 258, 106, 126, 162,
+ 170, 198, 222, 246, 270,
+ 13, AMR_OF(0, pul_il[1]), 122, 194, 218, 242, 266, 118, 134, 174,
+ 182, 206, 230, 254, 278,
+ 9, AMR_OF(0, pul_il[2]), 130, 138, 146, 154, 178, 202, 226, 250,
+ 274,
+ 9, AMR_OF(0, pul_il[3]), 142, 150, 158, 166, 190, 214, 238, 262,
+ 286,
+ 7, AMR_OF(0, vq_gain), 4, 19, 45, 27, 39, 33, 31,
+ 6, AMR_OF(1, adap), 35, 54, 79, 74, 84, 101,
+ 1, AMR_OF(1, ltp), 109,
+ 13, AMR_OF(1, pul_il[0]), 113, 185, 209, 233, 257, 105, 125, 161,
+ 169, 197, 221, 245, 269,
+ 13, AMR_OF(1, pul_il[1]), 121, 193, 217, 241, 265, 117, 133, 173,
+ 181, 205, 229, 253, 277,
+ 9, AMR_OF(1, pul_il[2]), 129, 137, 145, 153, 177, 201, 225, 249,
+ 273,
+ 9, AMR_OF(1, pul_il[3]), 141, 149, 157, 165, 189, 213, 237, 261,
+ 285,
+ 7, AMR_OF(1, vq_gain), 3, 18, 44, 26, 38, 32, 30,
+ 9, AMR_OF(2, adap), 8, 23, 22, 21, 20, 52, 65, 87,
+ 91,
+ 1, AMR_OF(2, ltp), 108,
+ 13, AMR_OF(2, pul_il[0]), 112, 184, 208, 232, 256, 104, 124, 160,
+ 168, 196, 220, 244, 268,
+ 13, AMR_OF(2, pul_il[1]), 120, 192, 216, 240, 264, 116, 132, 172,
+ 180, 204, 228, 252, 276,
+ 9, AMR_OF(2, pul_il[2]), 128, 136, 144, 152, 176, 200, 224, 248,
+ 272,
+ 9, AMR_OF(2, pul_il[3]), 140, 148, 156, 164, 188, 212, 236, 260,
+ 284,
+ 7, AMR_OF(2, vq_gain), 2, 17, 43, 25, 37, 47, 29,
+ 6, AMR_OF(3, adap), 34, 53, 78, 73, 83, 100,
+ 1, AMR_OF(3, ltp), 107,
+ 13, AMR_OF(3, pul_il[0]), 127, 199, 223, 247, 271, 119, 123, 175,
+ 183, 195, 219, 243, 267,
+ 13, AMR_OF(3, pul_il[1]), 135, 207, 231, 255, 279, 115, 131, 171,
+ 179, 203, 227, 251, 275,
+ 9, AMR_OF(3, pul_il[2]), 143, 151, 159, 167, 191, 215, 239, 263,
+ 287,
+ 9, AMR_OF(3, pul_il[3]), 139, 147, 155, 163, 187, 211, 235, 259,
+ 283,
+ 7, AMR_OF(3, vq_gain), 1, 16, 42, 24, 36, 46, 28,
+ 0
+};
+
+static const uint16_t order_MODE_15k85[] = {
+ 1, AMR_BIT(vad), 7,
+ 8, AMR_BIT(isp_id[0]), 55, 40, 14, 6, 15, 5, 0, 68,
+ 8, AMR_BIT(isp_id[1]), 77, 58, 75, 49, 59, 57, 67, 61,
+ 6, AMR_BIT(isp_id[2]), 48, 63, 51, 62, 50, 70,
+ 7, AMR_BIT(isp_id[3]), 71, 56, 60, 69, 85, 86, 80,
+ 7, AMR_BIT(isp_id[4]), 93, 64, 94, 76, 82, 81, 89,
+ 5, AMR_BIT(isp_id[5]), 90, 103, 88, 102, 99,
+ 5, AMR_BIT(isp_id[6]), 98, 97, 96, 111, 95,
+ 9, AMR_OF(0, adap), 13, 12, 11, 10, 9, 41, 66, 72,
+ 92,
+ 1, AMR_OF(0, ltp), 110,
+ 13, AMR_OF(0, pul_il[0]), 122, 154, 170, 218, 266, 138, 106, 182,
+ 230, 278, 178, 226, 274,
+ 13, AMR_OF(0, pul_il[1]), 134, 166, 190, 238, 286, 150, 118, 186,
+ 234, 282, 198, 246, 294,
+ 13, AMR_OF(0, pul_il[2]), 130, 162, 194, 242, 290, 146, 114, 206,
+ 254, 302, 202, 250, 298,
+ 13, AMR_OF(0, pul_il[3]), 142, 174, 214, 262, 310, 158, 126, 210,
+ 258, 306, 222, 270, 318,
+ 7, AMR_OF(0, vq_gain), 4, 19, 45, 27, 39, 33, 31,
+ 6, AMR_OF(1, adap), 35, 54, 79, 74, 84, 101,
+ 1, AMR_OF(1, ltp), 109,
+ 13, AMR_OF(1, pul_il[0]), 121, 153, 169, 217, 265, 137, 105, 181,
+ 229, 277, 177, 225, 273,
+ 13, AMR_OF(1, pul_il[1]), 133, 165, 189, 237, 285, 149, 117, 185,
+ 233, 281, 197, 245, 293,
+ 13, AMR_OF(1, pul_il[2]), 129, 161, 193, 241, 289, 145, 113, 205,
+ 253, 301, 201, 249, 297,
+ 13, AMR_OF(1, pul_il[3]), 141, 173, 213, 261, 309, 157, 125, 209,
+ 257, 305, 221, 269, 317,
+ 7, AMR_OF(1, vq_gain), 3, 18, 44, 26, 38, 32, 30,
+ 9, AMR_OF(2, adap), 8, 23, 22, 21, 20, 52, 65, 87,
+ 91,
+ 1, AMR_OF(2, ltp), 108,
+ 13, AMR_OF(2, pul_il[0]), 120, 152, 168, 216, 264, 136, 104, 180,
+ 228, 276, 176, 224, 272,
+ 13, AMR_OF(2, pul_il[1]), 132, 164, 188, 236, 284, 148, 116, 184,
+ 232, 280, 196, 244, 292,
+ 13, AMR_OF(2, pul_il[2]), 128, 160, 192, 240, 288, 144, 112, 204,
+ 252, 300, 200, 248, 296,
+ 13, AMR_OF(2, pul_il[3]), 140, 172, 212, 260, 308, 156, 124, 208,
+ 256, 304, 220, 268, 316,
+ 7, AMR_OF(2, vq_gain), 2, 17, 43, 25, 37, 47, 29,
+ 6, AMR_OF(3, adap), 34, 53, 78, 73, 83, 100,
+ 1, AMR_OF(3, ltp), 107,
+ 13, AMR_OF(3, pul_il[0]), 135, 167, 183, 231, 279, 151, 119, 179,
+ 227, 275, 191, 239, 287,
+ 13, AMR_OF(3, pul_il[1]), 131, 163, 187, 235, 283, 147, 115, 199,
+ 247, 295, 195, 243, 291,
+ 13, AMR_OF(3, pul_il[2]), 143, 175, 207, 255, 303, 159, 127, 203,
+ 251, 299, 215, 263, 311,
+ 13, AMR_OF(3, pul_il[3]), 139, 171, 211, 259, 307, 155, 123, 223,
+ 271, 319, 219, 267, 315,
+ 7, AMR_OF(3, vq_gain), 1, 16, 42, 24, 36, 46, 28,
+ 0
+};
+
+static const uint16_t order_MODE_18k25[] = {
+ 1, AMR_BIT(vad), 7,
+ 8, AMR_BIT(isp_id[0]), 55, 40, 14, 6, 15, 5, 0, 68,
+ 8, AMR_BIT(isp_id[1]), 77, 58, 75, 49, 59, 57, 67, 61,
+ 6, AMR_BIT(isp_id[2]), 48, 63, 51, 62, 50, 70,
+ 7, AMR_BIT(isp_id[3]), 71, 56, 60, 69, 85, 86, 80,
+ 7, AMR_BIT(isp_id[4]), 93, 64, 94, 76, 82, 81, 89,
+ 5, AMR_BIT(isp_id[5]), 90, 103, 88, 102, 99,
+ 5, AMR_BIT(isp_id[6]), 98, 97, 96, 111, 95,
+ 9, AMR_OF(0, adap), 13, 12, 11, 10, 9, 41, 66, 72,
+ 92,
+ 1, AMR_OF(0, ltp), 110,
+ 2, AMR_OF(0, pul_ih[0]), 124, 115,
+ 2, AMR_OF(0, pul_ih[1]), 150, 117,
+ 2, AMR_OF(0, pul_ih[2]), 129, 114,
+ 2, AMR_OF(0, pul_ih[3]), 121, 131,
+ 14, AMR_OF(0, pul_il[0]), 161, 257, 343, 199, 177, 303, 204, 173,
+ 168, 260, 277, 307, 338, 128,
+ 14, AMR_OF(0, pul_il[1]), 194, 286, 347, 222, 214, 316, 236, 152,
+ 166, 242, 284, 308, 344, 142,
+ 14, AMR_OF(0, pul_il[2]), 169, 273, 353, 202, 189, 311, 240, 200,
+ 171, 261, 309, 296, 345, 130,
+ 14, AMR_OF(0, pul_il[3]), 198, 275, 349, 187, 163, 282, 193, 195,
+ 175, 234, 265, 289, 328, 119,
+ 7, AMR_OF(0, vq_gain), 4, 19, 45, 27, 39, 33, 31,
+ 6, AMR_OF(1, adap), 35, 54, 79, 74, 84, 101,
+ 1, AMR_OF(1, ltp), 109,
+ 2, AMR_OF(1, pul_ih[0]), 139, 104,
+ 2, AMR_OF(1, pul_ih[1]), 135, 118,
+ 2, AMR_OF(1, pul_ih[2]), 112, 127,
+ 2, AMR_OF(1, pul_ih[3]), 140, 141,
+ 14, AMR_OF(1, pul_il[0]), 179, 276, 340, 225, 223, 321, 235, 190,
+ 182, 271, 310, 315, 352, 125,
+ 14, AMR_OF(1, pul_il[1]), 153, 264, 329, 232, 209, 323, 231, 165,
+ 191, 279, 290, 312, 367, 134,
+ 14, AMR_OF(1, pul_il[2]), 167, 269, 341, 205, 197, 298, 224, 160,
+ 170, 259, 280, 317, 357, 148,
+ 14, AMR_OF(1, pul_il[3]), 203, 272, 342, 227, 192, 299, 233, 172,
+ 183, 256, 283, 326, 355, 106,
+ 7, AMR_OF(1, vq_gain), 3, 18, 44, 26, 38, 32, 30,
+ 9, AMR_OF(2, adap), 8, 23, 22, 21, 20, 52, 65, 87,
+ 91,
+ 1, AMR_OF(2, ltp), 108,
+ 2, AMR_OF(2, pul_ih[0]), 144, 120,
+ 2, AMR_OF(2, pul_ih[1]), 157, 123,
+ 2, AMR_OF(2, pul_ih[2]), 145, 138,
+ 2, AMR_OF(2, pul_ih[3]), 132, 154,
+ 14, AMR_OF(2, pul_il[0]), 241, 319, 365, 252, 253, 331, 254, 230,
+ 220, 263, 285, 314, 364, 156,
+ 14, AMR_OF(2, pul_il[1]), 247, 291, 339, 249, 250, 332, 267, 196,
+ 207, 268, 304, 324, 356, 158,
+ 14, AMR_OF(2, pul_il[2]), 210, 300, 348, 243, 237, 333, 246, 206,
+ 219, 266, 318, 335, 363, 159,
+ 14, AMR_OF(2, pul_il[3]), 239, 306, 366, 221, 226, 297, 251, 184,
+ 178, 258, 292, 305, 346, 116,
+ 7, AMR_OF(2, vq_gain), 2, 17, 43, 25, 37, 47, 29,
+ 6, AMR_OF(3, adap), 34, 53, 78, 73, 83, 100,
+ 1, AMR_OF(3, ltp), 107,
+ 2, AMR_OF(3, pul_ih[0]), 143, 126,
+ 2, AMR_OF(3, pul_ih[1]), 137, 122,
+ 2, AMR_OF(3, pul_ih[2]), 149, 105,
+ 2, AMR_OF(3, pul_ih[3]), 133, 136,
+ 14, AMR_OF(3, pul_il[0]), 162, 287, 337, 244, 229, 322, 218, 180,
+ 186, 262, 274, 288, 351, 146,
+ 14, AMR_OF(3, pul_il[1]), 212, 294, 358, 248, 228, 334, 215, 174,
+ 176, 270, 293, 301, 354, 147,
+ 14, AMR_OF(3, pul_il[2]), 185, 327, 336, 211, 213, 313, 245, 181,
+ 188, 255, 281, 325, 350, 151,
+ 14, AMR_OF(3, pul_il[3]), 201, 295, 359, 216, 208, 320, 238, 164,
+ 155, 217, 278, 302, 330, 113,
+ 7, AMR_OF(3, vq_gain), 1, 16, 42, 24, 36, 46, 28,
+ 0
+};
+
+static const uint16_t order_MODE_19k85[] = {
+ 1, AMR_BIT(vad), 7,
+ 8, AMR_BIT(isp_id[0]), 55, 40, 14, 6, 15, 5, 0, 68,
+ 8, AMR_BIT(isp_id[1]), 77, 58, 75, 49, 59, 57, 67, 61,
+ 6, AMR_BIT(isp_id[2]), 48, 63, 51, 62, 50, 70,
+ 7, AMR_BIT(isp_id[3]), 71, 56, 60, 69, 85, 86, 80,
+ 7, AMR_BIT(isp_id[4]), 93, 64, 94, 76, 82, 81, 89,
+ 5, AMR_BIT(isp_id[5]), 90, 103, 88, 102, 99,
+ 5, AMR_BIT(isp_id[6]), 98, 97, 96, 111, 95,
+ 9, AMR_OF(0, adap), 13, 12, 11, 10, 9, 41, 66, 72,
+ 92,
+ 1, AMR_OF(0, ltp), 110,
+ 10, AMR_OF(0, pul_ih[0]), 134, 153, 263, 342, 399, 154, 106, 177,
+ 317, 265,
+ 10, AMR_OF(0, pul_ih[1]), 128, 167, 270, 351, 385, 160, 105, 213,
+ 329, 259,
+ 2, AMR_OF(0, pul_ih[2]), 123, 147,
+ 2, AMR_OF(0, pul_ih[3]), 131, 143,
+ 10, AMR_OF(0, pul_il[0]), 346, 118, 170, 201, 296, 368, 250, 284,
+ 341, 391,
+ 10, AMR_OF(0, pul_il[1]), 345, 104, 166, 196, 281, 374, 242, 269,
+ 327, 390,
+ 14, AMR_OF(0, pul_il[2]), 141, 171, 291, 364, 229, 210, 308, 228,
+ 206, 200, 258, 295, 313, 361,
+ 14, AMR_OF(0, pul_il[3]), 144, 188, 282, 366, 217, 216, 309, 218,
+ 193, 182, 245, 287, 300, 367,
+ 7, AMR_OF(0, vq_gain), 4, 19, 45, 27, 39, 33, 31,
+ 6, AMR_OF(1, adap), 35, 54, 79, 74, 84, 101,
+ 1, AMR_OF(1, ltp), 109,
+ 10, AMR_OF(1, pul_ih[0]), 139, 169, 267, 348, 389, 163, 116, 189,
+ 343, 268,
+ 10, AMR_OF(1, pul_ih[1]), 120, 161, 249, 339, 397, 152, 114, 230,
+ 334, 303,
+ 2, AMR_OF(1, pul_ih[2]), 125, 138,
+ 2, AMR_OF(1, pul_ih[3]), 112, 129,
+ 10, AMR_OF(1, pul_il[0]), 349, 122, 162, 203, 288, 372, 278, 274,
+ 312, 377,
+ 10, AMR_OF(1, pul_il[1]), 357, 126, 165, 214, 298, 362, 252, 260,
+ 321, 378,
+ 14, AMR_OF(1, pul_il[2]), 150, 199, 266, 355, 211, 180, 285, 241,
+ 195, 198, 243, 275, 323, 375,
+ 14, AMR_OF(1, pul_il[3]), 142, 191, 256, 353, 208, 220, 314, 237,
+ 190, 212, 255, 304, 318, 371,
+ 7, AMR_OF(1, vq_gain), 3, 18, 44, 26, 38, 32, 30,
+ 9, AMR_OF(2, adap), 8, 23, 22, 21, 20, 52, 65, 87,
+ 91,
+ 1, AMR_OF(2, ltp), 108,
+ 10, AMR_OF(2, pul_ih[0]), 159, 168, 302, 356, 395, 178, 132, 185,
+ 330, 286,
+ 10, AMR_OF(2, pul_ih[1]), 158, 181, 292, 358, 396, 176, 133, 235,
+ 331, 276,
+ 2, AMR_OF(2, pul_ih[2]), 130, 157,
+ 2, AMR_OF(2, pul_ih[3]), 124, 136,
+ 10, AMR_OF(2, pul_il[0]), 354, 121, 194, 246, 322, 379, 272, 273,
+ 332, 398,
+ 10, AMR_OF(2, pul_il[1]), 359, 140, 186, 236, 333, 376, 290, 301,
+ 338, 387,
+ 14, AMR_OF(2, pul_il[2]), 155, 227, 319, 369, 253, 254, 350, 248,
+ 224, 239, 240, 293, 315, 383,
+ 14, AMR_OF(2, pul_il[3]), 156, 209, 297, 373, 225, 215, 326, 247,
+ 197, 184, 232, 289, 310, 365,
+ 7, AMR_OF(2, vq_gain), 2, 17, 43, 25, 37, 47, 29,
+ 6, AMR_OF(3, adap), 34, 53, 78, 73, 83, 100,
+ 1, AMR_OF(3, ltp), 107,
+ 10, AMR_OF(3, pul_ih[0]), 148, 164, 264, 340, 388, 183, 117, 205,
+ 336, 261,
+ 10, AMR_OF(3, pul_ih[1]), 146, 174, 257, 335, 384, 173, 113, 187,
+ 320, 279,
+ 2, AMR_OF(3, pul_ih[2]), 127, 151,
+ 2, AMR_OF(3, pul_ih[3]), 119, 137,
+ 10, AMR_OF(3, pul_il[0]), 352, 135, 172, 238, 306, 381, 262, 271,
+ 328, 382,
+ 10, AMR_OF(3, pul_il[1]), 347, 115, 179, 219, 305, 380, 277, 294,
+ 337, 386,
+ 14, AMR_OF(3, pul_il[2]), 145, 192, 307, 370, 234, 223, 324, 244,
+ 202, 204, 251, 299, 325, 360,
+ 14, AMR_OF(3, pul_il[3]), 149, 221, 311, 363, 226, 222, 316, 231,
+ 207, 175, 233, 280, 283, 344,
+ 7, AMR_OF(3, vq_gain), 1, 16, 42, 24, 36, 46, 28,
+ 0
+};
+
+static const uint16_t order_MODE_23k05[] = {
+ 1, AMR_BIT(vad), 7,
+ 8, AMR_BIT(isp_id[0]), 55, 40, 14, 6, 15, 5, 0, 68,
+ 8, AMR_BIT(isp_id[1]), 77, 58, 75, 49, 59, 57, 67, 61,
+ 6, AMR_BIT(isp_id[2]), 48, 63, 51, 62, 50, 70,
+ 7, AMR_BIT(isp_id[3]), 71, 56, 60, 69, 85, 86, 80,
+ 7, AMR_BIT(isp_id[4]), 93, 64, 94, 76, 82, 81, 89,
+ 5, AMR_BIT(isp_id[5]), 90, 103, 88, 102, 99,
+ 5, AMR_BIT(isp_id[6]), 98, 97, 96, 111, 95,
+ 9, AMR_OF(0, adap), 13, 12, 11, 10, 9, 41, 66, 72,
+ 92,
+ 1, AMR_OF(0, ltp), 110,
+ 11, AMR_OF(0, pul_ih[0]), 118, 129, 131, 153, 170, 282, 298, 210,
+ 191, 357, 317,
+ 11, AMR_OF(0, pul_ih[1]), 126, 146, 135, 165, 187, 273, 345, 295,
+ 172, 338, 340,
+ 11, AMR_OF(0, pul_ih[2]), 119, 137, 141, 167, 208, 304, 366, 256,
+ 177, 339, 328,
+ 11, AMR_OF(0, pul_ih[3]), 116, 130, 120, 166, 190, 252, 311, 239,
+ 173, 343, 318,
+ 11, AMR_OF(0, pul_il[0]), 245, 180, 342, 424, 259, 277, 266, 380,
+ 398, 423, 440,
+ 11, AMR_OF(0, pul_il[1]), 218, 207, 367, 434, 201, 240, 275, 363,
+ 399, 419, 452,
+ 11, AMR_OF(0, pul_il[2]), 274, 188, 348, 425, 242, 204, 262, 365,
+ 402, 431, 463,
+ 11, AMR_OF(0, pul_il[3]), 221, 183, 337, 439, 243, 216, 251, 354,
+ 390, 411, 462,
+ 7, AMR_OF(0, vq_gain), 4, 19, 45, 27, 39, 33, 31,
+ 6, AMR_OF(1, adap), 35, 54, 79, 74, 84, 101,
+ 1, AMR_OF(1, ltp), 109,
+ 11, AMR_OF(1, pul_ih[0]), 115, 140, 142, 161, 230, 291, 351, 235,
+ 181, 293, 310,
+ 11, AMR_OF(1, pul_ih[1]), 104, 138, 132, 162, 211, 315, 347, 233,
+ 176, 320, 329,
+ 11, AMR_OF(1, pul_ih[2]), 106, 134, 125, 154, 205, 267, 306, 220,
+ 185, 330, 297,
+ 11, AMR_OF(1, pul_ih[3]), 105, 148, 122, 152, 215, 302, 350, 254,
+ 178, 319, 313,
+ 11, AMR_OF(1, pul_il[0]), 269, 189, 382, 432, 272, 228, 263, 383,
+ 406, 422, 453,
+ 11, AMR_OF(1, pul_il[1]), 286, 206, 377, 446, 226, 222, 265, 368,
+ 404, 416, 454,
+ 11, AMR_OF(1, pul_il[2]), 247, 195, 358, 445, 224, 236, 309, 341,
+ 375, 408, 449,
+ 11, AMR_OF(1, pul_il[3]), 225, 192, 359, 436, 250, 258, 290, 389,
+ 400, 420, 448,
+ 7, AMR_OF(1, vq_gain), 3, 18, 44, 26, 38, 32, 30,
+ 9, AMR_OF(2, adap), 8, 23, 22, 21, 20, 52, 65, 87,
+ 91,
+ 1, AMR_OF(2, ltp), 108,
+ 11, AMR_OF(2, pul_ih[0]), 139, 144, 145, 169, 234, 327, 395, 299,
+ 244, 356, 379,
+ 11, AMR_OF(2, pul_ih[1]), 127, 156, 158, 171, 231, 308, 397, 355,
+ 261, 371, 335,
+ 11, AMR_OF(2, pul_ih[2]), 123, 155, 157, 193, 241, 362, 384, 323,
+ 238, 392, 361,
+ 11, AMR_OF(2, pul_ih[3]), 114, 147, 121, 175, 196, 333, 373, 303,
+ 184, 353, 322,
+ 11, AMR_OF(2, pul_il[0]), 271, 203, 385, 442, 307, 276, 334, 405,
+ 412, 427, 459,
+ 11, AMR_OF(2, pul_il[1]), 278, 200, 388, 447, 292, 288, 296, 403,
+ 415, 429, 460,
+ 11, AMR_OF(2, pul_il[2]), 312, 214, 393, 433, 279, 301, 314, 391,
+ 410, 426, 450,
+ 11, AMR_OF(2, pul_il[3]), 280, 186, 376, 437, 268, 260, 255, 364,
+ 414, 417, 441,
+ 7, AMR_OF(2, vq_gain), 2, 17, 43, 25, 37, 47, 29,
+ 6, AMR_OF(3, adap), 34, 53, 78, 73, 83, 100,
+ 1, AMR_OF(3, ltp), 107,
+ 11, AMR_OF(3, pul_ih[0]), 112, 159, 143, 164, 213, 281, 332, 284,
+ 168, 344, 325,
+ 11, AMR_OF(3, pul_ih[1]), 113, 150, 149, 179, 199, 316, 324, 285,
+ 237, 360, 336,
+ 11, AMR_OF(3, pul_ih[2]), 124, 136, 151, 174, 209, 326, 349, 248,
+ 198, 374, 331,
+ 11, AMR_OF(3, pul_ih[3]), 117, 128, 133, 163, 202, 300, 372, 305,
+ 194, 387, 321,
+ 11, AMR_OF(3, pul_il[0]), 249, 182, 352, 428, 253, 264, 289, 413,
+ 407, 418, 461,
+ 11, AMR_OF(3, pul_il[1]), 287, 212, 369, 444, 223, 246, 217, 346,
+ 394, 401, 451,
+ 11, AMR_OF(3, pul_il[2]), 219, 197, 378, 435, 229, 257, 283, 396,
+ 409, 430, 455,
+ 11, AMR_OF(3, pul_il[3]), 232, 160, 370, 438, 227, 270, 294, 381,
+ 386, 421, 443,
+ 7, AMR_OF(3, vq_gain), 1, 16, 42, 24, 36, 46, 28,
+ 0
+};
+
+static const uint16_t order_MODE_23k85[] = {
+ 1, AMR_BIT(vad), 7,
+ 8, AMR_BIT(isp_id[0]), 55, 40, 14, 6, 15, 5, 0, 68,
+ 8, AMR_BIT(isp_id[1]), 93, 58, 91, 49, 59, 57, 67, 61,
+ 6, AMR_BIT(isp_id[2]), 48, 63, 51, 62, 50, 70,
+ 7, AMR_BIT(isp_id[3]), 71, 56, 60, 69, 101, 102, 96,
+ 7, AMR_BIT(isp_id[4]), 109, 64, 110, 92, 98, 97, 105,
+ 5, AMR_BIT(isp_id[5]), 106, 119, 104, 118, 115,
+ 5, AMR_BIT(isp_id[6]), 114, 113, 112, 127, 111,
+ 9, AMR_OF(0, adap), 13, 12, 11, 10, 9, 41, 66, 88,
+ 108,
+ 1, AMR_OF(0, ltp), 126,
+ 11, AMR_OF(0, pul_ih[0]), 134, 145, 147, 169, 186, 298, 314, 226,
+ 207, 373, 333,
+ 11, AMR_OF(0, pul_ih[1]), 142, 162, 151, 181, 203, 289, 361, 311,
+ 188, 354, 356,
+ 11, AMR_OF(0, pul_ih[2]), 135, 153, 157, 183, 224, 320, 382, 272,
+ 193, 355, 344,
+ 11, AMR_OF(0, pul_ih[3]), 132, 146, 136, 182, 206, 268, 327, 255,
+ 189, 359, 334,
+ 11, AMR_OF(0, pul_il[0]), 261, 196, 358, 440, 275, 293, 282, 396,
+ 414, 439, 456,
+ 11, AMR_OF(0, pul_il[1]), 234, 223, 383, 450, 217, 256, 291, 379,
+ 415, 435, 468,
+ 11, AMR_OF(0, pul_il[2]), 290, 204, 364, 441, 258, 220, 278, 381,
+ 418, 447, 479,
+ 11, AMR_OF(0, pul_il[3]), 237, 199, 353, 455, 259, 232, 267, 370,
+ 406, 427, 478,
+ 7, AMR_OF(0, vq_gain), 4, 19, 45, 27, 39, 33, 31,
+ 4, AMR_OF(0, hb_gain), 79, 78, 77, 76,
+ 6, AMR_OF(1, adap), 35, 54, 95, 90, 100, 117,
+ 1, AMR_OF(1, ltp), 125,
+ 11, AMR_OF(1, pul_ih[0]), 131, 156, 158, 177, 246, 307, 367, 251,
+ 197, 309, 326,
+ 11, AMR_OF(1, pul_ih[1]), 120, 154, 148, 178, 227, 331, 363, 249,
+ 192, 336, 345,
+ 11, AMR_OF(1, pul_ih[2]), 122, 150, 141, 170, 221, 283, 322, 236,
+ 201, 346, 313,
+ 11, AMR_OF(1, pul_ih[3]), 121, 164, 138, 168, 231, 318, 366, 270,
+ 194, 335, 329,
+ 11, AMR_OF(1, pul_il[0]), 285, 205, 398, 448, 288, 244, 279, 399,
+ 422, 438, 469,
+ 11, AMR_OF(1, pul_il[1]), 302, 222, 393, 462, 242, 238, 281, 384,
+ 420, 432, 470,
+ 11, AMR_OF(1, pul_il[2]), 263, 211, 374, 461, 240, 252, 325, 357,
+ 391, 424, 465,
+ 11, AMR_OF(1, pul_il[3]), 241, 208, 375, 452, 266, 274, 306, 405,
+ 416, 436, 464,
+ 7, AMR_OF(1, vq_gain), 3, 18, 44, 26, 38, 32, 30,
+ 4, AMR_OF(1, hb_gain), 75, 74, 73, 72,
+ 9, AMR_OF(2, adap), 8, 23, 22, 21, 20, 52, 65, 103,
+ 107,
+ 1, AMR_OF(2, ltp), 124,
+ 11, AMR_OF(2, pul_ih[0]), 155, 160, 161, 185, 250, 343, 411, 315,
+ 260, 372, 395,
+ 11, AMR_OF(2, pul_ih[1]), 143, 172, 174, 187, 247, 324, 413, 371,
+ 277, 387, 351,
+ 11, AMR_OF(2, pul_ih[2]), 139, 171, 173, 209, 257, 378, 400, 339,
+ 254, 408, 377,
+ 11, AMR_OF(2, pul_ih[3]), 130, 163, 137, 191, 212, 349, 389, 319,
+ 200, 369, 338,
+ 11, AMR_OF(2, pul_il[0]), 287, 219, 401, 458, 323, 292, 350, 421,
+ 428, 443, 475,
+ 11, AMR_OF(2, pul_il[1]), 294, 216, 404, 463, 308, 304, 312, 419,
+ 431, 445, 476,
+ 11, AMR_OF(2, pul_il[2]), 328, 230, 409, 449, 295, 317, 330, 407,
+ 426, 442, 466,
+ 11, AMR_OF(2, pul_il[3]), 296, 202, 392, 453, 284, 276, 271, 380,
+ 430, 433, 457,
+ 7, AMR_OF(2, vq_gain), 2, 17, 43, 25, 37, 47, 29,
+ 4, AMR_OF(2, hb_gain), 87, 86, 85, 84,
+ 6, AMR_OF(3, adap), 34, 53, 94, 89, 99, 116,
+ 1, AMR_OF(3, ltp), 123,
+ 11, AMR_OF(3, pul_ih[0]), 128, 175, 159, 180, 229, 297, 348, 300,
+ 184, 360, 341,
+ 11, AMR_OF(3, pul_ih[1]), 129, 166, 165, 195, 215, 332, 340, 301,
+ 253, 376, 352,
+ 11, AMR_OF(3, pul_ih[2]), 140, 152, 167, 190, 225, 342, 365, 264,
+ 214, 390, 347,
+ 11, AMR_OF(3, pul_ih[3]), 133, 144, 149, 179, 218, 316, 388, 321,
+ 210, 403, 337,
+ 11, AMR_OF(3, pul_il[0]), 265, 198, 368, 444, 269, 280, 305, 429,
+ 423, 434, 477,
+ 11, AMR_OF(3, pul_il[1]), 303, 228, 385, 460, 239, 262, 233, 362,
+ 410, 417, 467,
+ 11, AMR_OF(3, pul_il[2]), 235, 213, 394, 451, 245, 273, 299, 412,
+ 425, 446, 471,
+ 11, AMR_OF(3, pul_il[3]), 248, 176, 386, 454, 243, 286, 310, 397,
+ 402, 437, 459,
+ 7, AMR_OF(3, vq_gain), 1, 16, 42, 24, 36, 46, 28,
+ 4, AMR_OF(3, hb_gain), 83, 82, 81, 80,
+ 0
+};
+
+/** Reordering array addresses for each mode */
+static const uint16_t* amr_bit_orderings_by_mode[] = {
+ order_MODE_6k60,
+ order_MODE_8k85,
+ order_MODE_12k65,
+ order_MODE_14k25,
+ order_MODE_15k85,
+ order_MODE_18k25,
+ order_MODE_19k85,
+ order_MODE_23k05,
+ order_MODE_23k85
+};
+
+// Extracted from 3GPP TS 26.173 V9.0.0 (qpisf_2s.tab)
+// The *_36b tables are used in 6k60 mode
+// Stored in fixed-point to save some space
+/** Indexed tables for retrieval of quantized ISF vectors in Q15 */
+static const int16_t dico1_isf[256][9] = {
+ { 579, 1081, 1035, 390, 3, -263, -198, -82, 38},
+ { 18, -68, -12, 313, 761, 405, 249, 111, -76},
+ { 740, 1263, 1292, 1006, 997, 1019, 1017, 976, 923},
+ { -91, 827, 948, 648, 613, 535, 522, 490, 421},
+ { 41, -44, -281, -472, 652, 534, 193, 135, -90},
+ { 41, -121, -356, -60, 663, 307, 61, -48, -344},
+ { 557, 946, 1049, 867, 846, 990, 1112, 1262, 1241},
+ { -118, -204, 328, 512, 870, 793, 610, 402, 186},
+ { 156, 293, 74, -338, -475, -897, -594, -161, -497},
+ { 226, 131, -138, 307, 169, -271, -164, -387, -624},
+ { 62, -32, -61, -252, -541, -828, -1027, -523, -662},
+ { 102, -61, 141, 112, -270, -251, -541, 25, -150},
+ { 6, -132, -356, -686, -96, -322, -522, -31, -326},
+ { -36, -209, -521, -229, 307, -132, -5, -99, -384},
+ { 60, -51, -237, -668, -973, -407, -708, -75, -172},
+ { 26, -138, -266, 111, -302, 43, -278, -356, -359},
+ { 570, 822, 496, -154, -312, -92, 137, 279, 371},
+ { -146, 368, 409, 68, 6, 77, 167, 202, 162},
+ { 633, 898, 996, 756, 662, 683, 783, 909, 996},
+ { -103, 294, 607, 415, 483, 462, 480, 431, 408},
+ { -120, -338, -612, -524, 584, 331, 92, 433, 276},
+ { -178, -293, -154, -41, 269, 100, -9, 213, 160},
+ { 830, 736, 278, 820, 1254, 686, 712, 1039, 473},
+ { -218, -304, 463, 454, 397, 273, 202, 286, 273},
+ { -232, 7, 6, -388, -472, -427, -378, -167, -100},
+ { -294, -183, 134, -47, 101, -88, -84, -117, -3},
+ { 57, 17, -202, -634, -989, -1119, -533, 176, -36},
+ { 120, -28, 23, 111, -319, 318, -22, -77, 266},
+ { -271, -464, -434, -658, -640, -385, -385, -99, -69},
+ { -198, -259, -266, -44, -39, -139, -137, 171, 66},
+ { 9, -145, -377, -846, -1000, -111, -325, 342, 135},
+ { -81, -286, -380, 192, -57, 307, 76, -24, -140},
+ { 677, 702, 247, 56, 249, 141, -105, -236, -99},
+ { 36, -39, -69, 348, 198, -93, 322, 91, -72},
+ { 503, 885, 1508, 1307, 1282, 1172, 1119, 1209, 1061},
+ { 416, 719, 989, 1227, 1001, 1052, 954, 741, 1044},
+ { -127, -376, -657, 139, 623, 223, 501, 306, 220},
+ { -113, -384, -796, 504, 438, 85, 213, -83, -194},
+ { 585, 1132, 1233, 1091, 1247, 1433, 1512, 1448, 1314},
+ { -174, -422, 7, 1155, 1089, 1182, 1003, 945, 806},
+ { 8, -126, -317, -103, -351, -695, -98, -268, -537},
+ { 33, -103, -290, 167, -39, -407, 44, -208, -375},
+ { 104, -23, -64, -291, -637, -851, -1084, -61, -112},
+ { -75, -306, -434, 218, -148, -354, -680, -133, -216},
+ { -121, -377, -718, -97, -130, -361, -156, -379, -599},
+ { -56, -254, -586, 235, 157, -214, 11, -260, -149},
+ { -124, -267, -397, -580, -593, -527, -805, -385, 346},
+ { -193, -440, -708, -351, -141, -255, -499, -147, -185},
+ { 448, 660, 494, 208, 509, 461, 338, 291, 149},
+ { -223, 88, 335, 159, 212, 191, 286, 308, 205},
+ { -31, 469, 803, 659, 619, 658, 843, 987, 1113},
+ { -171, -242, 514, 362, 295, 524, 552, 694, 585},
+ { -64, -308, -448, -21, 284, 786, 446, 289, 92},
+ { -218, -390, -7, 169, 206, 330, 352, 408, 358},
+ { -36, 702, 959, 859, 861, 1115, 1269, 1357, 1305},
+ { -133, -341, -65, 678, 417, 440, 486, 518, 780},
+ { 33, -44, -191, -344, -461, -755, -201, 217, -31},
+ { -353, -547, -44, 123, -61, -68, -79, 29, 60},
+ { 73, -57, -406, -766, -1243, -1203, 240, 400, 165},
+ { -73, -282, -601, -213, -171, -375, 332, 35, -103},
+ { -29, -207, -553, -476, -638, -908, 172, -22, -135},
+ { -192, -239, -164, -103, -111, -47, 153, 125, 110},
+ { -1, -203, -570, -1030, -1424, -535, 155, 1, 147},
+ { -333, -653, -865, -197, -158, -21, -44, 95, 108},
+ { 389, 588, 490, 33, -237, -524, -628, -136, -260},
+ { 40, -177, -462, 453, 862, 380, 131, -130, -405},
+ { 842, 1678, 1841, 1549, 1474, 1256, 1082, 905, 742},
+ { 370, 1216, 1768, 1633, 1212, 636, 22, -330, 71},
+ { -76, -281, -741, -742, 898, 619, 277, 71, -222},
+ { -32, -265, -556, -25, 994, 682, 305, 126, -165},
+ { 73, 738, 893, 968, 993, 1768, 2273, 1840, 1391},
+ { -69, -349, -585, 234, 1158, 903, 626, 510, 251},
+ { -1, -99, -272, -210, -603, -351, -540, -811, -383},
+ { -16, -230, -504, 410, 149, -205, -343, -651, -639},
+ { 103, -9, -227, -205, -562, -781, -1079, -1208, -156},
+ { 143, 63, -135, -67, -317, -602, -784, -1154, -640},
+ { -144, -391, -674, -622, -200, -254, -660, -947, -395},
+ { -40, -250, -625, 27, 543, 94, -131, -386, -673},
+ { -123, -371, -757, -451, -564, -614, -415, -711, -35},
+ { -116, -309, -593, -268, 239, -33, -338, -650, -135},
+ { 94, 251, 554, 57, -312, -423, -154, -57, 235},
+ { -268, -71, 381, 114, -44, -87, 125, 173, 133},
+ { 1513, 1714, 1238, 534, 276, 315, 461, 459, 508},
+ { -131, -19, 1149, 670, 486, 356, 309, 369, 296},
+ { -223, -501, -899, -722, -70, 6, 131, 310, 394},
+ { -99, -303, -517, 249, 64, -53, 135, -11, 453},
+ { -147, -399, -730, -401, 817, 738, 802, 749, 575},
+ { -154, -435, -739, 800, 593, 366, 529, 318, 326},
+ { -224, 45, -39, -387, -515, -518, -608, -384, -321},
+ { -315, -377, 143, -101, -113, -377, -177, -144, -12},
+ { 117, 40, -239, -651, -1051, -581, -737, -990, -328},
+ { 26, -50, -157, -23, -453, -283, -531, -546, 192},
+ { -252, -501, -743, -589, -627, -499, -328, -118, -72},
+ { -324, -494, -244, -306, -144, -177, -262, -135, -78},
+ { -36, -234, -519, -961, -1290, -314, -479, -371, -45},
+ { -95, -292, -535, -8, -300, 112, -164, -277, 198},
+ { -99, -128, 880, 836, 579, 351, 23, -95, -217},
+ { -27, -258, 124, 1011, 597, 425, 144, 7, -73},
+ { 421, 1293, 1640, 1623, 1742, 1617, 1499, 1284, 1006},
+ { -95, 752, 1680, 1569, 1618, 1436, 1200, 980, 712},
+ { -69, -300, -683, -435, 1132, 899, 504, 332, 109},
+ { -74, -323, -637, 563, 1074, 608, 371, 105, -49},
+ { -78, 831, 1194, 1110, 1378, 1481, 1492, 1365, 1217},
+ { -259, -121, 1440, 1334, 1628, 1490, 1438, 1223, 933},
+ { -82, -306, -613, -222, -378, -675, -545, -671, -845},
+ { 53, -124, -347, 422, 52, -125, -270, -529, 9},
+ { 79, -89, -320, -662, -999, -1199, -1243, -676, -297},
+ { -68, -273, -611, 137, -146, -397, -627, -845, -220},
+ { -112, -346, -797, -826, 234, -132, -188, -278, -522},
+ { -159, -405, -734, -419, 293, 74, -167, -167, 184},
+ { -153, -437, -833, -1080, -336, -472, -561, -340, -253},
+ { -169, -423, -820, -904, -131, -19, -346, -604, 31},
+ { 33, -31, 312, 62, -148, 49, -59, 564, 486},
+ { -306, -333, 194, -44, 67, 72, 147, 205, 243},
+ { -207, -49, 1360, 983, 969, 991, 1014, 1110, 973},
+ { -211, -172, 883, 627, 711, 674, 705, 798, 746},
+ { -88, -325, -763, -974, 687, 908, 514, 382, 172},
+ { -292, -612, -805, 63, 131, 270, 259, 352, 348},
+ { -235, -84, 955, 818, 1120, 1289, 1559, 1480, 1285},
+ { -180, -461, -614, 657, 691, 745, 854, 783, 713},
+ { -97, -309, -477, -614, -777, -734, -768, -526, -472},
+ { -344, -476, -35, -169, 49, -77, -150, -240, -141},
+ { -52, -268, -639, -919, -1278, -1113, -342, -333, -151},
+ { -68, -242, -585, -73, -209, -478, -159, -429, 133},
+ { -197, -499, -1005, -1268, -272, -224, -105, -67, 17},
+ { -363, -618, -414, -116, -62, 20, 10, 116, 108},
+ { -195, -475, -906, -1260, -891, -441, -277, -142, -28},
+ { -226, -519, -950, -700, -275, -266, -116, -105, 82},
+ { 404, 511, 520, 327, 17, -194, -333, -536, -586},
+ { -114, -130, 276, 237, 204, 342, 135, -16, -111},
+ { 670, 1208, 1168, 860, 742, 601, 528, 403, 309},
+ { 397, 621, 966, 752, 579, 398, 400, 329, 252},
+ { 191, 180, -137, -467, 272, 106, -95, 17, -192},
+ { -80, -290, -626, 194, 598, 196, 21, -281, 77},
+ { 510, 864, 1108, 807, 939, 902, 925, 717, 481},
+ { 137, 367, 534, 764, 670, 382, 296, 153, 84},
+ { 303, 497, 144, -85, -125, -539, -482, -464, -764},
+ { 233, 347, 68, -147, 169, -210, -242, -226, -482},
+ { 307, 422, 154, -175, -386, -722, -724, -904, -1015},
+ { 309, 308, 160, -60, -470, -420, -598, -791, -219},
+ { 68, 121, -137, -560, -146, -446, -515, -494, -729},
+ { 130, 53, -227, 46, 474, 32, -161, -192, -490},
+ { 213, 164, -71, -465, -876, -161, -456, -587, -48},
+ { 218, 117, 39, 177, -194, -88, -226, -418, 50},
+ { 210, 547, 569, 279, 121, -44, -50, 10, -84},
+ { 58, 140, 182, -5, 267, 117, 106, 211, 198},
+ { 539, 835, 913, 719, 617, 544, 591, 565, 642},
+ { 153, 559, 872, 460, 222, 108, 188, 180, 183},
+ { 158, 119, 284, -153, -271, 229, 87, 110, -57},
+ { -183, 82, 118, 21, 13, 40, 118, 191, 185},
+ { 162, 889, 654, 108, -34, 244, 488, 561, 532},
+ { 163, 56, 609, 341, 50, 329, 68, 266, 218},
+ { 100, 206, 18, -304, -107, -436, -487, -65, -306},
+ { -86, 154, 134, -30, -45, -73, -104, -80, -96},
+ { 245, 330, 10, -440, -849, -1082, 79, 40, -265},
+ { 196, 372, 272, -181, -493, -389, 275, 80, -59},
+ { 2, -12, -246, -505, -100, -436, 21, -187, -431},
+ { -221, -48, 36, -271, -186, -147, -109, 26, 71},
+ { 213, 140, 72, -351, -620, -84, -363, 69, 46},
+ { 91, 167, -3, -95, -99, -105, -48, 114, 147},
+ { 259, 249, 172, 607, 406, 52, 59, -189, -320},
+ { 115, -85, -54, 574, 128, 226, -59, -253, 130},
+ { -62, 1033, 1308, 1035, 1127, 1098, 1029, 961, 823},
+ { 39, 364, 757, 940, 728, 660, 659, 583, 770},
+ { -115, -338, -760, -471, 394, 37, 441, 178, 6},
+ { -57, -305, -525, 796, 453, 188, -4, -114, 248},
+ { 71, 444, 797, 731, 1096, 1157, 1222, 1029, 811},
+ { 135, 359, 551, 425, 749, 815, 874, 704, 502},
+ { 132, 247, 0, -206, -449, -750, -258, -514, -633},
+ { 248, 249, 91, 121, -195, -499, -90, -282, -435},
+ { 78, 20, -277, -623, -983, -1224, -415, -458, -639},
+ { 347, 509, 208, -179, -464, -728, -76, -237, -486},
+ { -103, -343, -756, -713, -265, -609, -191, -398, -636},
+ { -121, -383, -749, 567, 252, -36, -354, -417, -50},
+ { 204, 100, -149, -650, -1081, -47, -7, -263, 111},
+ { -46, -180, -267, -324, -562, -394, -692, 398, 292},
+ { 482, 670, 683, 624, 442, 165, 116, 36, -149},
+ { 108, 247, 291, 247, 355, 122, 109, 224, 296},
+ { -14, 945, 990, 801, 755, 815, 847, 913, 892},
+ { 292, 349, 725, 482, 388, 329, 429, 620, 667},
+ { -34, 197, 213, -127, 84, 494, 620, 575, 375},
+ { 126, 207, 172, 167, 362, 202, 296, 395, 455},
+ { -6, 250, 539, 467, 636, 801, 1149, 1287, 1118},
+ { 27, 240, 369, 280, 440, 411, 634, 892, 953},
+ { 159, 170, -58, -395, -797, -690, 77, -211, -334},
+ { -5, -28, -13, -74, -335, -603, 300, 88, -205},
+ { 82, -33, -364, -698, -1203, -1153, 110, -146, -289},
+ { 113, 1, -243, -588, -994, -496, 414, 160, 42},
+ { -56, -247, -440, -693, -996, -479, 11, -178, -357},
+ { -151, -353, -327, -211, -340, 141, 65, 425, 453},
+ { 34, -169, -455, -932, -1215, 138, 499, 256, 324},
+ { 68, 139, -15, -547, -478, 17, 306, 502, 481},
+ { -32, -134, 445, 129, -143, -244, -503, -507, -599},
+ { 61, -140, -345, 496, 458, -2, 20, -227, -514},
+ { 394, 1765, 1666, 1339, 1117, 806, 642, 479, 380},
+ { 215, 519, 920, 1053, 1090, 791, 528, 290, 155},
+ { -54, -233, -647, -602, 639, 294, -2, -167, -442},
+ { -78, -315, -791, -113, 820, 403, 158, -116, -356},
+ { 529, 1851, 2003, 1228, 622, -41, -416, 344, 819},
+ { -105, -379, -236, 1224, 893, 749, 568, 356, 214},
+ { -17, -199, -144, 50, -283, -247, -578, -846, -1087},
+ { 69, -11, -381, -206, 209, -284, -387, -416, -716},
+ { 39, -5, -145, -374, -682, -909, -1074, -1169, -1066},
+ { 287, 226, 67, -221, -662, -171, -421, -642, -707},
+ { -132, -348, -538, -448, -20, -4, -354, -748, -933},
+ { 4, -75, -289, -598, 317, 52, -208, -297, -559},
+ { -88, -264, -358, -589, -631, -248, -523, -822, -1071},
+ { 70, -8, 54, -314, -515, 92, -146, -274, -493},
+ { 199, 62, 391, 158, -141, 71, -219, -203, -207},
+ { 152, 40, 329, 162, -29, 48, -149, 108, 127},
+ { 635, 1058, 883, 492, 372, 312, 317, 274, 241},
+ { 267, 722, 1256, 882, 625, 248, 8, -81, -60},
+ { -58, -138, -291, -600, -12, -2, -39, 147, 117},
+ { -107, -345, -513, 459, 76, 92, -272, 388, 262},
+ { 362, 516, 203, -409, -716, -831, -331, 185, 209},
+ { -117, -391, -298, 671, 292, 538, 257, 166, -38},
+ { -102, -319, -194, -283, -573, -262, -579, -219, -444},
+ { -235, 78, 11, -168, -101, -229, -263, -321, -123},
+ { 70, 50, -170, -599, -996, -588, -263, -516, -455},
+ { 394, 363, 229, -136, -538, 21, -183, -348, -201},
+ { -124, -368, -640, -879, -847, -209, -409, -494, -515},
+ { -127, -341, -541, -425, -510, -10, -252, -473, -291},
+ { 84, -69, -201, -676, -868, 103, -311, -132, -320},
+ { 5, -173, -188, -297, -628, 197, -57, 7, -11},
+ { 49, -160, 56, 558, 111, 33, -311, -440, -463},
+ { -1, -246, -307, 862, 453, 139, -170, -355, -232},
+ { 279, 966, 1642, 1478, 1463, 1123, 795, 525, 339},
+ { -197, -38, 1702, 1331, 1252, 950, 692, 504, 426},
+ { -108, -344, -861, -1172, 444, 354, 88, -46, -220},
+ { -53, -321, -494, 1113, 744, 364, 198, -34, -75},
+ { 457, 955, 1177, 1214, 1427, 1457, 1345, 917, 539},
+ { -69, 199, 897, 1140, 1343, 1183, 977, 742, 522},
+ { 122, 44, -269, 27, -155, -562, -307, -590, -773},
+ { 154, 42, -160, 252, -129, -305, -471, -733, -371},
+ { 135, 185, -82, -416, -722, -913, -504, -743, -880},
+ { 149, 214, -84, -329, -680, -835, -426, -661, -81},
+ { -128, -380, -735, -998, -337, 17, -182, -467, -697},
+ { -84, -290, -510, -592, 13, 440, 154, -38, -279},
+ { 70, -61, -246, -727, -1047, -80, -381, -535, -704},
+ { 178, -2, -146, -670, -938, 482, 138, 63, 65},
+ { -11, 15, 772, 443, 142, -20, -209, -126, -161},
+ { -32, -249, 95, 552, 124, 30, -343, 82, -86},
+ { 148, 751, 1515, 1105, 867, 606, 474, 448, 399},
+ { -163, -257, 899, 1097, 906, 751, 502, 390, 294},
+ { -51, -258, -447, -806, -368, 763, 464, 364, 183},
+ { -166, -374, -367, 87, 35, 399, 418, 856, 833},
+ { -205, -310, 588, 778, 785, 1065, 1118, 1245, 1157},
+ { -173, -312, 107, 345, 400, 790, 870, 1113, 1001},
+ { -7, -120, -387, -410, -614, -943, -226, -384, -491},
+ { -203, -288, -51, -331, -90, -178, -408, -573, -338},
+ { 56, -29, -273, -627, -1041, -798, -247, -467, 148},
+ { 66, -2, -205, -205, -575, -349, -57, -352, -58},
+ { -45, -225, -471, -924, -497, 77, -32, 44, -135},
+ { -277, -491, -497, -502, -424, -202, -137, 77, 96},
+ { 26, -179, -469, -1008, -1260, 262, -35, -132, -259},
+ { -66, -232, -447, -533, -789, -191, -100, -267, 364}
+};
+
+static const int16_t dico2_isf[256][7] = {
+ { 1357, 1313, 1136, 784, 438, 181, 145},
+ { 636, 648, 667, 568, 442, 217, 362},
+ { 427, 440, 674, 524, 332, 117, -417},
+ { 121, 295, 468, 465, 230, 44, -221},
+ { -147, -240, 149, 80, 390, 278, 106},
+ { -418, -556, 552, 511, 235, 144, -95},
+ { 43, 193, 274, 150, 67, 34, -273},
+ { -43, -126, 171, 416, 282, 63, -354},
+ { -372, -86, -344, -108, -94, -182, -89},
+ { -600, -840, -200, 465, 258, -11, -253},
+ { -48, 329, 97, -290, -543, -795, -354},
+ { -570, -117, 187, 10, -133, -416, -76},
+ { -618, -129, -247, -371, 45, -76, 277},
+ { -1022, -1079, 126, 474, 254, 127, 52},
+ { -281, 76, -167, -361, -283, -551, -283},
+ { -119, -52, -1, 134, -32, -204, -415},
+ { 1064, 827, 637, 684, 464, 209, 12},
+ { 482, 416, 449, 371, 335, 294, 194},
+ { 719, 576, 365, 135, 113, 91, -199},
+ { 298, 176, 493, 366, 194, 163, 36},
+ { -35, -236, -259, -36, -4, 99, 152},
+ { -98, -306, -27, 228, 90, 111, -86},
+ { 91, 13, -211, -258, -106, 86, -64},
+ { 73, -35, -57, -31, 162, 35, -192},
+ { -109, -335, -629, -66, -61, -128, 322},
+ { -495, -669, -728, 193, 31, -220, 122},
+ { 324, 95, -89, -91, -409, -710, -154},
+ { 0, -234, 92, 33, -343, -609, -220},
+ { -343, -408, -476, -655, -153, 82, 222},
+ { -490, -745, -255, 49, -48, 135, -127},
+ { 119, -67, -328, -390, -272, -545, -56},
+ { -57, -130, -10, -7, -164, -47, -22},
+ { 984, 1064, 961, 568, 210, -27, 16},
+ { 811, 691, 754, 514, 224, -35, 166},
+ { 662, 704, 618, 386, 57, -211, -257},
+ { 510, 359, 418, 393, 91, -144, -18},
+ { -193, -31, -27, 223, 89, -143, 24},
+ { -112, -98, 471, 319, 185, 3, 175},
+ { 252, 146, -47, 272, 48, -211, -234},
+ { 146, 69, 203, 364, 68, -52, 51},
+ { -259, -478, -697, -349, -758, -501, 63},
+ { -501, -769, -289, 79, -311, -497, -106},
+ { 251, 53, -235, -469, -895, -884, 145},
+ { -416, -551, 140, -133, -523, -775, 44},
+ { -326, -423, -713, -497, -86, -431, 99},
+ { -757, -772, -160, -76, -46, -32, 379},
+ { 85, -35, -200, -401, -663, -1040, -247},
+ { -180, -330, -92, -376, 27, -183, -110},
+ { 1279, 1086, 781, 502, 324, 164, 157},
+ { 682, 466, 449, 277, 146, 28, 409},
+ { 635, 472, 390, 107, -232, -538, -139},
+ { 196, 396, 332, 213, 209, -29, -81},
+ { 150, -95, -312, 76, -77, -320, -50},
+ { 46, 9, 47, 175, 139, 30, 384},
+ { 218, 206, -24, -250, -96, -276, -183},
+ { 26, 119, 38, 14, -4, -133, -52},
+ { -477, -614, -987, -715, -631, -813, 200},
+ { -744, -1009, -1065, -745, -631, -171, 18},
+ { -137, -251, -483, -613, -980, -1203, 12},
+ { -605, -767, -562, -686, -1088, -515, 58},
+ { -202, -428, -782, -1072, -96, -234, -179},
+ { -480, -709, -1070, -897, -131, -92, 321},
+ { -145, -193, -512, -729, -572, -765, -210},
+ { -331, -585, -525, -631, -281, -208, -303},
+ { 1165, 1104, 939, 828, 716, 426, 155},
+ { 6, -109, 820, 778, 415, 113, -27},
+ { 381, 339, 314, 265, 121, -9, -474},
+ { -373, 47, 584, 442, 99, -231, -113},
+ { -496, -38, -285, 262, 305, 170, 4},
+ { -587, -556, 69, 66, 471, 354, 13},
+ { -138, 70, -18, 106, 67, 167, -302},
+ { -445, -141, 185, 191, 151, 83, -133},
+ { -257, -521, -720, -198, 134, -46, -182},
+ { -819, -1168, -777, 512, 359, 95, -113},
+ { 137, -2, -74, -138, -401, -114, -371},
+ { -242, -466, 204, 223, -31, -212, -192},
+ { -532, -637, -466, -686, 256, 277, -139},
+ { -1141, -1244, -381, -75, -54, 14, 88},
+ { -311, 115, -143, -499, -343, 124, -416},
+ { -616, -147, -135, 43, -4, 121, -369},
+ { 835, 783, 641, 390, 355, 350, 64},
+ { 72, 194, 443, 467, 436, 219, 372},
+ { 464, 369, 192, 4, -156, -72, -226},
+ { 57, 206, 303, 205, 188, 101, 265},
+ { -40, -205, -488, -184, 276, 64, -26},
+ { -217, -433, -297, 137, 328, 308, -289},
+ { 378, 81, -308, -465, 57, -37, 227},
+ { -100, 24, -36, -151, 199, 8, 143},
+ { -426, -697, -1059, -133, 388, 161, 321},
+ { -644, -1023, -1271, 39, 66, -123, 70},
+ { 372, 177, -173, -556, -553, -304, -189},
+ { -117, -369, -425, -122, -462, -152, -73},
+ { -649, -850, -1189, -767, 497, 360, 222},
+ { -798, -1139, -1455, -190, 430, 234, 179},
+ { 42, -94, -405, -692, 38, -202, -246},
+ { -169, -366, -290, -88, -64, 32, -292},
+ { 1010, 923, 938, 710, 465, 230, 342},
+ { 217, 300, 1054, 675, 68, -458, -179},
+ { 78, 453, 316, 18, -237, -496, -243},
+ { 167, 21, 424, 215, -91, -303, -170},
+ { -290, -81, -70, -67, 40, 54, -59},
+ { -353, -427, -90, 53, 94, 9, 54},
+ { -28, 318, 283, 15, -240, -58, 79},
+ { -75, -121, 229, 35, 58, 6, -133},
+ { -351, -514, -744, -834, -705, -137, 164},
+ { -1124, -1388, -1055, -230, -73, 40, 36},
+ { -163, -233, -532, -785, -1170, -697, 96},
+ { -788, -959, -246, -430, -624, -165, -8},
+ { -856, -540, -630, -907, -337, -70, 76},
+ { -937, -1042, -659, -733, -208, 199, -26},
+ { -523, 78, -98, -501, -869, -890, -81},
+ { -624, -703, -45, -348, -25, 87, -186},
+ { 1005, 823, 546, 249, 90, -22, 207},
+ { 298, 397, 381, 319, 200, 62, 303},
+ { 473, 379, 133, -247, -632, -441, 75},
+ { 284, 208, 391, 115, -25, 44, 95},
+ { -72, 79, -95, -63, -129, -293, 203},
+ { -164, -349, 115, 122, 69, -1, 378},
+ { 348, 170, 99, 58, -179, -302, 188},
+ { -190, -2, 150, 23, -51, -11, 216},
+ { -615, -863, -1090, -1427, -802, -48, -6},
+ { -961, -1276, -1548, -727, -58, 56, 223},
+ { -124, -255, -561, -988, -1277, -148, -82},
+ { -480, -660, -891, -1191, -1339, -325, 20},
+ { -621, -917, -1296, -1350, 264, 289, 50},
+ { -844, -1022, -1345, -1329, -293, 46, 278},
+ { -260, -468, -829, -1176, -533, -560, -78},
+ { -215, -484, -822, -1233, -791, 15, -138},
+ { 1301, 1317, 1262, 1048, 716, 357, -64},
+ { 578, 824, 925, 802, 630, 362, 102},
+ { 470, 925, 767, 514, 327, 190, -112},
+ { 225, 492, 495, 437, 598, 384, -45},
+ { 43, 82, -42, 175, 519, 342, -64},
+ { -304, -154, 159, 576, 403, 221, 327},
+ { 214, 244, 122, -62, 312, 92, -160},
+ { 218, 208, 310, 268, 306, 323, -199},
+ { -285, -269, -79, -124, -143, -153, 236},
+ { -205, -384, -426, 344, 59, -185, -184},
+ { -272, 247, 126, -210, -518, -468, 78},
+ { -99, -120, 502, 160, -280, -557, 304},
+ { -423, -17, -283, -443, 215, 212, -140},
+ { -564, -684, -228, 510, 361, 130, 323},
+ { -428, 335, 98, -65, 36, -215, -246},
+ { -362, 51, 364, -16, -234, 150, -165},
+ { 914, 883, 751, 653, 676, 464, -153},
+ { 631, 545, 535, 720, 596, 360, -81},
+ { 783, 712, 512, 439, 341, 251, -391},
+ { 497, 417, 249, 372, 295, 173, -193},
+ { 128, -110, -385, 93, 39, 173, -231},
+ { 216, -59, -253, 462, 389, 154, 69},
+ { 455, 270, -4, -337, -49, 233, -322},
+ { 307, 143, 53, 218, 128, 236, -156},
+ { -37, -186, -240, -411, -110, 9, 399},
+ { -140, -365, -628, 258, 380, 214, 277},
+ { 131, 454, 177, -285, -520, 108, -214},
+ { 77, -141, 201, -123, -490, -131, 60},
+ { -14, -194, -521, -741, 273, 362, -33},
+ { -362, -566, -287, -228, 161, 237, 317},
+ { -269, 195, -75, -375, -204, 11, 77},
+ { -128, -264, -156, -223, -475, 265, 27},
+ { 1238, 1147, 916, 689, 432, 210, -280},
+ { 800, 664, 879, 726, 411, 160, -164},
+ { 454, 686, 536, 275, 147, 46, 111},
+ { 303, 486, 512, 355, 241, 181, -69},
+ { 79, 92, 29, 147, 233, 52, 17},
+ { -171, 289, 131, 439, 271, 3, -10},
+ { 413, 241, 144, 174, 155, -2, 14},
+ { 58, 217, 247, 219, 149, 175, -18},
+ { 228, -8, -240, -206, -513, -191, 202},
+ { -96, -272, -454, 33, -300, -575, 46},
+ { -10, -108, -246, -347, -770, -535, 9},
+ { -326, -430, -61, -321, -704, -299, 201},
+ { -1, -280, -603, -419, -185, 18, -36},
+ { -516, -522, -379, -291, -181, -97, 27},
+ { -159, -313, -525, -224, -510, -831, -197},
+ { -292, -459, -59, -310, -562, -143, -351},
+ { 1066, 912, 631, 389, 207, 86, -224},
+ { 596, 512, 596, 505, 314, 122, -48},
+ { 787, 861, 441, -93, -303, 33, -190},
+ { 257, 469, 337, 51, 15, 298, -93},
+ { 295, 73, -119, 25, 36, 23, 108},
+ { -28, -3, -32, 114, 21, 185, 107},
+ { 482, 305, 15, -279, -319, 52, 96},
+ { 226, 46, 115, 72, -136, 133, -125},
+ { 18, -207, -559, -590, -503, -482, 321},
+ { -571, -789, -951, -172, -441, -538, 113},
+ { 181, 14, -310, -641, -1001, -202, 159},
+ { -136, -393, -433, -513, -911, -144, -22},
+ { 72, -265, -706, -954, -159, 53, 332},
+ { -338, -591, -852, -383, -395, 56, 44},
+ { 43, -158, -464, -897, -631, -157, -294},
+ { -161, -128, -328, -573, -483, -125, 11},
+ { 1017, 906, 1051, 1005, 679, 341, -102},
+ { 359, 334, 1567, 1314, 723, 105, 10},
+ { -65, 726, 529, 301, 220, 43, -273},
+ { -510, 436, 719, 566, 358, 179, 114},
+ { -560, 298, 133, -120, 342, 225, 14},
+ { -899, -101, 217, 617, 400, 146, -58},
+ { -41, 352, 82, -196, 39, 121, -167},
+ { -212, 59, 447, 284, 423, 250, -169},
+ { -371, -484, -596, 30, -41, 249, 22},
+ { -372, -650, -794, 477, 445, 216, -79},
+ { -352, 275, 17, -443, -929, 92, 19},
+ { -699, -696, 431, 264, -49, -310, 182},
+ { -978, -217, -430, -400, 101, 261, 72},
+ { -929, -889, -357, -13, 463, 378, 236},
+ { -826, 56, 30, -299, -360, -128, -51},
+ { -878, -299, -111, 75, 65, 36, 3},
+ { 817, 368, -25, 354, 697, 591, -173},
+ { 309, 212, 222, 751, 484, 140, -56},
+ { 593, 379, 70, -8, 258, 180, 110},
+ { 165, -46, 255, 297, 219, 273, 105},
+ { 160, -70, -358, -181, 379, 330, 319},
+ { -238, -369, -198, 740, 580, 319, -143},
+ { 201, 109, -202, -456, 328, 276, -141},
+ { 203, 170, 111, 42, 207, 360, 188},
+ { -345, -399, -513, -233, 650, 422, 81},
+ { -635, -961, -1220, 463, 539, 204, 209},
+ { 202, -25, -194, -498, -787, 193, -143},
+ { -449, -538, 195, -106, -331, 68, 62},
+ { -228, -477, -840, -576, 317, 128, 283},
+ { -671, -937, -807, -114, 391, 335, -62},
+ { 246, 2, -314, -679, -303, 180, -88},
+ { -107, -272, 90, -198, -28, 290, -112},
+ { 885, 1149, 1021, 712, 496, 281, -83},
+ { 269, 492, 787, 643, 347, 70, 124},
+ { 336, 636, 499, 92, -229, -179, 191},
+ { 26, 402, 564, 340, 149, -11, 135},
+ { -440, 561, 470, 204, -72, -186, 140},
+ { -720, 14, 355, 229, 68, -133, 465},
+ { 110, 310, 103, 12, 106, 29, 158},
+ { -178, 113, 161, 142, 121, 115, 27},
+ { -651, -414, -645, -152, -164, -13, -429},
+ { -639, -944, -681, -104, -81, 52, -189},
+ { -663, -164, -316, -683, -954, -205, -83},
+ { -609, -669, -172, -517, -694, 283, -80},
+ { -646, -152, -383, -678, -246, -40, -143},
+ { -747, -796, -745, -390, -98, 43, 275},
+ { -599, -199, -398, -433, -436, -538, 31},
+ { -1107, -568, -376, -265, -126, -21, 1},
+ { 847, 573, 308, 392, 305, 101, 55},
+ { 273, 293, 201, 267, 346, 201, 123},
+ { 727, 480, 226, 2, -65, -138, 164},
+ { 273, 208, 173, 292, 12, 253, 174},
+ { 340, 207, 180, 88, 116, 46, 475},
+ { -460, -166, -30, 13, 110, 173, 396},
+ { 137, 88, 43, -137, -94, 34, 284},
+ { 96, -14, 226, 40, 63, 70, 130},
+ { -467, -735, -1012, -1174, -307, 305, -67},
+ { -612, -920, -1146, -567, -8, 92, -25},
+ { -182, -271, -492, -754, -857, 287, -75},
+ { -494, -787, -689, -683, -709, 137, -326},
+ { -288, -550, -903, -1105, 334, 321, -62},
+ { -354, -653, -834, -445, 1, 377, -152},
+ { -162, -306, -608, -937, -297, 247, -192},
+ { -234, -477, -244, -488, -266, 342, -332}
+};
+
+static const int16_t dico21_isf[64][3] = {
+ { 329, 409, 249}, { -33, 505, 160},
+ { -29, -14, 582}, { -262, 127, 354},
+ { 145, 237, 175}, { -152, 245, 122},
+ { 27, 42, 340}, { -84, -93, 311},
+ { 285, 222, -156}, { 47, -43, -504},
+ { 234, 121, 385}, { 104, -317, 45},
+ { 176, 195, 8}, { 104, -59, -94},
+ { 177, 53, 192}, { -34, -127, 152},
+ { 570, 277, -34}, { -67, -329, -639},
+ { -157, -272, 462}, { -177, -462, 198},
+ { 322, 179, 115}, { -386, 171, 19},
+ { 19, -12, 195}, { -120, -252, 201},
+ { 304, 36, -336}, { -128, -221, -380},
+ { 171, -185, 296}, { -242, -312, 23},
+ { 198, 39, 16}, { -3, -177, -111},
+ { 111, -93, 76}, { -92, -223, 4},
+ { 177, 406, -44}, { -168, 380, -149},
+ { -4, 273, 331}, { -420, 513, 277},
+ { 21, 247, 47}, { -58, 131, -2},
+ { -3, 134, 180}, { -145, 40, 175},
+ { 189, 74, -145}, { -27, -45, -325},
+ { 370, -114, -21}, { -83, -415, -173},
+ { 77, 95, -51}, { -40, -30, -67},
+ { 71, 88, 86}, { -35, -98, 14},
+ { 69, 197, -334}, { -196, 79, -231},
+ { -348, -137, 218}, { -352, -89, -85},
+ { 47, 201, -130}, { -165, 37, -15},
+ { -43, 3, 86}, { -161, -108, 79},
+ { 83, 21, -237}, { -81, -149, -238},
+ { 150, -186, -251}, { -186, -249, -162},
+ { -19, 66, -139}, { -26, -50, -181},
+ { 24, 11, 0}, { -130, -105, -98}
+};
+
+static const int16_t dico22_isf[128][3] = {
+ { -127, 310, 42}, { -242, 197, 5},
+ { -151, 84, -17}, { -214, 127, -149},
+ { -247, -131, 159}, { -268, -267, -95},
+ { -217, 1, -79}, { -271, -80, -185},
+ { -45, 436, 159}, { 165, 199, 391},
+ { -33, 81, 187}, { -66, -42, 355},
+ { -298, -57, 343}, { -108, -537, 226},
+ { -144, -23, 193}, { 176, -402, 87},
+ { 53, 296, 25}, { -84, 253, -104},
+ { -58, 105, -126}, { -169, 174, -314},
+ { -48, 44, -294}, { -164, -417, -242},
+ { -139, 3, -194}, { -155, -207, -211},
+ { 119, 322, 213}, { 333, 50, 380},
+ { 237, 247, -2}, { 466, -16, 201},
+ { 238, -255, -107}, { 67, -440, -149},
+ { 122, -88, -139}, { 88, -247, -73},
+ { -41, 231, 167}, { -62, 155, 16},
+ { -65, 16, 77}, { -68, -2, -63},
+ { -151, -300, 160}, { -18, -333, 54},
+ { -56, -94, 5}, { 2, -190, 14},
+ { 92, 148, 209}, { 108, 9, 272},
+ { 108, 35, 110}, { 142, -85, 145},
+ { 47, -157, 279}, { 3, -320, 246},
+ { 43, -72, 68}, { 86, -217, 135},
+ { 36, 140, 79}, { 56, 175, -49},
+ { 26, 45, 3}, { 73, 55, -101},
+ { 109, -183, -242}, { -4, -283, -242},
+ { 48, -68, -48}, { -6, -153, -122},
+ { 161, 196, 96}, { 232, 80, 190},
+ { 165, 97, 11}, { 258, -31, 71},
+ { 267, -77, -91}, { 311, -209, 87},
+ { 152, -14, -22}, { 150, -149, 9},
+ { -324, 557, 187}, { -384, 307, 46},
+ { -251, 27, 77}, { -365, 77, -52},
+ { -482, -84, 160}, { -424, -515, -64},
+ { -294, -120, -4}, { -476, -116, -109},
+ { -97, 318, 365}, { 106, 627, 445},
+ { -190, 120, 287}, { -146, 65, 619},
+ { -427, 242, 363}, { -361, -371, 432},
+ { -347, 102, 168}, { -629, 195, -14},
+ { -65, 476, -47}, { -297, 320, -168},
+ { -55, 356, -264}, { -391, 82, -286},
+ { -51, -31, -556}, { -178, -399, -586},
+ { -205, -49, -360}, { -343, -238, -337},
+ { 220, 457, 58}, { 561, 467, 259},
+ { 340, 270, -168}, { 450, 77, -280},
+ { 60, 167, -413}, { 133, -252, -492},
+ { 216, 157, -290}, { 282, 0, -495},
+ { -226, 293, 183}, { -157, 135, 122},
+ { -158, -59, 39}, { -133, -118, -97},
+ { -332, -309, 113}, { -160, -425, -6},
+ { -149, -211, 24}, { -80, -277, -90},
+ { -11, 125, 338}, { 130, -71, 465},
+ { 5, -45, 184}, { 237, -95, 253},
+ { -139, -197, 297}, { -19, -300, 511},
+ { -63, -152, 139}, { 250, -289, 336},
+ { 124, 339, -150}, { 34, 176, -208},
+ { 171, 166, -116}, { 94, 38, -229},
+ { 75, -65, -339}, { -78, -205, -385},
+ { 0, -30, -163}, { -56, -110, -242},
+ { 321, 244, 194}, { 505, 238, -1},
+ { 317, 116, 65}, { 309, 88, -74},
+ { 452, -51, -50}, { 334, -217, -290},
+ { 211, 41, -152}, { 238, -55, -260}
+};
+
+static const int16_t dico23_isf[128][3] = {
+ { -10, 151, 359}, { 136, 298, 223},
+ { 255, -104, 290}, { 423, 6, 183},
+ { -270, -269, -98}, { -52, -82, 13},
+ { -82, -274, -97}, { 90, -246, -72},
+ { -299, -70, 421}, { -88, 365, 430},
+ { 187, -318, 381}, { 380, 37, 488},
+ { -373, -316, 79}, { -308, -101, 5},
+ { -135, -451, 8}, { 72, -421, -154},
+ { 180, 170, -121}, { 62, 177, -40},
+ { 326, 80, -105}, { 248, 263, -5},
+ { -168, -181, -221}, { -2, -23, -158},
+ { -14, -149, -121}, { 119, -91, -147},
+ { 119, 332, -153}, { 49, 303, 34},
+ { 442, -55, -69}, { 217, 454, 58},
+ { -359, -187, -375}, { -42, 50, -274},
+ { -8, -267, -249}, { 85, -86, -346},
+ { -77, -40, 345}, { 89, 134, 219},
+ { 156, -80, 160}, { 108, 40, 116},
+ { -158, -206, 29}, { 5, -32, 175},
+ { -65, -158, 146}, { 55, -78, 73},
+ { -114, -222, 353}, { -47, 81, 211},
+ { 49, -151, 268}, { 105, 4, 302},
+ { -263, -132, 183}, { -151, -28, 201},
+ { -177, -307, 166}, { 101, -221, 130},
+ { 74, 58, -98}, { 32, 44, 13},
+ { 194, 30, -142}, { 170, 96, 8},
+ { -136, -119, -91}, { -65, 8, -55},
+ { 3, -188, 12}, { 45, -63, -49},
+ { 149, -21, -19}, { 24, 144, 95},
+ { 254, -22, 60}, { 161, 196, 96},
+ { -158, -61, 48}, { -70, 33, 82},
+ { -23, -321, 58}, { 155, -147, 5},
+ { -364, 328, 77}, { -21, 453, 173},
+ { -108, 82, 630}, { 367, 263, 208},
+ { -300, -62, -176}, { -205, 143, -158},
+ { -169, -410, -264}, { 257, -269, -100},
+ { -636, 289, -2}, { -292, 627, 173},
+ { -382, -363, 387}, { 248, 524, 447},
+ { -521, -111, -107}, { -395, 118, -274},
+ { -343, -680, -125}, { -172, -447, -663},
+ { 75, 148, -367}, { -79, 263, -94},
+ { 249, 148, -286}, { 380, 271, -162},
+ { -142, -4, -186}, { -57, 111, -125},
+ { -35, -108, -254}, { 100, 29, -242},
+ { -80, 303, -264}, { -78, 464, -57},
+ { 248, -22, -494}, { 661, 662, 44},
+ { -193, -40, -330}, { -178, 145, -337},
+ { -90, -199, -400}, { -40, -23, -498},
+ { -192, 114, 315}, { -41, 244, 190},
+ { 88, -97, 485}, { 241, 80, 212},
+ { -246, 40, 87}, { -156, 147, 134},
+ { -2, -334, 239}, { 308, -203, 110},
+ { -459, 251, 422}, { -218, 310, 228},
+ { -86, -346, 654}, { 184, 175, 425},
+ { -481, -63, 169}, { -349, 117, 188},
+ { -125, -560, 310}, { 158, -416, 94},
+ { 46, 171, -192}, { -63, 157, 14},
+ { 256, -35, -271}, { 322, 123, 53},
+ { -214, 4, -76}, { -156, 86, -18},
+ { 128, -197, -232}, { 265, -90, -98},
+ { -308, 332, -145}, { -131, 308, 58},
+ { 509, 59, -339}, { 562, 196, -14},
+ { -378, 100, -47}, { -234, 202, 1},
+ { 104, -270, -493}, { 319, -210, -325}
+};
+
+static const int16_t dico24_isf[32][3] = {
+ { -79, -89, -4}, { -171, 77, -211},
+ { 160, -193, 98}, { 120, -103, 323},
+ { 32, -22, -129}, { 72, 78, -268},
+ { 182, -76, -66}, { 309, 99, -145},
+ { -229, -157, -84}, { -383, 98, -71},
+ { -90, -352, 12}, { -284, -178, 178},
+ { -65, -125, -166}, { -87, -175, -351},
+ { 42, -198, -48}, { 154, -140, -243},
+ { -77, 18, 108}, { -39, 355, 91},
+ { 87, 8, 155}, { -4, 158, 239},
+ { 128, 95, -54}, { 7, 246, -124},
+ { 258, 15, 89}, { 206, 216, 98},
+ { -201, 9, 18}, { -312, 233, 204},
+ { -39, -174, 155}, { -144, -9, 284},
+ { -57, 70, -69}, { -157, 187, 18},
+ { 54, -30, 23}, { 24, 135, 55}
+};
+
+static const int16_t dico25_isf[32][4] = {
+ { 169, 142, -119, 115}, { 206, -20, 94, 226},
+ { -106, 313, -21, 16}, { -62, 161, 71, 255},
+ { -89, 101, -185, 125}, { 72, -30, -201, 344},
+ { -258, 33, -8, 81}, { -104, -154, 72, 296},
+ { 144, -68, -268, -25}, { 81, -78, -87, 106},
+ { 22, 155, -186, -119}, { -46, -28, 27, 91},
+ { -114, -37, -175, -33}, { -94, -222, -189, 122},
+ { -132, -119, -191, -270}, { -172, -173, 18, -43},
+ { 279, 135, -42, -128}, { 187, -86, 229, -138},
+ { 159, 240, 140, 46}, { 69, 25, 227, 77},
+ { 21, 115, 13, 8}, { 68, -248, 126, 81},
+ { -150, 137, 207, -9}, { -154, -133, 289, 67},
+ { 143, -37, -86, -326}, { 180, -32, 19, -23},
+ { 26, 168, 116, -233}, { -32, -26, 118, -78},
+ { 3, -8, -45, -115}, { 57, -215, -54, -83},
+ { -209, 112, -22, -167}, { -91, -151, 168, -262}
+};
+
+static const int16_t dico21_isf_36b[128][5] = {
+ { -52, -96, 212, 315, -73}, { 82, -204, 363, 136, -197},
+ { -126, -331, 183, 218, 143}, { -49, -41, 557, 230, 72},
+ { 2, -73, 163, 377, 221}, { 133, 111, 278, 215, -110},
+ { -102, -20, 284, 113, 273}, { 84, 319, 290, 18, 85},
+ { -25, -5, 125, 132, -204}, { -38, -5, 286, -9, -356},
+ { -140, -256, 92, 117, -189}, { -144, 191, 313, 51, -98},
+ { 167, -10, 44, 247, 36}, { 381, 197, 238, 74, 6},
+ { 38, -408, 29, -3, -85}, { 92, 266, 157, -25, -200},
+ { 161, -121, 70, 84, -140}, { -16, -86, 112, -94, -189},
+ { -269, -270, 351, 107, -24}, { -68, -67, 492, -103, -155},
+ { -53, -131, 62, 122, 10}, { 135, 84, 283, -55, -120},
+ { -12, -219, 331, -81, 167}, { 220, -136, 147, -172, -42},
+ { 140, -95, -109, -88, -194}, { 0, -2, -4, -33, -381},
+ { -66, -217, 152, -186, -402}, { 244, 108, 156, -140, -395},
+ { 113, -136, -196, 110, -24}, { 214, 118, 11, -64, -131},
+ { -110, -286, -6, -332, 16}, { 94, 97, 79, -291, -205},
+ { -5, -39, -20, 252, -96}, { 76, 174, 101, 163, 61},
+ { -69, -239, -55, 399, 6}, { -115, 319, 164, 275, 196},
+ { -15, 36, -47, 331, 121}, { 226, 209, 271, 325, 184},
+ { 13, -80, -218, 471, 353}, { 288, 378, 16, -51, 251},
+ { 174, 116, 52, 149, -279}, { 235, 276, 39, 120, -48},
+ { 0, -108, -108, 241, -339}, { -93, 534, 45, 33, -87},
+ { 194, 149, -71, 405, -44}, { 409, 370, 81, -186, -154},
+ { 25, -102, -448, 124, -173}, { 22, 408, -110, -310, -214},
+ { -26, 23, -83, 114, 14}, { -110, 164, 52, 223, -82},
+ { 37, -25, -263, 306, -15}, { -466, 415, 292, 165, -18},
+ { 29, -19, -171, 155, 182}, { 179, 144, -27, 231, 258},
+ { -103, -247, -396, 238, 113}, { 375, -154, -109, -4, 156},
+ { 98, 85, -292, -5, -124}, { 116, 139, -116, -98, -294},
+ { -14, -83, -278, -117, -378}, { 106, 33, -106, -344, -484},
+ { 119, 17, -412, 138, 166}, { 384, 101, -204, 88, -156},
+ { -121, -284, -300, -1, -166}, { 280, 33, -152, -313, -81},
+ { -37, 22, 229, 153, 37}, { -60, -83, 236, -8, -41},
+ { -169, -228, 126, -20, 363}, { -235, 17, 364, -156, 156},
+ { -25, -30, 72, 144, 156}, { 153, -26, 256, 97, 144},
+ { -21, -37, 48, -65, 250}, { 63, 77, 273, -128, 124},
+ { -129, -26, 40, 9, -115}, { -6, 82, 38, -90, -182},
+ { -336, -13, 28, 158, 91}, { -30, 241, 137, -170, -17},
+ { 146, 14, -11, 33, 61}, { 192, 197, 54, -84, 85},
+ { 23, -200, -78, -29, 140}, { 122, 237, 106, -341, 136},
+ { -57, -142, -85, -16, -74}, { -59, -90, -8, -187, -20},
+ { -211, -267, 216, -179, -110}, { -50, -7, 220, -267, -70},
+ { -57, -42, -17, -15, 71}, { 32, 21, 63, -137, 33},
+ { -137, -175, 104, -68, 97}, { -67, -43, 133, -301, 221},
+ { -116, -200, -81, -92, -272}, { -64, -41, -54, -244, -220},
+ { -287, -242, -50, -87, -89}, { -245, 236, 102, -166, -295},
+ { 66, 24, -162, -71, 95}, { 66, 136, -90, -220, -36},
+ { -98, -161, -222, -188, 29}, { -18, 18, -19, -415, 9},
+ { 49, 61, 100, 39, -56}, { -111, 82, 135, -31, 52},
+ { -90, -153, -93, 189, 182}, { -214, 295, 119, -74, 284},
+ { 2, 137, 37, 47, 182}, { 92, 117, 184, -53, 373},
+ { -21, -14, -35, 136, 391}, { 146, 129, -164, -28, 333},
+ { 92, 80, -84, 100, -134}, { -8, 217, -32, 3, -47},
+ { -151, 251, -215, 142, 92}, { -224, 310, -172, -275, 98},
+ { 159, 155, -177, 112, 53}, { 205, 27, 8, -240, 192},
+ { 169, 120, -319, -201, 106}, { 11, 36, -86, -237, 455},
+ { -109, -154, -163, 174, -55}, { -38, 32, -101, -78, -59},
+ { -205, -321, -97, 69, 79}, { -310, 44, 18, -185, 34},
+ { -115, -20, -148, -39, 203}, { -29, 154, -30, -158, 166},
+ { -45, -131, -317, -24, 363}, { -165, -205, -112, -222, 265},
+ { -32, -44, -150, 54, -193}, { -6, -38, -255, -169, -115},
+ { -266, 87, -189, -36, -169}, { -60, -87, -266, -436, -170},
+ { -68, -81, -278, 24, 38}, { -23, -19, -155, -256, 141},
+ { -61, -226, -565, -175, 71}, { 9, -29, -237, -515, 263}
+};
+
+static const int16_t dico22_isf_36b[128][4] = {
+ { -298, -6, 95, 31}, { -213, -87, -122, 261},
+ { 4, -49, 208, 14}, { -129, -110, 30, 118},
+ { -214, 258, 110, -235}, { -41, -18, -126, 120},
+ { 103, 65, 127, -37}, { 126, -36, -24, 25},
+ { -138, -67, -278, -186}, { -164, -194, -201, 78},
+ { -211, -87, -51, -221}, { -174, -79, -94, -39},
+ { 23, -6, -157, -240}, { 22, -110, -153, -68},
+ { 148, -5, -2, -149}, { -1, -135, -39, -179},
+ { 68, 360, -117, -15}, { 137, 47, -278, 146},
+ { 136, 260, 135, 65}, { 61, 116, -45, 97},
+ { 231, 379, 87, -120}, { 338, 177, -272, 3},
+ { 266, 156, 28, -69}, { 260, 84, -85, 86},
+ { -266, 154, -256, -182}, { -17, -65, -304, -6},
+ { -40, 175, -151, -180}, { -27, 27, -87, -63},
+ { 121, 114, -166, -469}, { 159, -66, -323, -231},
+ { 214, 152, -141, -212}, { 137, 36, -184, -51},
+ { -282, -237, 40, 10}, { -48, -235, -37, 251},
+ { -54, -323, 136, 29}, { -88, -174, 213, 198},
+ { -390, 99, -63, -375}, { 107, -169, -164, 424},
+ { 69, -111, 141, -167}, { 74, -129, 65, 144},
+ { -353, -207, -205, -109}, { -160, -386, -355, 98},
+ { -176, -493, -20, -143}, { -252, -432, -2, 216},
+ { -90, -174, -168, -411}, { 13, -284, -229, -160},
+ { -87, -279, 34, -251}, { -75, -263, -58, -42},
+ { 420, 53, -211, -358}, { 384, -35, -374, 396},
+ { 68, -228, 323, -2}, { 167, -307, 192, 194},
+ { 459, 329, -5, -332}, { 375, 79, -7, 313},
+ { 282, -124, 200, -92}, { 271, -162, -70, 180},
+ { -157, -298, -514, -309}, { 58, -163, -546, 18},
+ { 124, -364, 167, -238}, { 83, -411, -117, 96},
+ { 140, -112, -388, -624}, { 259, -133, -317, 41},
+ { 163, -130, -64, -334}, { 226, -165, -124, -110},
+ { -466, -61, 6, 229}, { -153, 205, -145, 242},
+ { -159, 48, 195, 148}, { -58, 28, 31, 279},
+ { -303, 185, 279, -4}, { -61, 197, 59, 86},
+ { -114, 123, 168, -52}, { 35, 36, 100, 126},
+ { -407, 102, -77, -40}, { -338, -1, -342, 156},
+ { -179, 105, -34, -97}, { -185, 84, -35, 108},
+ { -133, 107, -91, -357}, { -180, 54, -229, 24},
+ { -44, 47, 47, -182}, { -66, 13, 45, 4},
+ { -339, 251, 64, 226}, { -42, 101, -350, 275},
+ { -99, 398, 142, 121}, { 111, 12, -102, 260},
+ { 0, 505, 260, -94}, { 161, 285, -96, 224},
+ { -4, 206, 314, 33}, { 167, 139, 88, 204},
+ { -235, 316, -60, -25}, { -8, -150, -312, 201},
+ { -36, 292, 61, -104}, { -40, 174, -162, 42},
+ { -21, 402, -29, -351}, { 21, 152, -360, -93},
+ { 57, 191, 212, -196}, { 76, 158, -21, -69},
+ { -328, -185, 331, 119}, { -53, 285, 56, 337},
+ { -107, -24, 405, 29}, { -18, 137, 272, 277},
+ { -255, 22, 173, -191}, { 295, 322, 325, 302},
+ { 21, -27, 332, -178}, { 119, 13, 271, 129},
+ { -455, -180, 116, -191}, { -227, 62, -148, 524},
+ { -176, -287, 282, -157}, { -243, 13, 199, 430},
+ { -59, -49, 115, -365}, { 72, -172, -137, 93},
+ { -138, -126, 141, -84}, { 5, -124, 38, -20},
+ { -258, 311, 601, 213}, { 94, 130, -61, 502},
+ { -1, -157, 485, 313}, { 146, -74, 158, 345},
+ { 276, 135, 280, -57}, { 490, 252, 99, 43},
+ { 267, -74, 429, 105}, { 278, -23, 119, 94},
+ { -542, 488, 257, -115}, { -84, -244, -438, 478},
+ { -113, -545, 387, 101}, { -95, -306, 111, 498},
+ { 95, 166, 22, -301}, { 420, -15, -58, -78},
+ { 270, 29, 122, -282}, { 160, -240, 50, -38}
+};
+
+static const int16_t dico23_isf_36b[64][7] = {
+ { 81, -18, 68, -27, -122, -280, -4},
+ { 45, -177, 209, -30, -136, -74, 131},
+ { -44, 101, -75, -88, -48, -137, -54},
+ { -245, -28, 63, -18, -112, -103, 58},
+ { -79, -6, 220, -65, 114, -35, -50},
+ { 109, -65, 143, -114, 129, 76, 125},
+ { 166, 90, -61, -242, 186, -74, -43},
+ { -46, -92, 49, -227, 24, -155, 39},
+ { 67, 85, 99, -42, 53, -184, -281},
+ { 142, -122, 0, 21, -142, -15, -17},
+ { 223, 92, -21, -48, -82, -14, -167},
+ { 51, -37, -243, -30, -90, 18, -56},
+ { 54, 105, 74, 86, 69, 13, -101},
+ { 196, 72, -89, 43, 65, 19, 39},
+ { 121, 34, 131, -82, 25, 213, -156},
+ { 101, -102, -136, -21, 57, 214, 22},
+ { 36, -124, 205, 204, 58, -156, -83},
+ { 83, -117, 137, 137, 85, 116, 44},
+ { -92, -148, -68, 11, -102, -197, -220},
+ { -76, -185, -58, 132, -26, -183, 85},
+ { -7, -31, -2, 23, 205, -151, 10},
+ { -27, -37, -5, -18, 292, 131, 1},
+ { 117, -168, 9, -93, 80, -59, -125},
+ { -182, -244, 98, -24, 135, -22, 94},
+ { 221, 97, 106, 42, 43, -160, 83},
+ { 25, -64, -21, 6, 14, -15, 154},
+ { 126, 15, -140, 150, -10, -207, -114},
+ { 79, -63, -211, -70, -28, -217, 165},
+ { 46, 38, -22, 281, 132, -62, 109},
+ { 112, 54, -112, -93, 208, 27, 296},
+ { 115, 10, -147, 41, 216, 42, -276},
+ { 50, -115, -254, 167, 117, -2, 61},
+ { 17, 144, 34, -72, -186, -150, 272},
+ { -29, -66, -89, -95, -149, 129, 251},
+ { 122, 0, -50, -234, -91, 36, 26},
+ { -105, -102, -88, -121, -236, -7, -11},
+ { -204, 109, 5, -191, 105, -15, 163},
+ { -80, 32, -24, -209, 41, 294, 70},
+ { -106, -94, -204, -118, 120, -50, -37},
+ { -82, -241, 46, -131, -29, 150, -55},
+ { 33, 155, 120, -89, -8, 7, 62},
+ { 213, 82, 61, 18, -161, 144, 152},
+ { 30, 131, 65, -87, -255, -17, -107},
+ { -8, 85, -64, 51, -162, 223, -53},
+ { -134, 261, 69, -56, 218, 72, -111},
+ { 2, 155, -113, -87, 49, 85, -28},
+ { -163, 42, -1, -196, 7, 39, -245},
+ { 14, -137, -79, 11, -160, 202, -293},
+ { -94, 33, 208, 100, 56, -44, 326},
+ { -78, -41, 232, 13, -142, 227, 80},
+ { -16, -87, 201, 33, -133, 15, -183},
+ { -58, -192, -47, 184, -128, 133, 99},
+ { -205, 11, -155, 78, 52, 72, 141},
+ { -246, 26, 99, 151, 59, 115, -64},
+ { -79, -47, -16, -14, 6, 47, -43},
+ { -72, -178, -27, 162, 112, 43, -174},
+ { -175, 238, 186, 71, -54, -188, -76},
+ { -225, 233, 39, -39, -158, 122, 44},
+ { -26, 43, 84, 130, -93, -51, 22},
+ { 3, 92, -150, 136, -182, -57, 97},
+ { -131, 179, -78, 80, 91, -165, 90},
+ { -2, 148, 15, 130, 65, 175, 117},
+ { -138, 114, -137, 132, 3, -10, -186},
+ { 140, -4, -37, 254, -62, 92, -109}
+};
+
+/** Means of ISF vectors in Q15 */
+static const int16_t isf_mean[LP_ORDER] = {
+ 738, 1326, 2336, 3578, 4596, 5662, 6711, 7730,
+ 8750, 9753, 10705, 11728, 12833, 13971, 15043, 4037
+};
+
+/** Initialization tables for the processed ISF vector in Q15 */
+static const int16_t isf_init[LP_ORDER] = {
+ 1024, 2048, 3072, 4096, 5120, 6144, 7168, 8192,
+ 9216, 10240, 11264, 12288, 13312, 14336, 15360, 3840
+};
+
+/** ISF/ISP interpolation coefficients for each subframe */
+static const float isfp_inter[4] = { 0.45, 0.8, 0.96, 1.0 };
+
+/** Coefficients for FIR interpolation of excitation vector
+ * at pitch lag resulting the adaptive codebook vector */
+static const float ac_inter[65] = {
+ 9.400024e-01,
+ 8.563843e-01, 6.322632e-01, 3.375854e-01, 5.908203e-02,
+ -1.310425e-01, -1.994019e-01, -1.585693e-01, -5.633545e-02,
+ 4.760742e-02, 1.067505e-01, 1.036987e-01, 5.206299e-02,
+ -1.519775e-02, -6.372070e-02, -7.366943e-02, -4.650879e-02,
+ -9.765625e-04, 3.820801e-02, 5.316162e-02, 4.003906e-02,
+ 9.338379e-03, -2.166748e-02, -3.778076e-02, -3.320312e-02,
+ -1.300049e-02, 1.068115e-02, 2.587891e-02, 2.630615e-02,
+ 1.379395e-02, -3.662109e-03, -1.678467e-02, -1.983643e-02,
+ -1.275635e-02, -5.493164e-04, 1.007080e-02, 1.409912e-02,
+ 1.068115e-02, 2.624512e-03, -5.371094e-03, -9.338379e-03,
+ -8.117676e-03, -3.173828e-03, 2.319336e-03, 5.615234e-03,
+ 5.554199e-03, 2.868652e-03, -6.103516e-04, -2.990723e-03,
+ -3.356934e-03, -2.014160e-03, -1.220703e-04, 1.342773e-03,
+ 1.708984e-03, 1.159668e-03, 2.441406e-04, -4.272461e-04,
+ -6.103516e-04, -4.272461e-04, -1.220703e-04, 6.103516e-05,
+ 1.220703e-04, 6.103516e-05, 0.000000e+00, 0.000000e+00
+};
+
+/** [i][j] is the number of pulses present in track j at mode i */
+static const uint8_t pulses_nb_per_mode_tr[][4] = {
+ {1, 1, 0, 0}, {1, 1, 1, 1}, {2, 2, 2, 2},
+ {3, 3, 2, 2}, {3, 3, 3, 3}, {4, 4, 4, 4},
+ {5, 5, 4, 4}, {6, 6, 6, 6}, {6, 6, 6, 6}
+};
+
+/** Tables for decoding quantized gains { pitch (Q14), fixed factor (Q11) } */
+static const int16_t qua_gain_6b[64][2] = {
+ { 1566, 1332}, { 1577, 3557},
+ { 3071, 6490}, { 4193, 10163},
+ { 4496, 2534}, { 5019, 4488},
+ { 5586, 15614}, { 5725, 1422},
+ { 6453, 580}, { 6724, 6831},
+ { 7657, 3527}, { 8072, 2099},
+ { 8232, 5319}, { 8827, 8775},
+ { 9740, 2868}, { 9856, 1465},
+ { 10087, 12488}, { 10241, 4453},
+ { 10859, 6618}, { 11321, 3587},
+ { 11417, 1800}, { 11643, 2428},
+ { 11718, 988}, { 12312, 5093},
+ { 12523, 8413}, { 12574, 26214},
+ { 12601, 3396}, { 13172, 1623},
+ { 13285, 2423}, { 13418, 6087},
+ { 13459, 12810}, { 13656, 3607},
+ { 14111, 4521}, { 14144, 1229},
+ { 14425, 1871}, { 14431, 7234},
+ { 14445, 2834}, { 14628, 10036},
+ { 14860, 17496}, { 15161, 3629},
+ { 15209, 5819}, { 15299, 2256},
+ { 15518, 4722}, { 15663, 1060},
+ { 15759, 7972}, { 15939, 11964},
+ { 16020, 2996}, { 16086, 1707},
+ { 16521, 4254}, { 16576, 6224},
+ { 16894, 2380}, { 16906, 681},
+ { 17213, 8406}, { 17610, 3418},
+ { 17895, 5269}, { 18168, 11748},
+ { 18230, 1575}, { 18607, 32767},
+ { 18728, 21684}, { 19137, 2543},
+ { 19422, 6577}, { 19446, 4097},
+ { 19450, 9056}, { 20371, 14885}
+};
+
+static const int16_t qua_gain_7b[128][2] = {
+ { 204, 441}, { 464, 1977},
+ { 869, 1077}, { 1072, 3062},
+ { 1281, 4759}, { 1647, 1539},
+ { 1845, 7020}, { 1853, 634},
+ { 1995, 2336}, { 2351, 15400},
+ { 2661, 1165}, { 2702, 3900},
+ { 2710, 10133}, { 3195, 1752},
+ { 3498, 2624}, { 3663, 849},
+ { 3984, 5697}, { 4214, 3399},
+ { 4415, 1304}, { 4695, 2056},
+ { 5376, 4558}, { 5386, 676},
+ { 5518, 23554}, { 5567, 7794},
+ { 5644, 3061}, { 5672, 1513},
+ { 5957, 2338}, { 6533, 1060},
+ { 6804, 5998}, { 6820, 1767},
+ { 6937, 3837}, { 7277, 414},
+ { 7305, 2665}, { 7466, 11304},
+ { 7942, 794}, { 8007, 1982},
+ { 8007, 1366}, { 8326, 3105},
+ { 8336, 4810}, { 8708, 7954},
+ { 8989, 2279}, { 9031, 1055},
+ { 9247, 3568}, { 9283, 1631},
+ { 9654, 6311}, { 9811, 2605},
+ { 10120, 683}, { 10143, 4179},
+ { 10245, 1946}, { 10335, 1218},
+ { 10468, 9960}, { 10651, 3000},
+ { 10951, 1530}, { 10969, 5290},
+ { 11203, 2305}, { 11325, 3562},
+ { 11771, 6754}, { 11839, 1849},
+ { 11941, 4495}, { 11954, 1298},
+ { 11975, 15223}, { 11977, 883},
+ { 11986, 2842}, { 12438, 2141},
+ { 12593, 3665}, { 12636, 8367},
+ { 12658, 1594}, { 12886, 2628},
+ { 12984, 4942}, { 13146, 1115},
+ { 13224, 524}, { 13341, 3163},
+ { 13399, 1923}, { 13549, 5961},
+ { 13606, 1401}, { 13655, 2399},
+ { 13782, 3909}, { 13868, 10923},
+ { 14226, 1723}, { 14232, 2939},
+ { 14278, 7528}, { 14439, 4598},
+ { 14451, 984}, { 14458, 2265},
+ { 14792, 1403}, { 14818, 3445},
+ { 14899, 5709}, { 15017, 15362},
+ { 15048, 1946}, { 15069, 2655},
+ { 15405, 9591}, { 15405, 4079},
+ { 15570, 7183}, { 15687, 2286},
+ { 15691, 1624}, { 15699, 3068},
+ { 15772, 5149}, { 15868, 1205},
+ { 15970, 696}, { 16249, 3584},
+ { 16338, 1917}, { 16424, 2560},
+ { 16483, 4438}, { 16529, 6410},
+ { 16620, 11966}, { 16839, 8780},
+ { 17030, 3050}, { 17033, 18325},
+ { 17092, 1568}, { 17123, 5197},
+ { 17351, 2113}, { 17374, 980},
+ { 17566, 26214}, { 17609, 3912},
+ { 17639, 32767}, { 18151, 7871},
+ { 18197, 2516}, { 18202, 5649},
+ { 18679, 3283}, { 18930, 1370},
+ { 19271, 13757}, { 19317, 4120},
+ { 19460, 1973}, { 19654, 10018},
+ { 19764, 6792}, { 19912, 5135},
+ { 20040, 2841}, { 21234, 19833}
+};
+
+/** 4-tap moving average prediction coefficients in reverse order */
+static const float energy_pred_fac[4] = { 0.2, 0.3, 0.4, 0.5 };
+
+/** impulse response filter tables converted to float from Q15
+ * used for anti-sparseness processing */
+static const float ir_filter_str[64] = {
+ 6.159058e-01, 2.958069e-01, 9.979248e-02, -1.048889e-01,
+ 8.740234e-02, -1.599121e-01, 4.849243e-02, -4.141235e-02,
+ 1.831055e-02, 1.188049e-01, -4.568481e-02, -2.130127e-02,
+ 3.671265e-02, -1.601868e-01, 3.659058e-02, 1.639099e-01,
+ -4.541016e-02, -2.151489e-02, -8.810425e-02, 6.030273e-02,
+ 2.740479e-02, 2.200317e-02, -1.182861e-01, 1.289978e-01,
+ -1.560059e-01, 1.953125e-01, -3.149414e-02, -1.441956e-01,
+ 1.249084e-01, -1.328125e-01, 9.780884e-02, 6.500244e-02,
+ -6.091309e-02, -5.599976e-02, 8.081055e-02, -5.450439e-02,
+ -1.239014e-02, 1.748657e-02, 7.580566e-02, -1.101074e-01,
+ 9.579468e-02, -4.159546e-02, -7.830811e-02, 1.162109e-01,
+ -1.950073e-02, -6.259155e-02, -1.651001e-02, 7.250977e-02,
+ 1.199951e-01, -1.911011e-01, 4.370117e-02, -1.098938e-01,
+ 1.492004e-01, 1.129150e-02, 1.730347e-02, -3.549194e-02,
+ -8.709717e-02, 5.841064e-02, 1.190186e-03, -7.379150e-02,
+ 1.054077e-01, 9.078979e-02, -1.227112e-01, 1.047058e-01
+};
+
+static const float ir_filter_mid[64] = {
+ 7.354126e-01, 3.192139e-01, -1.606140e-01, -2.328491e-02,
+ 6.250000e-02, -2.828979e-02, 5.349731e-02, -1.014099e-01,
+ 6.750488e-02, 1.989746e-02, -6.549072e-02, 7.589722e-02,
+ -1.080017e-01, 1.253967e-01, -6.430054e-02, -1.141357e-02,
+ -1.910400e-02, 1.303101e-01, -1.673889e-01, 6.820679e-02,
+ 5.670166e-02, -8.450317e-02, 2.270508e-02, 3.479004e-02,
+ -2.328491e-02, -4.928589e-02, 1.239014e-01, -1.395874e-01,
+ 9.100342e-02, -3.549194e-02, 2.230835e-02, -3.350830e-02,
+ 2.450562e-02, 5.096436e-03, -2.178955e-02, 1.849365e-02,
+ -1.708984e-02, 1.950073e-02, 1.312256e-03, -5.389404e-02,
+ 9.851074e-02, -8.489990e-02, 2.029419e-02, 2.328491e-02,
+ 7.110596e-03, -6.109619e-02, 3.939819e-02, 5.709839e-02,
+ -1.058960e-01, 3.149414e-02, 8.270264e-02, -1.232910e-01,
+ 1.105957e-01, -1.286011e-01, 1.614990e-01, -1.303101e-01,
+ 4.769897e-02, 3.295898e-03, -1.770020e-02, 5.010986e-02,
+ -7.501221e-02, 2.920532e-02, 1.660156e-02, 7.751465e-02
+};
+
+static const float *ir_filters_lookup[2] = {
+ ir_filter_str, ir_filter_mid
+};
+
+/** High-pass filters coefficients for 31 Hz and 400 Hz cutoff */
+static const float hpf_zeros[2] = { -2.0, 1.0 };
+static const float hpf_31_poles[2] = { -1.978881836, 0.979125977 };
+static const float hpf_31_gain = 0.989501953;
+
+static const float hpf_400_poles[2] = { -1.787109375, 0.864257812 };
+static const float hpf_400_gain = 0.893554687;
+
+/** Interpolation coefficients for 5/4 signal upsampling
+ * Table from the reference source was reordered for efficiency */
+static const float upsample_fir[4][24] = {
+ { -6.103516e-05, 7.324219e-04, -2.014160e-03, 4.150391e-03,
+ -7.263184e-03, 1.165771e-02, -1.776123e-02, 2.624512e-02,
+ -3.869629e-02, 5.877686e-02, -9.863281e-02, 2.314453e-01,
+ 9.348755e-01, -1.523438e-01, 7.861328e-02, -4.937744e-02,
+ 3.308105e-02, -2.252197e-02, 1.507568e-02, -9.765625e-03,
+ 5.859375e-03, -3.173828e-03, 1.403809e-03, -3.662109e-04 },
+ { -2.441406e-04, 1.464844e-03, -3.784180e-03, 7.568359e-03,
+ -1.300049e-02, 2.062988e-02, -3.112793e-02, 4.589844e-02,
+ -6.781006e-02, 1.042480e-01, -1.815186e-01, 5.016479e-01,
+ 7.548828e-01, -2.094727e-01, 1.148071e-01, -7.348633e-02,
+ 4.956055e-02, -3.369141e-02, 2.246094e-02, -1.434326e-02,
+ 8.483887e-03, -4.455566e-03, 1.831055e-03, -4.272461e-04 },
+ { -4.272461e-04, 1.831055e-03, -4.455566e-03, 8.483887e-03,
+ -1.434326e-02, 2.246094e-02, -3.369141e-02, 4.956055e-02,
+ -7.348633e-02, 1.148071e-01, -2.094727e-01, 7.548828e-01,
+ 5.016479e-01, -1.815186e-01, 1.042480e-01, -6.781006e-02,
+ 4.589844e-02, -3.112793e-02, 2.062988e-02, -1.300049e-02,
+ 7.568359e-03, -3.784180e-03, 1.464844e-03, -2.441406e-04 },
+ { -3.662109e-04, 1.403809e-03, -3.173828e-03, 5.859375e-03,
+ -9.765625e-03, 1.507568e-02, -2.252197e-02, 3.308105e-02,
+ -4.937744e-02, 7.861328e-02, -1.523438e-01, 9.348755e-01,
+ 2.314453e-01, -9.863281e-02, 5.877686e-02, -3.869629e-02,
+ 2.624512e-02, -1.776123e-02, 1.165771e-02, -7.263184e-03,
+ 4.150391e-03, -2.014160e-03, 7.324219e-04, -6.103516e-05 }
+};
+
+/** High band quantized gains for 23k85 in Q14 */
+static const uint16_t qua_hb_gain[16] = {
+ 3624, 4673, 5597, 6479, 7425, 8378, 9324, 10264,
+ 11210, 12206, 13391, 14844, 16770, 19655, 24289, 32728
+};
+
+/** High-band post-processing FIR filters coefficients from Q15 */
+static const float bpf_6_7_coef[31] = { // band pass, 6kHz and 7kHz cutoffs
+ -2.441406e-04, 3.585815e-04, 2.441406e-04,
+ -2.059937e-04, -2.815248e-03, 8.560180e-03,
+ -1.084137e-02, 0.000000e+00, 2.897645e-02,
+ -6.774902e-02, 9.421540e-02, -8.380128e-02,
+ 2.706910e-02, 5.924987e-02, -1.373367e-01,
+ 1.687469e-01,
+ -1.373367e-01, 5.924987e-02, 2.706910e-02,
+ -8.380128e-02, 9.421540e-02, -6.774902e-02,
+ 2.897645e-02, 0.000000e+00, -1.084137e-02,
+ 8.560180e-03, -2.815248e-03, -2.059937e-04,
+ 2.441406e-04, 3.585815e-04, -2.441406e-04
+};
+
+static const float lpf_7_coef[31] = { // low pass, 7kHz cutoff
+ -6.408691e-04, 1.434326e-03, -2.716064e-03,
+ 4.455566e-03, -6.195068e-03, 6.988525e-03,
+ -5.401611e-03, 0.000000e+00, 1.022339e-02,
+ -2.560425e-02, 4.531860e-02, -6.747437e-02,
+ 8.944702e-02, -1.080933e-01, 1.206360e-01,
+ 8.753052e-01,
+ 1.206360e-01, -1.080933e-01, 8.944702e-02,
+ -6.747437e-02, 4.531860e-02, -2.560425e-02,
+ 1.022339e-02, 0.000000e+00, -5.401611e-03,
+ 6.988525e-03, -6.195068e-03, 4.455566e-03,
+ -2.716064e-03, 1.434326e-03, -6.408691e-04
+};
+
+/** Core frame sizes in each mode */
+static const uint16_t cf_sizes_wb[] = {
+ 132, 177, 253, 285, 317, 365, 397, 461, 477,
+ 40 /// SID/comfort noise frame
+};
+
+#endif
Added: trunk/libavcodec/amrwbdec.c
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ trunk/libavcodec/amrwbdec.c Sat Dec 18 04:03:18 2010 (r26051)
@@ -0,0 +1,1237 @@
+/*
+ * AMR wideband decoder
+ * Copyright (c) 2010 Marcelo Galvao Povoa
+ *
+ * This file is part of FFmpeg.
+ *
+ * FFmpeg is free software; you can redistribute it and/or
+ * modify it under the terms of the GNU Lesser General Public
+ * License as published by the Free Software Foundation; either
+ * version 2.1 of the License, or (at your option) any later version.
+ *
+ * FFmpeg is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A particular PURPOSE. See the GNU
+ * Lesser General Public License for more details.
+ *
+ * You should have received a copy of the GNU Lesser General Public
+ * License along with FFmpeg; if not, write to the Free Software
+ * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
+ */
+
+/**
+ * @file
+ * AMR wideband decoder
+ */
+
+#include "libavutil/lfg.h"
+
+#include "avcodec.h"
+#include "get_bits.h"
+#include "lsp.h"
+#include "celp_math.h"
+#include "celp_filters.h"
+#include "acelp_filters.h"
+#include "acelp_vectors.h"
+#include "acelp_pitch_delay.h"
+
+#define AMR_USE_16BIT_TABLES
+#include "amr.h"
+
+#include "amrwbdata.h"
+
+typedef struct {
+ AMRWBFrame frame; ///< AMRWB parameters decoded from bitstream
+ enum Mode fr_cur_mode; ///< mode index of current frame
+ uint8_t fr_quality; ///< frame quality index (FQI)
+ float isf_cur[LP_ORDER]; ///< working ISF vector from current frame
+ float isf_q_past[LP_ORDER]; ///< quantized ISF vector of the previous frame
+ float isf_past_final[LP_ORDER]; ///< final processed ISF vector of the previous frame
+ double isp[4][LP_ORDER]; ///< ISP vectors from current frame
+ double isp_sub4_past[LP_ORDER]; ///< ISP vector for the 4th subframe of the previous frame
+
+ float lp_coef[4][LP_ORDER]; ///< Linear Prediction Coefficients from ISP vector
+
+ uint8_t base_pitch_lag; ///< integer part of pitch lag for the next relative subframe
+ uint8_t pitch_lag_int; ///< integer part of pitch lag of the previous subframe
+
+ float excitation_buf[AMRWB_P_DELAY_MAX + LP_ORDER + 2 + AMRWB_SFR_SIZE]; ///< current excitation and all necessary excitation history
+ float *excitation; ///< points to current excitation in excitation_buf[]
+
+ float pitch_vector[AMRWB_SFR_SIZE]; ///< adaptive codebook (pitch) vector for current subframe
+ float fixed_vector[AMRWB_SFR_SIZE]; ///< algebraic codebook (fixed) vector for current subframe
+
+ float prediction_error[4]; ///< quantified prediction errors {20log10(^gamma_gc)} for previous four subframes
+ float pitch_gain[6]; ///< quantified pitch gains for the current and previous five subframes
+ float fixed_gain[2]; ///< quantified fixed gains for the current and previous subframes
+
+ float tilt_coef; ///< {beta_1} related to the voicing of the previous subframe
+
+ float prev_sparse_fixed_gain; ///< previous fixed gain; used by anti-sparseness to determine "onset"
+ uint8_t prev_ir_filter_nr; ///< previous impulse response filter "impNr": 0 - strong, 1 - medium, 2 - none
+ float prev_tr_gain; ///< previous initial gain used by noise enhancer for threshold
+
+ float samples_az[LP_ORDER + AMRWB_SFR_SIZE]; ///< low-band samples and memory from synthesis at 12.8kHz
+ float samples_up[UPS_MEM_SIZE + AMRWB_SFR_SIZE]; ///< low-band samples and memory processed for upsampling
+ float samples_hb[LP_ORDER_16k + AMRWB_SFR_SIZE_16k]; ///< high-band samples and memory from synthesis at 16kHz
+
+ float hpf_31_mem[2], hpf_400_mem[2]; ///< previous values in the high pass filters
+ float demph_mem[1]; ///< previous value in the de-emphasis filter
+ float bpf_6_7_mem[HB_FIR_SIZE]; ///< previous values in the high-band band pass filter
+ float lpf_7_mem[HB_FIR_SIZE]; ///< previous values in the high-band low pass filter
+
+ AVLFG prng; ///< random number generator for white noise excitation
+ uint8_t first_frame; ///< flag active during decoding of the first frame
+} AMRWBContext;
+
+static av_cold int amrwb_decode_init(AVCodecContext *avctx)
+{
+ AMRWBContext *ctx = avctx->priv_data;
+ int i;
+
+ avctx->sample_fmt = SAMPLE_FMT_FLT;
+
+ av_lfg_init(&ctx->prng, 1);
+
+ ctx->excitation = &ctx->excitation_buf[AMRWB_P_DELAY_MAX + LP_ORDER + 1];
+ ctx->first_frame = 1;
+
+ for (i = 0; i < LP_ORDER; i++)
+ ctx->isf_past_final[i] = isf_init[i] * (1.0f / (1 << 15));
+
+ for (i = 0; i < 4; i++)
+ ctx->prediction_error[i] = MIN_ENERGY;
+
+ return 0;
+}
+
+/**
+ * Decode the frame header in the "MIME/storage" format. This format
+ * is simpler and does not carry the auxiliary information of the frame
+ *
+ * @param[in] ctx The Context
+ * @param[in] buf Pointer to the input buffer
+ *
+ * @return The decoded header length in bytes
+ */
+static int decode_mime_header(AMRWBContext *ctx, const uint8_t *buf)
+{
+ GetBitContext gb;
+ init_get_bits(&gb, buf, 8);
+
+ /* Decode frame header (1st octet) */
+ skip_bits(&gb, 1); // padding bit
+ ctx->fr_cur_mode = get_bits(&gb, 4);
+ ctx->fr_quality = get_bits1(&gb);
+ skip_bits(&gb, 2); // padding bits
+
+ return 1;
+}
+
+/**
+ * Decodes quantized ISF vectors using 36-bit indexes (6K60 mode only)
+ *
+ * @param[in] ind Array of 5 indexes
+ * @param[out] isf_q Buffer for isf_q[LP_ORDER]
+ *
+ */
+static void decode_isf_indices_36b(uint16_t *ind, float *isf_q)
+{
+ int i;
+
+ for (i = 0; i < 9; i++)
+ isf_q[i] = dico1_isf[ind[0]][i] * (1.0f / (1 << 15));
+
+ for (i = 0; i < 7; i++)
+ isf_q[i + 9] = dico2_isf[ind[1]][i] * (1.0f / (1 << 15));
+
+ for (i = 0; i < 5; i++)
+ isf_q[i] += dico21_isf_36b[ind[2]][i] * (1.0f / (1 << 15));
+
+ for (i = 0; i < 4; i++)
+ isf_q[i + 5] += dico22_isf_36b[ind[3]][i] * (1.0f / (1 << 15));
+
+ for (i = 0; i < 7; i++)
+ isf_q[i + 9] += dico23_isf_36b[ind[4]][i] * (1.0f / (1 << 15));
+}
+
+/**
+ * Decodes quantized ISF vectors using 46-bit indexes (except 6K60 mode)
+ *
+ * @param[in] ind Array of 7 indexes
+ * @param[out] isf_q Buffer for isf_q[LP_ORDER]
+ *
+ */
+static void decode_isf_indices_46b(uint16_t *ind, float *isf_q)
+{
+ int i;
+
+ for (i = 0; i < 9; i++)
+ isf_q[i] = dico1_isf[ind[0]][i] * (1.0f / (1 << 15));
+
+ for (i = 0; i < 7; i++)
+ isf_q[i + 9] = dico2_isf[ind[1]][i] * (1.0f / (1 << 15));
+
+ for (i = 0; i < 3; i++)
+ isf_q[i] += dico21_isf[ind[2]][i] * (1.0f / (1 << 15));
+
+ for (i = 0; i < 3; i++)
+ isf_q[i + 3] += dico22_isf[ind[3]][i] * (1.0f / (1 << 15));
+
+ for (i = 0; i < 3; i++)
+ isf_q[i + 6] += dico23_isf[ind[4]][i] * (1.0f / (1 << 15));
+
+ for (i = 0; i < 3; i++)
+ isf_q[i + 9] += dico24_isf[ind[5]][i] * (1.0f / (1 << 15));
+
+ for (i = 0; i < 4; i++)
+ isf_q[i + 12] += dico25_isf[ind[6]][i] * (1.0f / (1 << 15));
+}
+
+/**
+ * Apply mean and past ISF values using the prediction factor
+ * Updates past ISF vector
+ *
+ * @param[in,out] isf_q Current quantized ISF
+ * @param[in,out] isf_past Past quantized ISF
+ *
+ */
+static void isf_add_mean_and_past(float *isf_q, float *isf_past)
+{
+ int i;
+ float tmp;
+
+ for (i = 0; i < LP_ORDER; i++) {
+ tmp = isf_q[i];
+ isf_q[i] += isf_mean[i] * (1.0f / (1 << 15));
+ isf_q[i] += PRED_FACTOR * isf_past[i];
+ isf_past[i] = tmp;
+ }
+}
+
+/**
+ * Interpolate the fourth ISP vector from current and past frames
+ * to obtain a ISP vector for each subframe
+ *
+ * @param[in,out] isp_q ISPs for each subframe
+ * @param[in] isp4_past Past ISP for subframe 4
+ */
+static void interpolate_isp(double isp_q[4][LP_ORDER], const double *isp4_past)
+{
+ int i, k;
+
+ for (k = 0; k < 3; k++) {
+ float c = isfp_inter[k];
+ for (i = 0; i < LP_ORDER; i++)
+ isp_q[k][i] = (1.0 - c) * isp4_past[i] + c * isp_q[3][i];
+ }
+}
+
+/**
+ * Decode an adaptive codebook index into pitch lag (except 6k60, 8k85 modes)
+ * Calculate integer lag and fractional lag always using 1/4 resolution
+ * In 1st and 3rd subframes the index is relative to last subframe integer lag
+ *
+ * @param[out] lag_int Decoded integer pitch lag
+ * @param[out] lag_frac Decoded fractional pitch lag
+ * @param[in] pitch_index Adaptive codebook pitch index
+ * @param[in,out] base_lag_int Base integer lag used in relative subframes
+ * @param[in] subframe Current subframe index (0 to 3)
+ */
+static void decode_pitch_lag_high(int *lag_int, int *lag_frac, int pitch_index,
+ uint8_t *base_lag_int, int subframe)
+{
+ if (subframe == 0 || subframe == 2) {
+ if (pitch_index < 376) {
+ *lag_int = (pitch_index + 137) >> 2;
+ *lag_frac = pitch_index - (*lag_int << 2) + 136;
+ } else if (pitch_index < 440) {
+ *lag_int = (pitch_index + 257 - 376) >> 1;
+ *lag_frac = (pitch_index - (*lag_int << 1) + 256 - 376) << 1;
+ /* the actual resolution is 1/2 but expressed as 1/4 */
+ } else {
+ *lag_int = pitch_index - 280;
+ *lag_frac = 0;
+ }
+ /* minimum lag for next subframe */
+ *base_lag_int = av_clip(*lag_int - 8 - (*lag_frac < 0),
+ AMRWB_P_DELAY_MIN, AMRWB_P_DELAY_MAX - 15);
+ // XXX: the spec states clearly that *base_lag_int should be
+ // the nearest integer to *lag_int (minus 8), but the ref code
+ // actually always uses its floor, I'm following the latter
+ } else {
+ *lag_int = (pitch_index + 1) >> 2;
+ *lag_frac = pitch_index - (*lag_int << 2);
+ *lag_int += *base_lag_int;
+ }
+}
+
+/**
+ * Decode a adaptive codebook index into pitch lag for 8k85 and 6k60 modes
+ * Description is analogous to decode_pitch_lag_high, but in 6k60 relative
+ * index is used for all subframes except the first
+ */
+static void decode_pitch_lag_low(int *lag_int, int *lag_frac, int pitch_index,
+ uint8_t *base_lag_int, int subframe, enum Mode mode)
+{
+ if (subframe == 0 || (subframe == 2 && mode != MODE_6k60)) {
+ if (pitch_index < 116) {
+ *lag_int = (pitch_index + 69) >> 1;
+ *lag_frac = (pitch_index - (*lag_int << 1) + 68) << 1;
+ } else {
+ *lag_int = pitch_index - 24;
+ *lag_frac = 0;
+ }
+ // XXX: same problem as before
+ *base_lag_int = av_clip(*lag_int - 8 - (*lag_frac < 0),
+ AMRWB_P_DELAY_MIN, AMRWB_P_DELAY_MAX - 15);
+ } else {
+ *lag_int = (pitch_index + 1) >> 1;
+ *lag_frac = (pitch_index - (*lag_int << 1)) << 1;
+ *lag_int += *base_lag_int;
+ }
+}
+
+/**
+ * Find the pitch vector by interpolating the past excitation at the
+ * pitch delay, which is obtained in this function
+ *
+ * @param[in,out] ctx The context
+ * @param[in] amr_subframe Current subframe data
+ * @param[in] subframe Current subframe index (0 to 3)
+ */
+static void decode_pitch_vector(AMRWBContext *ctx,
+ const AMRWBSubFrame *amr_subframe,
+ const int subframe)
+{
+ int pitch_lag_int, pitch_lag_frac;
+ int i;
+ float *exc = ctx->excitation;
+ enum Mode mode = ctx->fr_cur_mode;
+
+ if (mode <= MODE_8k85) {
+ decode_pitch_lag_low(&pitch_lag_int, &pitch_lag_frac, amr_subframe->adap,
+ &ctx->base_pitch_lag, subframe, mode);
+ } else
+ decode_pitch_lag_high(&pitch_lag_int, &pitch_lag_frac, amr_subframe->adap,
+ &ctx->base_pitch_lag, subframe);
+
+ ctx->pitch_lag_int = pitch_lag_int;
+ pitch_lag_int += pitch_lag_frac > 0;
+
+ /* Calculate the pitch vector by interpolating the past excitation at the
+ pitch lag using a hamming windowed sinc function */
+ ff_acelp_interpolatef(exc, exc + 1 - pitch_lag_int,
+ ac_inter, 4,
+ pitch_lag_frac + (pitch_lag_frac > 0 ? 0 : 4),
+ LP_ORDER, AMRWB_SFR_SIZE + 1);
+
+ /* Check which pitch signal path should be used
+ * 6k60 and 8k85 modes have the ltp flag set to 0 */
+ if (amr_subframe->ltp) {
+ memcpy(ctx->pitch_vector, exc, AMRWB_SFR_SIZE * sizeof(float));
+ } else {
+ for (i = 0; i < AMRWB_SFR_SIZE; i++)
+ ctx->pitch_vector[i] = 0.18 * exc[i - 1] + 0.64 * exc[i] +
+ 0.18 * exc[i + 1];
+ memcpy(exc, ctx->pitch_vector, AMRWB_SFR_SIZE * sizeof(float));
+ }
+}
+
+/** Get x bits in the index interval [lsb,lsb+len-1] inclusive */
+#define BIT_STR(x,lsb,len) (((x) >> (lsb)) & ((1 << (len)) - 1))
+
+/** Get the bit at specified position */
+#define BIT_POS(x, p) (((x) >> (p)) & 1)
+
+/**
+ * The next six functions decode_[i]p_track decode exactly i pulses
+ * positions and amplitudes (-1 or 1) in a subframe track using
+ * an encoded pulse indexing (TS 26.190 section 5.8.2)
+ *
+ * The results are given in out[], in which a negative number means
+ * amplitude -1 and vice versa (i.e., ampl(x) = x / abs(x) )
+ *
+ * @param[out] out Output buffer (writes i elements)
+ * @param[in] code Pulse index (no. of bits varies, see below)
+ * @param[in] m (log2) Number of potential positions
+ * @param[in] off Offset for decoded positions
+ */
+static inline void decode_1p_track(int *out, int code, int m, int off)
+{
+ int pos = BIT_STR(code, 0, m) + off; ///code: m+1 bits
+
+ out[0] = BIT_POS(code, m) ? -pos : pos;
+}
+
+static inline void decode_2p_track(int *out, int code, int m, int off) ///code: 2m+1 bits
+{
+ int pos0 = BIT_STR(code, m, m) + off;
+ int pos1 = BIT_STR(code, 0, m) + off;
+
+ out[0] = BIT_POS(code, 2*m) ? -pos0 : pos0;
+ out[1] = BIT_POS(code, 2*m) ? -pos1 : pos1;
+ out[1] = pos0 > pos1 ? -out[1] : out[1];
+}
+
+static void decode_3p_track(int *out, int code, int m, int off) ///code: 3m+1 bits
+{
+ int half_2p = BIT_POS(code, 2*m - 1) << (m - 1);
+
+ decode_2p_track(out, BIT_STR(code, 0, 2*m - 1),
+ m - 1, off + half_2p);
+ decode_1p_track(out + 2, BIT_STR(code, 2*m, m + 1), m, off);
+}
+
+static void decode_4p_track(int *out, int code, int m, int off) ///code: 4m bits
+{
+ int half_4p, subhalf_2p;
+ int b_offset = 1 << (m - 1);
+
+ switch (BIT_STR(code, 4*m - 2, 2)) { /* case ID (2 bits) */
+ case 0: /* 0 pulses in A, 4 pulses in B or vice versa */
+ half_4p = BIT_POS(code, 4*m - 3) << (m - 1); // which has 4 pulses
+ subhalf_2p = BIT_POS(code, 2*m - 3) << (m - 2);
+
+ decode_2p_track(out, BIT_STR(code, 0, 2*m - 3),
+ m - 2, off + half_4p + subhalf_2p);
+ decode_2p_track(out + 2, BIT_STR(code, 2*m - 2, 2*m - 1),
+ m - 1, off + half_4p);
+ break;
+ case 1: /* 1 pulse in A, 3 pulses in B */
+ decode_1p_track(out, BIT_STR(code, 3*m - 2, m),
+ m - 1, off);
+ decode_3p_track(out + 1, BIT_STR(code, 0, 3*m - 2),
+ m - 1, off + b_offset);
+ break;
+ case 2: /* 2 pulses in each half */
+ decode_2p_track(out, BIT_STR(code, 2*m - 1, 2*m - 1),
+ m - 1, off);
+ decode_2p_track(out + 2, BIT_STR(code, 0, 2*m - 1),
+ m - 1, off + b_offset);
+ break;
+ case 3: /* 3 pulses in A, 1 pulse in B */
+ decode_3p_track(out, BIT_STR(code, m, 3*m - 2),
+ m - 1, off);
+ decode_1p_track(out + 3, BIT_STR(code, 0, m),
+ m - 1, off + b_offset);
+ break;
+ }
+}
+
+static void decode_5p_track(int *out, int code, int m, int off) ///code: 5m bits
+{
+ int half_3p = BIT_POS(code, 5*m - 1) << (m - 1);
+
+ decode_3p_track(out, BIT_STR(code, 2*m + 1, 3*m - 2),
+ m - 1, off + half_3p);
+
+ decode_2p_track(out + 3, BIT_STR(code, 0, 2*m + 1), m, off);
+}
+
+static void decode_6p_track(int *out, int code, int m, int off) ///code: 6m-2 bits
+{
+ int b_offset = 1 << (m - 1);
+ /* which half has more pulses in cases 0 to 2 */
+ int half_more = BIT_POS(code, 6*m - 5) << (m - 1);
+ int half_other = b_offset - half_more;
+
+ switch (BIT_STR(code, 6*m - 4, 2)) { /* case ID (2 bits) */
+ case 0: /* 0 pulses in A, 6 pulses in B or vice versa */
+ decode_1p_track(out, BIT_STR(code, 0, m),
+ m - 1, off + half_more);
+ decode_5p_track(out + 1, BIT_STR(code, m, 5*m - 5),
+ m - 1, off + half_more);
+ break;
+ case 1: /* 1 pulse in A, 5 pulses in B or vice versa */
+ decode_1p_track(out, BIT_STR(code, 0, m),
+ m - 1, off + half_other);
+ decode_5p_track(out + 1, BIT_STR(code, m, 5*m - 5),
+ m - 1, off + half_more);
+ break;
+ case 2: /* 2 pulses in A, 4 pulses in B or vice versa */
+ decode_2p_track(out, BIT_STR(code, 0, 2*m - 1),
+ m - 1, off + half_other);
+ decode_4p_track(out + 2, BIT_STR(code, 2*m - 1, 4*m - 4),
+ m - 1, off + half_more);
+ break;
+ case 3: /* 3 pulses in A, 3 pulses in B */
+ decode_3p_track(out, BIT_STR(code, 3*m - 2, 3*m - 2),
+ m - 1, off);
+ decode_3p_track(out + 3, BIT_STR(code, 0, 3*m - 2),
+ m - 1, off + b_offset);
+ break;
+ }
+}
+
+/**
+ * Decode the algebraic codebook index to pulse positions and signs,
+ * then construct the algebraic codebook vector
+ *
+ * @param[out] fixed_vector Buffer for the fixed codebook excitation
+ * @param[in] pulse_hi MSBs part of the pulse index array (higher modes only)
+ * @param[in] pulse_lo LSBs part of the pulse index array
+ * @param[in] mode Mode of the current frame
+ */
+static void decode_fixed_vector(float *fixed_vector, const uint16_t *pulse_hi,
+ const uint16_t *pulse_lo, const enum Mode mode)
+{
+ /* sig_pos stores for each track the decoded pulse position indexes
+ * (1-based) multiplied by its corresponding amplitude (+1 or -1) */
+ int sig_pos[4][6];
+ int spacing = (mode == MODE_6k60) ? 2 : 4;
+ int i, j;
+
+ switch (mode) {
+ case MODE_6k60:
+ for (i = 0; i < 2; i++)
+ decode_1p_track(sig_pos[i], pulse_lo[i], 5, 1);
+ break;
+ case MODE_8k85:
+ for (i = 0; i < 4; i++)
+ decode_1p_track(sig_pos[i], pulse_lo[i], 4, 1);
+ break;
+ case MODE_12k65:
+ for (i = 0; i < 4; i++)
+ decode_2p_track(sig_pos[i], pulse_lo[i], 4, 1);
+ break;
+ case MODE_14k25:
+ for (i = 0; i < 2; i++)
+ decode_3p_track(sig_pos[i], pulse_lo[i], 4, 1);
+ for (i = 2; i < 4; i++)
+ decode_2p_track(sig_pos[i], pulse_lo[i], 4, 1);
+ break;
+ case MODE_15k85:
+ for (i = 0; i < 4; i++)
+ decode_3p_track(sig_pos[i], pulse_lo[i], 4, 1);
+ break;
+ case MODE_18k25:
+ for (i = 0; i < 4; i++)
+ decode_4p_track(sig_pos[i], (int) pulse_lo[i] +
+ ((int) pulse_hi[i] << 14), 4, 1);
+ break;
+ case MODE_19k85:
+ for (i = 0; i < 2; i++)
+ decode_5p_track(sig_pos[i], (int) pulse_lo[i] +
+ ((int) pulse_hi[i] << 10), 4, 1);
+ for (i = 2; i < 4; i++)
+ decode_4p_track(sig_pos[i], (int) pulse_lo[i] +
+ ((int) pulse_hi[i] << 14), 4, 1);
+ break;
+ case MODE_23k05:
+ case MODE_23k85:
+ for (i = 0; i < 4; i++)
+ decode_6p_track(sig_pos[i], (int) pulse_lo[i] +
+ ((int) pulse_hi[i] << 11), 4, 1);
+ break;
+ }
+
+ memset(fixed_vector, 0, sizeof(float) * AMRWB_SFR_SIZE);
+
+ for (i = 0; i < 4; i++)
+ for (j = 0; j < pulses_nb_per_mode_tr[mode][i]; j++) {
+ int pos = (FFABS(sig_pos[i][j]) - 1) * spacing + i;
+
+ fixed_vector[pos] += sig_pos[i][j] < 0 ? -1.0 : 1.0;
+ }
+}
+
+/**
+ * Decode pitch gain and fixed gain correction factor
+ *
+ * @param[in] vq_gain Vector-quantized index for gains
+ * @param[in] mode Mode of the current frame
+ * @param[out] fixed_gain_factor Decoded fixed gain correction factor
+ * @param[out] pitch_gain Decoded pitch gain
+ */
+static void decode_gains(const uint8_t vq_gain, const enum Mode mode,
+ float *fixed_gain_factor, float *pitch_gain)
+{
+ const int16_t *gains = (mode <= MODE_8k85 ? qua_gain_6b[vq_gain] :
+ qua_gain_7b[vq_gain]);
+
+ *pitch_gain = gains[0] * (1.0f / (1 << 14));
+ *fixed_gain_factor = gains[1] * (1.0f / (1 << 11));
+}
+
+/**
+ * Apply pitch sharpening filters to the fixed codebook vector
+ *
+ * @param[in] ctx The context
+ * @param[in,out] fixed_vector Fixed codebook excitation
+ */
+// XXX: Spec states this procedure should be applied when the pitch
+// lag is less than 64, but this checking seems absent in reference and AMR-NB
+static void pitch_sharpening(AMRWBContext *ctx, float *fixed_vector)
+{
+ int i;
+
+ /* Tilt part */
+ for (i = AMRWB_SFR_SIZE - 1; i != 0; i--)
+ fixed_vector[i] -= fixed_vector[i - 1] * ctx->tilt_coef;
+
+ /* Periodicity enhancement part */
+ for (i = ctx->pitch_lag_int; i < AMRWB_SFR_SIZE; i++)
+ fixed_vector[i] += fixed_vector[i - ctx->pitch_lag_int] * 0.85;
+}
+
+/**
+ * Calculate the voicing factor (-1.0 = unvoiced to 1.0 = voiced)
+ *
+ * @param[in] p_vector, f_vector Pitch and fixed excitation vectors
+ * @param[in] p_gain, f_gain Pitch and fixed gains
+ */
+// XXX: There is something wrong with the precision here! The magnitudes
+// of the energies are not correct. Please check the reference code carefully
+static float voice_factor(float *p_vector, float p_gain,
+ float *f_vector, float f_gain)
+{
+ double p_ener = (double) ff_dot_productf(p_vector, p_vector,
+ AMRWB_SFR_SIZE) * p_gain * p_gain;
+ double f_ener = (double) ff_dot_productf(f_vector, f_vector,
+ AMRWB_SFR_SIZE) * f_gain * f_gain;
+
+ return (p_ener - f_ener) / (p_ener + f_ener);
+}
+
+/**
+ * Reduce fixed vector sparseness by smoothing with one of three IR filters
+ * Also known as "adaptive phase dispersion"
+ *
+ * @param[in] ctx The context
+ * @param[in,out] fixed_vector Unfiltered fixed vector
+ * @param[out] buf Space for modified vector if necessary
+ *
+ * @return The potentially overwritten filtered fixed vector address
+ */
+static float *anti_sparseness(AMRWBContext *ctx,
+ float *fixed_vector, float *buf)
+{
+ int ir_filter_nr;
+
+ if (ctx->fr_cur_mode > MODE_8k85) // no filtering in higher modes
+ return fixed_vector;
+
+ if (ctx->pitch_gain[0] < 0.6) {
+ ir_filter_nr = 0; // strong filtering
+ } else if (ctx->pitch_gain[0] < 0.9) {
+ ir_filter_nr = 1; // medium filtering
+ } else
+ ir_filter_nr = 2; // no filtering
+
+ /* detect 'onset' */
+ if (ctx->fixed_gain[0] > 3.0 * ctx->fixed_gain[1]) {
+ if (ir_filter_nr < 2)
+ ir_filter_nr++;
+ } else {
+ int i, count = 0;
+
+ for (i = 0; i < 6; i++)
+ if (ctx->pitch_gain[i] < 0.6)
+ count++;
+
+ if (count > 2)
+ ir_filter_nr = 0;
+
+ if (ir_filter_nr > ctx->prev_ir_filter_nr + 1)
+ ir_filter_nr--;
+ }
+
+ /* update ir filter strength history */
+ ctx->prev_ir_filter_nr = ir_filter_nr;
+
+ ir_filter_nr += (ctx->fr_cur_mode == MODE_8k85);
+
+ if (ir_filter_nr < 2) {
+ int i;
+ const float *coef = ir_filters_lookup[ir_filter_nr];
+
+ /* Circular convolution code in the reference
+ * decoder was modified to avoid using one
+ * extra array. The filtered vector is given by:
+ *
+ * c2(n) = sum(i,0,len-1){ c(i) * coef( (n - i + len) % len ) }
+ */
+
+ memset(buf, 0, sizeof(float) * AMRWB_SFR_SIZE);
+ for (i = 0; i < AMRWB_SFR_SIZE; i++)
+ if (fixed_vector[i])
+ ff_celp_circ_addf(buf, buf, coef, i, fixed_vector[i],
+ AMRWB_SFR_SIZE);
+ fixed_vector = buf;
+ }
+
+ return fixed_vector;
+}
+
+/**
+ * Calculate a stability factor {teta} based on distance between
+ * current and past isf. A value of 1 shows maximum signal stability
+ */
+static float stability_factor(const float *isf, const float *isf_past)
+{
+ int i;
+ float acc = 0.0;
+
+ for (i = 0; i < LP_ORDER - 1; i++)
+ acc += (isf[i] - isf_past[i]) * (isf[i] - isf_past[i]);
+
+ // XXX: This part is not so clear from the reference code
+ // the result is more accurate changing the "/ 256" to "* 512"
+ return FFMAX(0.0, 1.25 - acc * 0.8 * 512);
+}
+
+/**
+ * Apply a non-linear fixed gain smoothing in order to reduce
+ * fluctuation in the energy of excitation
+ *
+ * @param[in] fixed_gain Unsmoothed fixed gain
+ * @param[in,out] prev_tr_gain Previous threshold gain (updated)
+ * @param[in] voice_fac Frame voicing factor
+ * @param[in] stab_fac Frame stability factor
+ *
+ * @return The smoothed gain
+ */
+static float noise_enhancer(float fixed_gain, float *prev_tr_gain,
+ float voice_fac, float stab_fac)
+{
+ float sm_fac = 0.5 * (1 - voice_fac) * stab_fac;
+ float g0;
+
+ // XXX: the following fixed-point constants used to in(de)crement
+ // gain by 1.5dB were taken from the reference code, maybe it could
+ // be simpler
+ if (fixed_gain < *prev_tr_gain) {
+ g0 = FFMIN(*prev_tr_gain, fixed_gain + fixed_gain *
+ (6226 * (1.0f / (1 << 15)))); // +1.5 dB
+ } else
+ g0 = FFMAX(*prev_tr_gain, fixed_gain *
+ (27536 * (1.0f / (1 << 15)))); // -1.5 dB
+
+ *prev_tr_gain = g0; // update next frame threshold
+
+ return sm_fac * g0 + (1 - sm_fac) * fixed_gain;
+}
+
+/**
+ * Filter the fixed_vector to emphasize the higher frequencies
+ *
+ * @param[in,out] fixed_vector Fixed codebook vector
+ * @param[in] voice_fac Frame voicing factor
+ */
+static void pitch_enhancer(float *fixed_vector, float voice_fac)
+{
+ int i;
+ float cpe = 0.125 * (1 + voice_fac);
+ float last = fixed_vector[0]; // holds c(i - 1)
+
+ fixed_vector[0] -= cpe * fixed_vector[1];
+
+ for (i = 1; i < AMRWB_SFR_SIZE - 1; i++) {
+ float cur = fixed_vector[i];
+
+ fixed_vector[i] -= cpe * (last + fixed_vector[i + 1]);
+ last = cur;
+ }
+
+ fixed_vector[AMRWB_SFR_SIZE - 1] -= cpe * last;
+}
+
+/**
+ * Conduct 16th order linear predictive coding synthesis from excitation
+ *
+ * @param[in] ctx Pointer to the AMRWBContext
+ * @param[in] lpc Pointer to the LPC coefficients
+ * @param[out] excitation Buffer for synthesis final excitation
+ * @param[in] fixed_gain Fixed codebook gain for synthesis
+ * @param[in] fixed_vector Algebraic codebook vector
+ * @param[in,out] samples Pointer to the output samples and memory
+ */
+static void synthesis(AMRWBContext *ctx, float *lpc, float *excitation,
+ float fixed_gain, const float *fixed_vector,
+ float *samples)
+{
+ ff_weighted_vector_sumf(excitation, ctx->pitch_vector, fixed_vector,
+ ctx->pitch_gain[0], fixed_gain, AMRWB_SFR_SIZE);
+
+ /* emphasize pitch vector contribution in low bitrate modes */
+ if (ctx->pitch_gain[0] > 0.5 && ctx->fr_cur_mode <= MODE_8k85) {
+ int i;
+ float energy = ff_dot_productf(excitation, excitation,
+ AMRWB_SFR_SIZE);
+
+ // XXX: Weird part in both ref code and spec. A unknown parameter
+ // {beta} seems to be identical to the current pitch gain
+ float pitch_factor = 0.25 * ctx->pitch_gain[0] * ctx->pitch_gain[0];
+
+ for (i = 0; i < AMRWB_SFR_SIZE; i++)
+ excitation[i] += pitch_factor * ctx->pitch_vector[i];
+
+ ff_scale_vector_to_given_sum_of_squares(excitation, excitation,
+ energy, AMRWB_SFR_SIZE);
+ }
+
+ ff_celp_lp_synthesis_filterf(samples, lpc, excitation,
+ AMRWB_SFR_SIZE, LP_ORDER);
+}
+
+/**
+ * Apply to synthesis a de-emphasis filter of the form:
+ * H(z) = 1 / (1 - m * z^-1)
+ *
+ * @param[out] out Output buffer
+ * @param[in] in Input samples array with in[-1]
+ * @param[in] m Filter coefficient
+ * @param[in,out] mem State from last filtering
+ */
+static void de_emphasis(float *out, float *in, float m, float mem[1])
+{
+ int i;
+
+ out[0] = in[0] + m * mem[0];
+
+ for (i = 1; i < AMRWB_SFR_SIZE; i++)
+ out[i] = in[i] + out[i - 1] * m;
+
+ mem[0] = out[AMRWB_SFR_SIZE - 1];
+}
+
+/**
+ * Upsample a signal by 5/4 ratio (from 12.8kHz to 16kHz) using
+ * a FIR interpolation filter. Uses past data from before *in address
+ *
+ * @param[out] out Buffer for interpolated signal
+ * @param[in] in Current signal data (length 0.8*o_size)
+ * @param[in] o_size Output signal length
+ */
+static void upsample_5_4(float *out, const float *in, int o_size)
+{
+ const float *in0 = in - UPS_FIR_SIZE + 1;
+ int i, j, k;
+ int int_part = 0, frac_part;
+
+ i = 0;
+ for (j = 0; j < o_size / 5; j++) {
+ out[i] = in[int_part];
+ frac_part = 4;
+ i++;
+
+ for (k = 1; k < 5; k++) {
+ out[i] = ff_dot_productf(in0 + int_part, upsample_fir[4 - frac_part],
+ UPS_MEM_SIZE);
+ int_part++;
+ frac_part--;
+ i++;
+ }
+ }
+}
+
+/**
+ * Calculate the high-band gain based on encoded index (23k85 mode) or
+ * on the low-band speech signal and the Voice Activity Detection flag
+ *
+ * @param[in] ctx The context
+ * @param[in] synth LB speech synthesis at 12.8k
+ * @param[in] hb_idx Gain index for mode 23k85 only
+ * @param[in] vad VAD flag for the frame
+ */
+static float find_hb_gain(AMRWBContext *ctx, const float *synth,
+ uint16_t hb_idx, uint8_t vad)
+{
+ int wsp = (vad > 0);
+ float tilt;
+
+ if (ctx->fr_cur_mode == MODE_23k85)
+ return qua_hb_gain[hb_idx] * (1.0f / (1 << 14));
+
+ tilt = ff_dot_productf(synth, synth + 1, AMRWB_SFR_SIZE - 1) /
+ ff_dot_productf(synth, synth, AMRWB_SFR_SIZE);
+
+ /* return gain bounded by [0.1, 1.0] */
+ return av_clipf((1.0 - FFMAX(0.0, tilt)) * (1.25 - 0.25 * wsp), 0.1, 1.0);
+}
+
+/**
+ * Generate the high-band excitation with the same energy from the lower
+ * one and scaled by the given gain
+ *
+ * @param[in] ctx The context
+ * @param[out] hb_exc Buffer for the excitation
+ * @param[in] synth_exc Low-band excitation used for synthesis
+ * @param[in] hb_gain Wanted excitation gain
+ */
+static void scaled_hb_excitation(AMRWBContext *ctx, float *hb_exc,
+ const float *synth_exc, float hb_gain)
+{
+ int i;
+ float energy = ff_dot_productf(synth_exc, synth_exc, AMRWB_SFR_SIZE);
+
+ /* Generate a white-noise excitation */
+ for (i = 0; i < AMRWB_SFR_SIZE_16k; i++)
+ hb_exc[i] = 32768.0 - (uint16_t) av_lfg_get(&ctx->prng);
+
+ ff_scale_vector_to_given_sum_of_squares(hb_exc, hb_exc,
+ energy * hb_gain * hb_gain,
+ AMRWB_SFR_SIZE_16k);
+}
+
+/**
+ * Calculate the auto-correlation for the ISF difference vector
+ */
+static float auto_correlation(float *diff_isf, float mean, int lag)
+{
+ int i;
+ float sum = 0.0;
+
+ for (i = 7; i < LP_ORDER - 2; i++) {
+ float prod = (diff_isf[i] - mean) * (diff_isf[i - lag] - mean);
+ sum += prod * prod;
+ }
+ return sum;
+}
+
+/**
+ * Extrapolate a ISF vector to the 16kHz range (20th order LP)
+ * used at mode 6k60 LP filter for the high frequency band
+ *
+ * @param[out] out Buffer for extrapolated isf
+ * @param[in] isf Input isf vector
+ */
+static void extrapolate_isf(float out[LP_ORDER_16k], float isf[LP_ORDER])
+{
+ float diff_isf[LP_ORDER - 2], diff_mean;
+ float *diff_hi = diff_isf - LP_ORDER + 1; // diff array for extrapolated indexes
+ float corr_lag[3];
+ float est, scale;
+ int i, i_max_corr;
+
+ memcpy(out, isf, (LP_ORDER - 1) * sizeof(float));
+ out[LP_ORDER_16k - 1] = isf[LP_ORDER - 1];
+
+ /* Calculate the difference vector */
+ for (i = 0; i < LP_ORDER - 2; i++)
+ diff_isf[i] = isf[i + 1] - isf[i];
+
+ diff_mean = 0.0;
+ for (i = 2; i < LP_ORDER - 2; i++)
+ diff_mean += diff_isf[i] * (1.0f / (LP_ORDER - 4));
+
+ /* Find which is the maximum autocorrelation */
+ i_max_corr = 0;
+ for (i = 0; i < 3; i++) {
+ corr_lag[i] = auto_correlation(diff_isf, diff_mean, i + 2);
+
+ if (corr_lag[i] > corr_lag[i_max_corr])
+ i_max_corr = i;
+ }
+ i_max_corr++;
+
+ for (i = LP_ORDER - 1; i < LP_ORDER_16k - 1; i++)
+ out[i] = isf[i - 1] + isf[i - 1 - i_max_corr]
+ - isf[i - 2 - i_max_corr];
+
+ /* Calculate an estimate for ISF(18) and scale ISF based on the error */
+ est = 7965 + (out[2] - out[3] - out[4]) / 6.0;
+ scale = 0.5 * (FFMIN(est, 7600) - out[LP_ORDER - 2]) /
+ (out[LP_ORDER_16k - 2] - out[LP_ORDER - 2]);
+
+ for (i = LP_ORDER - 1; i < LP_ORDER_16k - 1; i++)
+ diff_hi[i] = scale * (out[i] - out[i - 1]);
+
+ /* Stability insurance */
+ for (i = LP_ORDER; i < LP_ORDER_16k - 1; i++)
+ if (diff_hi[i] + diff_hi[i - 1] < 5.0) {
+ if (diff_hi[i] > diff_hi[i - 1]) {
+ diff_hi[i - 1] = 5.0 - diff_hi[i];
+ } else
+ diff_hi[i] = 5.0 - diff_hi[i - 1];
+ }
+
+ for (i = LP_ORDER - 1; i < LP_ORDER_16k - 1; i++)
+ out[i] = out[i - 1] + diff_hi[i] * (1.0f / (1 << 15));
+
+ /* Scale the ISF vector for 16000 Hz */
+ for (i = 0; i < LP_ORDER_16k - 1; i++)
+ out[i] *= 0.8;
+}
+
+/**
+ * Spectral expand the LP coefficients using the equation:
+ * y[i] = x[i] * (gamma ** i)
+ *
+ * @param[out] out Output buffer (may use input array)
+ * @param[in] lpc LP coefficients array
+ * @param[in] gamma Weighting factor
+ * @param[in] size LP array size
+ */
+static void lpc_weighting(float *out, const float *lpc, float gamma, int size)
+{
+ int i;
+ float fac = gamma;
+
+ for (i = 0; i < size; i++) {
+ out[i] = lpc[i] * fac;
+ fac *= gamma;
+ }
+}
+
+/**
+ * Conduct 20th order linear predictive coding synthesis for the high
+ * frequency band excitation at 16kHz
+ *
+ * @param[in] ctx The context
+ * @param[in] subframe Current subframe index (0 to 3)
+ * @param[in,out] samples Pointer to the output speech samples
+ * @param[in] exc Generated white-noise scaled excitation
+ * @param[in] isf Current frame isf vector
+ * @param[in] isf_past Past frame final isf vector
+ */
+static void hb_synthesis(AMRWBContext *ctx, int subframe, float *samples,
+ const float *exc, const float *isf, const float *isf_past)
+{
+ float hb_lpc[LP_ORDER_16k];
+ enum Mode mode = ctx->fr_cur_mode;
+
+ if (mode == MODE_6k60) {
+ float e_isf[LP_ORDER_16k]; // ISF vector for extrapolation
+ double e_isp[LP_ORDER_16k];
+
+ ff_weighted_vector_sumf(e_isf, isf_past, isf, isfp_inter[subframe],
+ 1.0 - isfp_inter[subframe], LP_ORDER);
+
+ extrapolate_isf(e_isf, e_isf);
+
+ e_isf[LP_ORDER_16k - 1] *= 2.0;
+ ff_acelp_lsf2lspd(e_isp, e_isf, LP_ORDER_16k);
+ ff_amrwb_lsp2lpc(e_isp, hb_lpc, LP_ORDER_16k);
+
+ lpc_weighting(hb_lpc, hb_lpc, 0.9, LP_ORDER_16k);
+ } else {
+ lpc_weighting(hb_lpc, ctx->lp_coef[subframe], 0.6, LP_ORDER);
+ }
+
+ ff_celp_lp_synthesis_filterf(samples, hb_lpc, exc, AMRWB_SFR_SIZE_16k,
+ (mode == MODE_6k60) ? LP_ORDER_16k : LP_ORDER);
+}
+
+/**
+ * Apply to high-band samples a 15th order filter
+ * The filter characteristic depends on the given coefficients
+ *
+ * @param[out] out Buffer for filtered output
+ * @param[in] fir_coef Filter coefficients
+ * @param[in,out] mem State from last filtering (updated)
+ * @param[in] in Input speech data (high-band)
+ *
+ * @remark It is safe to pass the same array in in and out parameters
+ */
+static void hb_fir_filter(float *out, const float fir_coef[HB_FIR_SIZE + 1],
+ float mem[HB_FIR_SIZE], const float *in)
+{
+ int i, j;
+ float data[AMRWB_SFR_SIZE_16k + HB_FIR_SIZE]; // past and current samples
+
+ memcpy(data, mem, HB_FIR_SIZE * sizeof(float));
+ memcpy(data + HB_FIR_SIZE, in, AMRWB_SFR_SIZE_16k * sizeof(float));
+
+ for (i = 0; i < AMRWB_SFR_SIZE_16k; i++) {
+ out[i] = 0.0;
+ for (j = 0; j <= HB_FIR_SIZE; j++)
+ out[i] += data[i + j] * fir_coef[j];
+ }
+
+ memcpy(mem, data + AMRWB_SFR_SIZE_16k, HB_FIR_SIZE * sizeof(float));
+}
+
+/**
+ * Update context state before the next subframe
+ */
+static void update_sub_state(AMRWBContext *ctx)
+{
+ memmove(&ctx->excitation_buf[0], &ctx->excitation_buf[AMRWB_SFR_SIZE],
+ (AMRWB_P_DELAY_MAX + LP_ORDER + 1) * sizeof(float));
+
+ memmove(&ctx->pitch_gain[1], &ctx->pitch_gain[0], 5 * sizeof(float));
+ memmove(&ctx->fixed_gain[1], &ctx->fixed_gain[0], 1 * sizeof(float));
+
+ memmove(&ctx->samples_az[0], &ctx->samples_az[AMRWB_SFR_SIZE],
+ LP_ORDER * sizeof(float));
+ memmove(&ctx->samples_up[0], &ctx->samples_up[AMRWB_SFR_SIZE],
+ UPS_MEM_SIZE * sizeof(float));
+ memmove(&ctx->samples_hb[0], &ctx->samples_hb[AMRWB_SFR_SIZE_16k],
+ LP_ORDER_16k * sizeof(float));
+}
+
+static int amrwb_decode_frame(AVCodecContext *avctx, void *data, int *data_size,
+ AVPacket *avpkt)
+{
+ AMRWBContext *ctx = avctx->priv_data;
+ AMRWBFrame *cf = &ctx->frame;
+ const uint8_t *buf = avpkt->data;
+ int buf_size = avpkt->size;
+ int expected_fr_size, header_size;
+ float *buf_out = data;
+ float spare_vector[AMRWB_SFR_SIZE]; // extra stack space to hold result from anti-sparseness processing
+ float fixed_gain_factor; // fixed gain correction factor (gamma)
+ float *synth_fixed_vector; // pointer to the fixed vector that synthesis should use
+ float synth_fixed_gain; // the fixed gain that synthesis should use
+ float voice_fac, stab_fac; // parameters used for gain smoothing
+ float synth_exc[AMRWB_SFR_SIZE]; // post-processed excitation for synthesis
+ float hb_exc[AMRWB_SFR_SIZE_16k]; // excitation for the high frequency band
+ float hb_samples[AMRWB_SFR_SIZE_16k]; // filtered high-band samples from synthesis
+ float hb_gain;
+ int sub, i;
+
+ header_size = decode_mime_header(ctx, buf);
+ expected_fr_size = ((cf_sizes_wb[ctx->fr_cur_mode] + 7) >> 3) + 1;
+
+ if (buf_size < expected_fr_size) {
+ av_log(avctx, AV_LOG_ERROR,
+ "Frame too small (%d bytes). Truncated file?\n", buf_size);
+ *data_size = 0;
+ return buf_size;
+ }
+
+ if (!ctx->fr_quality || ctx->fr_cur_mode > MODE_SID)
+ av_log(avctx, AV_LOG_ERROR, "Encountered a bad or corrupted frame\n");
+
+ if (ctx->fr_cur_mode == MODE_SID) /* Comfort noise frame */
+ av_log_missing_feature(avctx, "SID mode", 1);
+
+ if (ctx->fr_cur_mode >= MODE_SID)
+ return -1;
+
+ ff_amr_bit_reorder((uint16_t *) &ctx->frame, sizeof(AMRWBFrame),
+ buf + header_size, amr_bit_orderings_by_mode[ctx->fr_cur_mode]);
+
+ /* Decode the quantized ISF vector */
+ if (ctx->fr_cur_mode == MODE_6k60) {
+ decode_isf_indices_36b(cf->isp_id, ctx->isf_cur);
+ } else {
+ decode_isf_indices_46b(cf->isp_id, ctx->isf_cur);
+ }
+
+ isf_add_mean_and_past(ctx->isf_cur, ctx->isf_q_past);
+ ff_set_min_dist_lsf(ctx->isf_cur, MIN_ISF_SPACING, LP_ORDER - 1);
+
+ stab_fac = stability_factor(ctx->isf_cur, ctx->isf_past_final);
+
+ ctx->isf_cur[LP_ORDER - 1] *= 2.0;
+ ff_acelp_lsf2lspd(ctx->isp[3], ctx->isf_cur, LP_ORDER);
+
+ /* Generate a ISP vector for each subframe */
+ if (ctx->first_frame) {
+ ctx->first_frame = 0;
+ memcpy(ctx->isp_sub4_past, ctx->isp[3], LP_ORDER * sizeof(double));
+ }
+ interpolate_isp(ctx->isp, ctx->isp_sub4_past);
+
+ for (sub = 0; sub < 4; sub++)
+ ff_amrwb_lsp2lpc(ctx->isp[sub], ctx->lp_coef[sub], LP_ORDER);
+
+ for (sub = 0; sub < 4; sub++) {
+ const AMRWBSubFrame *cur_subframe = &cf->subframe[sub];
+ float *sub_buf = buf_out + sub * AMRWB_SFR_SIZE_16k;
+
+ /* Decode adaptive codebook (pitch vector) */
+ decode_pitch_vector(ctx, cur_subframe, sub);
+ /* Decode innovative codebook (fixed vector) */
+ decode_fixed_vector(ctx->fixed_vector, cur_subframe->pul_ih,
+ cur_subframe->pul_il, ctx->fr_cur_mode);
+
+ pitch_sharpening(ctx, ctx->fixed_vector);
+
+ decode_gains(cur_subframe->vq_gain, ctx->fr_cur_mode,
+ &fixed_gain_factor, &ctx->pitch_gain[0]);
+
+ ctx->fixed_gain[0] =
+ ff_amr_set_fixed_gain(fixed_gain_factor,
+ ff_dot_productf(ctx->fixed_vector, ctx->fixed_vector,
+ AMRWB_SFR_SIZE) / AMRWB_SFR_SIZE,
+ ctx->prediction_error,
+ ENERGY_MEAN, energy_pred_fac);
+
+ /* Calculate voice factor and store tilt for next subframe */
+ voice_fac = voice_factor(ctx->pitch_vector, ctx->pitch_gain[0],
+ ctx->fixed_vector, ctx->fixed_gain[0]);
+ ctx->tilt_coef = voice_fac * 0.25 + 0.25;
+
+ /* Construct current excitation */
+ for (i = 0; i < AMRWB_SFR_SIZE; i++) {
+ ctx->excitation[i] *= ctx->pitch_gain[0];
+ ctx->excitation[i] += ctx->fixed_gain[0] * ctx->fixed_vector[i];
+ ctx->excitation[i] = truncf(ctx->excitation[i]);
+ }
+
+ /* Post-processing of excitation elements */
+ synth_fixed_gain = noise_enhancer(ctx->fixed_gain[0], &ctx->prev_tr_gain,
+ voice_fac, stab_fac);
+
+ synth_fixed_vector = anti_sparseness(ctx, ctx->fixed_vector,
+ spare_vector);
+
+ pitch_enhancer(synth_fixed_vector, voice_fac);
+
+ synthesis(ctx, ctx->lp_coef[sub], synth_exc, synth_fixed_gain,
+ synth_fixed_vector, &ctx->samples_az[LP_ORDER]);
+
+ /* Synthesis speech post-processing */
+ de_emphasis(&ctx->samples_up[UPS_MEM_SIZE],
+ &ctx->samples_az[LP_ORDER], PREEMPH_FAC, ctx->demph_mem);
+
+ ff_acelp_apply_order_2_transfer_function(&ctx->samples_up[UPS_MEM_SIZE],
+ &ctx->samples_up[UPS_MEM_SIZE], hpf_zeros, hpf_31_poles,
+ hpf_31_gain, ctx->hpf_31_mem, AMRWB_SFR_SIZE);
+
+ upsample_5_4(sub_buf, &ctx->samples_up[UPS_FIR_SIZE],
+ AMRWB_SFR_SIZE_16k);
+
+ /* High frequency band (6.4 - 7.0 kHz) generation part */
+ ff_acelp_apply_order_2_transfer_function(hb_samples,
+ &ctx->samples_up[UPS_MEM_SIZE], hpf_zeros, hpf_400_poles,
+ hpf_400_gain, ctx->hpf_400_mem, AMRWB_SFR_SIZE);
+
+ hb_gain = find_hb_gain(ctx, hb_samples,
+ cur_subframe->hb_gain, cf->vad);
+
+ scaled_hb_excitation(ctx, hb_exc, synth_exc, hb_gain);
+
+ hb_synthesis(ctx, sub, &ctx->samples_hb[LP_ORDER_16k],
+ hb_exc, ctx->isf_cur, ctx->isf_past_final);
+
+ /* High-band post-processing filters */
+ hb_fir_filter(hb_samples, bpf_6_7_coef, ctx->bpf_6_7_mem,
+ &ctx->samples_hb[LP_ORDER_16k]);
+
+ if (ctx->fr_cur_mode == MODE_23k85)
+ hb_fir_filter(hb_samples, lpf_7_coef, ctx->lpf_7_mem,
+ hb_samples);
+
+ /* Add the low and high frequency bands */
+ for (i = 0; i < AMRWB_SFR_SIZE_16k; i++)
+ sub_buf[i] = (sub_buf[i] + hb_samples[i]) * (1.0f / (1 << 15));
+
+ /* Update buffers and history */
+ update_sub_state(ctx);
+ }
+
+ /* update state for next frame */
+ memcpy(ctx->isp_sub4_past, ctx->isp[3], LP_ORDER * sizeof(ctx->isp[3][0]));
+ memcpy(ctx->isf_past_final, ctx->isf_cur, LP_ORDER * sizeof(float));
+
+ /* report how many samples we got */
+ *data_size = 4 * AMRWB_SFR_SIZE_16k * sizeof(float);
+
+ return expected_fr_size;
+}
+
+AVCodec amrwb_decoder = {
+ .name = "amrwb",
+ .type = CODEC_TYPE_AUDIO,
+ .id = CODEC_ID_AMR_WB,
+ .priv_data_size = sizeof(AMRWBContext),
+ .init = amrwb_decode_init,
+ .decode = amrwb_decode_frame,
+ .long_name = NULL_IF_CONFIG_SMALL("Adaptive Multi-Rate WideBand"),
+ .sample_fmts = (enum SampleFormat[]){SAMPLE_FMT_FLT,SAMPLE_FMT_NONE},
+};
More information about the ffmpeg-cvslog
mailing list