summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorGravatar Alexey Yakovenko <wakeroid@gmail.com>2010-06-02 20:16:51 +0200
committerGravatar Alexey Yakovenko <wakeroid@gmail.com>2010-06-02 20:17:25 +0200
commita1d0a56df66331a842c261ec71303343ce48f26c (patch)
treed135b919691dfbb2b1e230e95309394d18b13acc
parent15af1ef566b2aa2ddb820b34c5ea65e6ef3ee65b (diff)
ported musepack plugin to libmpcdev SV8 (r435)
Diffstat
-rw-r--r--configure.ac7
-rw-r--r--plugins/musepack/Makefile.am25
-rw-r--r--plugins/musepack/crc32.c56
-rw-r--r--plugins/musepack/decoder.h101
-rw-r--r--plugins/musepack/huffman.c358
-rw-r--r--plugins/musepack/huffman.h83
-rw-r--r--plugins/musepack/internal.h94
-rw-r--r--plugins/musepack/mpc_bits_reader.c179
-rw-r--r--plugins/musepack/mpc_bits_reader.h149
-rw-r--r--plugins/musepack/mpc_decoder.c681
-rw-r--r--plugins/musepack/mpc_demux.c661
-rw-r--r--plugins/musepack/mpc_reader.c144
-rw-r--r--plugins/musepack/mpcdec.h148
-rw-r--r--plugins/musepack/mpcdec_math.h135
-rw-r--r--plugins/musepack/musepack.c104
-rw-r--r--plugins/musepack/reader.h98
-rw-r--r--plugins/musepack/requant.c124
-rw-r--r--plugins/musepack/requant.h61
-rw-r--r--plugins/musepack/streaminfo.c244
-rw-r--r--plugins/musepack/streaminfo.h109
-rw-r--r--plugins/musepack/synth_filter.c430
21 files changed, 3941 insertions, 50 deletions
diff --git a/configure.ac b/configure.ac
index 3b676a15..3f8455c2 100644
--- a/configure.ac
+++ b/configure.ac
@@ -355,12 +355,7 @@ if test "x$enable_emidi" != "xno" ; then
fi
if test "x$enable_musepack" != "xno" ; then
- AC_CHECK_LIB([mpcdec], [main], [HAVE_MPCDEC=yes])
- if test "x$HAVE_MPCDEC" = "xyes" ; then
- HAVE_MUSEPACK=yes
- MUSEPACK_LIBS="-lmpcdec"
- AC_SUBST(MUSEPACK_LIBS)
- fi
+ HAVE_MUSEPACK=yes
fi
if test "x$enable_wildmidi" != "xno" ; then
diff --git a/plugins/musepack/Makefile.am b/plugins/musepack/Makefile.am
index d2202fcb..496f1468 100644
--- a/plugins/musepack/Makefile.am
+++ b/plugins/musepack/Makefile.am
@@ -2,9 +2,28 @@ if HAVE_MUSEPACK
pkglib_LTLIBRARIES = musepack.la
-musepack_la_SOURCES = musepack.c
-musepack_la_LDFLAGS = -module $(MUSEPACK_LIBS)
+musepack_la_SOURCES = musepack.c\
+huffman.c\
+mpc_bits_reader.c\
+mpc_decoder.c\
+mpc_demux.c\
+mpc_reader.c\
+requant.c\
+streaminfo.c\
+synth_filter.c\
+crc32.c\
+decoder.h\
+huffman.h\
+internal.h\
+mpc_bits_reader.h\
+mpcdec.h\
+mpcdec_math.h\
+reader.h\
+requant.h\
+streaminfo.h
-AM_CFLAGS = $(CFLAGS) -std=c99 -fPIC
+musepack_la_LDFLAGS = -module
+
+AM_CFLAGS = $(CFLAGS) -fPIC
endif
diff --git a/plugins/musepack/crc32.c b/plugins/musepack/crc32.c
new file mode 100644
index 00000000..3347d3cf
--- /dev/null
+++ b/plugins/musepack/crc32.c
@@ -0,0 +1,56 @@
+/*
+* C Implementation: crc32
+*
+* code from http://www.w3.org/TR/PNG/#D-CRCAppendix
+*
+*/
+
+/* Table of CRCs of all 8-bit messages. */
+static unsigned long crc_table[256];
+
+/* Flag: has the table been computed? Initially false. */
+static int crc_table_computed = 0;
+
+/* Make the table for a fast CRC. */
+static void make_crc_table(void)
+{
+ unsigned long c;
+ int n, k;
+
+ for (n = 0; n < 256; n++) {
+ c = (unsigned long) n;
+ for (k = 0; k < 8; k++) {
+ if (c & 1)
+ c = 0xedb88320L ^ (c >> 1);
+ else
+ c = c >> 1;
+ }
+ crc_table[n] = c;
+ }
+ crc_table_computed = 1;
+}
+
+
+/* Update a running CRC with the bytes buf[0..len-1]--the CRC
+ should be initialized to all 1's, and the transmitted value
+ is the 1's complement of the final running CRC (see the
+ crc() routine below). */
+
+static unsigned long update_crc(unsigned long crc, unsigned char *buf, int len)
+{
+ unsigned long c = crc;
+ int n;
+
+ if (!crc_table_computed)
+ make_crc_table();
+ for (n = 0; n < len; n++) {
+ c = crc_table[(c ^ buf[n]) & 0xff] ^ (c >> 8);
+ }
+ return c;
+}
+
+/* Return the CRC of the bytes buf[0..len-1]. */
+unsigned long crc32(unsigned char *buf, int len)
+{
+ return update_crc(0xffffffffL, buf, len) ^ 0xffffffffL;
+}
diff --git a/plugins/musepack/decoder.h b/plugins/musepack/decoder.h
new file mode 100644
index 00000000..3813fc02
--- /dev/null
+++ b/plugins/musepack/decoder.h
@@ -0,0 +1,101 @@
+/*
+ Copyright (c) 2005-2009, The Musepack Development Team
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are
+ met:
+
+ * Redistributions of source code must retain the above copyright
+ notice, this list of conditions and the following disclaimer.
+
+ * Redistributions in binary form must reproduce the above
+ copyright notice, this list of conditions and the following
+ disclaimer in the documentation and/or other materials provided
+ with the distribution.
+
+ * Neither the name of the The Musepack Development Team nor the
+ names of its contributors may be used to endorse or promote
+ products derived from this software without specific prior
+ written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+ "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+ LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+ A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
+ OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
+ SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
+ LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
+ DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
+ THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
+ (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+ OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+*/
+/// \file decoder.h
+#ifndef _MPCDEC_DECODER_H_
+#define _MPCDEC_DECODER_H_
+#ifdef WIN32
+#pragma once
+#endif
+
+#include <mpc/reader.h>
+
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#define SEEKING_TABLE_SIZE 256u
+// set it to SLOW_SEEKING_WINDOW to not use fast seeking
+#define FAST_SEEKING_WINDOW 32
+// set it to FAST_SEEKING_WINDOW to only use fast seeking
+#define SLOW_SEEKING_WINDOW 0x80000000
+
+enum {
+ MPC_V_MEM = 2304,
+ MPC_DECODER_MEMSIZE = 16384, // overall buffer size
+};
+
+struct mpc_decoder_t {
+ /// @name internal state variables
+ //@{
+ mpc_uint32_t stream_version; ///< Streamversion of stream
+ mpc_int32_t max_band; ///< Maximum band-index used in stream (0...31)
+ mpc_uint32_t ms; ///< Mid/side stereo (0: off, 1: on)
+ mpc_uint32_t channels; ///< Number of channels in stream
+
+ mpc_uint64_t samples; ///< Number of samples in stream
+
+ mpc_uint64_t decoded_samples; ///< Number of samples decoded from file begining
+ mpc_uint32_t samples_to_skip; ///< Number samples to skip (used for seeking)
+ mpc_int32_t last_max_band; ///< number of bands used in the last frame
+
+ // randomizer state variables
+ mpc_uint32_t __r1;
+ mpc_uint32_t __r2;
+
+ mpc_int32_t SCF_Index_L [32] [3];
+ mpc_int32_t SCF_Index_R [32] [3]; // holds scalefactor-indices
+ mpc_quantizer Q [32]; // holds quantized samples
+ mpc_int32_t Res_L [32];
+ mpc_int32_t Res_R [32]; // holds the chosen quantizer for each subband
+ mpc_bool_t DSCF_Flag_L [32];
+ mpc_bool_t DSCF_Flag_R [32]; // differential SCF used?
+ mpc_int32_t SCFI_L [32];
+ mpc_int32_t SCFI_R [32]; // describes order of transmitted SCF
+ mpc_bool_t MS_Flag[32]; // MS used?
+#ifdef MPC_FIXED_POINT
+ mpc_uint8_t SCF_shift[256];
+#endif
+
+ MPC_SAMPLE_FORMAT V_L[MPC_V_MEM + 960];
+ MPC_SAMPLE_FORMAT V_R[MPC_V_MEM + 960];
+ MPC_SAMPLE_FORMAT Y_L[36][32];
+ MPC_SAMPLE_FORMAT Y_R[36][32];
+ MPC_SAMPLE_FORMAT SCF[256]; ///< holds adapted scalefactors (for clipping prevention)
+ //@}
+};
+
+#ifdef __cplusplus
+}
+#endif
+#endif
diff --git a/plugins/musepack/huffman.c b/plugins/musepack/huffman.c
new file mode 100644
index 00000000..52cdfd92
--- /dev/null
+++ b/plugins/musepack/huffman.c
@@ -0,0 +1,358 @@
+/*
+ Copyright (c) 2005-2009, The Musepack Development Team
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are
+ met:
+
+ * Redistributions of source code must retain the above copyright
+ notice, this list of conditions and the following disclaimer.
+
+ * Redistributions in binary form must reproduce the above
+ copyright notice, this list of conditions and the following
+ disclaimer in the documentation and/or other materials provided
+ with the distribution.
+
+ * Neither the name of the The Musepack Development Team nor the
+ names of its contributors may be used to endorse or promote
+ products derived from this software without specific prior
+ written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+ "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+ LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+ A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
+ OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
+ SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
+ LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
+ DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
+ THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
+ (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+ OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+*/
+/// \file huffman.c
+/// Implementations of sv7/sv8 huffman decoding functions.
+#include "huffman.h"
+
+
+// sv7 huffman tables
+static const mpc_huffman mpc_table_HuffHdr [10] = {
+{0x8000, 1, 0}, {0x6000, 3, 1}, {0x5e00, 7, -4}, {0x5d80, 9, 3}, {0x5d00, 9, 4}, {0x5c00, 8, -5}, {0x5800, 6, 2}, {0x5000, 5, -3}, {0x4000, 4, -2}, {0x0, 2, -1}
+};
+mpc_lut_data mpc_HuffHdr = {mpc_table_HuffHdr};
+
+const mpc_huffman mpc_table_HuffSCFI [4] = {
+ {0x8000, 1, 1}, {0x6000, 3, 2}, {0x4000, 3, 0}, {0x0, 2, 3}
+};
+
+static const mpc_huffman mpc_table_HuffDSCF [16] = {
+ {0xf800, 5, 5}, {0xf000, 5, -4}, {0xe000, 4, 3}, {0xd000, 4, -3}, {0xc000, 4, 8}, {0xa000, 3, 1}, {0x9000, 4, 0}, {0x8800, 5, -5}, {0x8400, 6, 7}, {0x8000, 6, -7}, {0x6000, 3, -1}, {0x4000, 3, 2}, {0x3000, 4, 4}, {0x2800, 5, 6}, {0x2000, 5, -6}, {0x0, 3, -2}
+};
+mpc_lut_data mpc_HuffDSCF = {mpc_table_HuffDSCF};
+
+static const mpc_huffman mpc_table_HuffQ1 [2] [27] = {
+ {
+ {0xe000, 3, 13}, {0xdc00, 6, 26}, {0xd800, 6, 0}, {0xd400, 6, 20}, {0xd000, 6, 6}, {0xc000, 4, 14}, {0xb000, 4, 12}, {0xa000, 4, 4}, {0x9000, 4, 22}, {0x8c00, 6, 8}, {0x8800, 6, 18}, {0x8400, 6, 24}, {0x8000, 6, 2}, {0x7000, 4, 16}, {0x6000, 4, 10}, {0x5800, 5, 17}, {0x5000, 5, 9}, {0x4800, 5, 1}, {0x4000, 5, 25}, {0x3800, 5, 5}, {0x3000, 5, 21}, {0x2800, 5, 3}, {0x2000, 5, 11}, {0x1800, 5, 15}, {0x1000, 5, 23}, {0x800, 5, 19}, {0x0, 5, 7}
+ }, {
+ {0x8000, 1, 13}, {0x7e00, 7, 15}, {0x7c00, 7, 1}, {0x7a00, 7, 11}, {0x7800, 7, 7}, {0x7600, 7, 17}, {0x7400, 7, 25}, {0x7200, 7, 19}, {0x7180, 9, 8}, {0x7100, 9, 18}, {0x7080, 9, 2}, {0x7000, 9, 24}, {0x6e00, 7, 3}, {0x6c00, 7, 23}, {0x6a00, 7, 21}, {0x6800, 7, 5}, {0x6700, 8, 0}, {0x6600, 8, 26}, {0x6500, 8, 6}, {0x6400, 8, 20}, {0x6000, 6, 9}, {0x5000, 4, 14}, {0x4000, 4, 12}, {0x3000, 4, 4}, {0x2000, 4, 22}, {0x1000, 4, 16}, {0x0, 4, 10}
+ }
+};
+
+static const mpc_huffman mpc_table_HuffQ2 [2] [25] = {
+ {
+ {0xf000, 4, 13}, {0xe000, 4, 17}, {0xd000, 4, 7}, {0xc000, 4, 11}, {0xbc00, 6, 1}, {0xb800, 6, 23}, {0xb600, 7, 4}, {0xb400, 7, 20}, {0xb200, 7, 0}, {0xb000, 7, 24}, {0xa800, 5, 22}, {0xa000, 5, 10}, {0x8000, 3, 12}, {0x7800, 5, 2}, {0x7000, 5, 14}, {0x6000, 4, 6}, {0x5000, 4, 18}, {0x4000, 4, 8}, {0x3000, 4, 16}, {0x2800, 5, 9}, {0x2000, 5, 5}, {0x1800, 5, 15}, {0x1000, 5, 21}, {0x800, 5, 19}, {0x0, 5, 3}
+ }, {
+ {0xf800, 5, 18}, {0xf000, 5, 6}, {0xe800, 5, 8}, {0xe700, 8, 3}, {0xe6c0, 10, 24}, {0xe680, 10, 4}, {0xe640, 10, 0}, {0xe600, 10, 20}, {0xe400, 7, 23}, {0xe200, 7, 1}, {0xe000, 7, 19}, {0xd800, 5, 16}, {0xd600, 7, 15}, {0xd400, 7, 21}, {0xd200, 7, 9}, {0xd000, 7, 5}, {0xcc00, 6, 2}, {0xc800, 6, 10}, {0xc400, 6, 14}, {0xc000, 6, 22}, {0x8000, 2, 12}, {0x6000, 3, 13}, {0x4000, 3, 17}, {0x2000, 3, 11}, {0x0, 3, 7}
+ }
+};
+
+static const mpc_huffman mpc_table_HuffQ3 [2] [7] = {
+ {
+ {0xe000, 3, 1}, {0xd000, 4, 3}, {0xc000, 4, -3}, {0xa000, 3, 2}, {0x8000, 3, -2}, {0x4000, 2, 0}, {0x0, 2, -1}
+ }, {
+ {0xc000, 2, 0}, {0x8000, 2, -1}, {0x4000, 2, 1}, {0x3000, 4, -2}, {0x2800, 5, 3}, {0x2000, 5, -3}, {0x0, 3, 2}
+ }
+};
+
+static const mpc_huffman mpc_table_HuffQ4 [2] [9] = {
+ {
+ {0xe000, 3, 0}, {0xc000, 3, -1}, {0xa000, 3, 1}, {0x8000, 3, -2}, {0x6000, 3, 2}, {0x5000, 4, -4}, {0x4000, 4, 4}, {0x2000, 3, 3}, {0x0, 3, -3}
+ }, {
+ {0xe000, 3, 1}, {0xd000, 4, 2}, {0xc000, 4, -3}, {0x8000, 2, 0}, {0x6000, 3, -2}, {0x5000, 4, 3}, {0x4800, 5, -4}, {0x4000, 5, 4}, {0x0, 2, -1}
+ }
+};
+
+static const mpc_huffman mpc_table_HuffQ5 [2] [15] = {
+ {
+ {0xf000, 4, 2}, {0xe800, 5, 5}, {0xe400, 6, -7}, {0xe000, 6, 7}, {0xd000, 4, -3}, {0xc000, 4, 3}, {0xb800, 5, -6}, {0xb000, 5, 6}, {0xa000, 4, -4}, {0x9000, 4, 4}, {0x8000, 4, -5}, {0x6000, 3, 0}, {0x4000, 3, -1}, {0x2000, 3, 1}, {0x0, 3, -2}
+ }, {
+ {0xf000, 4, 3}, {0xe800, 5, 4}, {0xe600, 7, 6}, {0xe500, 8, -7}, {0xe400, 8, 7}, {0xe000, 6, -6}, {0xc000, 3, 0}, {0xa000, 3, -1}, {0x8000, 3, 1}, {0x6000, 3, -2}, {0x4000, 3, 2}, {0x3800, 5, -5}, {0x3000, 5, 5}, {0x2000, 4, -4}, {0x0, 3, -3}
+ }
+};
+
+static const mpc_huffman mpc_table_HuffQ6 [2] [31] = {
+ {
+ {0xf800, 5, 3}, {0xf000, 5, -4}, {0xec00, 6, -11}, {0xe800, 6, 12}, {0xe000, 5, 4}, {0xd800, 5, 6}, {0xd000, 5, -5}, {0xc800, 5, 5}, {0xc000, 5, 7}, {0xb800, 5, -7}, {0xb400, 6, -12}, {0xb000, 6, -13}, {0xa800, 5, -6}, {0xa000, 5, 8}, {0x9800, 5, -8}, {0x9000, 5, 9}, {0x8800, 5, -9}, {0x8400, 6, 13}, {0x8200, 7, -15}, {0x8000, 7, 15}, {0x7000, 4, 0}, {0x6800, 5, -10}, {0x6000, 5, 10}, {0x5000, 4, -1}, {0x4000, 4, 2}, {0x3000, 4, 1}, {0x2000, 4, -2}, {0x1c00, 6, 14}, {0x1800, 6, -14}, {0x1000, 5, 11}, {0x0, 4, -3}
+ }, {
+ {0xf800, 5, -6}, {0xf000, 5, 6}, {0xe000, 4, 1}, {0xd000, 4, -1}, {0xce00, 7, 10}, {0xcc00, 7, -10}, {0xcb00, 8, -11}, {0xca80, 9, -12}, {0xca60, 11, 13}, {0xca58, 13, 15}, {0xca50, 13, -14}, {0xca48, 13, 14}, {0xca40, 13, -15}, {0xca00, 10, -13}, {0xc900, 8, 11}, {0xc800, 8, 12}, {0xc400, 6, -9}, {0xc000, 6, 9}, {0xb000, 4, -2}, {0xa000, 4, 2}, {0x9000, 4, 3}, {0x8000, 4, -3}, {0x7800, 5, -7}, {0x7000, 5, 7}, {0x6000, 4, -4}, {0x5000, 4, 4}, {0x4800, 5, -8}, {0x4000, 5, 8}, {0x3000, 4, 5}, {0x2000, 4, -5}, {0x0, 3, 0}
+ }
+};
+
+static const mpc_huffman mpc_table_HuffQ7 [2] [63] = {
+ {
+ {0xfc00, 6, 7}, {0xf800, 6, 8}, {0xf400, 6, 9}, {0xf000, 6, -8}, {0xec00, 6, 11}, {0xea00, 7, 21}, {0xe900, 8, -28}, {0xe800, 8, 28}, {0xe400, 6, -9}, {0xe200, 7, -22}, {0xe000, 7, -21}, {0xdc00, 6, -10}, {0xd800, 6, -11}, {0xd400, 6, 10}, {0xd000, 6, 12}, {0xcc00, 6, -13}, {0xca00, 7, 22}, {0xc800, 7, 23}, {0xc400, 6, -12}, {0xc000, 6, 13}, {0xbc00, 6, 14}, {0xb800, 6, -14}, {0xb600, 7, -23}, {0xb500, 8, -29}, {0xb400, 8, 29}, {0xb000, 6, -15}, {0xac00, 6, 15}, {0xa800, 6, 16}, {0xa400, 6, -16}, {0xa200, 7, -24}, {0xa000, 7, 24}, {0x9c00, 6, 17}, {0x9a00, 7, -25}, {0x9900, 8, -30}, {0x9800, 8, 30}, {0x9400, 6, -17}, {0x9000, 6, 18}, {0x8c00, 6, -18}, {0x8a00, 7, 25}, {0x8800, 7, 26}, {0x8400, 6, 19}, {0x8200, 7, -26}, {0x8000, 7, -27}, {0x7800, 5, 2}, {0x7400, 6, -19}, {0x7000, 6, 20}, {0x6800, 5, -1}, {0x6700, 8, -31}, {0x6600, 8, 31}, {0x6400, 7, 27}, {0x6000, 6, -20}, {0x5800, 5, 1}, {0x5000, 5, -5}, {0x4800, 5, -3}, {0x4000, 5, 3}, {0x3800, 5, 0}, {0x3000, 5, -2}, {0x2800, 5, -4}, {0x2000, 5, 4}, {0x1800, 5, 5}, {0x1000, 5, -6}, {0x800, 5, 6}, {0x0, 5, -7}
+ }, {
+ {0xf800, 5, -1}, {0xf000, 5, 2}, {0xe800, 5, -2}, {0xe000, 5, 3}, {0xdf00, 8, -20}, {0xdec0, 10, 24}, {0xdebc, 14, 28}, {0xdeb8, 14, -28}, {0xdeb4, 14, -30}, {0xdeb0, 14, 30}, {0xdea0, 12, -27}, {0xde9c, 14, 29}, {0xde98, 14, -29}, {0xde94, 14, 31}, {0xde90, 14, -31}, {0xde80, 12, 27}, {0xde00, 9, -22}, {0xdc00, 7, -17}, {0xd800, 6, -11}, {0xd000, 5, -3}, {0xc800, 5, 4}, {0xc000, 5, -4}, {0xbe00, 7, 17}, {0xbd00, 8, 20}, {0xbc80, 9, 22}, {0xbc40, 10, -25}, {0xbc00, 10, -26}, {0xb800, 6, 12}, {0xb000, 5, 5}, {0xa800, 5, -5}, {0xa000, 5, 6}, {0x9800, 5, -6}, {0x9400, 6, -12}, {0x9200, 7, -18}, {0x9000, 7, 18}, {0x8c00, 6, 13}, {0x8800, 6, -13}, {0x8000, 5, -7}, {0x7c00, 6, 14}, {0x7b00, 8, 21}, {0x7a00, 8, -21}, {0x7800, 7, -19}, {0x7000, 5, 7}, {0x6800, 5, 8}, {0x6400, 6, -14}, {0x6000, 6, -15}, {0x5800, 5, -8}, {0x5400, 6, 15}, {0x5200, 7, 19}, {0x51c0, 10, 25}, {0x5180, 10, 26}, {0x5100, 9, -23}, {0x5080, 9, 23}, {0x5000, 9, -24}, {0x4800, 5, -9}, {0x4000, 5, 9}, {0x3c00, 6, 16}, {0x3800, 6, -16}, {0x3000, 5, 10}, {0x2000, 4, 0}, {0x1800, 5, -10}, {0x1000, 5, 11}, {0x0, 4, 1}
+ }
+};
+
+mpc_lut_data mpc_HuffQ [7] [2] = {
+ {{mpc_table_HuffQ1[0]}, {mpc_table_HuffQ1[1]}},
+ {{mpc_table_HuffQ2[0]}, {mpc_table_HuffQ2[1]}},
+ {{mpc_table_HuffQ3[0]}, {mpc_table_HuffQ3[1]}},
+ {{mpc_table_HuffQ4[0]}, {mpc_table_HuffQ4[1]}},
+ {{mpc_table_HuffQ5[0]}, {mpc_table_HuffQ5[1]}},
+ {{mpc_table_HuffQ6[0]}, {mpc_table_HuffQ6[1]}},
+ {{mpc_table_HuffQ7[0]}, {mpc_table_HuffQ7[1]}}
+};
+
+
+// sv8 huffman tables
+static const mpc_huffman mpc_huff_SCFI_1 [3] = {
+ {0x8000, 1, 1}, {0x4000, 2, 2}, {0x0, 3, 3}
+}; // 3
+static const mpc_int8_t mpc_sym_SCFI_1 [4] = {
+ 2, 3, 1, 0
+};
+static const mpc_huffman mpc_huff_SCFI_2 [5] = {
+ {0x8000, 2, 3}, {0x4000, 3, 5}, {0x1800, 5, 11}, {0x400, 6, 14}, {0x0, 7, 15}
+}; // 5
+static const mpc_int8_t mpc_sym_SCFI_2 [16] = {
+ 15, 10, 14, 11, 13, 9, 7, 6, 5, 12, 8, 3, 2, 0, 4, 1
+};
+mpc_can_data mpc_can_SCFI[2] = {{mpc_huff_SCFI_1, mpc_sym_SCFI_1}, {mpc_huff_SCFI_2, mpc_sym_SCFI_2}};
+
+static const mpc_huffman mpc_huff_DSCF_1 [12] = {
+ {0xa000, 3, 7}, {0x4000, 4, 12}, {0x2800, 5, 16}, {0x1800, 6, 21}, {0xe00, 7, 27}, {0x700, 8, 34}, {0x380, 9, 41}, {0x140, 10, 48}, {0x80, 11, 53}, {0x30, 12, 57}, {0x18, 13, 60}, {0x0, 14, 63}
+}; // 12
+static const mpc_int8_t mpc_sym_DSCF_1 [64] = {
+ 35, 34, 33, 36, 32, 30, 29, 27, 26, 37, 28, 25, 39, 38, 24, 23, 40, 22, 21, 20, 19, 43, 42, 41, 18, 17, 16, 15, 46, 45, 44, 14, 13, 12, 11, 49, 48, 47, 31, 10, 9, 8, 7, 6, 52, 51, 50, 5, 4, 3, 54, 53, 2, 1, 0, 57, 56, 55, 63, 62, 61, 60, 59, 58
+};
+static const mpc_huffman mpc_huff_DSCF_2 [13] = {
+ {0x6000, 3, 7}, {0x3000, 4, 10}, {0x1800, 5, 13}, {0x1000, 6, 16}, {0xa00, 7, 20}, {0x600, 8, 25}, {0x380, 9, 31}, {0x1c0, 10, 38}, {0xe0, 11, 45}, {0x50, 12, 52}, {0x20, 13, 57}, {0xc, 14, 61}, {0x0, 15, 64}
+}; // 13
+static const mpc_int8_t mpc_sym_DSCF_2 [65] = {
+ 33, 32, 31, 30, 29, 34, 28, 27, 36, 35, 26, 37, 25, 38, 24, 23, 40, 39, 22, 21, 42, 41, 20, 19, 18, 45, 44, 43, 17, 16, 15, 14, 48, 47, 46, 13, 12, 11, 10, 64, 52, 51, 50, 49, 9, 8, 7, 6, 55, 54, 53, 5, 4, 3, 58, 57, 56, 2, 1, 63, 62, 61, 60, 59, 0
+};
+mpc_can_data mpc_can_DSCF[2] = {{mpc_huff_DSCF_1, mpc_sym_DSCF_1}, {mpc_huff_DSCF_2, mpc_sym_DSCF_2}};
+
+static const mpc_huffman mpc_huff_Bands [12] = {
+ {0x8000, 1, 1}, {0x4000, 2, 2}, {0x2000, 3, 3}, {0x1000, 5, 6}, {0x800, 6, 8}, {0x600, 7, 10}, {0x300, 8, 13}, {0x200, 9, 16}, {0x140, 10, 20}, {0xc0, 11, 25}, {0x10, 12, 31}, {0x0, 13, 32}
+}; // 12
+static const mpc_int8_t mpc_sym_Bands [33] = {
+ 0, 32, 1, 31, 2, 30, 3, 4, 29, 6, 5, 28, 7, 27, 26, 8, 25, 24, 23, 9, 22, 21, 20, 18, 17, 16, 15, 14, 12, 11, 10, 19, 13
+};
+mpc_can_data mpc_can_Bands = {mpc_huff_Bands, mpc_sym_Bands};
+
+static const mpc_huffman mpc_huff_Res_1 [16] = {
+ {0x8000, 1, 1}, {0x4000, 2, 2}, {0x2000, 3, 3}, {0x1000, 4, 4}, {0x800, 5, 5}, {0x400, 6, 6}, {0x200, 7, 7}, {0x100, 8, 8}, {0x80, 9, 9}, {0x40, 10, 10}, {0x20, 11, 11}, {0x10, 12, 12}, {0x8, 13, 13}, {0x4, 14, 14}, {0x2, 15, 15}, {0x0, 16, 16}
+}; // 16
+static const mpc_int8_t mpc_sym_Res_1 [17] = {
+ 0, 1, 16, 2, 3, 4, 5, 15, 6, 7, 8, 9, 10, 11, 12, 14, 13
+};
+static const mpc_huffman mpc_huff_Res_2 [12] = {
+ {0x4000, 2, 3}, {0x2000, 3, 4}, {0x1000, 4, 5}, {0x800, 5, 6}, {0x400, 6, 7}, {0x200, 7, 8}, {0x100, 8, 9}, {0x80, 9, 10}, {0x40, 10, 11}, {0x20, 11, 12}, {0x10, 12, 13}, {0x0, 14, 16}
+}; // 12
+static const mpc_int8_t mpc_sym_Res_2 [17] = {
+ 16, 1, 0, 2, 15, 3, 14, 4, 5, 13, 6, 12, 7, 11, 10, 9, 8
+};
+mpc_can_data mpc_can_Res[2] = {{mpc_huff_Res_1, mpc_sym_Res_1}, {mpc_huff_Res_2, mpc_sym_Res_2}};
+
+static const mpc_huffman mpc_huff_Q1 [10] = {
+ {0x6000, 3, 7}, {0x1000, 4, 10}, {0x800, 5, 11}, {0x400, 6, 12}, {0x200, 7, 13}, {0x100, 8, 14}, {0x80, 9, 15}, {0x40, 10, 16}, {0x20, 11, 17}, {0x0, 12, 18}
+}; // 10
+static const mpc_int8_t mpc_sym_Q1 [19] = {
+ 7, 6, 5, 4, 3, 10, 9, 8, 2, 1, 11, 0, 12, 13, 14, 15, 16, 18, 17
+};
+mpc_can_data mpc_can_Q1 = {mpc_huff_Q1, mpc_sym_Q1};
+
+static const mpc_huffman mpc_huff_Q2_1 [10] = {
+ {0xe000, 3, 7}, {0x8000, 4, 14}, {0x3c00, 6, 38}, {0x2a00, 7, 53}, {0x1200, 8, 74}, {0x600, 9, 92}, {0x3c0, 10, 104}, {0x60, 11, 119}, {0x20, 12, 122}, {0x0, 13, 124}
+}; // 10
+static const mpc_int8_t mpc_sym_Q2_1 [125] = {
+ 62, 87, 67, 63, 61, 57, 37, 93, 92, 88, 86, 83, 82, 81, 68, 66, 58, 56, 42, 41, 38, 36, 32, 31, 112, 91, 72, 64, 60, 52, 43, 33, 12, 117, 113, 111, 107, 97, 89, 85, 77, 73, 71, 69, 65, 59, 55, 53, 51, 47, 39, 35, 27, 17, 13, 11, 7, 118, 116, 108, 106, 98, 96, 94, 90, 84, 80, 78, 76, 48, 46, 44, 40, 34, 30, 28, 26, 18, 16, 8, 6, 122, 110, 102, 74, 70, 54, 50, 22, 2, 123, 121, 119, 115, 114, 109, 105, 103, 101, 99, 95, 79, 75, 49, 45, 29, 25, 23, 21, 19, 15, 14, 10, 9, 5, 3, 1, 124, 104, 20, 0, 120, 100, 24, 4
+};
+static const mpc_huffman mpc_huff_Q2_2 [9] = {
+ {0xf000, 4, 15}, {0x7000, 5, 30}, {0x4800, 6, 44}, {0x3c00, 7, 62}, {0xc00, 8, 92}, {0x780, 9, 104}, {0xc0, 10, 119}, {0x40, 11, 122}, {0x0, 12, 124}
+}; // 9
+static const mpc_int8_t mpc_sym_Q2_2 [125] = {
+ 62, 92, 87, 86, 82, 68, 67, 66, 63, 61, 58, 57, 56, 42, 38, 37, 32, 93, 91, 88, 83, 81, 43, 41, 36, 33, 31, 112, 72, 64, 60, 52, 12, 118, 117, 116, 113, 111, 108, 107, 106, 98, 97, 96, 94, 90, 89, 85, 84, 80, 78, 77, 76, 73, 71, 69, 65, 59, 55, 53, 51, 48, 47, 46, 44, 40, 39, 35, 34, 30, 28, 27, 26, 18, 17, 16, 13, 11, 8, 7, 6, 122, 110, 74, 70, 54, 50, 22, 14, 2, 123, 121, 119, 115, 114, 109, 105, 103, 102, 101, 99, 95, 79, 75, 49, 45, 29, 25, 23, 21, 19, 15, 10, 9, 5, 3, 1, 124, 104, 20, 0, 120, 100, 24, 4
+};
+
+static const mpc_huffman mpc_huff_Q3 [7] = {
+ {0xe000, 3, 7}, {0x8000, 4, 14}, {0x5000, 5, 22}, {0x2400, 6, 32}, {0xa00, 7, 41}, {0x200, 8, 46}, {0x0, 9, 48}
+}; // 7
+static const mpc_int8_t mpc_sym_Q3 [49] = {
+ 0, 17, 16, 1, 15, -16, -1, 32, 31, 2, 14, -15, -32, 34, 33, 47, 46, 18, 30, -14, -2, -31, -17, -18, 49, 48, 63, 19, 29, 3, 13, -13, -3, -30, -47, -48, -33, 50, 62, 35, 45, -29, -19, -46, -34, 51, 61, -45, -35
+};
+
+static const mpc_huffman mpc_huff_Q4 [8] = {
+ {0xf000, 4, 15}, {0x9000, 5, 30}, {0x3400, 6, 48}, {0x1800, 7, 61}, {0x500, 8, 73}, {0x100, 9, 78}, {0x0, 10, 80}, {0x0, 0, 90}
+}; // 8
+static const mpc_int8_t mpc_sym_Q4 [91] = {
+ 0, 32, 17, 16, 31, 2, 1, 15, 14, -15, -16, -1, -32, 49, 48, 34, 33, 47, 46, 19, 18, 30, 29, 3, 13, -13, -14, -2, -3, -30, -31, -17, -18, -47, -48, -33, 64, 50, 63, 62, 35, 45, 4, 12, -29, -19, -46, -34, -64, -49, 66, 65, 79, 78, 51, 61, 36, 44, 20, 28, -12, -4, -28, -20, -45, -35, -62, -63, -50, 67, 77, 52, 60, -44, -36, -61, -51, 68, 76, -60, -52, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
+};
+
+static const mpc_huffman mpc_huff_Q5_1 [6] = {
+ {0xc000, 2, 3}, {0x4000, 3, 6}, {0x2000, 4, 8}, {0x1000, 5, 10}, {0x800, 6, 12}, {0x0, 7, 14}
+}; // 6
+static const mpc_int8_t mpc_sym_Q5_1 [15] = {
+ 0, 2, 1, -1, -2, 3, -3, 4, -4, 5, -5, 7, 6, -6, -7
+};
+static const mpc_huffman mpc_huff_Q5_2 [4] = {
+ {0x6000, 3, 7}, {0x2000, 4, 10}, {0x1000, 5, 12}, {0x0, 6, 14}
+}; // 4
+static const mpc_int8_t mpc_sym_Q5_2 [15] = {
+ 2, 1, 0, -1, -2, 4, 3, -3, -4, 5, -5, 7, 6, -6, -7
+};
+
+static const mpc_huffman mpc_huff_Q6_1 [8] = {
+ {0xc000, 2, 3}, {0x8000, 3, 6}, {0x4000, 4, 10}, {0x2800, 5, 14}, {0xc00, 6, 19}, {0x800, 7, 22}, {0x400, 8, 26}, {0x0, 9, 30}
+}; // 8
+static const mpc_int8_t mpc_sym_Q6_1 [31] = {
+ 0, 1, -1, 3, 2, -2, -3, 4, -4, -5, 8, 7, 6, 5, -6, -7, -8, 9, -9, 11, 10, -10, -11, 15, 14, 13, 12, -12, -13, -14, -15
+};
+static const mpc_huffman mpc_huff_Q6_2 [5] = {
+ {0x5000, 4, 15}, {0x2000, 5, 20}, {0x1000, 6, 24}, {0x400, 7, 28}, {0x0, 8, 30}
+}; // 5
+static const mpc_int8_t mpc_sym_Q6_2 [31] = {
+ 5, 4, 3, 2, 1, 0, -1, -2, -3, -4, -5, 8, 7, 6, -6, -7, -8, 10, 9, -9, -10, 13, 12, 11, -11, -12, -13, 15, 14, -14, -15
+};
+
+static const mpc_huffman mpc_huff_Q7_1 [9] = {
+ {0xc000, 2, 3}, {0x8000, 3, 6}, {0x6000, 4, 10}, {0x4000, 5, 16}, {0x2800, 6, 24}, {0x1400, 7, 34}, {0xa00, 8, 44}, {0x400, 9, 54}, {0x0, 10, 62}
+}; // 9
+static const mpc_int8_t mpc_sym_Q7_1 [63] = {
+ 0, 1, -1, 2, -2, 4, 3, -3, -4, 7, 6, 5, -5, -6, -7, 13, 11, 10, 9, 8, -8, -9, -10, -11, -12, 17, 16, 15, 14, 12, -13, -14, -15, -16, -17, 28, 27, 21, 20, 19, 18, -18, -19, -20, -21, -27, -28, 31, 30, 29, 26, 25, 24, 23, 22, -22, -23, -24, -25, -26, -29, -30, -31
+};
+static const mpc_huffman mpc_huff_Q7_2 [5] = {
+ {0x6000, 5, 31}, {0x2400, 6, 43}, {0x1000, 7, 52}, {0x200, 8, 60}, {0x0, 9, 62}
+}; // 5
+static const mpc_int8_t mpc_sym_Q7_2 [63] = {
+ 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, -1, -2, -3, -4, -5, -6, -7, -8, -9, 17, 16, 15, 14, 13, 12, 11, -10, -11, -12, -13, -14, -15, -16, -17, 22, 21, 20, 19, 18, -18, -19, -20, -21, -22, 29, 28, 27, 26, 25, 24, 23, -23, -24, -25, -26, -27, -28, -29, 31, 30, -30, -31
+};
+
+static const mpc_huffman mpc_huff_Q8_1 [11] = {
+ {0xc000, 2, 3}, {0x8000, 3, 6}, {0x7000, 4, 10}, {0x5800, 5, 17}, {0x3800, 6, 28}, {0x2800, 7, 42}, {0x1900, 8, 62}, {0xd00, 9, 87}, {0x280, 10, 113}, {0x60, 11, 123}, {0x0, 12, 126}
+}; // 11
+static const mpc_int8_t mpc_sym_Q8_1 [127] = {
+ 0, 1, -1, -2, 3, 2, -3, 7, 6, 5, 4, -4, -5, -6, -7, 11, 10, 9, 8, -8, -9, -10, -11, 19, 18, 17, 16, 15, 14, 13, 12, -12, -13, -14, -15, -16, -17, -19, 56, 55, 31, 28, 27, 26, 25, 24, 23, 22, 21, 20, -18, -20, -21, -22, -23, -24, -25, -26, -27, -33, -54, -56, 63, 62, 61, 60, 59, 58, 57, 54, 53, 43, 40, 39, 38, 37, 36, 35, 34, 33, 32, 30, 29, -28, -29, -30, -31, -32, -34, -35, -36, -37, -38, -39, -40, -41, -43, -53, -55, -57, -58, -59, -60, -61, 49, 47, 46, 45, 44, 42, 41, -42, -44, -45, -46, -47, -48, -49, -50, -62, -63, 52, 51, 50, 48, -51, -52
+};
+static const mpc_huffman mpc_huff_Q8_2 [4] = {
+ {0x9800, 6, 63}, {0x2a00, 7, 101}, {0x400, 8, 122}, {0x0, 9, 126}
+}; // 4
+static const mpc_int8_t mpc_sym_Q8_2 [127] = {
+ 13, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, -1, -2, -3, -4, -5, -6, -7, -8, -9, -10, -11, -12, -13, 39, 38, 37, 36, 35, 34, 33, 32, 31, 30, 29, 28, 27, 26, 25, 24, 23, 22, 21, 20, 19, 18, 17, 16, 15, 14, 12, -14, -15, -16, -17, -18, -19, -20, -21, -22, -23, -24, -25, -26, -27, -28, -29, -30, -31, -32, -33, -34, -35, -36, -37, -38, -39, -40, -41, 59, 58, 57, 56, 55, 54, 53, 52, 51, 50, 49, 48, 47, 46, 45, 44, 43, 42, 41, 40, -42, -43, -44, -45, -46, -47, -48, -49, -50, -51, -52, -53, -54, -55, -56, -57, -58, -59, 63, 62, 61, 60, -60, -61, -62, -63
+};
+
+static const mpc_huffman mpc_huff_Q9up [6] = {
+ {0xf800, 6, 63}, {0xac00, 7, 125}, {0x2600, 8, -45}, {0x280, 9, -7}, {0x40, 10, -2}, {0x0, 11, -1}
+}; // 6
+static const mpc_int8_t mpc_sym_Q9up [256] = {
+ -128, 127, -108, -110, -111, -112, -113, -114, -115, -116, -117, -118, -119, -120, -121, -122, -123, -124, -125, -126, -127, 126, 125, 124, 123, 122, 121, 120, 119, 118, 117, 116, 115, 114, 113, 112, 111, 110, 109, 108, -44, -45, -46, -47, -48, -49, -50, -51, -52, -53, -54, -55, -56, -57, -58, -59, -60, -61, -62, -63, -64, -65, -66, -67, -68, -69, -70, -71, -72, -73, -74, -75, -76, -77, -78, -79, -80, -81, -82, -83, -84, -85, -86, -87, -88, -89, -90, -91, -92, -93, -94, -95, -96, -97, -98, -99, -100, -101, -102, -103, -104, -105, -106, -107, -109, 107, 106, 105, 104, 103, 102, 101, 100, 99, 98, 97, 96, 95, 94, 93, 92, 91, 90, 89, 88, 87, 86, 85, 84, 83, 82, 81, 80, 79, 78, 77, 76, 75, 74, 73, 72, 71, 70, 69, 68, 67, 66, 65, 64, 63, 62, 61, 60, 59, 58, 57, 56, 55, 54, 53, 52, 51, 50, 49, 48, 47, 46, 45, 44, 43, 42, 40, 20, 19, -7, -8, -9, -10, -11, -12, -13, -14, -15, -16, -17, -18, -19, -20, -21, -22, -23, -24, -25, -26, -27, -28, -29, -30, -31, -32, -33, -34, -35, -36, -37, -38, -39, -40, -41, -42, -43, 41, 39, 38, 37, 36, 35, 34, 33, 32, 31, 30, 29, 28, 27, 26, 25, 24, 23, 22, 21, 18, 17, 16, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, -3, -4, -5, -6, 4, 3, 2, 1, 0, -1, -2
+};
+mpc_can_data mpc_can_Q9up = {mpc_huff_Q9up, mpc_sym_Q9up};
+
+
+mpc_can_data mpc_can_Q [6][2] = {
+ {{mpc_huff_Q2_1, mpc_sym_Q2_1}, {mpc_huff_Q2_2, mpc_sym_Q2_2}},
+ {{mpc_huff_Q3, mpc_sym_Q3}, {mpc_huff_Q4, mpc_sym_Q4}},
+ {{mpc_huff_Q5_1, mpc_sym_Q5_1}, {mpc_huff_Q5_2, mpc_sym_Q5_2}},
+ {{mpc_huff_Q6_1, mpc_sym_Q6_1}, {mpc_huff_Q6_2, mpc_sym_Q6_2}},
+ {{mpc_huff_Q7_1, mpc_sym_Q7_1}, {mpc_huff_Q7_2, mpc_sym_Q7_2}},
+ {{mpc_huff_Q8_1, mpc_sym_Q8_1}, {mpc_huff_Q8_2, mpc_sym_Q8_2}}
+};
+
+static void huff_fill_lut(const mpc_huffman * table, mpc_huff_lut * lut, const int bits)
+{
+ int i, idx = 0;
+ const int shift = 16 - bits;
+ for (i = (1 << bits) - 1; i >= 0 ; i--) {
+ if ((table[idx].Code >> shift) < i) {
+ lut[i].Length = table[idx].Length;
+ lut[i].Value = table[idx].Value;
+ } else {
+ if (table[idx].Length <= bits) {
+ lut[i].Length = table[idx].Length;
+ lut[i].Value = table[idx].Value;
+ } else {
+ lut[i].Length = 0;
+ lut[i].Value = idx;
+ }
+ if (i != 0)
+ do {
+ idx++;
+ } while ((table[idx].Code >> shift) == i);
+ }
+ }
+}
+
+static void can_fill_lut(mpc_can_data * data, const int bits)
+{
+ int i, idx = 0;
+ const int shift = 16 - bits;
+ const mpc_huffman * table = data->table;
+ const mpc_int8_t * sym = data->sym;
+ mpc_huff_lut * lut = data->lut;
+ for (i = (1 << bits) - 1; i >= 0 ; i--) {
+ if ((table[idx].Code >> shift) < i) {
+ if (table[idx].Length <= bits) {
+ lut[i].Length = table[idx].Length;
+ lut[i].Value = sym[(table[idx].Value - (i >> (bits - table[idx].Length))) & 0xFF];
+ } else {
+ lut[i].Length = 0;
+ lut[i].Value = idx;
+ }
+ } else {
+ if (table[idx].Length <= bits) {
+ lut[i].Length = table[idx].Length;
+ lut[i].Value = sym[(table[idx].Value - (i >> (bits - table[idx].Length))) & 0xFF];
+ } else {
+ lut[i].Length = 0;
+ lut[i].Value = idx;
+ }
+ if (i != 0)
+ do {
+ idx++;
+ } while ((table[idx].Code >> shift) == i);
+ }
+ }
+}
+
+void huff_init_lut(const int bits)
+{
+ int i, j;
+
+ huff_fill_lut(mpc_HuffDSCF.table, mpc_HuffDSCF.lut, bits);
+ huff_fill_lut(mpc_HuffHdr.table, mpc_HuffHdr.lut, bits);
+
+ can_fill_lut(&mpc_can_SCFI[0], bits);
+ can_fill_lut(&mpc_can_SCFI[1], bits);
+ can_fill_lut(&mpc_can_DSCF[0], bits);
+ can_fill_lut(&mpc_can_DSCF[1], bits);
+ can_fill_lut(&mpc_can_Res[0], bits);
+ can_fill_lut(&mpc_can_Res[1], bits);
+ can_fill_lut(&mpc_can_Q1, bits);
+ can_fill_lut(&mpc_can_Q9up, bits);
+
+
+ for( i = 0; i < 7; i++){
+ for( j = 0; j < 2; j++){
+ if (i != 6) can_fill_lut(&mpc_can_Q[i][j], bits);
+ huff_fill_lut(mpc_HuffQ[i][j].table, mpc_HuffQ[i][j].lut, bits);
+ }
+ }
+}
+
+
diff --git a/plugins/musepack/huffman.h b/plugins/musepack/huffman.h
new file mode 100644
index 00000000..5b12bff9
--- /dev/null
+++ b/plugins/musepack/huffman.h
@@ -0,0 +1,83 @@
+/*
+ Copyright (c) 2005-2009, The Musepack Development Team
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are
+ met:
+
+ * Redistributions of source code must retain the above copyright
+ notice, this list of conditions and the following disclaimer.
+
+ * Redistributions in binary form must reproduce the above
+ copyright notice, this list of conditions and the following
+ disclaimer in the documentation and/or other materials provided
+ with the distribution.
+
+ * Neither the name of the The Musepack Development Team nor the
+ names of its contributors may be used to endorse or promote
+ products derived from this software without specific prior
+ written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+ "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+ LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+ A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
+ OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
+ SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
+ LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
+ DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
+ THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
+ (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+ OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+*/
+/// \file huffman.h
+/// Data structures and functions for huffman coding.
+
+#ifndef _MPCDEC_HUFFMAN_H_
+#define _MPCDEC_HUFFMAN_H_
+#ifdef WIN32
+#pragma once
+#endif
+
+#include <mpc/mpc_types.h>
+
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+// LUT size parameter, LUT size is 1 << LUT_DEPTH
+#define LUT_DEPTH 6
+
+/// Huffman table entry.
+typedef struct mpc_huffman_t {
+ mpc_uint16_t Code;
+ mpc_uint8_t Length;
+ mpc_int8_t Value;
+} mpc_huffman;
+
+/// Huffman LUT entry.
+typedef struct mpc_huff_lut_t {
+ mpc_uint8_t Length;
+ mpc_int8_t Value;
+} mpc_huff_lut;
+
+/// Type used for huffman LUT decoding
+typedef struct mpc_lut_data_t {
+ mpc_huffman const * const table;
+ mpc_huff_lut lut[1 << LUT_DEPTH];
+} mpc_lut_data;
+
+/// Type used for canonical huffman decoding
+typedef struct mpc_can_data_t {
+ mpc_huffman const * const table;
+ mpc_int8_t const * const sym;
+ mpc_huff_lut lut[1 << LUT_DEPTH];
+} mpc_can_data;
+
+void huff_init_lut(const int bits);
+
+#ifdef __cplusplus
+}
+#endif
+#endif
diff --git a/plugins/musepack/internal.h b/plugins/musepack/internal.h
new file mode 100644
index 00000000..7810bcc2
--- /dev/null
+++ b/plugins/musepack/internal.h
@@ -0,0 +1,94 @@
+/*
+ Copyright (c) 2005-2009, The Musepack Development Team
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are
+ met:
+
+ * Redistributions of source code must retain the above copyright
+ notice, this list of conditions and the following disclaimer.
+
+ * Redistributions in binary form must reproduce the above
+ copyright notice, this list of conditions and the following
+ disclaimer in the documentation and/or other materials provided
+ with the distribution.
+
+ * Neither the name of the The Musepack Development Team nor the
+ names of its contributors may be used to endorse or promote
+ products derived from this software without specific prior
+ written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+ "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+ LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+ A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
+ OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
+ SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
+ LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
+ DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
+ THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
+ (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+ OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+*/
+/// \file internal.h
+/// Definitions and structures used only internally by the libmpcdec.
+#ifndef _MPCDEC_INTERNAL_H_
+#define _MPCDEC_INTERNAL_H_
+#ifdef WIN32
+#pragma once
+#endif
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include <mpc/mpcdec.h>
+
+/// Big/little endian 32 bit byte swapping routine.
+static mpc_inline
+mpc_uint32_t mpc_swap32(mpc_uint32_t val) {
+ return (((val & 0xFF000000) >> 24) | ((val & 0x00FF0000) >> 8)
+ | ((val & 0x0000FF00) << 8) | ((val & 0x000000FF) << 24));
+}
+
+typedef struct mpc_block_t {
+ char key[2]; // block key
+ mpc_uint64_t size; // block size minus the block header size
+} mpc_block;
+
+#define MAX_FRAME_SIZE 4352
+#define DEMUX_BUFFER_SIZE (65536 - MAX_FRAME_SIZE) // need some space as sand box
+
+struct mpc_demux_t {
+ mpc_reader * r;
+ mpc_decoder * d;
+ mpc_streaminfo si;
+
+ // buffer
+ mpc_uint8_t buffer[DEMUX_BUFFER_SIZE + MAX_FRAME_SIZE];
+ mpc_size_t bytes_total;
+ mpc_bits_reader bits_reader;
+ mpc_int32_t block_bits; /// bits remaining in current audio block
+ mpc_uint_t block_frames; /// frames remaining in current audio block
+
+ // seeking
+ mpc_seek_t * seek_table;
+ mpc_uint_t seek_pwr; /// distance between 2 frames in seek_table = 2^seek_pwr
+ mpc_uint32_t seek_table_size; /// used size in seek_table
+
+ // chapters
+ mpc_seek_t chap_pos; /// supposed position of the first chapter block
+ mpc_int_t chap_nb; /// number of chapters (-1 if unknown, 0 if no chapter)
+ mpc_chap_info * chap; /// chapters position and tag
+
+};
+
+/// helper functions used by multiple files
+mpc_uint32_t mpc_random_int(mpc_decoder *d); // in synth_filter.c
+void mpc_decoder_init_quant(mpc_decoder *d, double scale_factor);
+void mpc_decoder_synthese_filter_float(mpc_decoder *d, MPC_SAMPLE_FORMAT* OutData, mpc_int_t channels);
+
+#ifdef __cplusplus
+}
+#endif
+#endif
diff --git a/plugins/musepack/mpc_bits_reader.c b/plugins/musepack/mpc_bits_reader.c
new file mode 100644
index 00000000..5281288d
--- /dev/null
+++ b/plugins/musepack/mpc_bits_reader.c
@@ -0,0 +1,179 @@
+/*
+ Copyright (c) 2007-2009, The Musepack Development Team
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are
+ met:
+
+ * Redistributions of source code must retain the above copyright
+ notice, this list of conditions and the following disclaimer.
+
+ * Redistributions in binary form must reproduce the above
+ copyright notice, this list of conditions and the following
+ disclaimer in the documentation and/or other materials provided
+ with the distribution.
+
+ * Neither the name of the The Musepack Development Team nor the
+ names of its contributors may be used to endorse or promote
+ products derived from this software without specific prior
+ written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+ "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+ LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+ A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
+ OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
+ SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
+ LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
+ DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
+ THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
+ (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+ OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+*/
+
+#include <mpc/mpcdec.h>
+#include "internal.h"
+#include "huffman.h"
+#include "mpc_bits_reader.h"
+
+const mpc_uint32_t Cnk[MAX_ENUM / 2][MAX_ENUM] =
+{
+ {0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31},
+ {0, 0, 1, 3, 6, 10, 15, 21, 28, 36, 45, 55, 66, 78, 91, 105, 120, 136, 153, 171, 190, 210, 231, 253, 276, 300, 325, 351, 378, 406, 435, 465},
+ {0, 0, 0, 1, 4, 10, 20, 35, 56, 84, 120, 165, 220, 286, 364, 455, 560, 680, 816, 969, 1140, 1330, 1540, 1771, 2024, 2300, 2600, 2925, 3276, 3654, 4060, 4495},
+ {0, 0, 0, 0, 1, 5, 15, 35, 70, 126, 210, 330, 495, 715, 1001, 1365, 1820, 2380, 3060, 3876, 4845, 5985, 7315, 8855, 10626, 12650, 14950, 17550, 20475, 23751, 27405, 31465},
+ {0, 0, 0, 0, 0, 1, 6, 21, 56, 126, 252, 462, 792, 1287, 2002, 3003, 4368, 6188, 8568, 11628, 15504, 20349, 26334, 33649, 42504, 53130, 65780, 80730, 98280, 118755, 142506, 169911},
+ {0, 0, 0, 0, 0, 0, 1, 7, 28, 84, 210, 462, 924, 1716, 3003, 5005, 8008, 12376, 18564, 27132, 38760, 54264, 74613, 100947, 134596, 177100, 230230, 296010, 376740, 475020, 593775, 736281},
+ {0, 0, 0, 0, 0, 0, 0, 1, 8, 36, 120, 330, 792, 1716, 3432, 6435, 11440, 19448, 31824, 50388, 77520, 116280, 170544, 245157, 346104, 480700, 657800, 888030, 1184040, 1560780, 2035800, 2629575},
+ {0, 0, 0, 0, 0, 0, 0, 0, 1, 9, 45, 165, 495, 1287, 3003, 6435, 12870, 24310, 43758, 75582, 125970, 203490, 319770, 490314, 735471, 1081575, 1562275, 2220075, 3108105, 4292145, 5852925, 7888725},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 10, 55, 220, 715, 2002, 5005, 11440, 24310, 48620, 92378, 167960, 293930, 497420, 817190, 1307504, 2042975, 3124550, 4686825, 6906900, 10015005, 14307150, 20160075},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 11, 66, 286, 1001, 3003, 8008, 19448, 43758, 92378, 184756, 352716, 646646, 1144066, 1961256, 3268760, 5311735, 8436285, 13123110, 20030010, 30045015, 44352165},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 12, 78, 364, 1365, 4368, 12376, 31824, 75582, 167960, 352716, 705432, 1352078, 2496144, 4457400, 7726160, 13037895, 21474180, 34597290, 54627300, 84672315},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 13, 91, 455, 1820, 6188, 18564, 50388, 125970, 293930, 646646, 1352078, 2704156, 5200300, 9657700, 17383860, 30421755, 51895935, 86493225, 141120525},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 14, 105, 560, 2380, 8568, 27132, 77520, 203490, 497420, 1144066, 2496144, 5200300, 10400600, 20058300, 37442160, 67863915, 119759850, 206253075},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 15, 120, 680, 3060, 11628, 38760, 116280, 319770, 817190, 1961256, 4457400, 9657700, 20058300, 40116600, 77558760, 145422675, 265182525},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 16, 136, 816, 3876, 15504, 54264, 170544, 490314, 1307504, 3268760, 7726160, 17383860, 37442160, 77558760, 155117520, 300540195},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 17, 153, 969, 4845, 20349, 74613, 245157, 735471, 2042975, 5311735, 13037895, 30421755, 67863915, 145422675, 300540195}
+};
+
+const mpc_uint8_t Cnk_len[MAX_ENUM / 2][MAX_ENUM] =
+{
+ {0, 1, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5},
+ {0, 0, 2, 3, 4, 4, 5, 5, 6, 6, 6, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 9, 9, 9, 9},
+ {0, 0, 0, 2, 4, 5, 6, 6, 7, 7, 8, 8, 9, 9, 9, 10, 10, 10, 10, 11, 11, 11, 11, 11, 12, 12, 12, 12, 12, 12, 13, 13},
+ {0, 0, 0, 0, 3, 4, 6, 7, 7, 8, 9, 9, 10, 10, 11, 11, 12, 12, 12, 13, 13, 13, 14, 14, 14, 14, 15, 15, 15, 15, 15, 16},
+ {0, 0, 0, 0, 0, 3, 5, 6, 7, 8, 9, 10, 11, 11, 12, 13, 13, 14, 14, 14, 15, 15, 16, 16, 16, 17, 17, 17, 17, 18, 18, 18},
+ {0, 0, 0, 0, 0, 0, 3, 5, 7, 8, 9, 10, 11, 12, 13, 13, 14, 15, 15, 16, 16, 17, 17, 18, 18, 18, 19, 19, 19, 20, 20, 20},
+ {0, 0, 0, 0, 0, 0, 0, 3, 6, 7, 9, 10, 11, 12, 13, 14, 15, 15, 16, 17, 17, 18, 18, 19, 19, 20, 20, 21, 21, 21, 22, 22},
+ {0, 0, 0, 0, 0, 0, 0, 0, 4, 6, 8, 9, 11, 12, 13, 14, 15, 16, 17, 17, 18, 19, 19, 20, 21, 21, 22, 22, 23, 23, 23, 24},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 6, 8, 10, 11, 13, 14, 15, 16, 17, 18, 19, 19, 20, 21, 21, 22, 23, 23, 24, 24, 25, 25},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 7, 9, 10, 12, 13, 15, 16, 17, 18, 19, 20, 21, 21, 22, 23, 24, 24, 25, 25, 26, 26},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 7, 9, 11, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 23, 24, 25, 26, 26, 27, 27},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 7, 9, 11, 13, 15, 16, 17, 19, 20, 21, 22, 23, 24, 25, 25, 26, 27, 28, 28},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 7, 10, 12, 14, 15, 17, 18, 19, 21, 22, 23, 24, 25, 26, 27, 27, 28, 29},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 7, 10, 12, 14, 16, 17, 19, 20, 21, 23, 24, 25, 26, 27, 28, 28, 29},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 8, 10, 12, 14, 16, 18, 19, 21, 22, 23, 25, 26, 27, 28, 29, 30},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 5, 8, 10, 13, 15, 17, 18, 20, 21, 23, 24, 25, 27, 28, 29, 30}
+
+};
+
+const mpc_uint32_t Cnk_lost[MAX_ENUM / 2][MAX_ENUM] =
+{
+ {0, 0, 1, 0, 3, 2, 1, 0, 7, 6, 5, 4, 3, 2, 1, 0, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0},
+ {0, 0, 1, 2, 6, 1, 11, 4, 28, 19, 9, 62, 50, 37, 23, 8, 120, 103, 85, 66, 46, 25, 3, 236, 212, 187, 161, 134, 106, 77, 47, 16},
+ {0, 0, 0, 0, 6, 12, 29, 8, 44, 8, 91, 36, 226, 148, 57, 464, 344, 208, 55, 908, 718, 508, 277, 24, 1796, 1496, 1171, 820, 442, 36, 3697, 3232},
+ {0, 0, 0, 0, 3, 1, 29, 58, 2, 46, 182, 17, 309, 23, 683, 228, 1716, 1036, 220, 3347, 2207, 877, 7529, 5758, 3734, 1434, 15218, 12293, 9017, 5363, 1303, 29576},
+ {0, 0, 0, 0, 0, 2, 11, 8, 2, 4, 50, 232, 761, 46, 1093, 3824, 2004, 7816, 4756, 880, 12419, 6434, 31887, 23032, 12406, 65292, 50342, 32792, 12317, 119638, 92233, 60768},
+ {0, 0, 0, 0, 0, 0, 1, 4, 44, 46, 50, 100, 332, 1093, 3187, 184, 4008, 14204, 5636, 26776, 11272, 56459, 30125, 127548, 85044, 31914, 228278, 147548, 49268, 454801, 312295, 142384},
+ {0, 0, 0, 0, 0, 0, 0, 0, 28, 8, 182, 232, 332, 664, 1757, 4944, 13320, 944, 15148, 53552, 14792, 91600, 16987, 178184, 43588, 390776, 160546, 913112, 536372, 61352, 1564729, 828448},
+ {0, 0, 0, 0, 0, 0, 0, 0, 7, 19, 91, 17, 761, 1093, 1757, 3514, 8458, 21778, 55490, 5102, 58654, 204518, 33974, 313105, 1015577, 534877, 1974229, 1086199, 4096463, 2535683, 499883, 6258916},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 6, 9, 36, 309, 46, 3187, 4944, 8458, 16916, 38694, 94184, 230358, 26868, 231386, 789648, 54177, 1069754, 3701783, 1481708, 6762211, 2470066, 13394357, 5505632},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 5, 62, 226, 23, 1093, 184, 13320, 21778, 38694, 77388, 171572, 401930, 953086, 135896, 925544, 3076873, 8340931, 3654106, 13524422, 3509417, 22756699, 2596624},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 50, 148, 683, 3824, 4008, 944, 55490, 94184, 171572, 343144, 745074, 1698160, 3931208, 662448, 3739321, 12080252, 32511574, 12481564, 49545413, 5193248},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 37, 57, 228, 2004, 14204, 15148, 5102, 230358, 401930, 745074, 1490148, 3188308, 7119516, 16170572, 3132677, 15212929, 47724503, 127314931, 42642616},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 23, 464, 1716, 7816, 5636, 53552, 58654, 26868, 953086, 1698160, 3188308, 6376616, 13496132, 29666704, 66353813, 14457878, 62182381, 189497312},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 8, 344, 1036, 4756, 26776, 14792, 204518, 231386, 135896, 3931208, 7119516, 13496132, 26992264, 56658968, 123012781, 3252931, 65435312},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 120, 208, 220, 880, 11272, 91600, 33974, 789648, 925544, 662448, 16170572, 29666704, 56658968, 113317936, 236330717, 508019104},
+ {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 15, 103, 55, 3347, 12419, 56459, 16987, 313105, 54177, 3076873, 3739321, 3132677, 66353813, 123012781, 236330717}
+};
+
+static const mpc_uint8_t log2[32] =
+{ 1, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 6};
+
+static const mpc_uint8_t log2_lost[32] =
+{ 0, 1, 0, 3, 2, 1, 0, 7, 6, 5, 4, 3, 2, 1, 0, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, 31};
+
+mpc_int32_t mpc_bits_golomb_dec(mpc_bits_reader * r, const mpc_uint_t k)
+{
+ unsigned int l = 0;
+ unsigned int code = r->buff[0] & ((1 << r->count) - 1);
+
+ while( code == 0 ) {
+ l += r->count;
+ r->buff++;
+ code = r->buff[0];
+ r->count = 8;
+ }
+
+ while( ((1 << (r->count - 1)) & code) == 0 ) {
+ l++;
+ r->count--;
+ }
+ r->count--;
+
+ while( r->count < k ) {
+ r->buff++;
+ r->count += 8;
+ code = (code << 8) | r->buff[0];
+ }
+
+ r->count -= k;
+
+ return (l << k) | ((code >> r->count) & ((1 << k) - 1));
+}
+
+mpc_uint32_t mpc_bits_log_dec(mpc_bits_reader * r, mpc_uint_t max)
+{
+ mpc_uint32_t value = 0;
+ if (log2[max - 1] > 1)
+ value = mpc_bits_read(r, log2[max - 1] - 1);
+ if (value >= log2_lost[max - 1])
+ value = ((value << 1) | mpc_bits_read(r, 1)) - log2_lost[max - 1];
+ return value;
+}
+
+unsigned int mpc_bits_get_size(mpc_bits_reader * r, mpc_uint64_t * p_size)
+{
+ unsigned char tmp;
+ mpc_uint64_t size = 0;
+ unsigned int ret = 0;
+
+ do {
+ tmp = mpc_bits_read(r, 8);
+ size = (size << 7) | (tmp & 0x7F);
+ ret++;
+ } while((tmp & 0x80));
+
+ *p_size = size;
+ return ret;
+}
+
+int mpc_bits_get_block(mpc_bits_reader * r, mpc_block * p_block)
+{
+ int size = 2;
+
+ p_block->size = 0;
+ p_block->key[0] = mpc_bits_read(r, 8);
+ p_block->key[1] = mpc_bits_read(r, 8);
+
+ size += mpc_bits_get_size(r, &(p_block->size));
+
+ if (p_block->size >= size) // check if the block size doesn't conflict with the header size
+ p_block->size -= size;
+
+ return size;
+}
+
+
+
diff --git a/plugins/musepack/mpc_bits_reader.h b/plugins/musepack/mpc_bits_reader.h
new file mode 100644
index 00000000..a7d8d669
--- /dev/null
+++ b/plugins/musepack/mpc_bits_reader.h
@@ -0,0 +1,149 @@
+/*
+ Copyright (c) 2007-2009, The Musepack Development Team
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are
+ met:
+
+ * Redistributions of source code must retain the above copyright
+ notice, this list of conditions and the following disclaimer.
+
+ * Redistributions in binary form must reproduce the above
+ copyright notice, this list of conditions and the following
+ disclaimer in the documentation and/or other materials provided
+ with the distribution.
+
+ * Neither the name of the The Musepack Development Team nor the
+ names of its contributors may be used to endorse or promote
+ products derived from this software without specific prior
+ written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+ "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+ LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+ A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
+ OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
+ SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
+ LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
+ DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
+ THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
+ (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+ OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+*/
+
+#define MAX_ENUM 32
+
+MPC_API int mpc_bits_get_block(mpc_bits_reader * r, mpc_block * p_block);
+mpc_int32_t mpc_bits_golomb_dec(mpc_bits_reader * r, const mpc_uint_t k);
+MPC_API unsigned int mpc_bits_get_size(mpc_bits_reader * r, mpc_uint64_t * p_size);
+mpc_uint32_t mpc_bits_log_dec(mpc_bits_reader * r, mpc_uint_t max);
+
+extern const mpc_uint32_t Cnk[MAX_ENUM / 2][MAX_ENUM];
+extern const mpc_uint8_t Cnk_len[MAX_ENUM / 2][MAX_ENUM];
+extern const mpc_uint32_t Cnk_lost[MAX_ENUM / 2][MAX_ENUM];
+
+// can read up to 31 bits
+static mpc_inline mpc_uint32_t mpc_bits_read(mpc_bits_reader * r, const unsigned int nb_bits)
+{
+ mpc_uint32_t ret;
+
+ r->buff -= (int)(r->count - nb_bits) >> 3;
+ r->count = (r->count - nb_bits) & 0x07;
+
+ ret = (r->buff[0] | (r->buff[-1] << 8)) >> r->count;
+ if (nb_bits > (16 - r->count)) {
+ ret |= (mpc_uint32_t)((r->buff[-2] << 16) | (r->buff[-3] << 24)) >> r->count;
+ if (nb_bits > 24 && r->count != 0)
+ ret |= r->buff[-4] << (32 - r->count);
+ }
+
+ return ret & ((1 << nb_bits) - 1);
+}
+
+// basic huffman decoding routine
+// works with maximum lengths up to 16
+static mpc_inline mpc_int32_t mpc_bits_huff_dec(mpc_bits_reader * r, const mpc_huffman *Table)
+{
+ mpc_uint16_t code;
+ code = (mpc_uint16_t)((((r->buff[0] << 16) | (r->buff[1] << 8) | r->buff[2]) >> r->count) & 0xFFFF);
+
+ while (code < Table->Code) Table++;
+
+ r->buff -= (int)(r->count - Table->Length) >> 3;
+ r->count = (r->count - Table->Length) & 0x07;
+
+ return Table->Value;
+}
+
+static mpc_inline mpc_int32_t mpc_bits_can_dec(mpc_bits_reader * r, const mpc_can_data *can)
+{
+ mpc_uint16_t code;
+ mpc_huff_lut tmp;
+ const mpc_huffman * Table;
+ code = (mpc_uint16_t)((((r->buff[0] << 16) | (r->buff[1] << 8) | r->buff[2]) >> r->count) & 0xFFFF);
+
+ tmp = can->lut[code >> (16 - LUT_DEPTH)];
+ if (tmp.Length != 0) {
+ r->buff -= (int)(r->count - tmp.Length) >> 3;
+ r->count = (r->count - tmp.Length) & 0x07;
+ return tmp.Value;
+ }
+
+ Table = can->table + (unsigned char)tmp.Value;
+ while (code < Table->Code) Table++;
+
+ r->buff -= (int)(r->count - Table->Length) >> 3;
+ r->count = (r->count - Table->Length) & 0x07;
+
+ return can->sym[(Table->Value - (code >> (16 - Table->Length))) & 0xFF] ;
+}
+
+// LUT-based huffman decoding routine
+// works with maximum lengths up to 16
+static mpc_inline mpc_int32_t mpc_bits_huff_lut(mpc_bits_reader * r, const mpc_lut_data *lut)
+{
+ mpc_uint16_t code;
+ mpc_huff_lut tmp;
+ const mpc_huffman * Table;
+ code = (mpc_uint16_t)((((r->buff[0] << 16) | (r->buff[1] << 8) | r->buff[2]) >> r->count) & 0xFFFF);
+
+ tmp = lut->lut[code >> (16 - LUT_DEPTH)];
+ if (tmp.Length != 0) {
+ r->buff -= (int)(r->count - tmp.Length) >> 3;
+ r->count = (r->count - tmp.Length) & 0x07;
+ return tmp.Value;
+ }
+
+ Table = lut->table + (unsigned char)tmp.Value;
+ while (code < Table->Code) Table++;
+
+ r->buff -= (int)(r->count - Table->Length) >> 3;
+ r->count = (r->count - Table->Length) & 0x07;
+
+ return Table->Value;
+}
+
+static mpc_inline mpc_uint32_t mpc_bits_enum_dec(mpc_bits_reader * r, mpc_uint_t k, mpc_uint_t n)
+{
+ mpc_uint32_t bits = 0;
+ mpc_uint32_t code;
+ const mpc_uint32_t * C = Cnk[k-1];
+
+ code = mpc_bits_read(r, Cnk_len[k-1][n-1] - 1);
+
+ if (code >= Cnk_lost[k-1][n-1])
+ code = ((code << 1) | mpc_bits_read(r, 1)) - Cnk_lost[k-1][n-1];
+
+ do {
+ n--;
+ if (code >= C[n]) {
+ bits |= 1 << n;
+ code -= C[n];
+ C -= MAX_ENUM;
+ k--;
+ }
+ } while(k > 0);
+
+ return bits;
+}
diff --git a/plugins/musepack/mpc_decoder.c b/plugins/musepack/mpc_decoder.c
new file mode 100644
index 00000000..a7732bf2
--- /dev/null
+++ b/plugins/musepack/mpc_decoder.c
@@ -0,0 +1,681 @@
+/*
+ Copyright (c) 2005-2009, The Musepack Development Team
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are
+ met:
+
+ * Redistributions of source code must retain the above copyright
+ notice, this list of conditions and the following disclaimer.
+
+ * Redistributions in binary form must reproduce the above
+ copyright notice, this list of conditions and the following
+ disclaimer in the documentation and/or other materials provided
+ with the distribution.
+
+ * Neither the name of the The Musepack Development Team nor the
+ names of its contributors may be used to endorse or promote
+ products derived from this software without specific prior
+ written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+ "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+ LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+ A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
+ OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
+ SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
+ LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
+ DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
+ THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
+ (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+ OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+*/
+/// \file mpc_decoder.c
+/// Core decoding routines and logic.
+
+#include <string.h>
+#include <mpc/mpcdec.h>
+#include <mpc/minimax.h>
+#include "decoder.h"
+#include "huffman.h"
+#include "internal.h"
+#include "mpcdec_math.h"
+#include "requant.h"
+#include "mpc_bits_reader.h"
+
+//SV7 tables
+extern const mpc_lut_data mpc_HuffQ [7] [2];
+extern const mpc_lut_data mpc_HuffHdr;
+extern const mpc_huffman mpc_table_HuffSCFI [ 4];
+extern const mpc_lut_data mpc_HuffDSCF;
+
+//SV8 tables
+extern const mpc_can_data mpc_can_Bands;
+extern const mpc_can_data mpc_can_SCFI[2];
+extern const mpc_can_data mpc_can_DSCF[2];
+extern const mpc_can_data mpc_can_Res [2];
+extern const mpc_can_data mpc_can_Q [8][2];
+extern const mpc_can_data mpc_can_Q1;
+extern const mpc_can_data mpc_can_Q9up;
+
+//------------------------------------------------------------------------------
+// types
+//------------------------------------------------------------------------------
+enum
+{
+ MEMSIZE = MPC_DECODER_MEMSIZE, // overall buffer size
+ MEMSIZE2 = (MEMSIZE/2), // size of one buffer
+ MEMMASK = (MEMSIZE-1)
+};
+
+//------------------------------------------------------------------------------
+// forward declarations
+//------------------------------------------------------------------------------
+void mpc_decoder_read_bitstream_sv7(mpc_decoder * d, mpc_bits_reader * r);
+void mpc_decoder_read_bitstream_sv8(mpc_decoder * d, mpc_bits_reader * r,
+ mpc_bool_t is_key_frame);
+static void mpc_decoder_requantisierung(mpc_decoder *d);
+
+/**
+ * set the scf indexes for seeking use
+ * needed only for sv7 seeking
+ * @param d
+ */
+void mpc_decoder_reset_scf(mpc_decoder * d, int value)
+{
+ memset(d->SCF_Index_L, value, sizeof d->SCF_Index_L );
+ memset(d->SCF_Index_R, value, sizeof d->SCF_Index_R );
+}
+
+
+void mpc_decoder_setup(mpc_decoder *d)
+{
+ memset(d, 0, sizeof *d);
+
+ d->__r1 = 1;
+ d->__r2 = 1;
+
+ mpc_decoder_init_quant(d, 1.0f);
+}
+
+void mpc_decoder_set_streaminfo(mpc_decoder *d, mpc_streaminfo *si)
+{
+ d->stream_version = si->stream_version;
+ d->ms = si->ms;
+ d->max_band = si->max_band;
+ d->channels = si->channels;
+ d->samples_to_skip = MPC_DECODER_SYNTH_DELAY + si->beg_silence;
+
+ if (si->stream_version == 7 && si->is_true_gapless)
+ d->samples = ((si->samples + MPC_FRAME_LENGTH - 1) / MPC_FRAME_LENGTH) * MPC_FRAME_LENGTH;
+ else
+ d->samples = si->samples;
+}
+
+mpc_decoder * mpc_decoder_init(mpc_streaminfo *si)
+{
+ mpc_decoder* p_tmp = malloc(sizeof(mpc_decoder));
+
+ if (p_tmp != 0) {
+ mpc_decoder_setup(p_tmp);
+ mpc_decoder_set_streaminfo(p_tmp, si);
+ huff_init_lut(LUT_DEPTH); // FIXME : this needs to be called only once when the library is loaded
+ }
+
+ return p_tmp;
+}
+
+void mpc_decoder_exit(mpc_decoder *d)
+{
+ free(d);
+}
+
+void mpc_decoder_decode_frame(mpc_decoder * d,
+ mpc_bits_reader * r,
+ mpc_frame_info * i)
+{
+ mpc_bits_reader r_sav = *r;
+ mpc_int64_t samples_left;
+
+ samples_left = d->samples - d->decoded_samples + MPC_DECODER_SYNTH_DELAY;
+
+ if (samples_left <= 0 && d->samples != 0) {
+ i->samples = 0;
+ i->bits = -1;
+ return;
+ }
+
+ if (d->stream_version == 8)
+ mpc_decoder_read_bitstream_sv8(d, r, i->is_key_frame);
+ else
+ mpc_decoder_read_bitstream_sv7(d, r);
+
+ if (d->samples_to_skip < MPC_FRAME_LENGTH + MPC_DECODER_SYNTH_DELAY) {
+ mpc_decoder_requantisierung(d);
+ mpc_decoder_synthese_filter_float(d, i->buffer, d->channels);
+ }
+
+ d->decoded_samples += MPC_FRAME_LENGTH;
+
+ // reconstruct exact filelength
+ if (d->decoded_samples - d->samples < MPC_FRAME_LENGTH && d->stream_version == 7) {
+ int last_frame_samples = mpc_bits_read(r, 11);
+ if (d->decoded_samples == d->samples) {
+ if (last_frame_samples == 0) last_frame_samples = MPC_FRAME_LENGTH;
+ d->samples += last_frame_samples - MPC_FRAME_LENGTH;
+ samples_left += last_frame_samples - MPC_FRAME_LENGTH;
+ }
+ }
+
+ i->samples = samples_left > MPC_FRAME_LENGTH ? MPC_FRAME_LENGTH : samples_left < 0 ? 0 : (mpc_uint32_t) samples_left;
+ i->bits = (mpc_uint32_t) (((r->buff - r_sav.buff) << 3) + r_sav.count - r->count);
+
+ if (d->samples_to_skip) {
+ if (i->samples <= d->samples_to_skip) {
+ d->samples_to_skip -= i->samples;
+ i->samples = 0;
+ } else {
+ i->samples -= d->samples_to_skip;
+ memmove(i->buffer, i->buffer + d->samples_to_skip * d->channels,
+ i->samples * d->channels * sizeof (MPC_SAMPLE_FORMAT));
+ d->samples_to_skip = 0;
+ }
+ }
+}
+
+void
+mpc_decoder_requantisierung(mpc_decoder *d)
+{
+ mpc_int32_t Band;
+ mpc_int32_t n;
+ MPC_SAMPLE_FORMAT facL;
+ MPC_SAMPLE_FORMAT facR;
+ MPC_SAMPLE_FORMAT templ;
+ MPC_SAMPLE_FORMAT tempr;
+ MPC_SAMPLE_FORMAT* YL;
+ MPC_SAMPLE_FORMAT* YR;
+ mpc_int16_t* L;
+ mpc_int16_t* R;
+ const mpc_int32_t Last_Band = d->max_band;
+
+#ifdef MPC_FIXED_POINT
+#if MPC_FIXED_POINT_FRACTPART == 14
+#define MPC_MULTIPLY_SCF(CcVal, SCF_idx) \
+ MPC_MULTIPLY_EX(CcVal, d->SCF[SCF_idx], d->SCF_shift[SCF_idx])
+#else
+
+#error FIXME, Cc table is in 18.14 format
+
+#endif
+#else
+#define MPC_MULTIPLY_SCF(CcVal, SCF_idx) \
+ MPC_MULTIPLY(CcVal, d->SCF[SCF_idx])
+#endif
+ // requantization and scaling of subband-samples
+ for ( Band = 0; Band <= Last_Band; Band++ ) { // setting pointers
+ YL = d->Y_L[0] + Band;
+ YR = d->Y_R[0] + Band;
+ L = d->Q[Band].L;
+ R = d->Q[Band].R;
+ /************************** MS-coded **************************/
+ if ( d->MS_Flag [Band] ) {
+ if ( d->Res_L [Band] ) {
+ if ( d->Res_R [Band] ) { // M!=0, S!=0
+ facL = MPC_MULTIPLY_SCF( Cc[d->Res_L[Band]] , d->SCF_Index_L[Band][0] & 0xFF);
+ facR = MPC_MULTIPLY_SCF( Cc[d->Res_R[Band]] , d->SCF_Index_R[Band][0] & 0xFF);
+ for ( n = 0; n < 12; n++, YL += 32, YR += 32 ) {
+ *YL = (templ = MPC_MULTIPLY_FLOAT_INT(facL,*L++))+(tempr = MPC_MULTIPLY_FLOAT_INT(facR,*R++));
+ *YR = templ - tempr;
+ }
+ facL = MPC_MULTIPLY_SCF( Cc[d->Res_L[Band]] , d->SCF_Index_L[Band][1] & 0xFF);
+ facR = MPC_MULTIPLY_SCF( Cc[d->Res_R[Band]] , d->SCF_Index_R[Band][1] & 0xFF);
+ for ( ; n < 24; n++, YL += 32, YR += 32 ) {
+ *YL = (templ = MPC_MULTIPLY_FLOAT_INT(facL,*L++))+(tempr = MPC_MULTIPLY_FLOAT_INT(facR,*R++));
+ *YR = templ - tempr;
+ }
+ facL = MPC_MULTIPLY_SCF( Cc[d->Res_L[Band]] , d->SCF_Index_L[Band][2] & 0xFF);
+ facR = MPC_MULTIPLY_SCF( Cc[d->Res_R[Band]] , d->SCF_Index_R[Band][2] & 0xFF);
+ for ( ; n < 36; n++, YL += 32, YR += 32 ) {
+ *YL = (templ = MPC_MULTIPLY_FLOAT_INT(facL,*L++))+(tempr = MPC_MULTIPLY_FLOAT_INT(facR,*R++));
+ *YR = templ - tempr;
+ }
+ } else { // M!=0, S==0
+ facL = MPC_MULTIPLY_SCF( Cc[d->Res_L[Band]] , d->SCF_Index_L[Band][0] & 0xFF);
+ for ( n = 0; n < 12; n++, YL += 32, YR += 32 ) {
+ *YR = *YL = MPC_MULTIPLY_FLOAT_INT(facL,*L++);
+ }
+ facL = MPC_MULTIPLY_SCF( Cc[d->Res_L[Band]] , d->SCF_Index_L[Band][1] & 0xFF);
+ for ( ; n < 24; n++, YL += 32, YR += 32 ) {
+ *YR = *YL = MPC_MULTIPLY_FLOAT_INT(facL,*L++);
+ }
+ facL = MPC_MULTIPLY_SCF( Cc[d->Res_L[Band]] , d->SCF_Index_L[Band][2] & 0xFF);
+ for ( ; n < 36; n++, YL += 32, YR += 32 ) {
+ *YR = *YL = MPC_MULTIPLY_FLOAT_INT(facL,*L++);
+ }
+ }
+ } else {
+ if (d->Res_R[Band]) // M==0, S!=0
+ {
+ facR = MPC_MULTIPLY_SCF( Cc[d->Res_R[Band]] , d->SCF_Index_R[Band][0] & 0xFF);
+ for ( n = 0; n < 12; n++, YL += 32, YR += 32 ) {
+ *YR = - (*YL = MPC_MULTIPLY_FLOAT_INT(facR,*(R++)));
+ }
+ facR = MPC_MULTIPLY_SCF( Cc[d->Res_R[Band]] , d->SCF_Index_R[Band][1] & 0xFF);
+ for ( ; n < 24; n++, YL += 32, YR += 32 ) {
+ *YR = - (*YL = MPC_MULTIPLY_FLOAT_INT(facR,*(R++)));
+ }
+ facR = MPC_MULTIPLY_SCF( Cc[d->Res_R[Band]] , d->SCF_Index_R[Band][2] & 0xFF);
+ for ( ; n < 36; n++, YL += 32, YR += 32 ) {
+ *YR = - (*YL = MPC_MULTIPLY_FLOAT_INT(facR,*(R++)));
+ }
+ } else { // M==0, S==0
+ for ( n = 0; n < 36; n++, YL += 32, YR += 32 ) {
+ *YR = *YL = 0;
+ }
+ }
+ }
+ }
+ /************************** LR-coded **************************/
+ else {
+ if ( d->Res_L [Band] ) {
+ if ( d->Res_R [Band] ) { // L!=0, R!=0
+ facL = MPC_MULTIPLY_SCF( Cc[d->Res_L[Band]] , d->SCF_Index_L[Band][0] & 0xFF);
+ facR = MPC_MULTIPLY_SCF( Cc[d->Res_R[Band]] , d->SCF_Index_R[Band][0] & 0xFF);
+ for (n = 0; n < 12; n++, YL += 32, YR += 32 ) {
+ *YL = MPC_MULTIPLY_FLOAT_INT(facL,*L++);
+ *YR = MPC_MULTIPLY_FLOAT_INT(facR,*R++);
+ }
+ facL = MPC_MULTIPLY_SCF( Cc[d->Res_L[Band]] , d->SCF_Index_L[Band][1] & 0xFF);
+ facR = MPC_MULTIPLY_SCF( Cc[d->Res_R[Band]] , d->SCF_Index_R[Band][1] & 0xFF);
+ for (; n < 24; n++, YL += 32, YR += 32 ) {
+ *YL = MPC_MULTIPLY_FLOAT_INT(facL,*L++);
+ *YR = MPC_MULTIPLY_FLOAT_INT(facR,*R++);
+ }
+ facL = MPC_MULTIPLY_SCF( Cc[d->Res_L[Band]] , d->SCF_Index_L[Band][2] & 0xFF);
+ facR = MPC_MULTIPLY_SCF( Cc[d->Res_R[Band]] , d->SCF_Index_R[Band][2] & 0xFF);
+ for (; n < 36; n++, YL += 32, YR += 32 ) {
+ *YL = MPC_MULTIPLY_FLOAT_INT(facL,*L++);
+ *YR = MPC_MULTIPLY_FLOAT_INT(facR,*R++);
+ }
+ } else { // L!=0, R==0
+ facL = MPC_MULTIPLY_SCF( Cc[d->Res_L[Band]] , d->SCF_Index_L[Band][0] & 0xFF);
+ for ( n = 0; n < 12; n++, YL += 32, YR += 32 ) {
+ *YL = MPC_MULTIPLY_FLOAT_INT(facL,*L++);
+ *YR = 0;
+ }
+ facL = MPC_MULTIPLY_SCF( Cc[d->Res_L[Band]] , d->SCF_Index_L[Band][1] & 0xFF);
+ for ( ; n < 24; n++, YL += 32, YR += 32 ) {
+ *YL = MPC_MULTIPLY_FLOAT_INT(facL,*L++);
+ *YR = 0;
+ }
+ facL = MPC_MULTIPLY_SCF( Cc[d->Res_L[Band]] , d->SCF_Index_L[Band][2] & 0xFF);
+ for ( ; n < 36; n++, YL += 32, YR += 32 ) {
+ *YL = MPC_MULTIPLY_FLOAT_INT(facL,*L++);
+ *YR = 0;
+ }
+ }
+ }
+ else {
+ if ( d->Res_R [Band] ) { // L==0, R!=0
+ facR = MPC_MULTIPLY_SCF( Cc[d->Res_R[Band]] , d->SCF_Index_R[Band][0] & 0xFF);
+ for ( n = 0; n < 12; n++, YL += 32, YR += 32 ) {
+ *YL = 0;
+ *YR = MPC_MULTIPLY_FLOAT_INT(facR,*R++);
+ }
+ facR = MPC_MULTIPLY_SCF( Cc[d->Res_R[Band]] , d->SCF_Index_R[Band][1] & 0xFF);
+ for ( ; n < 24; n++, YL += 32, YR += 32 ) {
+ *YL = 0;
+ *YR = MPC_MULTIPLY_FLOAT_INT(facR,*R++);
+ }
+ facR = MPC_MULTIPLY_SCF( Cc[d->Res_R[Band]] , d->SCF_Index_R[Band][2] & 0xFF);
+ for ( ; n < 36; n++, YL += 32, YR += 32 ) {
+ *YL = 0;
+ *YR = MPC_MULTIPLY_FLOAT_INT(facR,*R++);
+ }
+ } else { // L==0, R==0
+ for ( n = 0; n < 36; n++, YL += 32, YR += 32 ) {
+ *YR = *YL = 0;
+ }
+ }
+ }
+ }
+ }
+}
+
+void mpc_decoder_read_bitstream_sv7(mpc_decoder * d, mpc_bits_reader * r)
+{
+ // these arrays hold decoding results for bundled quantizers (3- and 5-step)
+ static const mpc_int32_t idx30[] = { -1, 0, 1,-1, 0, 1,-1, 0, 1,-1, 0, 1,-1, 0, 1,-1, 0, 1,-1, 0, 1,-1, 0, 1,-1, 0, 1};
+ static const mpc_int32_t idx31[] = { -1,-1,-1, 0, 0, 0, 1, 1, 1,-1,-1,-1, 0, 0, 0, 1, 1, 1,-1,-1,-1, 0, 0, 0, 1, 1, 1};
+ static const mpc_int32_t idx32[] = { -1,-1,-1,-1,-1,-1,-1,-1,-1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1};
+ static const mpc_int32_t idx50[] = { -2,-1, 0, 1, 2,-2,-1, 0, 1, 2,-2,-1, 0, 1, 2,-2,-1, 0, 1, 2,-2,-1, 0, 1, 2};
+ static const mpc_int32_t idx51[] = { -2,-2,-2,-2,-2,-1,-1,-1,-1,-1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2};
+
+ mpc_int32_t n, idx, Max_used_Band = 0;
+
+ /***************************** Header *****************************/
+
+ // first subband
+ d->Res_L[0] = mpc_bits_read(r, 4);
+ d->Res_R[0] = mpc_bits_read(r, 4);
+ if (!(d->Res_L[0] == 0 && d->Res_R[0] == 0)) {
+ if (d->ms)
+ d->MS_Flag[0] = mpc_bits_read(r, 1);
+ Max_used_Band = 1;
+ }
+
+ // consecutive subbands
+ for ( n = 1; n <= d->max_band; n++ ) {
+ idx = mpc_bits_huff_lut(r, & mpc_HuffHdr);
+ d->Res_L[n] = (idx!=4) ? d->Res_L[n - 1] + idx : (int) mpc_bits_read(r, 4);
+
+ idx = mpc_bits_huff_lut(r, & mpc_HuffHdr);
+ d->Res_R[n] = (idx!=4) ? d->Res_R[n - 1] + idx : (int) mpc_bits_read(r, 4);
+
+ if (!(d->Res_L[n] == 0 && d->Res_R[n] == 0)) {
+ if (d->ms)
+ d->MS_Flag[n] = mpc_bits_read(r, 1);
+ Max_used_Band = n + 1;
+ }
+ }
+
+ /****************************** SCFI ******************************/
+ for ( n = 0; n < Max_used_Band; n++ ) {
+ if (d->Res_L[n])
+ d->SCFI_L[n] = mpc_bits_huff_dec(r, mpc_table_HuffSCFI);
+ if (d->Res_R[n])
+ d->SCFI_R[n] = mpc_bits_huff_dec(r, mpc_table_HuffSCFI);
+ }
+
+ /**************************** SCF/DSCF ****************************/
+ for ( n = 0; n < Max_used_Band; n++ ) {
+ mpc_int32_t * SCF = d->SCF_Index_L[n];
+ mpc_uint32_t Res = d->Res_L[n], SCFI = d->SCFI_L[n];
+ do {
+ if (Res) {
+ switch (SCFI) {
+ case 1:
+ idx = mpc_bits_huff_lut(r, & mpc_HuffDSCF);
+ SCF[0] = (idx!=8) ? SCF[2] + idx : (int) mpc_bits_read(r, 6);
+ idx = mpc_bits_huff_lut(r, & mpc_HuffDSCF);
+ SCF[1] = (idx!=8) ? SCF[0] + idx : (int) mpc_bits_read(r, 6);
+ SCF[2] = SCF[1];
+ break;
+ case 3:
+ idx = mpc_bits_huff_lut(r, & mpc_HuffDSCF);
+ SCF[0] = (idx!=8) ? SCF[2] + idx : (int) mpc_bits_read(r, 6);
+ SCF[1] = SCF[0];
+ SCF[2] = SCF[1];
+ break;
+ case 2:
+ idx = mpc_bits_huff_lut(r, & mpc_HuffDSCF);
+ SCF[0] = (idx!=8) ? SCF[2] + idx : (int) mpc_bits_read(r, 6);
+ SCF[1] = SCF[0];
+ idx = mpc_bits_huff_lut(r, & mpc_HuffDSCF);
+ SCF[2] = (idx!=8) ? SCF[1] + idx : (int) mpc_bits_read(r, 6);
+ break;
+ case 0:
+ idx = mpc_bits_huff_lut(r, & mpc_HuffDSCF);
+ SCF[0] = (idx!=8) ? SCF[2] + idx : (int) mpc_bits_read(r, 6);
+ idx = mpc_bits_huff_lut(r, & mpc_HuffDSCF);
+ SCF[1] = (idx!=8) ? SCF[0] + idx : (int) mpc_bits_read(r, 6);
+ idx = mpc_bits_huff_lut(r, & mpc_HuffDSCF);
+ SCF[2] = (idx!=8) ? SCF[1] + idx : (int) mpc_bits_read(r, 6);
+ break;
+ default:
+ return;
+ }
+ if (SCF[0] > 1024)
+ SCF[0] = 0x8080;
+ if (SCF[1] > 1024)
+ SCF[1] = 0x8080;
+ if (SCF[2] > 1024)
+ SCF[2] = 0x8080;
+ }
+ Res = d->Res_R[n];
+ SCFI = d->SCFI_R[n];
+ } while ( SCF == d->SCF_Index_L[n] && (SCF = d->SCF_Index_R[n]));
+ }
+
+// if (d->seeking == TRUE)
+// return;
+
+ /***************************** Samples ****************************/
+ for ( n = 0; n < Max_used_Band; n++ ) {
+ mpc_int16_t *q = d->Q[n].L, Res = d->Res_L[n];
+ do {
+ mpc_int32_t k;
+ const mpc_lut_data *Table;
+ switch (Res) {
+ case -2: case -3: case -4: case -5: case -6: case -7: case -8: case -9:
+ case -10: case -11: case -12: case -13: case -14: case -15: case -16: case -17: case 0:
+ break;
+ case -1:
+ for (k=0; k<36; k++ ) {
+ mpc_uint32_t tmp = mpc_random_int(d);
+ q[k] = ((tmp >> 24) & 0xFF) + ((tmp >> 16) & 0xFF) + ((tmp >> 8) & 0xFF) + ((tmp >> 0) & 0xFF) - 510;
+ }
+ break;
+ case 1:
+ Table = & mpc_HuffQ[0][mpc_bits_read(r, 1)];
+ for ( k = 0; k < 36; k += 3) {
+ idx = mpc_bits_huff_lut(r, Table);
+ q[k] = idx30[idx];
+ q[k + 1] = idx31[idx];
+ q[k + 2] = idx32[idx];
+ }
+ break;
+ case 2:
+ Table = & mpc_HuffQ[1][mpc_bits_read(r, 1)];
+ for ( k = 0; k < 36; k += 2) {
+ idx = mpc_bits_huff_lut(r, Table);
+ q[k] = idx50[idx];
+ q[k + 1] = idx51[idx];
+ }
+ break;
+ case 3:
+ case 4:
+ case 5:
+ case 6:
+ case 7:
+ Table = & mpc_HuffQ[Res - 1][mpc_bits_read(r, 1)];
+ for ( k = 0; k < 36; k++ )
+ q[k] = mpc_bits_huff_lut(r, Table);
+ break;
+ case 8: case 9: case 10: case 11: case 12: case 13: case 14: case 15: case 16: case 17:
+ for ( k = 0; k < 36; k++ )
+ q[k] = (mpc_int32_t)mpc_bits_read(r, Res_bit[Res]) - Dc[Res];
+ break;
+ default:
+ return;
+ }
+
+ Res = d->Res_R[n];
+ } while (q == d->Q[n].L && (q = d->Q[n].R));
+ }
+}
+
+void mpc_decoder_read_bitstream_sv8(mpc_decoder * d, mpc_bits_reader * r, mpc_bool_t is_key_frame)
+{
+ // these arrays hold decoding results for bundled quantizers (3- and 5-step)
+ static const mpc_int8_t idx50[125] = {-2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2, -2, -1, 0, 1, 2};
+ static const mpc_int8_t idx51[125] = {-2, -2, -2, -2, -2, -1, -1, -1, -1, -1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, -2, -2, -2, -2, -2, -1, -1, -1, -1, -1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, -2, -2, -2, -2, -2, -1, -1, -1, -1, -1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, -2, -2, -2, -2, -2, -1, -1, -1, -1, -1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, -2, -2, -2, -2, -2, -1, -1, -1, -1, -1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2};
+ static const mpc_int8_t idx52[125] = {-2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2};
+
+ mpc_int32_t n, Max_used_Band;
+ const mpc_can_data * Table, * Tables[2];
+
+ /***************************** Header *****************************/
+
+ if (is_key_frame == MPC_TRUE) {
+ Max_used_Band = mpc_bits_log_dec(r, d->max_band + 1);
+ } else {
+ Max_used_Band = d->last_max_band + mpc_bits_can_dec(r, & mpc_can_Bands);
+ if (Max_used_Band > 32) Max_used_Band -= 33;
+ }
+ d->last_max_band = Max_used_Band;
+
+ if (Max_used_Band) {
+ d->Res_L[Max_used_Band-1] = mpc_bits_can_dec(r, & mpc_can_Res[0]);
+ d->Res_R[Max_used_Band-1] = mpc_bits_can_dec(r, & mpc_can_Res[0]);
+ if (d->Res_L[Max_used_Band-1] > 15) d->Res_L[Max_used_Band-1] -= 17;
+ if (d->Res_R[Max_used_Band-1] > 15) d->Res_R[Max_used_Band-1] -= 17;
+ for ( n = Max_used_Band - 2; n >= 0; n--) {
+ d->Res_L[n] = mpc_bits_can_dec(r, & mpc_can_Res[d->Res_L[n + 1] > 2]) + d->Res_L[n + 1];
+ if (d->Res_L[n] > 15) d->Res_L[n] -= 17;
+ d->Res_R[n] = mpc_bits_can_dec(r, & mpc_can_Res[d->Res_R[n + 1] > 2]) + d->Res_R[n + 1];
+ if (d->Res_R[n] > 15) d->Res_R[n] -= 17;
+ }
+
+ if (d->ms) {
+ int cnt = 0, tot = 0;
+ mpc_uint32_t tmp = 0;
+ for( n = 0; n < Max_used_Band; n++)
+ if ( d->Res_L[n] != 0 || d->Res_R[n] != 0 )
+ tot++;
+ cnt = mpc_bits_log_dec(r, tot);
+ if (cnt != 0 && cnt != tot)
+ tmp = mpc_bits_enum_dec(r, mini(cnt, tot-cnt), tot);
+ if (cnt * 2 > tot) tmp = ~tmp;
+ for( n = Max_used_Band - 1; n >= 0; n--)
+ if ( d->Res_L[n] != 0 || d->Res_R[n] != 0 ) {
+ d->MS_Flag[n] = tmp & 1;
+ tmp >>= 1;
+ }
+ }
+ }
+
+ for( n = Max_used_Band; n <= d->max_band; n++)
+ d->Res_L[n] = d->Res_R[n] = 0;
+
+ /****************************** SCFI ******************************/
+ if (is_key_frame == MPC_TRUE){
+ for( n = 0; n < 32; n++)
+ d->DSCF_Flag_L[n] = d->DSCF_Flag_R[n] = 1; // new block -> force key frame
+ }
+
+ Tables[0] = & mpc_can_SCFI[0];
+ Tables[1] = & mpc_can_SCFI[1];
+ for ( n = 0; n < Max_used_Band; n++ ) {
+ int tmp = 0, cnt = -1;
+ if (d->Res_L[n]) cnt++;
+ if (d->Res_R[n]) cnt++;
+ if (cnt >= 0) {
+ tmp = mpc_bits_can_dec(r, Tables[cnt]);
+ if (d->Res_L[n]) d->SCFI_L[n] = tmp >> (2 * cnt);
+ if (d->Res_R[n]) d->SCFI_R[n] = tmp & 3;
+ }
+ }
+
+ /**************************** SCF/DSCF ****************************/
+
+ for ( n = 0; n < Max_used_Band; n++ ) {
+ mpc_int32_t * SCF = d->SCF_Index_L[n];
+ mpc_uint32_t Res = d->Res_L[n], SCFI = d->SCFI_L[n];
+ mpc_bool_t * DSCF_Flag = &d->DSCF_Flag_L[n];
+
+ do {
+ if ( Res ) {
+ int m;
+ if (*DSCF_Flag == 1) {
+ SCF[0] = (mpc_int32_t)mpc_bits_read(r, 7) - 6;
+ *DSCF_Flag = 0;
+ } else {
+ mpc_uint_t tmp = mpc_bits_can_dec(r, & mpc_can_DSCF[1]);
+ if (tmp == 64)
+ tmp += mpc_bits_read(r, 6);
+ SCF[0] = ((SCF[2] - 25 + tmp) & 127) - 6;
+ }
+ for( m = 0; m < 2; m++){
+ if (((SCFI << m) & 2) == 0) {
+ mpc_uint_t tmp = mpc_bits_can_dec(r, & mpc_can_DSCF[0]);
+ if (tmp == 31)
+ tmp = 64 + mpc_bits_read(r, 6);
+ SCF[m + 1] = ((SCF[m] - 25 + tmp) & 127) - 6;
+ } else
+ SCF[m + 1] = SCF[m];
+ }
+ }
+ Res = d->Res_R[n];
+ SCFI = d->SCFI_R[n];
+ DSCF_Flag = &d->DSCF_Flag_R[n];
+ } while ( SCF == d->SCF_Index_L[n] && (SCF = d->SCF_Index_R[n]));
+ }
+
+ /***************************** Samples ****************************/
+ for ( n = 0; n < Max_used_Band; n++ ) {
+ mpc_int16_t *q = d->Q[n].L, Res = d->Res_L[n];
+ static const unsigned int thres[] = {0, 0, 3, 0, 0, 1, 3, 4, 8};
+ static const mpc_int8_t HuffQ2_var[5*5*5] =
+ {6, 5, 4, 5, 6, 5, 4, 3, 4, 5, 4, 3, 2, 3, 4, 5, 4, 3, 4, 5, 6, 5, 4, 5, 6, 5, 4, 3, 4, 5, 4, 3, 2, 3, 4, 3, 2, 1, 2, 3, 4, 3, 2, 3, 4, 5, 4, 3, 4, 5, 4, 3, 2, 3, 4, 3, 2, 1, 2, 3, 2, 1, 0, 1, 2, 3, 2, 1, 2, 3, 4, 3, 2, 3, 4, 5, 4, 3, 4, 5, 4, 3, 2, 3, 4, 3, 2, 1, 2, 3, 4, 3, 2, 3, 4, 5, 4, 3, 4, 5, 6, 5, 4, 5, 6, 5, 4, 3, 4, 5, 4, 3, 2, 3, 4, 5, 4, 3, 4, 5, 6, 5, 4, 5, 6};
+
+ do {
+ mpc_int32_t k = 0, idx = 1;
+ if (Res != 0) {
+ if (Res == 2) {
+ Tables[0] = & mpc_can_Q [0][0];
+ Tables[1] = & mpc_can_Q [0][1];
+ idx = 2 * thres[Res];
+ for ( ; k < 36; k += 3) {
+ int tmp = mpc_bits_can_dec(r, Tables[idx > thres[Res]]);
+ q[k] = idx50[tmp];
+ q[k + 1] = idx51[tmp];
+ q[k + 2] = idx52[tmp];
+ idx = (idx >> 1) + HuffQ2_var[tmp];
+ }
+ } else if (Res == 1) {
+ Table = & mpc_can_Q1;
+ for( ; k < 36; ){
+ int kmax = k + 18;
+ mpc_uint_t cnt = mpc_bits_can_dec(r, Table);
+ idx = 0;
+ if (cnt > 0 && cnt < 18)
+ idx = mpc_bits_enum_dec(r, cnt <= 9 ? cnt : 18 - cnt, 18);
+ if (cnt > 9) idx = ~idx;
+ for ( ; k < kmax; k++) {
+ q[k] = 0;
+ if ( idx & (1 << 17) )
+ q[k] = (mpc_bits_read(r, 1) << 1) - 1;
+ idx <<= 1;
+ }
+ }
+ } else if (Res == -1) {
+ for ( ; k<36; k++ ) {
+ mpc_uint32_t tmp = mpc_random_int(d);
+ q[k] = ((tmp >> 24) & 0xFF) + ((tmp >> 16) & 0xFF) + ((tmp >> 8) & 0xFF) + ((tmp >> 0) & 0xFF) - 510;
+ }
+ } else if (Res <= 4) {
+ Table = & mpc_can_Q[1][Res - 3];
+ for ( ; k < 36; k += 2 ) {
+ union {
+ mpc_int8_t sym;
+ struct { mpc_int8_t s1:4, s2:4; };
+ } tmp;
+ tmp.sym = mpc_bits_can_dec(r, Table);
+ q[k] = tmp.s1;
+ q[k + 1] = tmp.s2;
+ }
+ } else if (Res <= 8) {
+ Tables[0] = & mpc_can_Q [Res - 3][0];
+ Tables[1] = & mpc_can_Q [Res - 3][1];
+ idx = 2 * thres[Res];
+ for ( ; k < 36; k++ ) {
+ q[k] = mpc_bits_can_dec(r, Tables[idx > thres[Res]]);
+ idx = (idx >> 1) + absi(q[k]);
+ }
+ } else {
+ for ( ; k < 36; k++ ) {
+ q[k] = (unsigned char) mpc_bits_can_dec(r, & mpc_can_Q9up);
+ if (Res != 9)
+ q[k] = (q[k] << (Res - 9)) | mpc_bits_read(r, Res - 9);
+ q[k] -= Dc[Res];
+ }
+ }
+ }
+
+ Res = d->Res_R[n];
+ } while (q == d->Q[n].L && (q = d->Q[n].R));
+ }
+}
+
diff --git a/plugins/musepack/mpc_demux.c b/plugins/musepack/mpc_demux.c
new file mode 100644
index 00000000..03bca9c3
--- /dev/null
+++ b/plugins/musepack/mpc_demux.c
@@ -0,0 +1,661 @@
+/*
+ Copyright (c) 2005-2009, The Musepack Development Team
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are
+ met:
+
+ * Redistributions of source code must retain the above copyright
+ notice, this list of conditions and the following disclaimer.
+
+ * Redistributions in binary form must reproduce the above
+ copyright notice, this list of conditions and the following
+ disclaimer in the documentation and/or other materials provided
+ with the distribution.
+
+ * Neither the name of the The Musepack Development Team nor the
+ names of its contributors may be used to endorse or promote
+ products derived from this software without specific prior
+ written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+ "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+ LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+ A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
+ OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
+ SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
+ LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
+ DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
+ THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
+ (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+ OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+*/
+
+#include <math.h>
+#include <string.h>
+#include <mpc/streaminfo.h>
+#include <mpc/mpcdec.h>
+#include "internal.h"
+#include "decoder.h"
+#include "huffman.h"
+#include "mpc_bits_reader.h"
+
+/// maximum number of seek points in the table. The distance between points will
+/// be adapted so this value is never exceeded.
+#define MAX_SEEK_TABLE_SIZE 65536
+
+// streaminfo.c
+mpc_status streaminfo_read_header_sv8(mpc_streaminfo* si,
+ const mpc_bits_reader * r_in,
+ mpc_size_t block_size);
+mpc_status streaminfo_read_header_sv7(mpc_streaminfo* si, mpc_bits_reader * r_in);
+void streaminfo_encoder_info(mpc_streaminfo* si, const mpc_bits_reader * r_in);
+void streaminfo_gain(mpc_streaminfo* si, const mpc_bits_reader * r_in);
+
+// mpc_decoder.c
+void mpc_decoder_reset_scf(mpc_decoder * d, int value);
+
+enum {
+ MPC_BUFFER_SWAP = 1,
+ MPC_BUFFER_FULL = 2,
+};
+
+static void mpc_demux_clear_buff(mpc_demux * d)
+{
+ d->bytes_total = 0;
+ d->bits_reader.buff = d->buffer;
+ d->bits_reader.count = 8;
+ d->block_bits = 0;
+ d->block_frames = 0;
+}
+
+static mpc_uint32_t
+mpc_demux_fill(mpc_demux * d, mpc_uint32_t min_bytes, int flags)
+{
+ mpc_uint32_t unread_bytes = d->bytes_total + d->buffer - d->bits_reader.buff
+ - ((8 - d->bits_reader.count) >> 3);
+ int offset = 0;
+
+ if (min_bytes == 0 || min_bytes > DEMUX_BUFFER_SIZE ||
+ (unread_bytes < min_bytes && flags & MPC_BUFFER_FULL))
+ min_bytes = DEMUX_BUFFER_SIZE;
+
+ if (unread_bytes < min_bytes) {
+ mpc_uint32_t bytes2read = min_bytes - unread_bytes;
+ mpc_uint32_t bytes_free = DEMUX_BUFFER_SIZE - d->bytes_total;
+
+ if (flags & MPC_BUFFER_SWAP) {
+ bytes2read &= -1 << 2;
+ offset = (unread_bytes + 3) & ( -1 << 2);
+ offset -= unread_bytes;
+ }
+
+ if (bytes2read > bytes_free) {
+ if (d->bits_reader.count == 0) {
+ d->bits_reader.count = 8;
+ d->bits_reader.buff++;
+ }
+ memmove(d->buffer + offset, d->bits_reader.buff, unread_bytes);
+ d->bits_reader.buff = d->buffer + offset;
+ d->bytes_total = unread_bytes + offset;
+ }
+ bytes2read = d->r->read(d->r, d->buffer + d->bytes_total, bytes2read);
+ if (flags & MPC_BUFFER_SWAP){
+ unsigned int i, * tmp = (unsigned int *) (d->buffer + d->bytes_total);
+ for(i = 0 ;i < (bytes2read >> 2); i++)
+ tmp[i] = mpc_swap32(tmp[i]);
+ }
+ d->bytes_total += bytes2read;
+ return bytes2read;
+ }
+
+ return (mpc_uint32_t) -1;
+}
+
+/**
+ * checks if a block key is valid
+ * @param key the two caracters key to check
+ * @return MPC_STATUS_INVALIDSV if the key is invalid, MPC_STATUS_OK else
+ */
+static mpc_inline mpc_status mpc_check_key(char * key)
+{
+ if (key[0] < 65 || key[0] > 90 || key[1] < 65 || key[1] > 90)
+ return MPC_STATUS_INVALIDSV;
+ return MPC_STATUS_OK;
+}
+
+/**
+ * seek to a bit position in the stream
+ * @param d demuxer context
+ * @param fpos position in the stream in bits from the beginning of mpc datas
+ * @param min_bytes number of bytes to load after seeking
+ */
+static void
+mpc_demux_seek(mpc_demux * d, mpc_seek_t fpos, mpc_uint32_t min_bytes) {
+ mpc_seek_t next_pos;
+ mpc_int_t bit_offset;
+
+ // FIXME : do not flush the buffer if fpos is in the current buffer
+
+ next_pos = fpos >> 3;
+ if (d->si.stream_version == 7)
+ next_pos = ((next_pos - d->si.header_position) & (-1 << 2)) + d->si.header_position;
+ bit_offset = (int) (fpos - (next_pos << 3));
+
+ d->r->seek(d->r, (mpc_int32_t) next_pos);
+ mpc_demux_clear_buff(d);
+ if (d->si.stream_version == 7)
+ mpc_demux_fill(d, (min_bytes + ((bit_offset + 7) >> 3) + 3) & (~3), MPC_BUFFER_SWAP);
+ else
+ mpc_demux_fill(d, min_bytes + ((bit_offset + 7) >> 3), 0);
+ d->bits_reader.buff += bit_offset >> 3;
+ d->bits_reader.count = 8 - (bit_offset & 7);
+}
+
+/**
+ * return the current position in the stream (in bits) from the beginning
+ * of the file
+ * @param d demuxer context
+ * @return current stream position in bits
+ */
+mpc_seek_t mpc_demux_pos(mpc_demux * d)
+{
+ return (((mpc_seek_t)(d->r->tell(d->r)) - d->bytes_total +
+ d->bits_reader.buff - d->buffer) << 3) + 8 - d->bits_reader.count;
+}
+
+/**
+ * Searches for a ID3v2-tag and reads the length (in bytes) of it.
+ *
+ * @param d demuxer context
+ * @return size of tag, in bytes
+ * @return MPC_STATUS_FILE on errors of any kind
+ */
+static mpc_int32_t mpc_demux_skip_id3v2(mpc_demux * d)
+{
+ mpc_uint8_t tmp [4];
+ mpc_bool_t footerPresent; // ID3v2.4-flag
+ mpc_int32_t size;
+
+ // we must be at the beginning of the stream
+ mpc_demux_fill(d, 3, 0);
+
+ // check id3-tag
+ if ( 0 != memcmp( d->bits_reader.buff, "ID3", 3 ) )
+ return 0;
+
+ mpc_demux_fill(d, 10, 0);
+
+ mpc_bits_read(&d->bits_reader, 24); // read ID3
+ mpc_bits_read(&d->bits_reader, 16); // read tag version
+
+ tmp[0] = mpc_bits_read(&d->bits_reader, 8); // read flags
+ footerPresent = tmp[0] & 0x10;
+ if ( tmp[0] & 0x0F )
+ return MPC_STATUS_FILE; // not (yet???) allowed
+
+ tmp[0] = mpc_bits_read(&d->bits_reader, 8); // read size
+ tmp[1] = mpc_bits_read(&d->bits_reader, 8); // read size
+ tmp[2] = mpc_bits_read(&d->bits_reader, 8); // read size
+ tmp[3] = mpc_bits_read(&d->bits_reader, 8); // read size
+
+ if ( (tmp[0] | tmp[1] | tmp[2] | tmp[3]) & 0x80 )
+ return MPC_STATUS_FILE; // not allowed
+
+ // read headerSize (syncsave: 4 * $0xxxxxxx = 28 significant bits)
+ size = tmp[0] << 21;
+ size |= tmp[1] << 14;
+ size |= tmp[2] << 7;
+ size |= tmp[3];
+
+ if ( footerPresent )
+ size += 10;
+
+ mpc_demux_fill(d, size, 0);
+ d->bits_reader.buff += size;
+
+ return size + 10;
+}
+
+static mpc_status mpc_demux_seek_init(mpc_demux * d)
+{
+ size_t seek_table_size;
+ if (d->seek_table != 0)
+ return MPC_STATUS_OK;
+
+ d->seek_pwr = 6;
+ if (d->si.block_pwr > d->seek_pwr)
+ d->seek_pwr = d->si.block_pwr;
+ seek_table_size = (2 + d->si.samples / (MPC_FRAME_LENGTH << d->seek_pwr));
+ while (seek_table_size > MAX_SEEK_TABLE_SIZE) {
+ d->seek_pwr++;
+ seek_table_size = (2 + d->si.samples / (MPC_FRAME_LENGTH << d->seek_pwr));
+ }
+ d->seek_table = malloc((size_t)(seek_table_size * sizeof(mpc_seek_t)));
+ if (d->seek_table == 0)
+ return MPC_STATUS_FILE;
+ d->seek_table[0] = (mpc_seek_t)mpc_demux_pos(d);
+ d->seek_table_size = 1;
+
+ return MPC_STATUS_OK;
+}
+
+static void mpc_demux_ST(mpc_demux * d)
+{
+ mpc_uint64_t tmp;
+ mpc_seek_t * table, last[2];
+ mpc_bits_reader r = d->bits_reader;
+ mpc_uint_t i, diff_pwr = 0, mask;
+ mpc_uint32_t file_table_size;
+
+ if (d->seek_table != 0)
+ return;
+
+ mpc_bits_get_size(&r, &tmp);
+ file_table_size = (mpc_seek_t) tmp;
+ d->seek_pwr = d->si.block_pwr + mpc_bits_read(&r, 4);
+
+ tmp = 2 + d->si.samples / (MPC_FRAME_LENGTH << d->seek_pwr);
+ while (tmp > MAX_SEEK_TABLE_SIZE) {
+ d->seek_pwr++;
+ diff_pwr++;
+ tmp = 2 + d->si.samples / (MPC_FRAME_LENGTH << d->seek_pwr);
+ }
+ if ((file_table_size >> diff_pwr) > tmp)
+ file_table_size = tmp << diff_pwr;
+ d->seek_table = malloc((size_t) (tmp * sizeof(mpc_seek_t)));
+ d->seek_table_size = (file_table_size + ((1 << diff_pwr) - 1)) >> diff_pwr;
+
+ table = d->seek_table;
+ mpc_bits_get_size(&r, &tmp);
+ table[0] = last[0] = (mpc_seek_t) (tmp + d->si.header_position) * 8;
+
+ if (d->seek_table_size == 1)
+ return;
+
+ mpc_bits_get_size(&r, &tmp);
+ last[1] = (mpc_seek_t) (tmp + d->si.header_position) * 8;
+ if (diff_pwr == 0) table[1] = last[1];
+
+ mask = (1 << diff_pwr) - 1;
+ for (i = 2; i < file_table_size; i++) {
+ int code = mpc_bits_golomb_dec(&r, 12);
+ if (code & 1)
+ code = -(code & (-1 << 1));
+ code <<= 2;
+ last[i & 1] = code + 2 * last[(i-1) & 1] - last[i & 1];
+ if ((i & mask) == 0)
+ table[i >> diff_pwr] = last[i & 1];
+ }
+}
+
+static void mpc_demux_SP(mpc_demux * d, int size, int block_size)
+{
+ mpc_seek_t cur;
+ mpc_uint64_t ptr;
+ mpc_block b;
+ int st_head_size;
+
+ cur = mpc_demux_pos(d);
+ mpc_bits_get_size(&d->bits_reader, &ptr);
+ mpc_demux_seek(d, (ptr - size) * 8 + cur, 11);
+ st_head_size = mpc_bits_get_block(&d->bits_reader, &b);
+ if (memcmp(b.key, "ST", 2) == 0) {
+ d->chap_pos = (ptr - size + b.size + st_head_size) * 8 + cur;
+ d->chap_nb = -1;
+ mpc_demux_fill(d, (mpc_uint32_t) b.size, 0);
+ mpc_demux_ST(d);
+ }
+ mpc_demux_seek(d, cur, 11 + block_size);
+}
+
+static void mpc_demux_chap_find(mpc_demux * d)
+{
+ mpc_block b;
+ int tag_size = 0, chap_size = 0, size, i = 0;
+
+ d->chap_nb = 0;
+
+ if (d->si.stream_version < 8)
+ return;
+
+ if (d->chap_pos == 0) {
+ mpc_uint64_t cur_pos = (d->si.header_position + 4) * 8;
+ mpc_demux_seek(d, cur_pos, 11); // seek to the beginning of the stream
+ size = mpc_bits_get_block(&d->bits_reader, &b);
+ while (memcmp(b.key, "SE", 2) != 0) {
+ if (mpc_check_key(b.key) != MPC_STATUS_OK)
+ return;
+ if (memcmp(b.key, "CT", 2) == 0) {
+ if (d->chap_pos == 0) d->chap_pos = cur_pos;
+ } else
+ d->chap_pos = 0;
+ cur_pos += (size + b.size) * 8;
+ mpc_demux_seek(d, cur_pos, 11);
+ size = mpc_bits_get_block(&d->bits_reader, &b);
+ }
+ if (d->chap_pos == 0)
+ d->chap_pos = cur_pos;
+ }
+
+ mpc_demux_seek(d, d->chap_pos, 20);
+ size = mpc_bits_get_block(&d->bits_reader, &b);
+ while (memcmp(b.key, "CT", 2) == 0) {
+ mpc_uint64_t chap_sample;
+ d->chap_nb++;
+ chap_size += size;
+ size = mpc_bits_get_size(&d->bits_reader, &chap_sample) + 4;
+ chap_size += size;
+ tag_size += b.size - size;
+ mpc_demux_seek(d, d->chap_pos + (chap_size + tag_size) * 8, 20);
+ size = mpc_bits_get_block(&d->bits_reader, &b);
+ }
+
+ if (d->chap_nb > 0) {
+ char * ptag;
+ d->chap = malloc(sizeof(mpc_chap_info) * d->chap_nb + tag_size);
+ ptag = (char*)(d->chap + d->chap_nb);
+
+ mpc_demux_seek(d, d->chap_pos, 11);
+ size = mpc_bits_get_block(&d->bits_reader, &b);
+ while (memcmp(b.key, "CT", 2) == 0) {
+ mpc_demux_fill(d, 11 + (mpc_uint32_t) b.size, 0);
+ size = mpc_bits_get_size(&d->bits_reader, &d->chap[i].sample) + 4;
+ d->chap[i].gain = (mpc_uint16_t) mpc_bits_read(&d->bits_reader, 16);
+ d->chap[i].peak = (mpc_uint16_t) mpc_bits_read(&d->bits_reader, 16);
+ memcpy(ptag, d->bits_reader.buff + ((8 - d->bits_reader.count) >> 3), b.size - size);
+ d->bits_reader.buff += b.size - size;
+ d->chap[i].tag_size = b.size - size;
+ d->chap[i].tag = ptag;
+ ptag += b.size - size;
+ i++;
+ size = mpc_bits_get_block(&d->bits_reader, &b);
+ }
+ }
+
+ d->bits_reader.buff -= size;
+}
+
+/**
+ * Gets the number of chapters in the stream
+ * @param d pointer to a musepack demuxer
+ * @return the number of chapters found in the stream
+ */
+mpc_int_t mpc_demux_chap_nb(mpc_demux * d)
+{
+ if (d->chap_nb == -1)
+ mpc_demux_chap_find(d);
+ return d->chap_nb;
+}
+
+/**
+ * Gets datas associated to a given chapter
+ * The chapter tag is an APEv2 tag without the preamble
+ * @param d pointer to a musepack demuxer
+ * @param chap_nb chapter number you want datas (from 0 to mpc_demux_chap_nb(d) - 1)
+ * @return the chapter information structure
+ */
+mpc_chap_info const * mpc_demux_chap(mpc_demux * d, int chap_nb)
+{
+ if (d->chap_nb == -1)
+ mpc_demux_chap_find(d);
+ if (chap_nb >= d->chap_nb || chap_nb < 0)
+ return 0;
+ return &d->chap[chap_nb];
+}
+
+static mpc_status mpc_demux_header(mpc_demux * d)
+{
+ char magic[4];
+
+ d->si.pns = 0xFF;
+ d->si.profile_name = "n.a.";
+
+ // get header position
+ d->si.header_position = mpc_demux_skip_id3v2(d);
+ if(d->si.header_position < 0) return MPC_STATUS_FILE;
+
+ d->si.tag_offset = d->si.total_file_length = d->r->get_size(d->r);
+
+ mpc_demux_fill(d, 4, 0);
+ magic[0] = mpc_bits_read(&d->bits_reader, 8);
+ magic[1] = mpc_bits_read(&d->bits_reader, 8);
+ magic[2] = mpc_bits_read(&d->bits_reader, 8);
+ magic[3] = mpc_bits_read(&d->bits_reader, 8);
+
+ if (memcmp(magic, "MP+", 3) == 0) {
+ d->si.stream_version = magic[3] & 15;
+ d->si.pns = magic[3] >> 4;
+ if (d->si.stream_version == 7) {
+ mpc_status ret;
+ mpc_demux_fill(d, 6 * 4, MPC_BUFFER_SWAP); // header block size + endian convertion
+ ret = streaminfo_read_header_sv7(&d->si, &d->bits_reader);
+ if (ret != MPC_STATUS_OK) return ret;
+ } else {
+ return MPC_STATUS_INVALIDSV;
+ }
+ } else if (memcmp(magic, "MPCK", 4) == 0) {
+ mpc_block b;
+ int size;
+ mpc_demux_fill(d, 11, 0); // max header block size
+ size = mpc_bits_get_block(&d->bits_reader, &b);
+ while( memcmp(b.key, "AP", 2) != 0 ){ // scan all blocks until audio
+ if (mpc_check_key(b.key) != MPC_STATUS_OK)
+ return MPC_STATUS_INVALIDSV;
+ if (b.size > (mpc_uint64_t) DEMUX_BUFFER_SIZE - 11)
+ return MPC_STATUS_INVALIDSV;
+ mpc_demux_fill(d, 11 + (mpc_uint32_t) b.size, 0);
+ if (memcmp(b.key, "SH", 2) == 0){
+ int ret = streaminfo_read_header_sv8(&d->si, &d->bits_reader, (mpc_uint32_t) b.size);
+ if (ret != MPC_STATUS_OK) return ret;
+ } else if (memcmp(b.key, "RG", 2) == 0)
+ streaminfo_gain(&d->si, &d->bits_reader);
+ else if (memcmp(b.key, "EI", 2) == 0)
+ streaminfo_encoder_info(&d->si, &d->bits_reader);
+ else if (memcmp(b.key, "SO", 2) == 0)
+ mpc_demux_SP(d, size, (mpc_uint32_t) b.size);
+ else if (memcmp(b.key, "ST", 2) == 0)
+ mpc_demux_ST(d);
+ d->bits_reader.buff += b.size;
+ size = mpc_bits_get_block(&d->bits_reader, &b);
+ }
+ d->bits_reader.buff -= size;
+ if (d->si.stream_version == 0) // si not initialized !!!
+ return MPC_STATUS_INVALIDSV;
+ } else
+ return MPC_STATUS_INVALIDSV;
+
+ return MPC_STATUS_OK;
+}
+
+mpc_demux * mpc_demux_init(mpc_reader * p_reader)
+{
+ mpc_demux* p_tmp = malloc(sizeof(mpc_demux));
+
+ if (p_tmp != 0) {
+ memset(p_tmp, 0, sizeof(mpc_demux));
+ p_tmp->r = p_reader;
+ p_tmp->chap_nb = -1;
+ mpc_demux_clear_buff(p_tmp);
+ if (mpc_demux_header(p_tmp) == MPC_STATUS_OK &&
+ mpc_demux_seek_init(p_tmp) == MPC_STATUS_OK) {
+ p_tmp->d = mpc_decoder_init(&p_tmp->si);
+ } else {
+ if (p_tmp->seek_table)
+ free(p_tmp->seek_table);
+ free(p_tmp);
+ p_tmp = 0;
+ }
+ }
+
+ return p_tmp;
+}
+
+void mpc_demux_exit(mpc_demux * d)
+{
+ mpc_decoder_exit(d->d);
+ free(d->seek_table);
+ free(d->chap);
+ free(d);
+}
+
+void mpc_demux_get_info(mpc_demux * d, mpc_streaminfo * i)
+{
+ memcpy(i, &d->si, sizeof d->si);
+}
+
+mpc_status mpc_demux_decode(mpc_demux * d, mpc_frame_info * i)
+{
+ mpc_bits_reader r;
+ if (d->si.stream_version >= 8) {
+ i->is_key_frame = MPC_FALSE;
+
+ if (d->block_frames == 0) {
+ mpc_block b = {{0,0},0};
+ d->bits_reader.count &= -8;
+ if (d->d->decoded_samples == (d->seek_table_size << d->seek_pwr) * MPC_FRAME_LENGTH) {
+ d->seek_table[d->seek_table_size] = (mpc_seek_t) mpc_demux_pos(d);
+ d->seek_table_size ++;
+ }
+ mpc_demux_fill(d, 11, 0); // max header block size
+ mpc_bits_get_block(&d->bits_reader, &b);
+ while( memcmp(b.key, "AP", 2) != 0 ) { // scan all blocks until audio
+ if (mpc_check_key(b.key) != MPC_STATUS_OK)
+ goto error;
+ if (memcmp(b.key, "SE", 2) == 0) { // end block
+ i->bits = -1;
+ return MPC_STATUS_OK;
+ }
+ if (mpc_demux_fill(d, 11 + (mpc_uint32_t) b.size, 0) == 0)
+ goto error;
+ d->bits_reader.buff += b.size;
+ mpc_bits_get_block(&d->bits_reader, &b);
+ }
+ d->block_bits = (mpc_uint32_t) b.size * 8;
+ d->block_frames = 1 << d->si.block_pwr;
+ i->is_key_frame = MPC_TRUE;
+ }
+ if (d->buffer + d->bytes_total - d->bits_reader.buff <= MAX_FRAME_SIZE)
+ mpc_demux_fill(d, (d->block_bits >> 3) + 1, 0);
+ r = d->bits_reader;
+ mpc_decoder_decode_frame(d->d, &d->bits_reader, i);
+ d->block_bits -= ((d->bits_reader.buff - r.buff) << 3) + r.count - d->bits_reader.count;
+ d->block_frames--;
+ if (d->block_bits < 0 || (d->block_frames == 0 && d->block_bits > 7))
+ goto error;
+ } else {
+ if (d->d->decoded_samples == (d->seek_table_size << d->seek_pwr) * MPC_FRAME_LENGTH) {
+ d->seek_table[d->seek_table_size] = (mpc_seek_t) mpc_demux_pos(d);
+ d->seek_table_size ++;
+ }
+ mpc_demux_fill(d, MAX_FRAME_SIZE, MPC_BUFFER_FULL | MPC_BUFFER_SWAP);
+ d->block_bits = (mpc_int_t) mpc_bits_read(&d->bits_reader, 20); // read frame size
+ if (MPC_FRAME_LENGTH > d->d->samples - d->d->decoded_samples - 1) d->block_bits += 11; // we will read last frame size
+ r = d->bits_reader;
+ mpc_decoder_decode_frame(d->d, &d->bits_reader, i);
+ if (i->bits != -1 && d->block_bits != ((d->bits_reader.buff - r.buff) << 3) + r.count - d->bits_reader.count)
+ goto error;
+ }
+ if (i->bits != -1 && d->buffer + d->bytes_total < d->bits_reader.buff + ((8 - d->bits_reader.count) >> 3))
+ goto error;
+
+ return MPC_STATUS_OK;
+error:
+ i->bits = -1; // we pretend it's end of file
+ return MPC_STATUS_INVALIDSV;
+}
+
+mpc_status mpc_demux_seek_second(mpc_demux * d, double seconds)
+{
+ return mpc_demux_seek_sample(d, (mpc_int64_t)(seconds * (double)d->si.sample_freq + 0.5));
+}
+
+mpc_status mpc_demux_seek_sample(mpc_demux * d, mpc_uint64_t destsample)
+{
+ mpc_uint32_t fwd, samples_to_skip, i;
+ mpc_uint32_t block_samples = MPC_FRAME_LENGTH << d->si.block_pwr;
+ mpc_seek_t fpos;
+
+ destsample += d->si.beg_silence;
+ if (destsample > d->si.samples) destsample = d->si.samples;
+ fwd = (mpc_uint32_t) (destsample / block_samples);
+ samples_to_skip = MPC_DECODER_SYNTH_DELAY +
+ (mpc_uint32_t) (destsample % block_samples);
+ if (d->si.stream_version == 7) {
+ if (fwd > 32) {
+ fwd -= 32;
+ samples_to_skip += MPC_FRAME_LENGTH * 32;
+ } else {
+ samples_to_skip += MPC_FRAME_LENGTH * fwd;
+ fwd = 0;
+ }
+ }
+
+ i = fwd >> (d->seek_pwr - d->si.block_pwr);
+ if (i >= d->seek_table_size)
+ i = d->seek_table_size - 1;
+ fpos = d->seek_table[i];
+ i <<= d->seek_pwr - d->si.block_pwr;
+ d->d->decoded_samples = i * block_samples;
+
+ if (d->si.stream_version >= 8) {
+ mpc_block b;
+ int size;
+ mpc_demux_seek(d, fpos, 11);
+ size = mpc_bits_get_block(&d->bits_reader, &b);
+ while(i < fwd) {
+ if (memcmp(b.key, "AP", 2) == 0) {
+ if (d->d->decoded_samples == (d->seek_table_size << d->seek_pwr) * MPC_FRAME_LENGTH) {
+ d->seek_table[d->seek_table_size] = (mpc_seek_t) mpc_demux_pos(d) - 8 * size;
+ d->seek_table_size ++;
+ }
+ d->d->decoded_samples += block_samples;
+ i++;
+ }
+ fpos += ((mpc_uint32_t)b.size + size) * 8;
+ mpc_demux_seek(d, fpos, 11);
+ size = mpc_bits_get_block(&d->bits_reader, &b);
+ }
+ d->bits_reader.buff -= size;
+ } else {
+ mpc_decoder_reset_scf(d->d, fwd != 0);
+ mpc_demux_seek(d, fpos, 4);
+ for( ; i < fwd; i++){
+ if (d->d->decoded_samples == (d->seek_table_size << d->seek_pwr) * MPC_FRAME_LENGTH) {
+ d->seek_table[d->seek_table_size] = (mpc_seek_t) mpc_demux_pos(d);
+ d->seek_table_size ++;
+ }
+ d->d->decoded_samples += block_samples;
+ fpos += mpc_bits_read(&d->bits_reader, 20) + 20;
+ mpc_demux_seek(d, fpos, 4);
+ }
+ }
+ d->d->samples_to_skip = samples_to_skip;
+ return MPC_STATUS_OK;
+}
+
+void mpc_set_replay_level(mpc_demux * d, float level, mpc_bool_t use_gain,
+ mpc_bool_t use_title, mpc_bool_t clip_prevention)
+{
+ float peak = use_title ? d->si.peak_title : d->si.peak_album;
+ float gain = use_title ? d->si.gain_title : d->si.gain_album;
+
+ if(!use_gain && !clip_prevention)
+ return;
+
+ if(!peak)
+ peak = 1.;
+ else
+ peak = (1 << 15) / pow(10, peak / (20 * 256));
+
+ if(!gain)
+ gain = 1.;
+ else
+ gain = pow(10, (level - gain / 256) / 20);
+
+ if(clip_prevention && (peak < gain || !use_gain))
+ gain = peak;
+
+ mpc_decoder_scale_output(d->d, gain);
+}
diff --git a/plugins/musepack/mpc_reader.c b/plugins/musepack/mpc_reader.c
new file mode 100644
index 00000000..e1d151fe
--- /dev/null
+++ b/plugins/musepack/mpc_reader.c
@@ -0,0 +1,144 @@
+/*
+ Copyright (c) 2005-2009, The Musepack Development Team
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are
+ met:
+
+ * Redistributions of source code must retain the above copyright
+ notice, this list of conditions and the following disclaimer.
+
+ * Redistributions in binary form must reproduce the above
+ copyright notice, this list of conditions and the following
+ disclaimer in the documentation and/or other materials provided
+ with the distribution.
+
+ * Neither the name of the The Musepack Development Team nor the
+ names of its contributors may be used to endorse or promote
+ products derived from this software without specific prior
+ written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+ "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+ LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+ A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
+ OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
+ SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
+ LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
+ DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
+ THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
+ (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+ OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+*/
+/// \file mpc_reader.c
+/// Contains implementations for simple file-based mpc_reader
+#include <mpc/reader.h>
+#include "internal.h"
+#include <stdio.h>
+
+#define STDIO_MAGIC 0xF34B963C ///< Just a random safe-check value...
+typedef struct mpc_reader_stdio_t {
+ FILE *p_file;
+ int file_size;
+ mpc_bool_t is_seekable;
+ mpc_int32_t magic;
+} mpc_reader_stdio;
+
+/// mpc_reader callback implementations
+static mpc_int32_t
+read_stdio(mpc_reader *p_reader, void *ptr, mpc_int32_t size)
+{
+ mpc_reader_stdio *p_stdio = (mpc_reader_stdio*) p_reader->data;
+ if(p_stdio->magic != STDIO_MAGIC) return MPC_STATUS_FILE;
+ return (mpc_int32_t) fread(ptr, 1, size, p_stdio->p_file);
+}
+
+static mpc_bool_t
+seek_stdio(mpc_reader *p_reader, mpc_int32_t offset)
+{
+ mpc_reader_stdio *p_stdio = (mpc_reader_stdio*) p_reader->data;
+ if(p_stdio->magic != STDIO_MAGIC) return MPC_FALSE;
+ return p_stdio->is_seekable ? fseek(p_stdio->p_file, offset, SEEK_SET) == 0 : MPC_FALSE;
+}
+
+static mpc_int32_t
+tell_stdio(mpc_reader *p_reader)
+{
+ mpc_reader_stdio *p_stdio = (mpc_reader_stdio*) p_reader->data;
+ if(p_stdio->magic != STDIO_MAGIC) return MPC_STATUS_FILE;
+ return ftell(p_stdio->p_file);
+}
+
+static mpc_int32_t
+get_size_stdio(mpc_reader *p_reader)
+{
+ mpc_reader_stdio *p_stdio = (mpc_reader_stdio*) p_reader->data;
+ if(p_stdio->magic != STDIO_MAGIC) return MPC_STATUS_FILE;
+ return p_stdio->file_size;
+}
+
+static mpc_bool_t
+canseek_stdio(mpc_reader *p_reader)
+{
+ mpc_reader_stdio *p_stdio = (mpc_reader_stdio*) p_reader->data;
+ if(p_stdio->magic != STDIO_MAGIC) return MPC_FALSE;
+ return p_stdio->is_seekable;
+}
+
+mpc_status
+mpc_reader_init_stdio_stream(mpc_reader * p_reader, FILE * p_file)
+{
+ mpc_reader tmp_reader; mpc_reader_stdio *p_stdio; int err;
+
+ p_stdio = NULL;
+ memset(&tmp_reader, 0, sizeof tmp_reader);
+ p_stdio = malloc(sizeof *p_stdio);
+ if(!p_stdio) return MPC_STATUS_FILE;
+ memset(p_stdio, 0, sizeof *p_stdio);
+
+ p_stdio->magic = STDIO_MAGIC;
+ p_stdio->p_file = p_file;
+ p_stdio->is_seekable = MPC_TRUE;
+ err = fseek(p_stdio->p_file, 0, SEEK_END);
+ if(err < 0) goto clean;
+ err = ftell(p_stdio->p_file);
+ if(err < 0) goto clean;
+ p_stdio->file_size = err;
+ err = fseek(p_stdio->p_file, 0, SEEK_SET);
+ if(err < 0) goto clean;
+
+ tmp_reader.data = p_stdio;
+ tmp_reader.canseek = canseek_stdio;
+ tmp_reader.get_size = get_size_stdio;
+ tmp_reader.read = read_stdio;
+ tmp_reader.seek = seek_stdio;
+ tmp_reader.tell = tell_stdio;
+
+ *p_reader = tmp_reader;
+ return MPC_STATUS_OK;
+clean:
+ if(p_stdio && p_stdio->p_file)
+ fclose(p_stdio->p_file);
+ free(p_stdio);
+ return MPC_STATUS_FILE;
+}
+
+mpc_status
+mpc_reader_init_stdio(mpc_reader *p_reader, const char *filename)
+{
+ FILE * stream = fopen(filename, "rb");
+ if (stream == NULL) return MPC_STATUS_FILE;
+ return mpc_reader_init_stdio_stream(p_reader,stream);
+}
+
+void
+mpc_reader_exit_stdio(mpc_reader *p_reader)
+{
+ mpc_reader_stdio *p_stdio = (mpc_reader_stdio*) p_reader->data;
+ if(p_stdio->magic != STDIO_MAGIC) return;
+ fclose(p_stdio->p_file);
+ free(p_stdio);
+ p_reader->data = NULL;
+}
+
diff --git a/plugins/musepack/mpcdec.h b/plugins/musepack/mpcdec.h
new file mode 100644
index 00000000..c7235951
--- /dev/null
+++ b/plugins/musepack/mpcdec.h
@@ -0,0 +1,148 @@
+/*
+ Copyright (c) 2005-2009, The Musepack Development Team
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are
+ met:
+
+ * Redistributions of source code must retain the above copyright
+ notice, this list of conditions and the following disclaimer.
+
+ * Redistributions in binary form must reproduce the above
+ copyright notice, this list of conditions and the following
+ disclaimer in the documentation and/or other materials provided
+ with the distribution.
+
+ * Neither the name of the The Musepack Development Team nor the
+ names of its contributors may be used to endorse or promote
+ products derived from this software without specific prior
+ written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+ "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+ LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+ A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
+ OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
+ SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
+ LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
+ DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
+ THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
+ (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+ OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+*/
+/// \file mpcdec.h
+/// Top level include file for libmpcdec.
+#ifndef _MPCDEC_H_
+#define _MPCDEC_H_
+#ifdef WIN32
+#pragma once
+#endif
+
+#include "reader.h"
+#include "streaminfo.h"
+
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+enum {
+ MPC_FRAME_LENGTH = (36 * 32), ///< Samples per mpc frame
+ MPC_DECODER_BUFFER_LENGTH = (MPC_FRAME_LENGTH * 4), ///< Required buffer size for decoder
+ MPC_DECODER_SYNTH_DELAY = 481
+};
+
+typedef struct mpc_decoder_t mpc_decoder;
+typedef struct mpc_demux_t mpc_demux;
+
+typedef struct mpc_bits_reader_t {
+ unsigned char * buff; /// pointer on current byte
+ unsigned int count; /// unread bits in current byte
+} mpc_bits_reader;
+
+typedef struct mpc_frame_info_t {
+ mpc_uint32_t samples; /// number of samples in the frame (counting once for multiple channels)
+ mpc_int32_t bits; /// number of bits consumed by this frame (-1) if end of stream
+ MPC_SAMPLE_FORMAT * buffer; /// frame samples buffer (size = samples * channels * sizeof(MPC_SAMPLE_FORMAT))
+ mpc_bool_t is_key_frame; /// 1 if this frame is a key frame (first in block) 0 else. Set by the demuxer.
+} mpc_frame_info;
+
+typedef struct mpc_chap_info_t {
+ mpc_uint64_t sample; /// sample where the chapter starts
+ mpc_uint16_t gain; /// replaygain chapter value
+ mpc_uint16_t peak; /// peak chapter loudness level
+ mpc_uint_t tag_size; /// size of the tag element (0 if no tag is present for this chapter)
+ char * tag; /// pointer to an APEv2 tag without the preamble
+} mpc_chap_info;
+
+/// Initializes mpc decoder with the supplied stream info parameters.
+/// \param si streaminfo structure indicating format of source stream
+/// \return pointer on the initialized decoder structure if successful, 0 if not
+MPC_API mpc_decoder * mpc_decoder_init(mpc_streaminfo *si);
+
+/// Releases input mpc decoder
+MPC_API void mpc_decoder_exit(mpc_decoder *p_dec);
+
+/**
+ * Sets decoder sample scaling factor. All decoded samples will be multiplied
+ * by this factor. Useful for applying replay gain.
+ * @param scale_factor multiplicative scaling factor
+ */
+MPC_API void mpc_decoder_scale_output(mpc_decoder *p_dec, double scale_factor);
+
+MPC_API void mpc_decoder_decode_frame(mpc_decoder * d, mpc_bits_reader * r, mpc_frame_info * i);
+
+// This is the gain reference used in old replaygain
+#define MPC_OLD_GAIN_REF 64.82
+
+/**
+ * init demuxer
+ * @param p_reader initialized mpc_reader pointer
+ * @return an initialized mpc_demux pointer
+ */
+MPC_API mpc_demux * mpc_demux_init(mpc_reader * p_reader);
+/// free demuxer
+MPC_API void mpc_demux_exit(mpc_demux * d);
+/**
+ * Calls mpc_decoder_scale_output to set the scaling factor according to the
+ * replay gain stream information and the supplied ouput level
+ * @param d pointer to a musepack demuxer
+ * @param level the desired ouput level (in db). Must be MPC_OLD_GAIN_REF (64.82 db) if you want to get the old replaygain behavior
+ * @param use_gain set it to MPC_TRUE if you want to set the scaling factor according to the stream gain
+ * @param use_title MPC_TRUE : uses the title gain, MPC_FALSE : uses the album gain
+ * @param clip_prevention MPC_TRUE : uses cliping prevention
+ */
+MPC_API void mpc_set_replay_level(mpc_demux * d, float level, mpc_bool_t use_gain,
+ mpc_bool_t use_title, mpc_bool_t clip_prevention);
+/// decode frame
+MPC_API mpc_status mpc_demux_decode(mpc_demux * d, mpc_frame_info * i);
+/// get streaminfo
+MPC_API void mpc_demux_get_info(mpc_demux * d, mpc_streaminfo * i);
+/// seeks to a given sample
+MPC_API mpc_status mpc_demux_seek_sample(mpc_demux * d, mpc_uint64_t destsample);
+/// seeks to a given second
+MPC_API mpc_status mpc_demux_seek_second(mpc_demux * d, double seconds);
+
+/// \return the current position in the stream (in bits) from the beginning of the file
+MPC_API mpc_seek_t mpc_demux_pos(mpc_demux * d);
+
+/// chapters : only for sv8 streams
+/**
+ * Gets the number of chapters in the stream
+ * @param d pointer to a musepack demuxer
+ * @return the number of chapters found in the stream
+ */
+MPC_API mpc_int_t mpc_demux_chap_nb(mpc_demux * d);
+/**
+ * Gets datas associated to a given chapter
+ * The chapter tag is an APEv2 tag without the preamble
+ * @param d pointer to a musepack demuxer
+ * @param chap_nb chapter number you want datas (from 0 to mpc_demux_chap_nb(d) - 1)
+ * @return the chapter information structure
+ */
+MPC_API mpc_chap_info const * mpc_demux_chap(mpc_demux * d, int chap_nb);
+
+#ifdef __cplusplus
+}
+#endif
+#endif
diff --git a/plugins/musepack/mpcdec_math.h b/plugins/musepack/mpcdec_math.h
new file mode 100644
index 00000000..27ae69c4
--- /dev/null
+++ b/plugins/musepack/mpcdec_math.h
@@ -0,0 +1,135 @@
+/*
+ Copyright (c) 2005-2009, The Musepack Development Team
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are
+ met:
+
+ * Redistributions of source code must retain the above copyright
+ notice, this list of conditions and the following disclaimer.
+
+ * Redistributions in binary form must reproduce the above
+ copyright notice, this list of conditions and the following
+ disclaimer in the documentation and/or other materials provided
+ with the distribution.
+
+ * Neither the name of the The Musepack Development Team nor the
+ names of its contributors may be used to endorse or promote
+ products derived from this software without specific prior
+ written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+ "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+ LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+ A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
+ OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
+ SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
+ LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
+ DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
+ THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
+ (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+ OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+*/
+/// \file math.h
+/// Libmpcdec internal math routines.
+#ifndef _MPCDEC_MATH_H_
+#define _MPCDEC_MATH_H_
+#ifdef WIN32
+#pragma once
+#endif
+
+#include <mpc/mpc_types.h>
+
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#ifdef MPC_FIXED_POINT
+
+#ifdef _WIN32_WCE
+#include <cmnintrin.h>
+#define MPC_HAVE_MULHIGH
+#endif
+
+//in fixedpoint mode, results in decode output buffer are in -MPC_FIXED_POINT_SCALE ... MPC_FIXED_POINT_SCALE range
+
+typedef mpc_int64_t MPC_SAMPLE_FORMAT_MULTIPLY;
+
+#define MAKE_MPC_SAMPLE(X) (MPC_SAMPLE_FORMAT)((double)(X) * (double)(((mpc_int64_t)1)<<MPC_FIXED_POINT_FRACTPART))
+#define MAKE_MPC_SAMPLE_EX(X,Y) (MPC_SAMPLE_FORMAT)((double)(X) * (double)(((mpc_int64_t)1)<<(Y)))
+
+#define MPC_MULTIPLY_NOTRUNCATE(X,Y) \
+ (((MPC_SAMPLE_FORMAT_MULTIPLY)(X) * (MPC_SAMPLE_FORMAT_MULTIPLY)(Y)) >> MPC_FIXED_POINT_FRACTPART)
+
+#define MPC_MULTIPLY_EX_NOTRUNCATE(X,Y,Z) \
+ (((MPC_SAMPLE_FORMAT_MULTIPLY)(X) * (MPC_SAMPLE_FORMAT_MULTIPLY)(Y)) >> (Z))
+
+#ifdef _DEBUG
+static mpc_inline MPC_SAMPLE_FORMAT MPC_MULTIPLY(MPC_SAMPLE_FORMAT item1,MPC_SAMPLE_FORMAT item2)
+{
+ MPC_SAMPLE_FORMAT_MULTIPLY temp = MPC_MULTIPLY_NOTRUNCATE(item1,item2);
+ assert(temp == (MPC_SAMPLE_FORMAT_MULTIPLY)(MPC_SAMPLE_FORMAT)temp);
+ return (MPC_SAMPLE_FORMAT)temp;
+}
+
+static mpc_inline MPC_SAMPLE_FORMAT MPC_MULTIPLY_EX(MPC_SAMPLE_FORMAT item1,MPC_SAMPLE_FORMAT item2,unsigned shift)
+{
+ MPC_SAMPLE_FORMAT_MULTIPLY temp = MPC_MULTIPLY_EX_NOTRUNCATE(item1,item2,shift);
+ assert(temp == (MPC_SAMPLE_FORMAT_MULTIPLY)(MPC_SAMPLE_FORMAT)temp);
+ return (MPC_SAMPLE_FORMAT)temp;
+}
+
+#else
+
+#define MPC_MULTIPLY(X,Y) ((MPC_SAMPLE_FORMAT)MPC_MULTIPLY_NOTRUNCATE(X,Y))
+#define MPC_MULTIPLY_EX(X,Y,Z) ((MPC_SAMPLE_FORMAT)MPC_MULTIPLY_EX_NOTRUNCATE(X,Y,Z))
+
+#endif
+
+#ifdef MPC_HAVE_MULHIGH
+#define MPC_MULTIPLY_FRACT(X,Y) _MulHigh(X,Y)
+#else
+#define MPC_MULTIPLY_FRACT(X,Y) MPC_MULTIPLY_EX(X,Y,32)
+#endif
+
+#define MPC_MAKE_FRACT_CONST(X) (MPC_SAMPLE_FORMAT)((X) * (double)(((mpc_int64_t)1)<<32) )
+#define MPC_MULTIPLY_FRACT_CONST(X,Y) MPC_MULTIPLY_FRACT(X,MPC_MAKE_FRACT_CONST(Y))
+#define MPC_MULTIPLY_FRACT_CONST_FIX(X,Y,Z) ( MPC_MULTIPLY_FRACT(X,MPC_MAKE_FRACT_CONST( Y / (1<<(Z)) )) << (Z) )
+#define MPC_MULTIPLY_FRACT_CONST_SHR(X,Y,Z) MPC_MULTIPLY_FRACT(X,MPC_MAKE_FRACT_CONST( Y / (1<<(Z)) ))
+
+#define MPC_MULTIPLY_FLOAT_INT(X,Y) ((X)*(Y))
+#define MPC_SCALE_CONST(X,Y,Z) MPC_MULTIPLY_EX(X,MAKE_MPC_SAMPLE_EX(Y,Z),(Z))
+#define MPC_SCALE_CONST_SHL(X,Y,Z,S) MPC_MULTIPLY_EX(X,MAKE_MPC_SAMPLE_EX(Y,Z),(Z)-(S))
+#define MPC_SCALE_CONST_SHR(X,Y,Z,S) MPC_MULTIPLY_EX(X,MAKE_MPC_SAMPLE_EX(Y,Z),(Z)+(S))
+#define MPC_SHR(X,Y) ((X)>>(Y))
+#define MPC_SHL(X,Y) ((X)<<(Y))
+
+#else
+
+//in floating-point mode, decoded samples are in -1...1 range
+
+#define MAKE_MPC_SAMPLE(X) ((MPC_SAMPLE_FORMAT)(X))
+#define MAKE_MPC_SAMPLE_EX(X,Y) ((MPC_SAMPLE_FORMAT)(X))
+
+#define MPC_MULTIPLY_FRACT(X,Y) ((X)*(Y))
+#define MPC_MAKE_FRACT_CONST(X) (X)
+#define MPC_MULTIPLY_FRACT_CONST(X,Y) MPC_MULTPLY_FRACT(X,MPC_MAKE_FRACT_CONST(Y))
+#define MPC_MULTIPLY_FRACT_CONST_SHR(X,Y,Z) MPC_MULTIPLY_FRACT(X,MPC_MAKE_FRACT_CONST( Y ))
+#define MPC_MULTIPLY_FRACT_CONST_FIX(X,Y,Z) MPC_MULTIPLY_FRACT(X,MPC_MAKE_FRACT_CONST( Y ))
+
+#define MPC_MULTIPLY_FLOAT_INT(X,Y) ((X)*(Y))
+#define MPC_MULTIPLY(X,Y) ((X)*(Y))
+#define MPC_MULTIPLY_EX(X,Y,Z) ((X)*(Y))
+#define MPC_SCALE_CONST(X,Y,Z) ((X)*(Y))
+#define MPC_SCALE_CONST_SHL(X,Y,Z,S) ((X)*(Y))
+#define MPC_SCALE_CONST_SHR(X,Y,Z,S) ((X)*(Y))
+#define MPC_SHR(X,Y) (X)
+#define MPC_SHL(X,Y) (X)
+
+#endif
+
+#ifdef __cplusplus
+}
+#endif
+#endif
diff --git a/plugins/musepack/musepack.c b/plugins/musepack/musepack.c
index 7bd363b7..44e364f4 100644
--- a/plugins/musepack/musepack.c
+++ b/plugins/musepack/musepack.c
@@ -20,8 +20,7 @@
#include <assert.h>
#include <limits.h>
#include <unistd.h>
-#include <mpcdec/mpcdec.h>
-#include <mpcdec/math.h>
+#include "mpcdec.h"
#ifdef HAVE_CONFIG_H
# include <config.h>
#endif
@@ -39,24 +38,25 @@ static DB_functions_t *deadbeef;
typedef struct {
DB_fileinfo_t info;
mpc_streaminfo si;
- mpc_decoder mpcdec;
+ mpc_demux *demux;
+// mpc_decoder *mpcdec;
mpc_reader reader;
int currentsample;
int startsample;
int endsample;
mpc_uint32_t vbr_update_acc;
mpc_uint32_t vbr_update_bits;
- float buffer[MPC_DECODER_BUFFER_LENGTH];
+ MPC_SAMPLE_FORMAT buffer[MPC_DECODER_BUFFER_LENGTH];
int remaining;
} musepack_info_t;
-mpc_int32_t musepack_vfs_read (void *t, void *ptr, mpc_int32_t size) {
- return deadbeef->fread(ptr, 1, size, (DB_FILE *)t);
+mpc_int32_t musepack_vfs_read (mpc_reader *r, void *ptr, mpc_int32_t size) {
+ return deadbeef->fread(ptr, 1, size, (DB_FILE *)r->data);
}
/// Seeks to byte position offset.
-mpc_bool_t musepack_vfs_seek (void *t, mpc_int32_t offset) {
- int res = deadbeef->fseek ((DB_FILE *)t, offset, SEEK_SET);
+mpc_bool_t musepack_vfs_seek (mpc_reader *r, mpc_int32_t offset) {
+ int res = deadbeef->fseek ((DB_FILE *)r->data, offset, SEEK_SET);
if (res == 0) {
return 1;
}
@@ -64,17 +64,17 @@ mpc_bool_t musepack_vfs_seek (void *t, mpc_int32_t offset) {
}
/// Returns the current byte offset in the stream.
-mpc_int32_t musepack_vfs_tell (void *t) {
- return deadbeef->ftell ((DB_FILE *)t);
+mpc_int32_t musepack_vfs_tell (mpc_reader *r) {
+ return deadbeef->ftell ((DB_FILE *)r->data);
}
/// Returns the total length of the source stream, in bytes.
-mpc_int32_t musepack_vfs_get_size (void *t) {
- return deadbeef->fgetlength ((DB_FILE *)t);
+mpc_int32_t musepack_vfs_get_size (mpc_reader *r) {
+ return deadbeef->fgetlength ((DB_FILE *)r->data);
}
/// True if the stream is a seekable stream.
-mpc_bool_t musepack_vfs_canseek (void *t) {
+mpc_bool_t musepack_vfs_canseek (mpc_reader *r) {
return 1;
}
@@ -102,20 +102,21 @@ musepack_init (DB_fileinfo_t *_info, DB_playItem_t *it) {
}
info->reader.data = fp;
- mpc_streaminfo_init (&info->si);
- mpc_int32_t err;
- if ((err = mpc_streaminfo_read (&info->si, &info->reader)) != ERROR_CODE_OK) {
- fprintf (stderr, "mpc_streaminfo_read failed (err=%d)\n", err);
- deadbeef->fclose ((DB_FILE *)info->reader.data);
- info->reader.data = NULL;
- return -1;
- }
- mpc_decoder_setup (&info->mpcdec, &info->reader);
- if (!mpc_decoder_initialize (&info->mpcdec, &info->si)) {
- deadbeef->fclose ((DB_FILE *)info->reader.data);
+ info->demux = mpc_demux_init (&info->reader);
+ if (!info->demux) {
+ fprintf (stderr, "mpc: mpc_demux_init failed\n");
+ deadbeef->fclose (fp);
info->reader.data = NULL;
return -1;
}
+ mpc_demux_get_info (info->demux, &info->si);
+
+// info->mpcdec = mpc_decoder_init (&info->si);
+// if (!info->mpcdec) {
+// deadbeef->fclose ((DB_FILE *)info->reader.data);
+// info->reader.data = NULL;
+// return -1;
+// }
info->vbr_update_acc = 0;
info->vbr_update_bits = 0;
info->remaining = 0;
@@ -143,6 +144,14 @@ static void
musepack_free (DB_fileinfo_t *_info) {
musepack_info_t *info = (musepack_info_t *)_info;
if (info) {
+// if (info->mpcdec) {
+// mpc_decoder_exit (info->mpcdec);
+// info->decoder = NULL;
+// }
+ if (info->demux) {
+ mpc_demux_exit (info->demux);
+ info->demux = NULL;
+ }
if (info->reader.data) {
deadbeef->fclose ((DB_FILE *)info->reader.data);
info->reader.data = NULL;
@@ -207,14 +216,17 @@ musepack_read_int16 (DB_fileinfo_t *_info, char *bytes, int size) {
}
if (size > 0 && !info->remaining) {
- mpc_uint32_t decoded = mpc_decoder_decode (&info->mpcdec, info->buffer, &info->vbr_update_acc, &info->vbr_update_bits);
- if (!decoded) {
+ mpc_frame_info frame;
+ frame.buffer = info->buffer;
+ mpc_status err = mpc_demux_decode (info->demux, &frame);
+ if (err != 0 || frame.bits == -1) {
break;
}
- assert (decoded <= MPC_DECODER_BUFFER_LENGTH);
- info->remaining = decoded;
+
+ info->remaining = frame.samples;
}
}
+ info->currentsample += (initsize-size) / sample_size;
return initsize-size;
}
@@ -259,26 +271,30 @@ musepack_read_float32 (DB_fileinfo_t *_info, char *bytes, int size) {
}
if (size > 0 && !info->remaining) {
- mpc_uint32_t decoded = mpc_decoder_decode (&info->mpcdec, info->buffer, &info->vbr_update_acc, &info->vbr_update_bits);
- if (!decoded) {
+ mpc_frame_info frame;
+ frame.buffer = info->buffer;
+ mpc_status err = mpc_demux_decode (info->demux, &frame);
+ if (err != 0 || frame.bits == -1) {
break;
}
- assert (decoded <= MPC_DECODER_BUFFER_LENGTH);
- info->remaining = decoded;
+
+ info->remaining = frame.samples;
}
}
+ info->currentsample += (initsize-size) / (4 * _info->channels);
return initsize-size;
}
static int
musepack_seek_sample (DB_fileinfo_t *_info, int sample) {
musepack_info_t *info = (musepack_info_t *)_info;
- mpc_bool_t res = mpc_decoder_seek_sample (&info->mpcdec, sample);
- if (!res) {
+ mpc_status err = mpc_demux_seek_sample (info->demux, sample + info->startsample);
+ if (err != 0) {
fprintf (stderr, "musepack: seek failed\n");
return -1;
}
- info->currentsample = sample;
- _info->readpos = (sample - info->startsample) / _info->samplerate;
+ info->currentsample = sample + info->startsample;
+ _info->readpos = (float)sample / _info->samplerate;
+ info->remaining = 0;
return 0;
}
@@ -300,22 +316,28 @@ musepack_insert (DB_playItem_t *after, const char *fname) {
DB_FILE *fp = deadbeef->fopen (fname);
if (!fp) {
+ fprintf (stderr, "mpc: insert failed to open %s\n", fname);
return NULL;
}
reader.data = fp;
- mpc_streaminfo si;
- mpc_streaminfo_init (&si);
- mpc_int32_t err;
- if ((err = mpc_streaminfo_read (&si, &reader)) != ERROR_CODE_OK) {
- fprintf (stderr, "mpc_streaminfo_read failed (err=%d)\n", err);
+ mpc_demux *demux = mpc_demux_init (&reader);
+ if (!demux) {
+ fprintf (stderr, "mpc: mpc_demux_init failed\n");
deadbeef->fclose (fp);
return NULL;
}
+ mpc_streaminfo si;
+ //mpc_streaminfo_init (&si);
+ mpc_demux_get_info (demux, &si);
+
int totalsamples = mpc_streaminfo_get_length_samples (&si);
double dur = mpc_streaminfo_get_length (&si);
+ mpc_demux_exit (demux);
+ demux = NULL;
+
DB_playItem_t *it = deadbeef->pl_item_alloc ();
it->decoder_id = deadbeef->plug_get_decoder_id (plugin.plugin.id);
it->fname = strdup (fname);
diff --git a/plugins/musepack/reader.h b/plugins/musepack/reader.h
new file mode 100644
index 00000000..1a93e067
--- /dev/null
+++ b/plugins/musepack/reader.h
@@ -0,0 +1,98 @@
+/*
+ Copyright (c) 2005-2009, The Musepack Development Team
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are
+ met:
+
+ * Redistributions of source code must retain the above copyright
+ notice, this list of conditions and the following disclaimer.
+
+ * Redistributions in binary form must reproduce the above
+ copyright notice, this list of conditions and the following
+ disclaimer in the documentation and/or other materials provided
+ with the distribution.
+
+ * Neither the name of the The Musepack Development Team nor the
+ names of its contributors may be used to endorse or promote
+ products derived from this software without specific prior
+ written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+ "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+ LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+ A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
+ OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
+ SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
+ LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
+ DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
+ THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
+ (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+ OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+*/
+/// \file reader.h
+#ifndef _MPCDEC_READER_H_
+#define _MPCDEC_READER_H_
+#ifdef WIN32
+#pragma once
+#endif
+
+#include <mpc/mpc_types.h>
+#include <stdio.h>
+
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+
+/// \brief Stream reader interface structure.
+///
+/// This is the structure you must supply to the musepack decoding library
+/// to feed it with raw data. Implement the five member functions to provide
+/// a functional reader.
+typedef struct mpc_reader_t mpc_reader;
+struct mpc_reader_t {
+ /// Reads size bytes of data into buffer at ptr.
+ mpc_int32_t (*read)(mpc_reader *p_reader, void *ptr, mpc_int32_t size);
+
+ /// Seeks to byte position offset.
+ mpc_bool_t (*seek)(mpc_reader *p_reader, mpc_int32_t offset);
+
+ /// Returns the current byte offset in the stream.
+ mpc_int32_t (*tell)(mpc_reader *p_reader);
+
+ /// Returns the total length of the source stream, in bytes.
+ mpc_int32_t (*get_size)(mpc_reader *p_reader);
+
+ /// True if the stream is a seekable stream.
+ mpc_bool_t (*canseek)(mpc_reader *p_reader);
+
+ /// Field that can be used to identify a particular instance of
+ /// reader or carry along data associated with that reader.
+ void *data;
+};
+
+/// Initializes reader with default stdio file reader implementation. Use
+/// this if you're just reading from a plain file.
+///
+/// \param r p_reader handle to initialize
+/// \param filename input filename to attach to the reader
+MPC_API mpc_status mpc_reader_init_stdio(mpc_reader *p_reader, const char *filename);
+
+/// Initializes reader with default stdio file reader implementation. Use
+/// this if you prefer to open the file yourself.
+///
+/// \param r p_reader handle to initialize
+/// \param p_file input file handle (already open)
+MPC_API mpc_status mpc_reader_init_stdio_stream(mpc_reader * p_reader, FILE * p_file);
+
+/// Release reader with default stdio file reader implementation.
+///
+/// \param r reader handle to release
+MPC_API void mpc_reader_exit_stdio(mpc_reader *p_reader);
+
+#ifdef __cplusplus
+}
+#endif
+#endif
diff --git a/plugins/musepack/requant.c b/plugins/musepack/requant.c
new file mode 100644
index 00000000..c74d145b
--- /dev/null
+++ b/plugins/musepack/requant.c
@@ -0,0 +1,124 @@
+/*
+ Copyright (c) 2005-2009, The Musepack Development Team
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are
+ met:
+
+ * Redistributions of source code must retain the above copyright
+ notice, this list of conditions and the following disclaimer.
+
+ * Redistributions in binary form must reproduce the above
+ copyright notice, this list of conditions and the following
+ disclaimer in the documentation and/or other materials provided
+ with the distribution.
+
+ * Neither the name of the The Musepack Development Team nor the
+ names of its contributors may be used to endorse or promote
+ products derived from this software without specific prior
+ written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+ "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+ LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+ A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
+ OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
+ SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
+ LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
+ DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
+ THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
+ (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+ OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+*/
+/// \file requant.c
+/// Requantization function implementations.
+/// \todo document me
+#include <mpc/mpcdec.h>
+
+#include "requant.h"
+#include "mpcdec_math.h"
+#include "decoder.h"
+
+/* C O N S T A N T S */
+// Bits per sample for chosen quantizer
+const mpc_uint8_t Res_bit [18] = {
+ 0, 0, 0, 0, 0, 0, 0, 0, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16
+};
+
+// Requantization coefficients
+// 65536/step bzw. 65536/(2*D+1)
+
+#define _(X) MAKE_MPC_SAMPLE_EX(X,14)
+
+const MPC_SAMPLE_FORMAT __Cc [1 + 18] = {
+ _(111.285962475327f), // 32768/2/255*sqrt(3)
+ _(65536.000000000000f), _(21845.333333333332f), _(13107.200000000001f), _(9362.285714285713f),
+ _(7281.777777777777f), _(4369.066666666666f), _(2114.064516129032f), _(1040.253968253968f),
+ _(516.031496062992f), _(257.003921568627f), _(128.250489236790f), _(64.062561094819f),
+ _(32.015632633121f), _(16.003907203907f), _(8.000976681723f), _(4.000244155527f),
+ _(2.000061037018f), _(1.000015259021f)
+};
+
+#undef _
+
+// Requantization offset
+// 2*D+1 = steps of quantizer
+const mpc_int16_t __Dc [1 + 18] = {
+ 2,
+ 0, 1, 2, 3, 4, 7, 15, 31, 63,
+ 127, 255, 511, 1023, 2047, 4095, 8191, 16383, 32767
+};
+
+#ifdef MPC_FIXED_POINT
+static mpc_uint32_t find_shift(double fval)
+{
+ mpc_int64_t val = (mpc_int64_t) fval;
+ mpc_uint32_t ptr = 0;
+ if(val<0)
+ val = -val;
+ while(val)
+ {
+ val >>= 1;
+ ptr++;
+ }
+ return ptr > 31 ? 0 : 31 - ptr;
+}
+#endif
+
+/* F U N C T I O N S */
+
+#define SET_SCF(N,X) d->SCF[N] = MAKE_MPC_SAMPLE_EX(X,d->SCF_shift[N] = (mpc_uint8_t) find_shift(X));
+
+void
+mpc_decoder_scale_output(mpc_decoder *d, double factor)
+{
+ mpc_int32_t n; double f1, f2;
+
+#ifndef MPC_FIXED_POINT
+ factor *= 1.0 / (double) (1<<(MPC_FIXED_POINT_SHIFT-1));
+#else
+ factor *= 1.0 / (double) (1<<(16-MPC_FIXED_POINT_SHIFT));
+#endif
+ f1 = f2 = factor;
+
+ // handles +1.58...-98.41 dB, where's scf[n] / scf[n-1] = 1.20050805774840750476
+
+ SET_SCF(1,factor);
+
+ f1 *= 0.83298066476582673961;
+ f2 *= 1/0.83298066476582673961;
+
+ for ( n = 1; n <= 128; n++ ) {
+ SET_SCF((mpc_uint8_t)(1+n),f1);
+ SET_SCF((mpc_uint8_t)(1-n),f2);
+ f1 *= 0.83298066476582673961;
+ f2 *= 1/0.83298066476582673961;
+ }
+}
+
+void
+mpc_decoder_init_quant(mpc_decoder *d, double scale_factor)
+{
+ mpc_decoder_scale_output(d, scale_factor);
+}
diff --git a/plugins/musepack/requant.h b/plugins/musepack/requant.h
new file mode 100644
index 00000000..eb3124db
--- /dev/null
+++ b/plugins/musepack/requant.h
@@ -0,0 +1,61 @@
+/*
+ Copyright (c) 2005-2009, The Musepack Development Team
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are
+ met:
+
+ * Redistributions of source code must retain the above copyright
+ notice, this list of conditions and the following disclaimer.
+
+ * Redistributions in binary form must reproduce the above
+ copyright notice, this list of conditions and the following
+ disclaimer in the documentation and/or other materials provided
+ with the distribution.
+
+ * Neither the name of the The Musepack Development Team nor the
+ names of its contributors may be used to endorse or promote
+ products derived from this software without specific prior
+ written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+ "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+ LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+ A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
+ OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
+ SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
+ LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
+ DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
+ THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
+ (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+ OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+*/
+/// \file requant.h
+/// Requantization function definitions.
+#ifndef _MPCDEC_REQUANT_H_
+#define _MPCDEC_REQUANT_H_
+#ifdef WIN32
+#pragma once
+#endif
+
+#include <mpc/mpc_types.h>
+
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+
+/* C O N S T A N T S */
+const mpc_uint8_t Res_bit [18]; ///< Bits per sample for chosen quantizer
+const MPC_SAMPLE_FORMAT __Cc [1 + 18]; ///< Requantization coefficients
+const mpc_int16_t __Dc [1 + 18]; ///< Requantization offset
+
+#define Cc (__Cc + 1)
+#define Dc (__Dc + 1)
+
+
+#ifdef __cplusplus
+}
+#endif
+#endif
diff --git a/plugins/musepack/streaminfo.c b/plugins/musepack/streaminfo.c
new file mode 100644
index 00000000..a1fedb9d
--- /dev/null
+++ b/plugins/musepack/streaminfo.c
@@ -0,0 +1,244 @@
+/*
+ Copyright (c) 2005-2009, The Musepack Development Team
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are
+ met:
+
+ * Redistributions of source code must retain the above copyright
+ notice, this list of conditions and the following disclaimer.
+
+ * Redistributions in binary form must reproduce the above
+ copyright notice, this list of conditions and the following
+ disclaimer in the documentation and/or other materials provided
+ with the distribution.
+
+ * Neither the name of the The Musepack Development Team nor the
+ names of its contributors may be used to endorse or promote
+ products derived from this software without specific prior
+ written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+ "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+ LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+ A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
+ OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
+ SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
+ LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
+ DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
+ THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
+ (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+ OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+*/
+/// \file streaminfo.c
+/// Implementation of streaminfo reading functions.
+
+#include <math.h>
+#include <mpc/mpcdec.h>
+#include <mpc/streaminfo.h>
+#include <stdio.h>
+#include "internal.h"
+#include "huffman.h"
+#include "mpc_bits_reader.h"
+
+unsigned long crc32(unsigned char *buf, int len);
+
+static const char na[] = "n.a.";
+static char const * const versionNames[] = {
+ na, "'Unstable/Experimental'", na, na, na, "'quality 0'", "'quality 1'",
+ "'Telephone'", "'Thumb'", "'Radio'", "'Standard'", "'Extreme'", "'Insane'",
+ "'BrainDead'", "'quality 9'", "'quality 10'"
+};
+static const mpc_int32_t samplefreqs[8] = { 44100, 48000, 37800, 32000 };
+
+static const char *
+mpc_get_version_string(float profile) // profile is 0...15, where 7...13 is used
+{
+ return profile >= sizeof versionNames / sizeof *versionNames ? na : versionNames[(int)profile];
+}
+
+static void
+mpc_get_encoder_string(mpc_streaminfo* si)
+{
+ int ver = si->encoder_version;
+ if (si->stream_version >= 8)
+ ver = (si->encoder_version >> 24) * 100 + ((si->encoder_version >> 16) & 0xFF);
+ if (ver <= 116) {
+ if (ver == 0) {
+ sprintf(si->encoder, "Buschmann 1.7.0...9, Klemm 0.90...1.05");
+ } else {
+ switch (ver % 10) {
+ case 0:
+ sprintf(si->encoder, "Release %u.%u", ver / 100,
+ ver / 10 % 10);
+ break;
+ case 2: case 4: case 6: case 8:
+ sprintf(si->encoder, "Beta %u.%02u", ver / 100,
+ ver % 100);
+ break;
+ default:
+ sprintf(si->encoder, "--Alpha-- %u.%02u",
+ ver / 100, ver % 100);
+ break;
+ }
+ }
+ } else {
+ int major = si->encoder_version >> 24;
+ int minor = (si->encoder_version >> 16) & 0xFF;
+ int build = (si->encoder_version >> 8) & 0xFF;
+ char * tmp = "--Stable--";
+
+ if (minor & 1)
+ tmp = "--Unstable--";
+
+ sprintf(si->encoder, "%s %u.%u.%u", tmp, major, minor, build);
+ }
+}
+
+static mpc_status check_streaminfo(mpc_streaminfo * si)
+{
+ if (si->max_band == 0 || si->max_band >= 32
+ || si->channels > 2)
+ return MPC_STATUS_FILE;
+ return MPC_STATUS_OK;
+}
+
+/// Reads streaminfo from SV7 header.
+mpc_status
+streaminfo_read_header_sv7(mpc_streaminfo* si, mpc_bits_reader * r)
+{
+ mpc_uint16_t Estimatedpeak_title = 0;
+ mpc_uint32_t frames, last_frame_samples;
+
+ si->bitrate = 0;
+ frames = (mpc_bits_read(r, 16) << 16) | mpc_bits_read(r, 16);
+ mpc_bits_read(r, 1); // intensity stereo : should be 0
+ si->ms = mpc_bits_read(r, 1);
+ si->max_band = mpc_bits_read(r, 6);
+ si->profile = mpc_bits_read(r, 4);
+ si->profile_name = mpc_get_version_string(si->profile);
+ mpc_bits_read(r, 2); // Link ?
+ si->sample_freq = samplefreqs[mpc_bits_read(r, 2)];
+ Estimatedpeak_title = (mpc_uint16_t) mpc_bits_read(r, 16); // read the ReplayGain data
+ si->gain_title = (mpc_uint16_t) mpc_bits_read(r, 16);
+ si->peak_title = (mpc_uint16_t) mpc_bits_read(r, 16);
+ si->gain_album = (mpc_uint16_t) mpc_bits_read(r, 16);
+ si->peak_album = (mpc_uint16_t) mpc_bits_read(r, 16);
+ si->is_true_gapless = mpc_bits_read(r, 1); // true gapless: used?
+ last_frame_samples = mpc_bits_read(r, 11); // true gapless: valid samples for last frame
+ si->fast_seek = mpc_bits_read(r, 1); // fast seeking
+ mpc_bits_read(r, 19); // unused
+ si->encoder_version = mpc_bits_read(r, 8);
+ si->channels = 2;
+ si->block_pwr = 0;
+
+ // convert gain info
+ if (si->gain_title != 0) {
+ int tmp = (int)((MPC_OLD_GAIN_REF - (mpc_int16_t)si->gain_title / 100.) * 256. + .5);
+ if (tmp >= (1 << 16) || tmp < 0) tmp = 0;
+ si->gain_title = (mpc_int16_t) tmp;
+ }
+
+ if (si->gain_album != 0) {
+ int tmp = (int)((MPC_OLD_GAIN_REF - (mpc_int16_t)si->gain_album / 100.) * 256. + .5);
+ if (tmp >= (1 << 16) || tmp < 0) tmp = 0;
+ si->gain_album = (mpc_int16_t) tmp;
+ }
+
+ if (si->peak_title != 0)
+ si->peak_title = (mpc_uint16_t) (log10(si->peak_title) * 20 * 256 + .5);
+
+ if (si->peak_album != 0)
+ si->peak_album = (mpc_uint16_t) (log10(si->peak_album) * 20 * 256 + .5);
+
+ mpc_get_encoder_string(si);
+
+ if (last_frame_samples == 0) last_frame_samples = MPC_FRAME_LENGTH;
+ si->samples = (mpc_int64_t) frames * MPC_FRAME_LENGTH;
+ if (si->is_true_gapless)
+ si->samples -= (MPC_FRAME_LENGTH - last_frame_samples);
+ else
+ si->samples -= MPC_DECODER_SYNTH_DELAY;
+
+ si->average_bitrate = (si->tag_offset - si->header_position) * 8.0
+ * si->sample_freq / si->samples;
+
+ return check_streaminfo(si);
+}
+
+/// Reads replay gain datas
+void streaminfo_gain(mpc_streaminfo* si, const mpc_bits_reader * r_in)
+{
+ mpc_bits_reader r = *r_in;
+
+ int version = mpc_bits_read(&r, 8); // gain version
+ if (version != 1) // we only know ver 1
+ return;
+ si->gain_title = (mpc_uint16_t) mpc_bits_read(&r, 16);
+ si->peak_title = (mpc_uint16_t) mpc_bits_read(&r, 16);
+ si->gain_album = (mpc_uint16_t) mpc_bits_read(&r, 16);
+ si->peak_album = (mpc_uint16_t) mpc_bits_read(&r, 16);
+}
+
+/// Reads streaminfo from SV8 header.
+mpc_status
+streaminfo_read_header_sv8(mpc_streaminfo* si, const mpc_bits_reader * r_in,
+ mpc_size_t block_size)
+{
+ mpc_uint32_t CRC;
+ mpc_bits_reader r = *r_in;
+
+ CRC = (mpc_bits_read(&r, 16) << 16) | mpc_bits_read(&r, 16);
+ if (CRC != crc32(r.buff + 1 - (r.count >> 3), (int)block_size - 4))
+ return MPC_STATUS_FILE;
+
+ si->stream_version = mpc_bits_read(&r, 8);
+ if (si->stream_version != 8)
+ return MPC_STATUS_INVALIDSV;
+
+ mpc_bits_get_size(&r, &si->samples);
+ mpc_bits_get_size(&r, &si->beg_silence);
+
+ si->is_true_gapless = 1;
+ si->sample_freq = samplefreqs[mpc_bits_read(&r, 3)];
+ si->max_band = mpc_bits_read(&r, 5) + 1;
+ si->channels = mpc_bits_read(&r, 4) + 1;
+ si->ms = mpc_bits_read(&r, 1);
+ si->block_pwr = mpc_bits_read(&r, 3) * 2;
+
+ si->bitrate = 0;
+
+ if ((si->samples - si->beg_silence) != 0)
+ si->average_bitrate = (si->tag_offset - si->header_position) * 8.0
+ * si->sample_freq / (si->samples - si->beg_silence);
+
+ return check_streaminfo(si);
+}
+
+/// Reads encoder informations
+void streaminfo_encoder_info(mpc_streaminfo* si, const mpc_bits_reader * r_in)
+{
+ mpc_bits_reader r = *r_in;
+
+ si->profile = mpc_bits_read(&r, 7) / 8.;
+ si->profile_name = mpc_get_version_string(si->profile);
+ si->pns = mpc_bits_read(&r, 1);
+ si->encoder_version = mpc_bits_read(&r, 8) << 24; // major
+ si->encoder_version |= mpc_bits_read(&r, 8) << 16; // minor
+ si->encoder_version |= mpc_bits_read(&r, 8) << 8; // build
+
+
+ mpc_get_encoder_string(si);
+}
+
+double
+mpc_streaminfo_get_length(mpc_streaminfo * si)
+{
+ return (double) (si->samples - si->beg_silence) / si->sample_freq;
+}
+
+mpc_int64_t mpc_streaminfo_get_length_samples(mpc_streaminfo *si)
+{
+ return si->samples - si->beg_silence;
+}
diff --git a/plugins/musepack/streaminfo.h b/plugins/musepack/streaminfo.h
new file mode 100644
index 00000000..a0a9470b
--- /dev/null
+++ b/plugins/musepack/streaminfo.h
@@ -0,0 +1,109 @@
+/*
+ Copyright (c) 2005-2009, The Musepack Development Team
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are
+ met:
+
+ * Redistributions of source code must retain the above copyright
+ notice, this list of conditions and the following disclaimer.
+
+ * Redistributions in binary form must reproduce the above
+ copyright notice, this list of conditions and the following
+ disclaimer in the documentation and/or other materials provided
+ with the distribution.
+
+ * Neither the name of the The Musepack Development Team nor the
+ names of its contributors may be used to endorse or promote
+ products derived from this software without specific prior
+ written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+ "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+ LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+ A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
+ OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
+ SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
+ LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
+ DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
+ THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
+ (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+ OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+*/
+/// \file streaminfo.h
+#ifndef _MPCDEC_STREAMINFO_H_
+#define _MPCDEC_STREAMINFO_H_
+#ifdef WIN32
+#pragma once
+#endif
+
+#include <mpc/mpc_types.h>
+
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+
+typedef mpc_int32_t mpc_streaminfo_off_t;
+
+/// \brief mpc stream properties structure
+///
+/// Structure containing all the properties of an mpc stream. Populated
+/// by the streaminfo_read function.
+typedef struct mpc_streaminfo {
+ /// @name Core mpc stream properties
+ //@{
+ mpc_uint32_t sample_freq; ///< Sample frequency of stream
+ mpc_uint32_t channels; ///< Number of channels in stream
+ mpc_uint32_t stream_version; ///< Streamversion of stream
+ mpc_uint32_t bitrate; ///< Bitrate of stream file (in bps)
+ double average_bitrate; ///< Average bitrate of stream (in bits/sec)
+ mpc_uint32_t max_band; ///< Maximum band-index used in stream (0...31)
+ mpc_uint32_t ms; ///< Mid/side stereo (0: off, 1: on)
+ mpc_uint32_t fast_seek; ///< True if stream supports fast-seeking (sv7)
+ mpc_uint32_t block_pwr; ///< Number of frames in a block = 2^block_pwr (sv8)
+ //@}
+
+ /// @name Replaygain properties
+ //@{
+ mpc_uint16_t gain_title; ///< Replaygain title value
+ mpc_uint16_t gain_album; ///< Replaygain album value
+ mpc_uint16_t peak_album; ///< Peak album loudness level
+ mpc_uint16_t peak_title; ///< Peak title loudness level
+ //@}
+
+ /// @name True gapless properties
+ //@{
+ mpc_uint32_t is_true_gapless; ///< True gapless? (0: no, 1: yes)
+ mpc_uint64_t samples; ///< Number of samples in the stream
+ mpc_uint64_t beg_silence; ///< Number of samples that must not be played at the beginning of the stream
+ //@}
+
+ /// @name Encoder informations
+ //@{
+ mpc_uint32_t encoder_version; ///< Version of encoder used
+ char encoder[256]; ///< Encoder name
+ mpc_bool_t pns; ///< pns used
+ float profile; ///< Quality profile of stream
+ const char* profile_name; ///< Name of profile used by stream
+ //@}
+
+
+ mpc_streaminfo_off_t header_position; ///< Byte offset of position of header in stream
+ mpc_streaminfo_off_t tag_offset; ///< Offset to file tags
+ mpc_streaminfo_off_t total_file_length; ///< Total length of underlying file
+} mpc_streaminfo;
+
+/// Gets length of stream si, in seconds.
+/// \return length of stream in seconds
+MPC_API double mpc_streaminfo_get_length(mpc_streaminfo *si);
+
+/// Returns length of stream si, in samples.
+/// \return length of stream in samples
+MPC_API mpc_int64_t mpc_streaminfo_get_length_samples(mpc_streaminfo *si);
+
+#ifdef __cplusplus
+}
+#endif
+#endif
diff --git a/plugins/musepack/synth_filter.c b/plugins/musepack/synth_filter.c
new file mode 100644
index 00000000..0d0c3453
--- /dev/null
+++ b/plugins/musepack/synth_filter.c
@@ -0,0 +1,430 @@
+/*
+ Copyright (c) 2005-2009, The Musepack Development Team
+ All rights reserved.
+
+ Redistribution and use in source and binary forms, with or without
+ modification, are permitted provided that the following conditions are
+ met:
+
+ * Redistributions of source code must retain the above copyright
+ notice, this list of conditions and the following disclaimer.
+
+ * Redistributions in binary form must reproduce the above
+ copyright notice, this list of conditions and the following
+ disclaimer in the documentation and/or other materials provided
+ with the distribution.
+
+ * Neither the name of the The Musepack Development Team nor the
+ names of its contributors may be used to endorse or promote
+ products derived from this software without specific prior
+ written permission.
+
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+ "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+ LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+ A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
+ OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
+ SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
+ LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
+ DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
+ THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
+ (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+ OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+*/
+/// \file synth_filter.c
+/// Synthesis functions.
+/// \todo document me
+#include <string.h>
+#include <mpc/mpcdec.h>
+#include "decoder.h"
+#include "mpcdec_math.h"
+
+/* C O N S T A N T S */
+#define MPC_FIXED_POINT_SYNTH_FIX 2
+
+#undef _
+#ifdef MPC_FIXED_POINT
+#define _(value) MPC_MAKE_FRACT_CONST((double)value/(double)(0x40000))
+#else
+#define _(value) MAKE_MPC_SAMPLE((double)value/(double)(0x10000))
+#endif
+
+
+static const MPC_SAMPLE_FORMAT Di_opt [32] [16] = {
+ { _( 0), _( -29), _( 213), _( -459), _( 2037), _(-5153), _( 6574), _(-37489), _(75038), _(37489), _(6574), _( 5153), _(2037), _( 459), _(213), _(29) },
+ { _( -1), _( -31), _( 218), _( -519), _( 2000), _(-5517), _( 5959), _(-39336), _(74992), _(35640), _(7134), _( 4788), _(2063), _( 401), _(208), _(26) },
+ { _( -1), _( -35), _( 222), _( -581), _( 1952), _(-5879), _( 5288), _(-41176), _(74856), _(33791), _(7640), _( 4425), _(2080), _( 347), _(202), _(24) },
+ { _( -1), _( -38), _( 225), _( -645), _( 1893), _(-6237), _( 4561), _(-43006), _(74630), _(31947), _(8092), _( 4063), _(2087), _( 294), _(196), _(21) },
+ { _( -1), _( -41), _( 227), _( -711), _( 1822), _(-6589), _( 3776), _(-44821), _(74313), _(30112), _(8492), _( 3705), _(2085), _( 244), _(190), _(19) },
+ { _( -1), _( -45), _( 228), _( -779), _( 1739), _(-6935), _( 2935), _(-46617), _(73908), _(28289), _(8840), _( 3351), _(2075), _( 197), _(183), _(17) },
+ { _( -1), _( -49), _( 228), _( -848), _( 1644), _(-7271), _( 2037), _(-48390), _(73415), _(26482), _(9139), _( 3004), _(2057), _( 153), _(176), _(16) },
+ { _( -2), _( -53), _( 227), _( -919), _( 1535), _(-7597), _( 1082), _(-50137), _(72835), _(24694), _(9389), _( 2663), _(2032), _( 111), _(169), _(14) },
+ { _( -2), _( -58), _( 224), _( -991), _( 1414), _(-7910), _( 70), _(-51853), _(72169), _(22929), _(9592), _( 2330), _(2001), _( 72), _(161), _(13) },
+ { _( -2), _( -63), _( 221), _(-1064), _( 1280), _(-8209), _( -998), _(-53534), _(71420), _(21189), _(9750), _( 2006), _(1962), _( 36), _(154), _(11) },
+ { _( -2), _( -68), _( 215), _(-1137), _( 1131), _(-8491), _( -2122), _(-55178), _(70590), _(19478), _(9863), _( 1692), _(1919), _( 2), _(147), _(10) },
+ { _( -3), _( -73), _( 208), _(-1210), _( 970), _(-8755), _( -3300), _(-56778), _(69679), _(17799), _(9935), _( 1388), _(1870), _( -29), _(139), _( 9) },
+ { _( -3), _( -79), _( 200), _(-1283), _( 794), _(-8998), _( -4533), _(-58333), _(68692), _(16155), _(9966), _( 1095), _(1817), _( -57), _(132), _( 8) },
+ { _( -4), _( -85), _( 189), _(-1356), _( 605), _(-9219), _( -5818), _(-59838), _(67629), _(14548), _(9959), _( 814), _(1759), _( -83), _(125), _( 7) },
+ { _( -4), _( -91), _( 177), _(-1428), _( 402), _(-9416), _( -7154), _(-61289), _(66494), _(12980), _(9916), _( 545), _(1698), _(-106), _(117), _( 7) },
+ { _( -5), _( -97), _( 163), _(-1498), _( 185), _(-9585), _( -8540), _(-62684), _(65290), _(11455), _(9838), _( 288), _(1634), _(-127), _(111), _( 6) },
+ { _( -5), _(-104), _( 146), _(-1567), _( -45), _(-9727), _( -9975), _(-64019), _(64019), _( 9975), _(9727), _( 45), _(1567), _(-146), _(104), _( 5) },
+ { _( -6), _(-111), _( 127), _(-1634), _( -288), _(-9838), _(-11455), _(-65290), _(62684), _( 8540), _(9585), _( -185), _(1498), _(-163), _( 97), _( 5) },
+ { _( -7), _(-117), _( 106), _(-1698), _( -545), _(-9916), _(-12980), _(-66494), _(61289), _( 7154), _(9416), _( -402), _(1428), _(-177), _( 91), _( 4) },
+ { _( -7), _(-125), _( 83), _(-1759), _( -814), _(-9959), _(-14548), _(-67629), _(59838), _( 5818), _(9219), _( -605), _(1356), _(-189), _( 85), _( 4) },
+ { _( -8), _(-132), _( 57), _(-1817), _(-1095), _(-9966), _(-16155), _(-68692), _(58333), _( 4533), _(8998), _( -794), _(1283), _(-200), _( 79), _( 3) },
+ { _( -9), _(-139), _( 29), _(-1870), _(-1388), _(-9935), _(-17799), _(-69679), _(56778), _( 3300), _(8755), _( -970), _(1210), _(-208), _( 73), _( 3) },
+ { _(-10), _(-147), _( -2), _(-1919), _(-1692), _(-9863), _(-19478), _(-70590), _(55178), _( 2122), _(8491), _(-1131), _(1137), _(-215), _( 68), _( 2) },
+ { _(-11), _(-154), _( -36), _(-1962), _(-2006), _(-9750), _(-21189), _(-71420), _(53534), _( 998), _(8209), _(-1280), _(1064), _(-221), _( 63), _( 2) },
+ { _(-13), _(-161), _( -72), _(-2001), _(-2330), _(-9592), _(-22929), _(-72169), _(51853), _( -70), _(7910), _(-1414), _( 991), _(-224), _( 58), _( 2) },
+ { _(-14), _(-169), _(-111), _(-2032), _(-2663), _(-9389), _(-24694), _(-72835), _(50137), _(-1082), _(7597), _(-1535), _( 919), _(-227), _( 53), _( 2) },
+ { _(-16), _(-176), _(-153), _(-2057), _(-3004), _(-9139), _(-26482), _(-73415), _(48390), _(-2037), _(7271), _(-1644), _( 848), _(-228), _( 49), _( 1) },
+ { _(-17), _(-183), _(-197), _(-2075), _(-3351), _(-8840), _(-28289), _(-73908), _(46617), _(-2935), _(6935), _(-1739), _( 779), _(-228), _( 45), _( 1) },
+ { _(-19), _(-190), _(-244), _(-2085), _(-3705), _(-8492), _(-30112), _(-74313), _(44821), _(-3776), _(6589), _(-1822), _( 711), _(-227), _( 41), _( 1) },
+ { _(-21), _(-196), _(-294), _(-2087), _(-4063), _(-8092), _(-31947), _(-74630), _(43006), _(-4561), _(6237), _(-1893), _( 645), _(-225), _( 38), _( 1) },
+ { _(-24), _(-202), _(-347), _(-2080), _(-4425), _(-7640), _(-33791), _(-74856), _(41176), _(-5288), _(5879), _(-1952), _( 581), _(-222), _( 35), _( 1) },
+ { _(-26), _(-208), _(-401), _(-2063), _(-4788), _(-7134), _(-35640), _(-74992), _(39336), _(-5959), _(5517), _(-2000), _( 519), _(-218), _( 31), _( 1) }
+};
+
+#undef _
+
+static void
+mpc_compute_new_V(const MPC_SAMPLE_FORMAT* p_sample, MPC_SAMPLE_FORMAT* pV)
+{
+ // Calculating new V-buffer values for left channel
+ // calculate new V-values (ISO-11172-3, p. 39)
+ // based upon fast-MDCT algorithm by Byeong Gi Lee
+ MPC_SAMPLE_FORMAT A00, A01, A02, A03, A04, A05, A06, A07, A08, A09, A10, A11, A12, A13, A14, A15;
+ MPC_SAMPLE_FORMAT B00, B01, B02, B03, B04, B05, B06, B07, B08, B09, B10, B11, B12, B13, B14, B15;
+ MPC_SAMPLE_FORMAT tmp;
+
+ A00 = p_sample[ 0] + p_sample[31];
+ A01 = p_sample[ 1] + p_sample[30];
+ A02 = p_sample[ 2] + p_sample[29];
+ A03 = p_sample[ 3] + p_sample[28];
+ A04 = p_sample[ 4] + p_sample[27];
+ A05 = p_sample[ 5] + p_sample[26];
+ A06 = p_sample[ 6] + p_sample[25];
+ A07 = p_sample[ 7] + p_sample[24];
+ A08 = p_sample[ 8] + p_sample[23];
+ A09 = p_sample[ 9] + p_sample[22];
+ A10 = p_sample[10] + p_sample[21];
+ A11 = p_sample[11] + p_sample[20];
+ A12 = p_sample[12] + p_sample[19];
+ A13 = p_sample[13] + p_sample[18];
+ A14 = p_sample[14] + p_sample[17];
+ A15 = p_sample[15] + p_sample[16];
+
+ B00 = A00 + A15;
+ B01 = A01 + A14;
+ B02 = A02 + A13;
+ B03 = A03 + A12;
+ B04 = A04 + A11;
+ B05 = A05 + A10;
+ B06 = A06 + A09;
+ B07 = A07 + A08;;
+ B08 = MPC_SCALE_CONST((A00 - A15) , 0.5024192929f , 31);
+ B09 = MPC_SCALE_CONST((A01 - A14) , 0.5224986076f , 31);
+ B10 = MPC_SCALE_CONST((A02 - A13) , 0.5669440627f , 31);
+ B11 = MPC_SCALE_CONST((A03 - A12) , 0.6468217969f , 31);
+ B12 = MPC_SCALE_CONST((A04 - A11) , 0.7881546021f , 31);
+ B13 = MPC_SCALE_CONST((A05 - A10) , 1.0606776476f , 30);
+ B14 = MPC_SCALE_CONST((A06 - A09) , 1.7224471569f , 30);
+ B15 = MPC_SCALE_CONST((A07 - A08) , 5.1011486053f , 28);
+
+ A00 = B00 + B07;
+ A01 = B01 + B06;
+ A02 = B02 + B05;
+ A03 = B03 + B04;
+ A04 = MPC_SCALE_CONST((B00 - B07) , 0.5097956061f , 31);
+ A05 = MPC_SCALE_CONST((B01 - B06) , 0.6013448834f , 31);
+ A06 = MPC_SCALE_CONST((B02 - B05) , 0.8999761939f , 31);
+ A07 = MPC_SCALE_CONST((B03 - B04) , 2.5629155636f , 29);
+ A08 = B08 + B15;
+ A09 = B09 + B14;
+ A10 = B10 + B13;
+ A11 = B11 + B12;
+ A12 = MPC_SCALE_CONST((B08 - B15) , 0.5097956061f , 31);
+ A13 = MPC_SCALE_CONST((B09 - B14) , 0.6013448834f , 31);
+ A14 = MPC_SCALE_CONST((B10 - B13) , 0.8999761939f , 31);
+ A15 = MPC_SCALE_CONST((B11 - B12) , 2.5629155636f , 29);
+
+ B00 = A00 + A03;
+ B01 = A01 + A02;
+ B02 = MPC_MULTIPLY_FRACT_CONST_FIX((A00 - A03) , 0.5411961079f , 1);
+ B03 = MPC_MULTIPLY_FRACT_CONST_FIX((A01 - A02) , 1.3065630198f , 2);
+ B04 = A04 + A07;
+ B05 = A05 + A06;
+ B06 = MPC_MULTIPLY_FRACT_CONST_FIX((A04 - A07) , 0.5411961079f , 1);
+ B07 = MPC_MULTIPLY_FRACT_CONST_FIX((A05 - A06) , 1.3065630198f , 2);
+ B08 = A08 + A11;
+ B09 = A09 + A10;
+ B10 = MPC_MULTIPLY_FRACT_CONST_FIX((A08 - A11) , 0.5411961079f , 1);
+ B11 = MPC_MULTIPLY_FRACT_CONST_FIX((A09 - A10) , 1.3065630198f , 2);
+ B12 = A12 + A15;
+ B13 = A13 + A14;
+ B14 = MPC_MULTIPLY_FRACT_CONST_FIX((A12 - A15) , 0.5411961079f , 1);
+ B15 = MPC_MULTIPLY_FRACT_CONST_FIX((A13 - A14) , 1.3065630198f , 2);
+
+ A00 = B00 + B01;
+ A01 = MPC_MULTIPLY_FRACT_CONST_FIX((B00 - B01) , 0.7071067691f , 1);
+ A02 = B02 + B03;
+ A03 = MPC_MULTIPLY_FRACT_CONST_FIX((B02 - B03) , 0.7071067691f , 1);
+ A04 = B04 + B05;
+ A05 = MPC_MULTIPLY_FRACT_CONST_FIX((B04 - B05) , 0.7071067691f , 1);
+ A06 = B06 + B07;
+ A07 = MPC_MULTIPLY_FRACT_CONST_FIX((B06 - B07) , 0.7071067691f , 1);
+ A08 = B08 + B09;
+ A09 = MPC_MULTIPLY_FRACT_CONST_FIX((B08 - B09) , 0.7071067691f , 1);
+ A10 = B10 + B11;
+ A11 = MPC_MULTIPLY_FRACT_CONST_FIX((B10 - B11) , 0.7071067691f , 1);
+ A12 = B12 + B13;
+ A13 = MPC_MULTIPLY_FRACT_CONST_FIX((B12 - B13) , 0.7071067691f , 1);
+ A14 = B14 + B15;
+ A15 = MPC_MULTIPLY_FRACT_CONST_FIX((B14 - B15) , 0.7071067691f , 1);
+
+ pV[48] = -A00;
+ pV[ 0] = A01;
+ pV[40] = -A02 - (pV[ 8] = A03);
+ pV[36] = -((pV[ 4] = A05 + (pV[12] = A07)) + A06);
+ pV[44] = - A04 - A06 - A07;
+ pV[ 6] = (pV[10] = A11 + (pV[14] = A15)) + A13;
+ pV[38] = (pV[34] = -(pV[ 2] = A09 + A13 + A15) - A14) + A09 - A10 - A11;
+ pV[46] = (tmp = -(A12 + A14 + A15)) - A08;
+ pV[42] = tmp - A10 - A11;
+
+ A00 = MPC_MULTIPLY_FRACT_CONST_SHR((p_sample[ 0] - p_sample[31]) , 0.5006030202f , MPC_FIXED_POINT_SYNTH_FIX);
+ A01 = MPC_MULTIPLY_FRACT_CONST_SHR((p_sample[ 1] - p_sample[30]) , 0.5054709315f , MPC_FIXED_POINT_SYNTH_FIX);
+ A02 = MPC_MULTIPLY_FRACT_CONST_SHR((p_sample[ 2] - p_sample[29]) , 0.5154473186f , MPC_FIXED_POINT_SYNTH_FIX);
+ A03 = MPC_MULTIPLY_FRACT_CONST_SHR((p_sample[ 3] - p_sample[28]) , 0.5310425758f , MPC_FIXED_POINT_SYNTH_FIX);
+ A04 = MPC_MULTIPLY_FRACT_CONST_SHR((p_sample[ 4] - p_sample[27]) , 0.5531039238f , MPC_FIXED_POINT_SYNTH_FIX);
+ A05 = MPC_MULTIPLY_FRACT_CONST_SHR((p_sample[ 5] - p_sample[26]) , 0.5829349756f , MPC_FIXED_POINT_SYNTH_FIX);
+ A06 = MPC_MULTIPLY_FRACT_CONST_SHR((p_sample[ 6] - p_sample[25]) , 0.6225041151f , MPC_FIXED_POINT_SYNTH_FIX);
+ A07 = MPC_MULTIPLY_FRACT_CONST_SHR((p_sample[ 7] - p_sample[24]) , 0.6748083234f , MPC_FIXED_POINT_SYNTH_FIX);
+ A08 = MPC_MULTIPLY_FRACT_CONST_SHR((p_sample[ 8] - p_sample[23]) , 0.7445362806f , MPC_FIXED_POINT_SYNTH_FIX);
+ A09 = MPC_MULTIPLY_FRACT_CONST_SHR((p_sample[ 9] - p_sample[22]) , 0.8393496275f , MPC_FIXED_POINT_SYNTH_FIX);
+ A10 = MPC_MULTIPLY_FRACT_CONST_SHR((p_sample[10] - p_sample[21]) , 0.9725682139f , MPC_FIXED_POINT_SYNTH_FIX);
+#if MPC_FIXED_POINT_SYNTH_FIX>=2
+ A11 = MPC_MULTIPLY_FRACT_CONST_SHR((p_sample[11] - p_sample[20]) , 1.1694399118f , MPC_FIXED_POINT_SYNTH_FIX);
+ A12 = MPC_MULTIPLY_FRACT_CONST_SHR((p_sample[12] - p_sample[19]) , 1.4841645956f , MPC_FIXED_POINT_SYNTH_FIX);
+#else
+ A11 = MPC_SCALE_CONST_SHR ((p_sample[11] - p_sample[20]) , 1.1694399118f , 30, MPC_FIXED_POINT_SYNTH_FIX);
+ A12 = MPC_SCALE_CONST_SHR ((p_sample[12] - p_sample[19]) , 1.4841645956f , 30, MPC_FIXED_POINT_SYNTH_FIX);
+#endif
+ A13 = MPC_SCALE_CONST_SHR ((p_sample[13] - p_sample[18]) , 2.0577809811f , 29, MPC_FIXED_POINT_SYNTH_FIX);
+ A14 = MPC_SCALE_CONST_SHR ((p_sample[14] - p_sample[17]) , 3.4076085091f , 29, MPC_FIXED_POINT_SYNTH_FIX);
+ A15 = MPC_SCALE_CONST_SHR ((p_sample[15] - p_sample[16]) , 10.1900081635f, 27 ,MPC_FIXED_POINT_SYNTH_FIX);
+
+ B00 = A00 + A15;
+ B01 = A01 + A14;
+ B02 = A02 + A13;
+ B03 = A03 + A12;
+ B04 = A04 + A11;
+ B05 = A05 + A10;
+ B06 = A06 + A09;
+ B07 = A07 + A08;
+ B08 = MPC_SCALE_CONST((A00 - A15) , 0.5024192929f , 31);
+ B09 = MPC_SCALE_CONST((A01 - A14) , 0.5224986076f , 31);
+ B10 = MPC_SCALE_CONST((A02 - A13) , 0.5669440627f , 31);
+ B11 = MPC_SCALE_CONST((A03 - A12) , 0.6468217969f , 31);
+ B12 = MPC_SCALE_CONST((A04 - A11) , 0.7881546021f , 31);
+ B13 = MPC_SCALE_CONST((A05 - A10) , 1.0606776476f , 30);
+ B14 = MPC_SCALE_CONST((A06 - A09) , 1.7224471569f , 30);
+ B15 = MPC_SCALE_CONST((A07 - A08) , 5.1011486053f , 28);
+
+ A00 = B00 + B07;
+ A01 = B01 + B06;
+ A02 = B02 + B05;
+ A03 = B03 + B04;
+ A04 = MPC_SCALE_CONST((B00 - B07) , 0.5097956061f , 31);
+ A05 = MPC_SCALE_CONST((B01 - B06) , 0.6013448834f , 31);
+ A06 = MPC_SCALE_CONST((B02 - B05) , 0.8999761939f , 31);
+ A07 = MPC_SCALE_CONST((B03 - B04) , 2.5629155636f , 29);
+ A08 = B08 + B15;
+ A09 = B09 + B14;
+ A10 = B10 + B13;
+ A11 = B11 + B12;
+ A12 = MPC_SCALE_CONST((B08 - B15) , 0.5097956061f , 31);
+ A13 = MPC_SCALE_CONST((B09 - B14) , 0.6013448834f , 31);
+ A14 = MPC_SCALE_CONST((B10 - B13) , 0.8999761939f , 31);
+ A15 = MPC_SCALE_CONST((B11 - B12) , 2.5629155636f , 29);
+
+ B00 = A00 + A03;
+ B01 = A01 + A02;
+ B02 = MPC_SCALE_CONST((A00 - A03) , 0.5411961079f , 31);
+ B03 = MPC_SCALE_CONST((A01 - A02) , 1.3065630198f , 30);
+ B04 = A04 + A07;
+ B05 = A05 + A06;
+ B06 = MPC_SCALE_CONST((A04 - A07) , 0.5411961079f , 31);
+ B07 = MPC_SCALE_CONST((A05 - A06) , 1.3065630198f , 30);
+ B08 = A08 + A11;
+ B09 = A09 + A10;
+ B10 = MPC_SCALE_CONST((A08 - A11) , 0.5411961079f , 31);
+ B11 = MPC_SCALE_CONST((A09 - A10) , 1.3065630198f , 30);
+ B12 = A12 + A15;
+ B13 = A13 + A14;
+ B14 = MPC_SCALE_CONST((A12 - A15) , 0.5411961079f , 31);
+ B15 = MPC_SCALE_CONST((A13 - A14) , 1.3065630198f , 30);
+
+ A00 = MPC_SHL(B00 + B01, MPC_FIXED_POINT_SYNTH_FIX);
+ A01 = MPC_SCALE_CONST_SHL((B00 - B01) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
+ A02 = MPC_SHL(B02 + B03, MPC_FIXED_POINT_SYNTH_FIX);
+ A03 = MPC_SCALE_CONST_SHL((B02 - B03) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
+ A04 = MPC_SHL(B04 + B05, MPC_FIXED_POINT_SYNTH_FIX);
+ A05 = MPC_SCALE_CONST_SHL((B04 - B05) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
+ A06 = MPC_SHL(B06 + B07, MPC_FIXED_POINT_SYNTH_FIX);
+ A07 = MPC_SCALE_CONST_SHL((B06 - B07) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
+ A08 = MPC_SHL(B08 + B09, MPC_FIXED_POINT_SYNTH_FIX);
+ A09 = MPC_SCALE_CONST_SHL((B08 - B09) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
+ A10 = MPC_SHL(B10 + B11, MPC_FIXED_POINT_SYNTH_FIX);
+ A11 = MPC_SCALE_CONST_SHL((B10 - B11) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
+ A12 = MPC_SHL(B12 + B13, MPC_FIXED_POINT_SYNTH_FIX);
+ A13 = MPC_SCALE_CONST_SHL((B12 - B13) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
+ A14 = MPC_SHL(B14 + B15, MPC_FIXED_POINT_SYNTH_FIX);
+ A15 = MPC_SCALE_CONST_SHL((B14 - B15) , 0.7071067691f , 31, MPC_FIXED_POINT_SYNTH_FIX);
+
+ // mehrfach verwendete Ausdrcke: A04+A06+A07, A09+A13+A15
+ pV[ 5] = (pV[11] = (pV[13] = A07 + (pV[15] = A15)) + A11) + A05 + A13;
+ pV[ 7] = (pV[ 9] = A03 + A11 + A15) + A13;
+ pV[33] = -(pV[ 1] = A01 + A09 + A13 + A15) - A14;
+ pV[35] = -(pV[ 3] = A05 + A07 + A09 + A13 + A15) - A06 - A14;
+ pV[37] = (tmp = -(A10 + A11 + A13 + A14 + A15)) - A05 - A06 - A07;
+ pV[39] = tmp - A02 - A03; // abh�gig vom Befehl drber
+ pV[41] = (tmp += A13 - A12) - A02 - A03; // abh�gig vom Befehl 2 drber
+ pV[43] = tmp - A04 - A06 - A07; // abh�gig von Befehlen 1 und 3 drber
+ pV[47] = (tmp = -(A08 + A12 + A14 + A15)) - A00;
+ pV[45] = tmp - A04 - A06 - A07; // abh�gig vom Befehl drber
+
+ pV[32] = -pV[ 0];
+ pV[31] = -pV[ 1];
+ pV[30] = -pV[ 2];
+ pV[29] = -pV[ 3];
+ pV[28] = -pV[ 4];
+ pV[27] = -pV[ 5];
+ pV[26] = -pV[ 6];
+ pV[25] = -pV[ 7];
+ pV[24] = -pV[ 8];
+ pV[23] = -pV[ 9];
+ pV[22] = -pV[10];
+ pV[21] = -pV[11];
+ pV[20] = -pV[12];
+ pV[19] = -pV[13];
+ pV[18] = -pV[14];
+ pV[17] = -pV[15];
+
+ pV[63] = pV[33];
+ pV[62] = pV[34];
+ pV[61] = pV[35];
+ pV[60] = pV[36];
+ pV[59] = pV[37];
+ pV[58] = pV[38];
+ pV[57] = pV[39];
+ pV[56] = pV[40];
+ pV[55] = pV[41];
+ pV[54] = pV[42];
+ pV[53] = pV[43];
+ pV[52] = pV[44];
+ pV[51] = pV[45];
+ pV[50] = pV[46];
+ pV[49] = pV[47];
+}
+
+static void
+mpc_synthese_filter_float_internal(MPC_SAMPLE_FORMAT* p_out, MPC_SAMPLE_FORMAT* pV, const MPC_SAMPLE_FORMAT* pY, mpc_int_t channels)
+{
+ mpc_uint32_t n;
+ for ( n = 0; n < 36; n++, pY += 32 )
+ {
+ MPC_SAMPLE_FORMAT* pData = p_out;
+ const MPC_SAMPLE_FORMAT* pD = (const MPC_SAMPLE_FORMAT*) &Di_opt;
+ mpc_int32_t k;
+ pV -= 64;
+ mpc_compute_new_V( pY, pV );
+ for ( k = 0; k < 32; k++, pD += 16, pV++ )
+ {
+ *pData = MPC_SHL(
+ MPC_MULTIPLY_FRACT(pV[ 0], pD[ 0]) + MPC_MULTIPLY_FRACT(pV[ 96], pD[ 1]) + MPC_MULTIPLY_FRACT(pV[128], pD[ 2]) + MPC_MULTIPLY_FRACT(pV[224], pD[ 3])
+ + MPC_MULTIPLY_FRACT(pV[256], pD[ 4]) + MPC_MULTIPLY_FRACT(pV[352], pD[ 5]) + MPC_MULTIPLY_FRACT(pV[384], pD[ 6]) + MPC_MULTIPLY_FRACT(pV[480], pD[ 7])
+ + MPC_MULTIPLY_FRACT(pV[512], pD[ 8]) + MPC_MULTIPLY_FRACT(pV[608], pD[ 9]) + MPC_MULTIPLY_FRACT(pV[640], pD[10]) + MPC_MULTIPLY_FRACT(pV[736], pD[11])
+ + MPC_MULTIPLY_FRACT(pV[768], pD[12]) + MPC_MULTIPLY_FRACT(pV[864], pD[13]) + MPC_MULTIPLY_FRACT(pV[896], pD[14]) + MPC_MULTIPLY_FRACT(pV[992], pD[15])
+ , 2);
+ pData += channels;
+ }
+ pV -= 32; //bleh
+ p_out += 32 * channels;
+ }
+}
+
+void
+mpc_decoder_synthese_filter_float(mpc_decoder* p_dec, MPC_SAMPLE_FORMAT* p_out, mpc_int_t channels)
+{
+ /********* left channel ********/
+ memmove(&p_dec->V_L[MPC_V_MEM], p_dec->V_L, 960 * sizeof *p_dec->V_L);
+ mpc_synthese_filter_float_internal(p_out, &p_dec->V_L[MPC_V_MEM], p_dec->Y_L[0], channels);
+
+ /******** right channel ********/
+ if (channels > 1) {
+ memmove(&p_dec->V_R[MPC_V_MEM], p_dec->V_R, 960 * sizeof *p_dec->V_R);
+ mpc_synthese_filter_float_internal(p_out + 1, &p_dec->V_R[MPC_V_MEM], p_dec->Y_R[0], channels);
+ }
+}
+
+/*******************************************/
+/* */
+/* dithered synthesis */
+/* */
+/*******************************************/
+
+static const mpc_uint8_t Parity [256] = { // parity
+ 0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,0,1,0,0,1,0,1,1,0,0,1,1,0,1,0,0,1,
+ 1,0,0,1,0,1,1,0,0,1,1,0,1,0,0,1,0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,0,
+ 1,0,0,1,0,1,1,0,0,1,1,0,1,0,0,1,0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,0,
+ 0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,0,1,0,0,1,0,1,1,0,0,1,1,0,1,0,0,1,
+ 1,0,0,1,0,1,1,0,0,1,1,0,1,0,0,1,0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,0,
+ 0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,0,1,0,0,1,0,1,1,0,0,1,1,0,1,0,0,1,
+ 0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,0,1,0,0,1,0,1,1,0,0,1,1,0,1,0,0,1,
+ 1,0,0,1,0,1,1,0,0,1,1,0,1,0,0,1,0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,0
+};
+
+/*
+ * This is a simple random number generator with good quality for audio purposes.
+ * It consists of two polycounters with opposite rotation direction and different
+ * periods. The periods are coprime, so the total period is the product of both.
+ *
+ * -------------------------------------------------------------------------------------------------
+ * +-> |31:30:29:28:27:26:25:24:23:22:21:20:19:18:17:16:15:14:13:12:11:10: 9: 8: 7: 6: 5: 4: 3: 2: 1: 0|
+ * | -------------------------------------------------------------------------------------------------
+ * | | | | | | |
+ * | +--+--+--+-XOR-+--------+
+ * | |
+ * +--------------------------------------------------------------------------------------+
+ *
+ * -------------------------------------------------------------------------------------------------
+ * |31:30:29:28:27:26:25:24:23:22:21:20:19:18:17:16:15:14:13:12:11:10: 9: 8: 7: 6: 5: 4: 3: 2: 1: 0| <-+
+ * ------------------------------------------------------------------------------------------------- |
+ * | | | | |
+ * +--+----XOR----+--+ |
+ * | |
+ * +----------------------------------------------------------------------------------------+
+ *
+ *
+ * The first has an period of 3*5*17*257*65537, the second of 7*47*73*178481,
+ * which gives a period of 18.410.713.077.675.721.215. The result is the
+ * XORed values of both generators.
+ */
+mpc_uint32_t
+mpc_random_int(mpc_decoder* p_dec)
+{
+#if 1
+ mpc_uint32_t t1, t2, t3, t4;
+
+ t3 = t1 = p_dec->__r1; t4 = t2 = p_dec->__r2; // Parity calculation is done via table lookup, this is also available
+ t1 &= 0xF5; t2 >>= 25; // on CPUs without parity, can be implemented in C and avoid unpredictable
+ t1 = Parity[t1]; t2 &= 0x63; // jumps and slow rotate through the carry flag operations.
+ t1 <<= 31; t2 = Parity[t2];
+
+ return (p_dec->__r1 = (t3 >> 1) | t1 ) ^ (p_dec->__r2 = (t4 + t4) | t2 );
+#else
+ return (p_dec->__r1 = (p_dec->__r1 >> 1) | ((mpc_uint32_t) Parity[ p_dec->__r1 & 0xF5] << 31))
+ ^ (p_dec->__r2 = (p_dec->__r2 << 1) | (mpc_uint32_t) Parity[(p_dec->__r2 >> 25) & 0x63]);
+#endif
+}
generated by cgit v1.2.3 (git 2.39.1) at 2024-06-26 03:41:04 +0000