From 878dd470f235a2c945f2410927175a45f052b234 Mon Sep 17 00:00:00 2001 From: markster Date: Tue, 7 Dec 1999 05:45:48 +0000 Subject: Version 0.1.1 from FTP git-svn-id: http://svn.digium.com/svn/asterisk/trunk@95 f38db490-d61c-443f-a65b-d21fe96a405b --- codecs/mp3/src/iupini.c | 389 ++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 389 insertions(+) create mode 100755 codecs/mp3/src/iupini.c (limited to 'codecs/mp3/src/iupini.c') diff --git a/codecs/mp3/src/iupini.c b/codecs/mp3/src/iupini.c new file mode 100755 index 000000000..83d5793b4 --- /dev/null +++ b/codecs/mp3/src/iupini.c @@ -0,0 +1,389 @@ +/*____________________________________________________________________________ + + 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 + +/* 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 */ +} +/*---------------------------------------------------------*/ -- cgit v1.2.3