summaryrefslogtreecommitdiff
path: root/codecs/mp3/src/cdct.c
diff options
context:
space:
mode:
Diffstat (limited to 'codecs/mp3/src/cdct.c')
-rwxr-xr-xcodecs/mp3/src/cdct.c349
1 files changed, 0 insertions, 349 deletions
diff --git a/codecs/mp3/src/cdct.c b/codecs/mp3/src/cdct.c
deleted file mode 100755
index f9ce28a7e..000000000
--- a/codecs/mp3/src/cdct.c
+++ /dev/null
@@ -1,349 +0,0 @@
-/*____________________________________________________________________________
-
- 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$
-____________________________________________________________________________*/
-
-/**** cdct.c ***************************************************
-
-mod 5/16/95 first stage in 8 pt dct does not drop last sb mono
-
-
-MPEG audio decoder, dct
-portable C
-
-******************************************************************/
-
-#include <stdlib.h>
-#include <stdio.h>
-#include <float.h>
-#include <math.h>
-#include "L3.h"
-#include "mhead.h"
-
-#ifdef ASM_X86
-extern void fdct32_asm(float*a, float*b);
-extern void fdct32_dual_asm(float*a, float*b);
-#endif /* ASM_X86 */
-
-/*------------------------------------------------------------*/
-float *dct_coef_addr(MPEG *m)
-{
- return m->cdct.coef32;
-}
-/*------------------------------------------------------------*/
-static void forward_bf(int m, int n, float x[], float f[], float coef[])
-{
- int i, j, n2;
- int p, q, p0, k;
-
- p0 = 0;
- n2 = n >> 1;
- for (i = 0; i < m; i++, p0 += n)
- {
- k = 0;
- p = p0;
- q = p + n - 1;
- for (j = 0; j < n2; j++, p++, q--, k++)
- {
- f[p] = x[p] + x[q];
- f[n2 + p] = coef[k] * (x[p] - x[q]);
- }
- }
-}
-/*------------------------------------------------------------*/
-static void back_bf(int m, int n, float x[], float f[])
-{
- int i, j, n2, n21;
- int p, q, p0;
-
- p0 = 0;
- n2 = n >> 1;
- n21 = n2 - 1;
- for (i = 0; i < m; i++, p0 += n)
- {
- p = p0;
- q = p0;
- for (j = 0; j < n2; j++, p += 2, q++)
- f[p] = x[q];
- p = p0 + 1;
- for (j = 0; j < n21; j++, p += 2, q++)
- f[p] = x[q] + x[q + 1];
- f[p] = x[q];
- }
-}
-/*------------------------------------------------------------*/
-
-#ifdef _EQUALIZER_ENABLE_
-extern float equalizer[32];
-extern int enableEQ;
-#endif
-
-void fdct32(MPEG *m, float x[], float c[])
-{
-#if (!defined(ASM_X86) && !defined(ASM_X86_OLD) || defined(_EQUALIZER_ENABLE_))
- float a[32]; /* ping pong buffers */
- float b[32];
- int p, q;
-#endif
-
- float *src = x;
-
-#ifdef _EQUALIZER_ENABLE_
- int i;
- float b[32];
- if (enableEQ) {
- for(i=0; i<32; i++)
- b[i] = x[i] * equalizer[i];
- src = b;
- }
-#endif /* _EQUALIZER_ENABLE_ */
-#undef _EQUALIZER_ENABLE_
-
-#ifdef ASM_X86
- fdct32_asm(src, c);
-#elif defined(ASM_X86_OLD)
- asm_fdct32(src, c);
-#else
-/* special first stage */
- for (p = 0, q = 31; p < 16; p++, q--)
- {
- a[p] = src[p] + src[q];
- a[16 + p] = m->cdct.coef32[p] * (src[p] - src[q]);
- }
- forward_bf(2, 16, a, b, m->cdct.coef32 + 16);
- forward_bf(4, 8, b, a, m->cdct.coef32 + 16 + 8);
- forward_bf(8, 4, a, b, m->cdct.coef32 + 16 + 8 + 4);
- forward_bf(16, 2, b, a, m->cdct.coef32 + 16 + 8 + 4 + 2);
- back_bf(8, 4, a, b);
- back_bf(4, 8, b, a);
- back_bf(2, 16, a, b);
- back_bf(1, 32, b, c);
-#endif
-}
-/*------------------------------------------------------------*/
-void fdct32_dual(MPEG *m, float x[], float c[])
-{
-#ifdef ASM_X86
- fdct32_dual_asm(x, c);
-#else
- float a[32]; /* ping pong buffers */
- float b[32];
- int p, pp, qq;
-
-/* special first stage for dual chan (interleaved x) */
- pp = 0;
- qq = 2 * 31;
- for (p = 0; p < 16; p++, pp += 2, qq -= 2)
- {
- a[p] = x[pp] + x[qq];
- a[16 + p] = m->cdct.coef32[p] * (x[pp] - x[qq]);
- }
- forward_bf(2, 16, a, b, m->cdct.coef32 + 16);
- forward_bf(4, 8, b, a, m->cdct.coef32 + 16 + 8);
- forward_bf(8, 4, a, b, m->cdct.coef32 + 16 + 8 + 4);
- forward_bf(16, 2, b, a, m->cdct.coef32 + 16 + 8 + 4 + 2);
- back_bf(8, 4, a, b);
- back_bf(4, 8, b, a);
- back_bf(2, 16, a, b);
- back_bf(1, 32, b, c);
-#endif
-}
-/*---------------convert dual to mono------------------------------*/
-void fdct32_dual_mono(MPEG *m, float x[], float c[])
-{
- float a[32]; /* ping pong buffers */
- float b[32];
- float t1, t2;
- int p, pp, qq;
-
-/* special first stage */
- pp = 0;
- qq = 2 * 31;
- for (p = 0; p < 16; p++, pp += 2, qq -= 2)
- {
- t1 = 0.5F * (x[pp] + x[pp + 1]);
- t2 = 0.5F * (x[qq] + x[qq + 1]);
- a[p] = t1 + t2;
- a[16 + p] = m->cdct.coef32[p] * (t1 - t2);
- }
- forward_bf(2, 16, a, b, m->cdct.coef32 + 16);
- forward_bf(4, 8, b, a, m->cdct.coef32 + 16 + 8);
- forward_bf(8, 4, a, b, m->cdct.coef32 + 16 + 8 + 4);
- forward_bf(16, 2, b, a, m->cdct.coef32 + 16 + 8 + 4 + 2);
- back_bf(8, 4, a, b);
- back_bf(4, 8, b, a);
- back_bf(2, 16, a, b);
- back_bf(1, 32, b, c);
-}
-/*------------------------------------------------------------*/
-/*---------------- 16 pt fdct -------------------------------*/
-void fdct16(MPEG *m, float x[], float c[])
-{
- float a[16]; /* ping pong buffers */
- float b[16];
- int p, q;
-
-/* special first stage (drop highest sb) */
- a[0] = x[0];
- a[8] = m->cdct.coef32[16] * x[0];
- for (p = 1, q = 14; p < 8; p++, q--)
- {
- a[p] = x[p] + x[q];
- a[8 + p] = m->cdct.coef32[16 + p] * (x[p] - x[q]);
- }
- forward_bf(2, 8, a, b, m->cdct.coef32 + 16 + 8);
- forward_bf(4, 4, b, a, m->cdct.coef32 + 16 + 8 + 4);
- forward_bf(8, 2, a, b, m->cdct.coef32 + 16 + 8 + 4 + 2);
- back_bf(4, 4, b, a);
- back_bf(2, 8, a, b);
- back_bf(1, 16, b, c);
-}
-/*------------------------------------------------------------*/
-/*---------------- 16 pt fdct dual chan---------------------*/
-void fdct16_dual(MPEG *m, float x[], float c[])
-{
- float a[16]; /* ping pong buffers */
- float b[16];
- int p, pp, qq;
-
-/* special first stage for interleaved input */
- a[0] = x[0];
- a[8] = m->cdct.coef32[16] * x[0];
- pp = 2;
- qq = 2 * 14;
- for (p = 1; p < 8; p++, pp += 2, qq -= 2)
- {
- a[p] = x[pp] + x[qq];
- a[8 + p] = m->cdct.coef32[16 + p] * (x[pp] - x[qq]);
- }
- forward_bf(2, 8, a, b, m->cdct.coef32 + 16 + 8);
- forward_bf(4, 4, b, a, m->cdct.coef32 + 16 + 8 + 4);
- forward_bf(8, 2, a, b, m->cdct.coef32 + 16 + 8 + 4 + 2);
- back_bf(4, 4, b, a);
- back_bf(2, 8, a, b);
- back_bf(1, 16, b, c);
-}
-/*------------------------------------------------------------*/
-/*---------------- 16 pt fdct dual to mono-------------------*/
-void fdct16_dual_mono(MPEG *m, float x[], float c[])
-{
- float a[16]; /* ping pong buffers */
- float b[16];
- float t1, t2;
- int p, pp, qq;
-
-/* special first stage */
- a[0] = 0.5F * (x[0] + x[1]);
- a[8] = m->cdct.coef32[16] * a[0];
- pp = 2;
- qq = 2 * 14;
- for (p = 1; p < 8; p++, pp += 2, qq -= 2)
- {
- t1 = 0.5F * (x[pp] + x[pp + 1]);
- t2 = 0.5F * (x[qq] + x[qq + 1]);
- a[p] = t1 + t2;
- a[8 + p] = m->cdct.coef32[16 + p] * (t1 - t2);
- }
- forward_bf(2, 8, a, b, m->cdct.coef32 + 16 + 8);
- forward_bf(4, 4, b, a, m->cdct.coef32 + 16 + 8 + 4);
- forward_bf(8, 2, a, b, m->cdct.coef32 + 16 + 8 + 4 + 2);
- back_bf(4, 4, b, a);
- back_bf(2, 8, a, b);
- back_bf(1, 16, b, c);
-}
-/*------------------------------------------------------------*/
-/*---------------- 8 pt fdct -------------------------------*/
-void fdct8(MPEG *m, float x[], float c[])
-{
- float a[8]; /* ping pong buffers */
- float b[8];
- int p, q;
-
-/* special first stage */
-
- b[0] = x[0] + x[7];
- b[4] = m->cdct.coef32[16 + 8] * (x[0] - x[7]);
- for (p = 1, q = 6; p < 4; p++, q--)
- {
- b[p] = x[p] + x[q];
- b[4 + p] = m->cdct.coef32[16 + 8 + p] * (x[p] - x[q]);
- }
-
- forward_bf(2, 4, b, a, m->cdct.coef32 + 16 + 8 + 4);
- forward_bf(4, 2, a, b, m->cdct.coef32 + 16 + 8 + 4 + 2);
- back_bf(2, 4, b, a);
- back_bf(1, 8, a, c);
-}
-/*------------------------------------------------------------*/
-/*---------------- 8 pt fdct dual chan---------------------*/
-void fdct8_dual(MPEG *m, float x[], float c[])
-{
- float a[8]; /* ping pong buffers */
- float b[8];
- int p, pp, qq;
-
-/* special first stage for interleaved input */
- b[0] = x[0] + x[14];
- b[4] = m->cdct.coef32[16 + 8] * (x[0] - x[14]);
- pp = 2;
- qq = 2 * 6;
- for (p = 1; p < 4; p++, pp += 2, qq -= 2)
- {
- b[p] = x[pp] + x[qq];
- b[4 + p] = m->cdct.coef32[16 + 8 + p] * (x[pp] - x[qq]);
- }
- forward_bf(2, 4, b, a, m->cdct.coef32 + 16 + 8 + 4);
- forward_bf(4, 2, a, b, m->cdct.coef32 + 16 + 8 + 4 + 2);
- back_bf(2, 4, b, a);
- back_bf(1, 8, a, c);
-}
-/*------------------------------------------------------------*/
-/*---------------- 8 pt fdct dual to mono---------------------*/
-void fdct8_dual_mono(MPEG *m, float x[], float c[])
-{
- float a[8]; /* ping pong buffers */
- float b[8];
- float t1, t2;
- int p, pp, qq;
-
-/* special first stage */
- t1 = 0.5F * (x[0] + x[1]);
- t2 = 0.5F * (x[14] + x[15]);
- b[0] = t1 + t2;
- b[4] = m->cdct.coef32[16 + 8] * (t1 - t2);
- pp = 2;
- qq = 2 * 6;
- for (p = 1; p < 4; p++, pp += 2, qq -= 2)
- {
- t1 = 0.5F * (x[pp] + x[pp + 1]);
- t2 = 0.5F * (x[qq] + x[qq + 1]);
- b[p] = t1 + t2;
- b[4 + p] = m->cdct.coef32[16 + 8 + p] * (t1 - t2);
- }
- forward_bf(2, 4, b, a, m->cdct.coef32 + 16 + 8 + 4);
- forward_bf(4, 2, a, b, m->cdct.coef32 + 16 + 8 + 4 + 2);
- back_bf(2, 4, b, a);
- back_bf(1, 8, a, c);
-}
-/*------------------------------------------------------------*/
generated by cgit v1.2.3 (git 2.39.1) at 2024-07-19 22:45:29 +0000