/*____________________________________________________________________________ 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 cup.c - include to cup.c mpeg audio decoder portable "c" mod 8/6/96 add 8 bit output mod 5/10/95 add quick (low precision) window mod 5/16/95 sb limit for reduced samprate output changed from 94% to 100% of Nyquist sb mod 11/15/95 for Layer I =========================================================*/ /*-- compiler bug, floating constant overflow w/ansi --*/ #ifdef _MSC_VER #pragma warning(disable:4056) #endif /* Read Only */ static long steps[18] = { 0, 3, 5, 7, 9, 15, 31, 63, 127, 255, 511, 1023, 2047, 4095, 8191, 16383, 32767, 65535}; /* 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 */ // mono extended beyond legal br index // 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,0,0,0,0,-1, /* 48ks */ // 1,3,3,0,0,0,1,1,1,1,1,1,1,1,1,-1, /* 32ks */ // legal mono {{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 sbt_mono(MPEG *m, float *sample, short *pcm, int n); void sbt_dual(MPEG *m, float *sample, short *pcm, int n); void sbt_dual_mono(MPEG *m, float *sample, short *pcm, int n); void sbt_dual_left(MPEG *m, float *sample, short *pcm, int n); void sbt_dual_right(MPEG *m, float *sample, short *pcm, int n); void sbt16_mono(MPEG *m, float *sample, short *pcm, int n); void sbt16_dual(MPEG *m, float *sample, short *pcm, int n); void sbt16_dual_mono(MPEG *m, float *sample, short *pcm, int n); void sbt16_dual_left(MPEG *m, float *sample, short *pcm, int n); void sbt16_dual_right(MPEG *m, float *sample, short *pcm, int n); void sbt8_mono(MPEG *m, float *sample, short *pcm, int n); void sbt8_dual(MPEG *m, float *sample, short *pcm, int n); void sbt8_dual_mono(MPEG *m, float *sample, short *pcm, int n); void sbt8_dual_left(MPEG *m, float *sample, short *pcm, int n); void sbt8_dual_right(MPEG *m, float *sample, short *pcm, int n); /*--- 8 bit output ---*/ void sbtB_mono(MPEG *m, float *sample, unsigned char *pcm, int n); void sbtB_dual(MPEG *m, float *sample, unsigned char *pcm, int n); void sbtB_dual_mono(MPEG *m, float *sample, unsigned char *pcm, int n); void sbtB_dual_left(MPEG *m, float *sample, unsigned char *pcm, int n); void sbtB_dual_right(MPEG *m, float *sample, unsigned char *pcm, int n); void sbtB16_mono(MPEG *m, float *sample, unsigned char *pcm, int n); void sbtB16_dual(MPEG *m, float *sample, unsigned char *pcm, int n); void sbtB16_dual_mono(MPEG *m, float *sample, unsigned char *pcm, int n); void sbtB16_dual_left(MPEG *m, float *sample, unsigned char *pcm, int n); void sbtB16_dual_right(MPEG *m, float *sample, unsigned char *pcm, int n); void sbtB8_mono(MPEG *m, float *sample, unsigned char *pcm, int n); void sbtB8_dual(MPEG *m, float *sample, unsigned char *pcm, int n); void sbtB8_dual_mono(MPEG *m, float *sample, unsigned char *pcm, int n); void sbtB8_dual_left(MPEG *m, float *sample, unsigned char *pcm, int n); void sbtB8_dual_right(MPEG *m, float *sample, unsigned char *pcm, int n); static SBT_FUNCTION_F sbt_table[2][3][5] = { {{sbt_mono, sbt_dual, sbt_dual_mono, sbt_dual_left, sbt_dual_right}, {sbt16_mono, sbt16_dual, sbt16_dual_mono, sbt16_dual_left, sbt16_dual_right}, {sbt8_mono, sbt8_dual, sbt8_dual_mono, sbt8_dual_left, sbt8_dual_right}}, {{(SBT_FUNCTION_F) sbtB_mono, (SBT_FUNCTION_F) sbtB_dual, (SBT_FUNCTION_F) sbtB_dual_mono, (SBT_FUNCTION_F) sbtB_dual_left, (SBT_FUNCTION_F) sbtB_dual_right}, {(SBT_FUNCTION_F) sbtB16_mono, (SBT_FUNCTION_F) sbtB16_dual, (SBT_FUNCTION_F) sbtB16_dual_mono, (SBT_FUNCTION_F) sbtB16_dual_left, (SBT_FUNCTION_F) sbtB16_dual_right}, {(SBT_FUNCTION_F) sbtB8_mono, (SBT_FUNCTION_F) sbtB8_dual, (SBT_FUNCTION_F) sbtB8_dual_mono, (SBT_FUNCTION_F) sbtB8_dual_left, (SBT_FUNCTION_F) sbtB8_dual_right}}, }; static int out_chans[5] = {1, 2, 1, 1, 1}; int audio_decode_initL1(MPEG_HEAD * h, int framebytes_arg, int reduction_code, int transform_code, int convert_code, int freq_limit); void sbt_init(); IN_OUT L1audio_decode(void *mv, unsigned char *bs, signed short *pcm); IN_OUT L2audio_decode(void *mv, unsigned char *bs, signed short *pcm); IN_OUT L3audio_decode(void *mv, unsigned char *bs, unsigned char *pcm); static AUDIO_DECODE_ROUTINE decode_routine_table[4] = { L2audio_decode, (AUDIO_DECODE_ROUTINE)L3audio_decode, L2audio_decode, L1audio_decode,}; extern void cup3_init(MPEG *m); void mpeg_init(MPEG *m) { memset(m, 0, sizeof(MPEG)); m->cup.nsb_limit = 6; m->cup.nbat[0] = 3; m->cup.nbat[1] = 8; m->cup.nbat[3] = 12; m->cup.nbat[4] = 7; m->cup.sbt = sbt_mono; m->cup.first_pass = 1; m->cup.first_pass_L1 = 1; m->cup.audio_decode_routine = L2audio_decode; m->cup.cs_factorL1 = m->cup.cs_factor[0]; m->cup.nbatL1 = 32; m->cupl.band_limit = 576; m->cupl.band_limit21 = 567; m->cupl.band_limit12 = 576; m->cupl.band_limit_nsb = 32; m->cupl.nsb_limit=32; m->cup.sample = (float *)&m->cupl.sample; m->csbt.first_pass = 1; cup3_init(m); } /*---------------------------------------------------------*/ static void table_init(MPEG *m) { int i, j; int code; /*-- c_values (dequant) --*/ for (i = 1; i < 18; i++) m->cup.look_c_value[i] = 2.0F / steps[i]; /*-- scale factor table, scale by 32768 for 16 pcm output --*/ for (i = 0; i < 64; i++) m->cup.sf_table[i] = (float) (32768.0 * 2.0 * pow(2.0, -i / 3.0)); /*-- grouped 3 level lookup table 5 bit token --*/ for (i = 0; i < 32; i++) { code = i; for (j = 0; j < 3; j++) { m->cup.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->cup.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->cup.group9_table[i][j] = (short) ((code % 9) - 4); code /= 9; } } } /*---------------------------------------------------------*/ int L1audio_decode_init(MPEG *m, MPEG_HEAD * h, int framebytes_arg, int reduction_code, int transform_code, int convert_code, int freq_limit); int L3audio_decode_init(MPEG *m, MPEG_HEAD * h, int framebytes_arg, int reduction_code, int transform_code, int convert_code, int freq_limit); /*---------------------------------------------------------*/ /* mpeg_head defined in mhead.h frame bytes is without pad */ int audio_decode_init(MPEG *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->cup.first_pass) { table_init(m); m->cup.first_pass = 0; } /* select decoder routine Layer I,II,III */ m->cup.audio_decode_routine = decode_routine_table[h->option & 3]; if (h->option == 3) /* layer I */ return L1audio_decode_init(m, h, framebytes_arg, reduction_code, transform_code, convert_code, freq_limit); if (h->option == 1) /* layer III */ return L3audio_decode_init(m, h, framebytes_arg, reduction_code, transform_code, convert_code, freq_limit); 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->cup.framebytes = framebytes_arg; /* check if code handles */ if (h->option != 2) return 0; /* layer II only */ if (h->sr_index == 3) return 0; /* reserved */ /* 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 */ if (abcd_index < 0) return 0; // fail invalid Layer II bit rate index for (i = 0; i < 4; i++) for (j = 0; j < 16; j++) m->cup.bat[i][j] = look_bat[abcd_index][i][j]; for (i = 0; i < 4; i++) m->cup.nbat[i] = look_nbat[abcd_index][i]; m->cup.max_sb = m->cup.nbat[0] + m->cup.nbat[1] + m->cup.nbat[2] + m->cup.nbat[3]; /*----- compute nsb_limit --------*/ samprate = sr_table[4 * h->id + h->sr_index]; m->cup.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->cup.nsb_limit > limit) m->cup.nsb_limit = limit; if (m->cup.nsb_limit > m->cup.max_sb) m->cup.nsb_limit = m->cup.max_sb; m->cup.outvalues = 1152 >> reduction_code; if (h->mode != 3) { /* adjust for 2 channel modes */ for (i = 0; i < 4; i++) m->cup.nbat[i] *= 2; m->cup.max_sb *= 2; m->cup.nsb_limit *= 2; } /* set sbt function */ k = 1 + convert_code; if (h->mode == 3) { k = 0; } m->cup.sbt = sbt_table[bit_code][reduction_code][k]; m->cup.outvalues *= out_chans[k]; if (bit_code) m->cup.outbytes = m->cup.outvalues; else m->cup.outbytes = sizeof(short) * m->cup.outvalues; m->cup.decinfo.channels = out_chans[k]; m->cup.decinfo.outvalues = m->cup.outvalues; m->cup.decinfo.samprate = samprate >> reduction_code; if (bit_code) m->cup.decinfo.bits = 8; else m->cup.decinfo.bits = sizeof(short) * 8; m->cup.decinfo.framebytes = m->cup.framebytes; m->cup.decinfo.type = 0; /* clear sample buffer, unused sub bands must be 0 */ for (i = 0; i < 2304; i++) m->cup.sample[i] = 0.0F; /* init sub-band transform */ sbt_init(); return 1; } /*---------------------------------------------------------*/ void audio_decode_info(MPEG *m, DEC_INFO * info) { *info = m->cup.decinfo; /* info return, call after init */ } /*---------------------------------------------------------*/ void decode_table_init(MPEG *m) { /* dummy for asm version compatability */ } /*---------------------------------------------------------*/