kenmirrors/asterisk
1
0
Fork 0
mirror of https://github.com/asterisk/asterisk.git synced 2026-05-04 20:33:53 +00:00
Code Issues Packages Projects Releases Wiki Activity
Files
c680adba74e2d47a9a9fc2fe4b4100cac66092b7
asterisk/codecs/mp3/src/iupini.c
Mark Spencer 0d803132c1 Version 0.1.7 from FTP 
git-svn-id: https://origsvn.digium.com/svn/asterisk/trunk@243 65c4cc65-6c06-0410-ace0-fbb531ad65f3
2001-03-20 20:11:20 +00:00

391 lines
12 KiB
C
Executable File
Raw Blame History

/*____________________________________________________________________________
FreeAmp - The Free MP3 Player
MP3 Decoder originally Copyright (C) 1995-1997 Xing Technology
Corp. http://www.xingtech.com
Portions Copyright (C) 1998-1999 EMusic.com
This program is free software; you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation; either version 2 of the License, or
(at your option) any later version.
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License
along with this program; if not, write to the Free Software
Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
$Id$
____________________________________________________________________________*/
/*=========================================================
initialization for iup.c - include to iup.c
mpeg audio decoder portable "c" integer
mods 11/15/95 for Layer I
mods 1/8/97 warnings
=========================================================*/
#include <limits.h>
#include <string.h>
/* Read only */
static long steps[18] =
{
0, 3, 5, 7, 9, 15, 31, 63, 127,
255, 511, 1023, 2047, 4095, 8191, 16383, 32767, 65535};
/* Read only */
static int stepbits[18] =
{
0, 2, 3, 3, 4, 4, 5, 6, 7,
8, 9, 10, 11, 12, 13, 14, 15, 16};
/* ABCD_INDEX = lookqt[mode][sr_index][br_index] */
/* -1 = invalid */
/* Read only */
static signed char lookqt[4][3][16] =
{
{{1, -1, -1, -1, 2, -1, 2, 0, 0, 0, 1, 1, 1, 1, 1, -1}, /* 44ks stereo */
{0, -1, -1, -1, 2, -1, 2, 0, 0, 0, 0, 0, 0, 0, 0, -1}, /* 48ks */
{1, -1, -1, -1, 3, -1, 3, 0, 0, 0, 1, 1, 1, 1, 1, -1}}, /* 32ks */
{{1, -1, -1, -1, 2, -1, 2, 0, 0, 0, 1, 1, 1, 1, 1, -1}, /* 44ks joint stereo */
{0, -1, -1, -1, 2, -1, 2, 0, 0, 0, 0, 0, 0, 0, 0, -1}, /* 48ks */
{1, -1, -1, -1, 3, -1, 3, 0, 0, 0, 1, 1, 1, 1, 1, -1}}, /* 32ks */
{{1, -1, -1, -1, 2, -1, 2, 0, 0, 0, 1, 1, 1, 1, 1, -1}, /* 44ks dual chan */
{0, -1, -1, -1, 2, -1, 2, 0, 0, 0, 0, 0, 0, 0, 0, -1}, /* 48ks */
{1, -1, -1, -1, 3, -1, 3, 0, 0, 0, 1, 1, 1, 1, 1, -1}}, /* 32ks */
{{1, 2, 2, 0, 0, 0, 1, 1, 1, 1, 1, -1, -1, -1, -1, -1}, /* 44ks single chan */
{0, 2, 2, 0, 0, 0, 0, 0, 0, 0, 0, -1, -1, -1, -1, -1}, /* 48ks */
{1, 3, 3, 0, 0, 0, 1, 1, 1, 1, 1, -1, -1, -1, -1, -1}}, /* 32ks */
};
/* Read only */
static long sr_table[8] =
{22050L, 24000L, 16000L, 1L,
44100L, 48000L, 32000L, 1L};
/* bit allocation table look up */
/* table per mpeg spec tables 3b2a/b/c/d /e is mpeg2 */
/* look_bat[abcd_index][4][16] */
/* Read only */
static unsigned char look_bat[5][4][16] =
{
/* LOOK_BATA */
{{0, 1, 3, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17},
{0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 17},
{0, 1, 2, 3, 4, 5, 6, 17, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 1, 2, 17, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
/* LOOK_BATB */
{{0, 1, 3, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17},
{0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 17},
{0, 1, 2, 3, 4, 5, 6, 17, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 1, 2, 17, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
/* LOOK_BATC */
{{0, 1, 2, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 1, 2, 4, 5, 6, 7, 8, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
/* LOOK_BATD */
{{0, 1, 2, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 1, 2, 4, 5, 6, 7, 8, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
/* LOOK_BATE */
{{0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15},
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 1, 2, 4, 5, 6, 7, 8, 0, 0, 0, 0, 0, 0, 0, 0},
{0, 1, 2, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}},
};
/* look_nbat[abcd_index]][4] */
/* Read only */
static unsigned char look_nbat[5][4] =
{
{3, 8, 12, 4},
{3, 8, 12, 7},
{2, 0, 6, 0},
{2, 0, 10, 0},
{4, 0, 7, 19},
};
void i_sbt_mono(SAMPLEINT * sample, short *pcm, int n);
void i_sbt_dual(SAMPLEINT * sample, short *pcm, int n);
void i_sbt_dual_mono(SAMPLEINT * sample, short *pcm, int n);
void i_sbt_dual_left(SAMPLEINT * sample, short *pcm, int n);
void i_sbt_dual_right(SAMPLEINT * sample, short *pcm, int n);
void i_sbt16_mono(SAMPLEINT * sample, short *pcm, int n);
void i_sbt16_dual(SAMPLEINT * sample, short *pcm, int n);
void i_sbt16_dual_mono(SAMPLEINT * sample, short *pcm, int n);
void i_sbt16_dual_left(SAMPLEINT * sample, short *pcm, int n);
void i_sbt16_dual_right(SAMPLEINT * sample, short *pcm, int n);
void i_sbt8_mono(SAMPLEINT * sample, short *pcm, int n);
void i_sbt8_dual(SAMPLEINT * sample, short *pcm, int n);
void i_sbt8_dual_mono(SAMPLEINT * sample, short *pcm, int n);
void i_sbt8_dual_left(SAMPLEINT * sample, short *pcm, int n);
void i_sbt8_dual_right(SAMPLEINT * sample, short *pcm, int n);
/*--- 8 bit output ---*/
void i_sbtB_mono(SAMPLEINT * sample, unsigned char *pcm, int n);
void i_sbtB_dual(SAMPLEINT * sample, unsigned char *pcm, int n);
void i_sbtB_dual_mono(SAMPLEINT * sample, unsigned char *pcm, int n);
void i_sbtB_dual_left(SAMPLEINT * sample, unsigned char *pcm, int n);
void i_sbtB_dual_right(SAMPLEINT * sample, unsigned char *pcm, int n);
void i_sbtB16_mono(SAMPLEINT * sample, unsigned char *pcm, int n);
void i_sbtB16_dual(SAMPLEINT * sample, unsigned char *pcm, int n);
void i_sbtB16_dual_mono(SAMPLEINT * sample, unsigned char *pcm, int n);
void i_sbtB16_dual_left(SAMPLEINT * sample, unsigned char *pcm, int n);
void i_sbtB16_dual_right(SAMPLEINT * sample, unsigned char *pcm, int n);
void i_sbtB8_mono(SAMPLEINT * sample, unsigned char *pcm, int n);
void i_sbtB8_dual(SAMPLEINT * sample, unsigned char *pcm, int n);
void i_sbtB8_dual_mono(SAMPLEINT * sample, unsigned char *pcm, int n);
void i_sbtB8_dual_left(SAMPLEINT * sample, unsigned char *pcm, int n);
void i_sbtB8_dual_right(SAMPLEINT * sample, unsigned char *pcm, int n);
/* Read Only */
static SBT_FUNCTION sbt_table[2][3][5] =
{
{{i_sbt_mono, i_sbt_dual, i_sbt_dual_mono, i_sbt_dual_left, i_sbt_dual_right},
{i_sbt16_mono, i_sbt16_dual, i_sbt16_dual_mono, i_sbt16_dual_left, i_sbt16_dual_right},
{i_sbt8_mono, i_sbt8_dual, i_sbt8_dual_mono, i_sbt8_dual_left, i_sbt8_dual_right}},
{{(SBT_FUNCTION) i_sbtB_mono,
(SBT_FUNCTION) i_sbtB_dual,
(SBT_FUNCTION) i_sbtB_dual_mono,
(SBT_FUNCTION) i_sbtB_dual_left,
(SBT_FUNCTION) i_sbtB_dual_right},
{(SBT_FUNCTION) i_sbtB16_mono,
(SBT_FUNCTION) i_sbtB16_dual,
(SBT_FUNCTION) i_sbtB16_dual_mono,
(SBT_FUNCTION) i_sbtB16_dual_left,
(SBT_FUNCTION) i_sbtB16_dual_right},
{(SBT_FUNCTION) i_sbtB8_mono,
(SBT_FUNCTION) i_sbtB8_dual,
(SBT_FUNCTION) i_sbtB8_dual_mono,
(SBT_FUNCTION) i_sbtB8_dual_left,
(SBT_FUNCTION) i_sbtB8_dual_right}},
};
/* Read Only */
static int out_chans[5] =
{1, 2, 1, 1, 1};
/*---------------------------------------------------------*/
#ifdef _MSC_VER
#pragma warning(disable: 4056)
#endif
void i_mpeg_init(MPEGI *m)
{
memset(&m->iup, 0, sizeof(m->iup));
m->iup.nsb_limit = 6;
m->iup.nbat[0] = 3;
m->iup.nbat[1] = 8;
m->iup.nbat[3] = 12;
m->iup.nbat[4] = 7;
m->iup.nsbt = 36;
m->iup.sbt = i_sbt_mono;
m->iup.unpack_routine = unpack;
m->iup.first_pass = 1;
m->iup.first_pass_L1 = 1;
m->iup.nbatL1 = 32;
m->iup.cs_factorL1 = m->iup.cs_factor[0];
}
static void table_init(MPEGI *m)
{
int i, j;
int code;
int bits;
long tmp, sfmax;
/*-- c_values (dequant) --*/
for (i = 1; i < 18; i++)
m->iup.look_c_value[i] = (int) (32768.0 * 2.0 / steps[i]);
for (i = 1; i < 18; i++)
m->iup.look_c_shift[i] = 16 - stepbits[i];
/*-- scale factor table, scale by 32768 for 16 pcm output --*/
bits = min(8 * sizeof(SAMPLEINT), 8 * sizeof(m->iup.sf_table[0]));
tmp = 1L << (bits - 2);
sfmax = tmp + (tmp - 1);
for (i = 0; i < 64; i++)
{
tmp = (long) (32768.0 * 2.0 * pow(2.0, -i / 3.0));
if (tmp > sfmax)
tmp = sfmax;
m->iup.sf_table[i] = tmp;
}
/*-- grouped 3 level lookup table 5 bit token --*/
for (i = 0; i < 32; i++)
{
code = i;
for (j = 0; j < 3; j++)
{
m->iup.group3_table[i][j] = (char) ((code % 3) - 1);
code /= 3;
}
}
/*-- grouped 5 level lookup table 7 bit token --*/
for (i = 0; i < 128; i++)
{
code = i;
for (j = 0; j < 3; j++)
{
m->iup.group5_table[i][j] = (char) ((code % 5) - 2);
code /= 5;
}
}
/*-- grouped 9 level lookup table 10 bit token --*/
for (i = 0; i < 1024; i++)
{
code = i;
for (j = 0; j < 3; j++)
{
m->iup.group9_table[i][j] = (short) ((code % 9) - 4);
code /= 9;
}
}
}
#ifdef _MSC_VER
#pragma warning(default: 4056)
#endif
/*---------------------------------------------------------*/
int i_audio_decode_initL1(MPEGI *m, MPEG_HEAD * h, int framebytes_arg,
int reduction_code, int transform_code, int convert_code,
int freq_limit);
void i_sbt_init();
/*---------------------------------------------------------*/
/* mpeg_head defined in mhead.h frame bytes is without pad */
int i_audio_decode_init(MPEGI *m, MPEG_HEAD * h, int framebytes_arg,
int reduction_code, int transform_code, int convert_code,
int freq_limit)
{
int i, j, k;
int abcd_index;
long samprate;
int limit;
int bit_code;
if (m->iup.first_pass)
{
table_init(m);
m->iup.first_pass = 0;
}
/* check if code handles */
if (h->option == 3) /* layer I */
return i_audio_decode_initL1(m, h, framebytes_arg,
reduction_code, transform_code, convert_code, freq_limit);
if (h->option != 2)
return 0; /* layer II only */
m->iup.unpack_routine = unpack;
transform_code = transform_code; /* not used, asm compatability */
bit_code = 0;
if (convert_code & 8)
bit_code = 1;
convert_code = convert_code & 3; /* higher bits used by dec8 freq cvt */
if (reduction_code < 0)
reduction_code = 0;
if (reduction_code > 2)
reduction_code = 2;
if (freq_limit < 1000)
freq_limit = 1000;
m->iup.framebytes = framebytes_arg;
/* compute abcd index for bit allo table selection */
if (h->id) /* mpeg 1 */
abcd_index = lookqt[h->mode][h->sr_index][h->br_index];
else
abcd_index = 4; /* mpeg 2 */
for (i = 0; i < 4; i++)
for (j = 0; j < 16; j++)
m->iup.bat[i][j] = look_bat[abcd_index][i][j];
for (i = 0; i < 4; i++)
m->iup.nbat[i] = look_nbat[abcd_index][i];
m->iup.max_sb = m->iup.nbat[0] + m->iup.nbat[1] + m->iup.nbat[2] + m->iup.nbat[3];
/*----- compute nsb_limit --------*/
samprate = sr_table[4 * h->id + h->sr_index];
m->iup.nsb_limit = (freq_limit * 64L + samprate / 2) / samprate;
/*- caller limit -*/
/*---- limit = 0.94*(32>>reduction_code); ----*/
limit = (32 >> reduction_code);
if (limit > 8)
limit--;
if (m->iup.nsb_limit > limit)
m->iup.nsb_limit = limit;
if (m->iup.nsb_limit > m->iup.max_sb)
m->iup.nsb_limit = m->iup.max_sb;
m->iup.outvalues = 1152 >> reduction_code;
if (h->mode != 3)
{ /* adjust for 2 channel modes */
for (i = 0; i < 4; i++)
m->iup.nbat[i] *= 2;
m->iup.max_sb *= 2;
m->iup.nsb_limit *= 2;
}
/* set sbt function */
m->iup.nsbt = 36;
k = 1 + convert_code;
if (h->mode == 3)
{
k = 0;
}
m->iup.sbt = sbt_table[bit_code][reduction_code][k];
m->iup.outvalues *= out_chans[k];
if (bit_code != 0)
m->iup.outbytes = m->iup.outvalues;
else
m->iup.outbytes = sizeof(short) * m->iup.outvalues;
m->iup.decinfo.channels = out_chans[k];
m->iup.decinfo.outvalues = m->iup.outvalues;
m->iup.decinfo.samprate = samprate >> reduction_code;
if (bit_code != 0)
m->iup.decinfo.bits = 8;
else
m->iup.decinfo.bits = sizeof(short) * 8;
m->iup.decinfo.framebytes = m->iup.framebytes;
m->iup.decinfo.type = 0;
/* clear sample buffer, unused sub bands must be 0 */
for (i = 0; i < 2304; i++)
m->iup.sample[i] = 0;
/* init sub-band transform */
i_sbt_init();
return 1;
}
/*---------------------------------------------------------*/
void i_audio_decode_info(MPEGI *m, DEC_INFO * info)
{
*info = m->iup.decinfo; /* info return, call after init */
}
/*---------------------------------------------------------*/
Reference in New Issue View Git Blame Copy Permalink