[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