cannam@90: /*
cannam@90: * libmad - MPEG audio decoder library
cannam@90: * Copyright (C) 2000-2004 Underbit Technologies, Inc.
cannam@90: *
cannam@90: * This program is free software; you can redistribute it and/or modify
cannam@90: * it under the terms of the GNU General Public License as published by
cannam@90: * the Free Software Foundation; either version 2 of the License, or
cannam@90: * (at your option) any later version.
cannam@90: *
cannam@90: * This program is distributed in the hope that it will be useful,
cannam@90: * but WITHOUT ANY WARRANTY; without even the implied warranty of
cannam@90: * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
cannam@90: * GNU General Public License for more details.
cannam@90: *
cannam@90: * You should have received a copy of the GNU General Public License
cannam@90: * along with this program; if not, write to the Free Software
cannam@90: * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
cannam@90: *
cannam@90: * If you would like to negotiate alternate licensing terms, you may do
cannam@90: * so by contacting: Underbit Technologies, Inc. <info@underbit.com>
cannam@90: */
cannam@90:
cannam@90: # ifdef __cplusplus
cannam@90: extern "C" {
cannam@90: # endif
cannam@90:
cannam@90: # define FPM_INTEL
cannam@90:
cannam@90:
cannam@90:
cannam@90: # define SIZEOF_INT 4
cannam@90: # define SIZEOF_LONG 4
cannam@90: # define SIZEOF_LONG_LONG 8
cannam@90:
cannam@90:
cannam@90: /* Id: version.h,v 1.26 2004/01/23 09:41:33 rob Exp */
cannam@90:
cannam@90: # ifndef LIBMAD_VERSION_H
cannam@90: # define LIBMAD_VERSION_H
cannam@90:
cannam@90: # define MAD_VERSION_MAJOR 0
cannam@90: # define MAD_VERSION_MINOR 15
cannam@90: # define MAD_VERSION_PATCH 1
cannam@90: # define MAD_VERSION_EXTRA " (beta)"
cannam@90:
cannam@90: # define MAD_VERSION_STRINGIZE(str) #str
cannam@90: # define MAD_VERSION_STRING(num) MAD_VERSION_STRINGIZE(num)
cannam@90:
cannam@90: # define MAD_VERSION MAD_VERSION_STRING(MAD_VERSION_MAJOR) "." \
cannam@90: MAD_VERSION_STRING(MAD_VERSION_MINOR) "." \
cannam@90: MAD_VERSION_STRING(MAD_VERSION_PATCH) \
cannam@90: MAD_VERSION_EXTRA
cannam@90:
cannam@90: # define MAD_PUBLISHYEAR "2000-2004"
cannam@90: # define MAD_AUTHOR "Underbit Technologies, Inc."
cannam@90: # define MAD_EMAIL "info@underbit.com"
cannam@90:
cannam@90: extern char const mad_version[];
cannam@90: extern char const mad_copyright[];
cannam@90: extern char const mad_author[];
cannam@90: extern char const mad_build[];
cannam@90:
cannam@90: # endif
cannam@90:
cannam@90: /* Id: fixed.h,v 1.38 2004/02/17 02:02:03 rob Exp */
cannam@90:
cannam@90: # ifndef LIBMAD_FIXED_H
cannam@90: # define LIBMAD_FIXED_H
cannam@90:
cannam@90: # if SIZEOF_INT >= 4
cannam@90: typedef signed int mad_fixed_t;
cannam@90:
cannam@90: typedef signed int mad_fixed64hi_t;
cannam@90: typedef unsigned int mad_fixed64lo_t;
cannam@90: # else
cannam@90: typedef signed long mad_fixed_t;
cannam@90:
cannam@90: typedef signed long mad_fixed64hi_t;
cannam@90: typedef unsigned long mad_fixed64lo_t;
cannam@90: # endif
cannam@90:
cannam@90: # if defined(_MSC_VER)
cannam@90: # define mad_fixed64_t signed __int64
cannam@90: # elif 1 || defined(__GNUC__)
cannam@90: # define mad_fixed64_t signed long long
cannam@90: # endif
cannam@90:
cannam@90: # if defined(FPM_FLOAT)
cannam@90: typedef double mad_sample_t;
cannam@90: # else
cannam@90: typedef mad_fixed_t mad_sample_t;
cannam@90: # endif
cannam@90:
cannam@90: /*
cannam@90: * Fixed-point format: 0xABBBBBBB
cannam@90: * A == whole part (sign + 3 bits)
cannam@90: * B == fractional part (28 bits)
cannam@90: *
cannam@90: * Values are signed two's complement, so the effective range is:
cannam@90: * 0x80000000 to 0x7fffffff
cannam@90: * -8.0 to +7.9999999962747097015380859375
cannam@90: *
cannam@90: * The smallest representable value is:
cannam@90: * 0x00000001 == 0.0000000037252902984619140625 (i.e. about 3.725e-9)
cannam@90: *
cannam@90: * 28 bits of fractional accuracy represent about
cannam@90: * 8.6 digits of decimal accuracy.
cannam@90: *
cannam@90: * Fixed-point numbers can be added or subtracted as normal
cannam@90: * integers, but multiplication requires shifting the 64-bit result
cannam@90: * from 56 fractional bits back to 28 (and rounding.)
cannam@90: *
cannam@90: * Changing the definition of MAD_F_FRACBITS is only partially
cannam@90: * supported, and must be done with care.
cannam@90: */
cannam@90:
cannam@90: # define MAD_F_FRACBITS 28
cannam@90:
cannam@90: # if MAD_F_FRACBITS == 28
cannam@90: # define MAD_F(x) ((mad_fixed_t) (x##L))
cannam@90: # else
cannam@90: # if MAD_F_FRACBITS < 28
cannam@90: # warning "MAD_F_FRACBITS < 28"
cannam@90: # define MAD_F(x) ((mad_fixed_t) \
cannam@90: (((x##L) + \
cannam@90: (1L << (28 - MAD_F_FRACBITS - 1))) >> \
cannam@90: (28 - MAD_F_FRACBITS)))
cannam@90: # elif MAD_F_FRACBITS > 28
cannam@90: # error "MAD_F_FRACBITS > 28 not currently supported"
cannam@90: # define MAD_F(x) ((mad_fixed_t) \
cannam@90: ((x##L) << (MAD_F_FRACBITS - 28)))
cannam@90: # endif
cannam@90: # endif
cannam@90:
cannam@90: # define MAD_F_MIN ((mad_fixed_t) -0x80000000L)
cannam@90: # define MAD_F_MAX ((mad_fixed_t) +0x7fffffffL)
cannam@90:
cannam@90: # define MAD_F_ONE MAD_F(0x10000000)
cannam@90:
cannam@90: # define mad_f_tofixed(x) ((mad_fixed_t) \
cannam@90: ((x) * (double) (1L << MAD_F_FRACBITS) + 0.5))
cannam@90: # define mad_f_todouble(x) ((double) \
cannam@90: ((x) / (double) (1L << MAD_F_FRACBITS)))
cannam@90:
cannam@90: # define mad_f_intpart(x) ((x) >> MAD_F_FRACBITS)
cannam@90: # define mad_f_fracpart(x) ((x) & ((1L << MAD_F_FRACBITS) - 1))
cannam@90: /* (x should be positive) */
cannam@90:
cannam@90: # define mad_f_fromint(x) ((x) << MAD_F_FRACBITS)
cannam@90:
cannam@90: # define mad_f_add(x, y) ((x) + (y))
cannam@90: # define mad_f_sub(x, y) ((x) - (y))
cannam@90:
cannam@90: # if defined(FPM_FLOAT)
cannam@90: # error "FPM_FLOAT not yet supported"
cannam@90:
cannam@90: # undef MAD_F
cannam@90: # define MAD_F(x) mad_f_todouble(x)
cannam@90:
cannam@90: # define mad_f_mul(x, y) ((x) * (y))
cannam@90: # define mad_f_scale64
cannam@90:
cannam@90: # undef ASO_ZEROCHECK
cannam@90:
cannam@90: # elif defined(FPM_64BIT)
cannam@90:
cannam@90: /*
cannam@90: * This version should be the most accurate if 64-bit types are supported by
cannam@90: * the compiler, although it may not be the most efficient.
cannam@90: */
cannam@90: # if defined(OPT_ACCURACY)
cannam@90: # define mad_f_mul(x, y) \
cannam@90: ((mad_fixed_t) \
cannam@90: ((((mad_fixed64_t) (x) * (y)) + \
cannam@90: (1L << (MAD_F_SCALEBITS - 1))) >> MAD_F_SCALEBITS))
cannam@90: # else
cannam@90: # define mad_f_mul(x, y) \
cannam@90: ((mad_fixed_t) (((mad_fixed64_t) (x) * (y)) >> MAD_F_SCALEBITS))
cannam@90: # endif
cannam@90:
cannam@90: # define MAD_F_SCALEBITS MAD_F_FRACBITS
cannam@90:
cannam@90: /* --- Intel --------------------------------------------------------------- */
cannam@90:
cannam@90: # elif defined(FPM_INTEL)
cannam@90:
cannam@90: # if defined(_MSC_VER)
cannam@90: # pragma warning(push)
cannam@90: # pragma warning(disable: 4035) /* no return value */
cannam@90: static __forceinline
cannam@90: mad_fixed_t mad_f_mul_inline(mad_fixed_t x, mad_fixed_t y)
cannam@90: {
cannam@90: enum {
cannam@90: fracbits = MAD_F_FRACBITS
cannam@90: };
cannam@90:
cannam@90: __asm {
cannam@90: mov eax, x
cannam@90: imul y
cannam@90: shrd eax, edx, fracbits
cannam@90: }
cannam@90:
cannam@90: /* implicit return of eax */
cannam@90: }
cannam@90: # pragma warning(pop)
cannam@90:
cannam@90: # define mad_f_mul mad_f_mul_inline
cannam@90: # define mad_f_scale64
cannam@90: # else
cannam@90: /*
cannam@90: * This Intel version is fast and accurate; the disposition of the least
cannam@90: * significant bit depends on OPT_ACCURACY via mad_f_scale64().
cannam@90: */
cannam@90: # define MAD_F_MLX(hi, lo, x, y) \
cannam@90: asm ("imull %3" \
cannam@90: : "=a" (lo), "=d" (hi) \
cannam@90: : "%a" (x), "rm" (y) \
cannam@90: : "cc")
cannam@90:
cannam@90: # if defined(OPT_ACCURACY)
cannam@90: /*
cannam@90: * This gives best accuracy but is not very fast.
cannam@90: */
cannam@90: # define MAD_F_MLA(hi, lo, x, y) \
cannam@90: ({ mad_fixed64hi_t __hi; \
cannam@90: mad_fixed64lo_t __lo; \
cannam@90: MAD_F_MLX(__hi, __lo, (x), (y)); \
cannam@90: asm ("addl %2,%0\n\t" \
cannam@90: "adcl %3,%1" \
cannam@90: : "=rm" (lo), "=rm" (hi) \
cannam@90: : "r" (__lo), "r" (__hi), "0" (lo), "1" (hi) \
cannam@90: : "cc"); \
cannam@90: })
cannam@90: # endif /* OPT_ACCURACY */
cannam@90:
cannam@90: # if defined(OPT_ACCURACY)
cannam@90: /*
cannam@90: * Surprisingly, this is faster than SHRD followed by ADC.
cannam@90: */
cannam@90: # define mad_f_scale64(hi, lo) \
cannam@90: ({ mad_fixed64hi_t __hi_; \
cannam@90: mad_fixed64lo_t __lo_; \
cannam@90: mad_fixed_t __result; \
cannam@90: asm ("addl %4,%2\n\t" \
cannam@90: "adcl %5,%3" \
cannam@90: : "=rm" (__lo_), "=rm" (__hi_) \
cannam@90: : "0" (lo), "1" (hi), \
cannam@90: "ir" (1L << (MAD_F_SCALEBITS - 1)), "ir" (0) \
cannam@90: : "cc"); \
cannam@90: asm ("shrdl %3,%2,%1" \
cannam@90: : "=rm" (__result) \
cannam@90: : "0" (__lo_), "r" (__hi_), "I" (MAD_F_SCALEBITS) \
cannam@90: : "cc"); \
cannam@90: __result; \
cannam@90: })
cannam@90: # elif defined(OPT_INTEL)
cannam@90: /*
cannam@90: * Alternate Intel scaling that may or may not perform better.
cannam@90: */
cannam@90: # define mad_f_scale64(hi, lo) \
cannam@90: ({ mad_fixed_t __result; \
cannam@90: asm ("shrl %3,%1\n\t" \
cannam@90: "shll %4,%2\n\t" \
cannam@90: "orl %2,%1" \
cannam@90: : "=rm" (__result) \
cannam@90: : "0" (lo), "r" (hi), \
cannam@90: "I" (MAD_F_SCALEBITS), "I" (32 - MAD_F_SCALEBITS) \
cannam@90: : "cc"); \
cannam@90: __result; \
cannam@90: })
cannam@90: # else
cannam@90: # define mad_f_scale64(hi, lo) \
cannam@90: ({ mad_fixed_t __result; \
cannam@90: asm ("shrdl %3,%2,%1" \
cannam@90: : "=rm" (__result) \
cannam@90: : "0" (lo), "r" (hi), "I" (MAD_F_SCALEBITS) \
cannam@90: : "cc"); \
cannam@90: __result; \
cannam@90: })
cannam@90: # endif /* OPT_ACCURACY */
cannam@90:
cannam@90: # define MAD_F_SCALEBITS MAD_F_FRACBITS
cannam@90: # endif
cannam@90:
cannam@90: /* --- ARM ----------------------------------------------------------------- */
cannam@90:
cannam@90: # elif defined(FPM_ARM)
cannam@90:
cannam@90: /*
cannam@90: * This ARM V4 version is as accurate as FPM_64BIT but much faster. The
cannam@90: * least significant bit is properly rounded at no CPU cycle cost!
cannam@90: */
cannam@90: # if 1
cannam@90: /*
cannam@90: * This is faster than the default implementation via MAD_F_MLX() and
cannam@90: * mad_f_scale64().
cannam@90: */
cannam@90: # define mad_f_mul(x, y) \
cannam@90: ({ mad_fixed64hi_t __hi; \
cannam@90: mad_fixed64lo_t __lo; \
cannam@90: mad_fixed_t __result; \
cannam@90: asm ("smull %0, %1, %3, %4\n\t" \
cannam@90: "movs %0, %0, lsr %5\n\t" \
cannam@90: "adc %2, %0, %1, lsl %6" \
cannam@90: : "=&r" (__lo), "=&r" (__hi), "=r" (__result) \
cannam@90: : "%r" (x), "r" (y), \
cannam@90: "M" (MAD_F_SCALEBITS), "M" (32 - MAD_F_SCALEBITS) \
cannam@90: : "cc"); \
cannam@90: __result; \
cannam@90: })
cannam@90: # endif
cannam@90:
cannam@90: # define MAD_F_MLX(hi, lo, x, y) \
cannam@90: asm ("smull %0, %1, %2, %3" \
cannam@90: : "=&r" (lo), "=&r" (hi) \
cannam@90: : "%r" (x), "r" (y))
cannam@90:
cannam@90: # define MAD_F_MLA(hi, lo, x, y) \
cannam@90: asm ("smlal %0, %1, %2, %3" \
cannam@90: : "+r" (lo), "+r" (hi) \
cannam@90: : "%r" (x), "r" (y))
cannam@90:
cannam@90: # define MAD_F_MLN(hi, lo) \
cannam@90: asm ("rsbs %0, %2, #0\n\t" \
cannam@90: "rsc %1, %3, #0" \
cannam@90: : "=r" (lo), "=r" (hi) \
cannam@90: : "0" (lo), "1" (hi) \
cannam@90: : "cc")
cannam@90:
cannam@90: # define mad_f_scale64(hi, lo) \
cannam@90: ({ mad_fixed_t __result; \
cannam@90: asm ("movs %0, %1, lsr %3\n\t" \
cannam@90: "adc %0, %0, %2, lsl %4" \
cannam@90: : "=&r" (__result) \
cannam@90: : "r" (lo), "r" (hi), \
cannam@90: "M" (MAD_F_SCALEBITS), "M" (32 - MAD_F_SCALEBITS) \
cannam@90: : "cc"); \
cannam@90: __result; \
cannam@90: })
cannam@90:
cannam@90: # define MAD_F_SCALEBITS MAD_F_FRACBITS
cannam@90:
cannam@90: /* --- MIPS ---------------------------------------------------------------- */
cannam@90:
cannam@90: # elif defined(FPM_MIPS)
cannam@90:
cannam@90: /*
cannam@90: * This MIPS version is fast and accurate; the disposition of the least
cannam@90: * significant bit depends on OPT_ACCURACY via mad_f_scale64().
cannam@90: */
cannam@90: # define MAD_F_MLX(hi, lo, x, y) \
cannam@90: asm ("mult %2,%3" \
cannam@90: : "=l" (lo), "=h" (hi) \
cannam@90: : "%r" (x), "r" (y))
cannam@90:
cannam@90: # if defined(HAVE_MADD_ASM)
cannam@90: # define MAD_F_MLA(hi, lo, x, y) \
cannam@90: asm ("madd %2,%3" \
cannam@90: : "+l" (lo), "+h" (hi) \
cannam@90: : "%r" (x), "r" (y))
cannam@90: # elif defined(HAVE_MADD16_ASM)
cannam@90: /*
cannam@90: * This loses significant accuracy due to the 16-bit integer limit in the
cannam@90: * multiply/accumulate instruction.
cannam@90: */
cannam@90: # define MAD_F_ML0(hi, lo, x, y) \
cannam@90: asm ("mult %2,%3" \
cannam@90: : "=l" (lo), "=h" (hi) \
cannam@90: : "%r" ((x) >> 12), "r" ((y) >> 16))
cannam@90: # define MAD_F_MLA(hi, lo, x, y) \
cannam@90: asm ("madd16 %2,%3" \
cannam@90: : "+l" (lo), "+h" (hi) \
cannam@90: : "%r" ((x) >> 12), "r" ((y) >> 16))
cannam@90: # define MAD_F_MLZ(hi, lo) ((mad_fixed_t) (lo))
cannam@90: # endif
cannam@90:
cannam@90: # if defined(OPT_SPEED)
cannam@90: # define mad_f_scale64(hi, lo) \
cannam@90: ((mad_fixed_t) ((hi) << (32 - MAD_F_SCALEBITS)))
cannam@90: # define MAD_F_SCALEBITS MAD_F_FRACBITS
cannam@90: # endif
cannam@90:
cannam@90: /* --- SPARC --------------------------------------------------------------- */
cannam@90:
cannam@90: # elif defined(FPM_SPARC)
cannam@90:
cannam@90: /*
cannam@90: * This SPARC V8 version is fast and accurate; the disposition of the least
cannam@90: * significant bit depends on OPT_ACCURACY via mad_f_scale64().
cannam@90: */
cannam@90: # define MAD_F_MLX(hi, lo, x, y) \
cannam@90: asm ("smul %2, %3, %0\n\t" \
cannam@90: "rd %%y, %1" \
cannam@90: : "=r" (lo), "=r" (hi) \
cannam@90: : "%r" (x), "rI" (y))
cannam@90:
cannam@90: /* --- PowerPC ------------------------------------------------------------- */
cannam@90:
cannam@90: # elif defined(FPM_PPC)
cannam@90:
cannam@90: /*
cannam@90: * This PowerPC version is fast and accurate; the disposition of the least
cannam@90: * significant bit depends on OPT_ACCURACY via mad_f_scale64().
cannam@90: */
cannam@90: # define MAD_F_MLX(hi, lo, x, y) \
cannam@90: do { \
cannam@90: asm ("mullw %0,%1,%2" \
cannam@90: : "=r" (lo) \
cannam@90: : "%r" (x), "r" (y)); \
cannam@90: asm ("mulhw %0,%1,%2" \
cannam@90: : "=r" (hi) \
cannam@90: : "%r" (x), "r" (y)); \
cannam@90: } \
cannam@90: while (0)
cannam@90:
cannam@90: # if defined(OPT_ACCURACY)
cannam@90: /*
cannam@90: * This gives best accuracy but is not very fast.
cannam@90: */
cannam@90: # define MAD_F_MLA(hi, lo, x, y) \
cannam@90: ({ mad_fixed64hi_t __hi; \
cannam@90: mad_fixed64lo_t __lo; \
cannam@90: MAD_F_MLX(__hi, __lo, (x), (y)); \
cannam@90: asm ("addc %0,%2,%3\n\t" \
cannam@90: "adde %1,%4,%5" \
cannam@90: : "=r" (lo), "=r" (hi) \
cannam@90: : "%r" (lo), "r" (__lo), \
cannam@90: "%r" (hi), "r" (__hi) \
cannam@90: : "xer"); \
cannam@90: })
cannam@90: # endif
cannam@90:
cannam@90: # if defined(OPT_ACCURACY)
cannam@90: /*
cannam@90: * This is slower than the truncating version below it.
cannam@90: */
cannam@90: # define mad_f_scale64(hi, lo) \
cannam@90: ({ mad_fixed_t __result, __round; \
cannam@90: asm ("rotrwi %0,%1,%2" \
cannam@90: : "=r" (__result) \
cannam@90: : "r" (lo), "i" (MAD_F_SCALEBITS)); \
cannam@90: asm ("extrwi %0,%1,1,0" \
cannam@90: : "=r" (__round) \
cannam@90: : "r" (__result)); \
cannam@90: asm ("insrwi %0,%1,%2,0" \
cannam@90: : "+r" (__result) \
cannam@90: : "r" (hi), "i" (MAD_F_SCALEBITS)); \
cannam@90: asm ("add %0,%1,%2" \
cannam@90: : "=r" (__result) \
cannam@90: : "%r" (__result), "r" (__round)); \
cannam@90: __result; \
cannam@90: })
cannam@90: # else
cannam@90: # define mad_f_scale64(hi, lo) \
cannam@90: ({ mad_fixed_t __result; \
cannam@90: asm ("rotrwi %0,%1,%2" \
cannam@90: : "=r" (__result) \
cannam@90: : "r" (lo), "i" (MAD_F_SCALEBITS)); \
cannam@90: asm ("insrwi %0,%1,%2,0" \
cannam@90: : "+r" (__result) \
cannam@90: : "r" (hi), "i" (MAD_F_SCALEBITS)); \
cannam@90: __result; \
cannam@90: })
cannam@90: # endif
cannam@90:
cannam@90: # define MAD_F_SCALEBITS MAD_F_FRACBITS
cannam@90:
cannam@90: /* --- Default ------------------------------------------------------------- */
cannam@90:
cannam@90: # elif defined(FPM_DEFAULT)
cannam@90:
cannam@90: /*
cannam@90: * This version is the most portable but it loses significant accuracy.
cannam@90: * Furthermore, accuracy is biased against the second argument, so care
cannam@90: * should be taken when ordering operands.
cannam@90: *
cannam@90: * The scale factors are constant as this is not used with SSO.
cannam@90: *
cannam@90: * Pre-rounding is required to stay within the limits of compliance.
cannam@90: */
cannam@90: # if defined(OPT_SPEED)
cannam@90: # define mad_f_mul(x, y) (((x) >> 12) * ((y) >> 16))
cannam@90: # else
cannam@90: # define mad_f_mul(x, y) ((((x) + (1L << 11)) >> 12) * \
cannam@90: (((y) + (1L << 15)) >> 16))
cannam@90: # endif
cannam@90:
cannam@90: /* ------------------------------------------------------------------------- */
cannam@90:
cannam@90: # else
cannam@90: # error "no FPM selected"
cannam@90: # endif
cannam@90:
cannam@90: /* default implementations */
cannam@90:
cannam@90: # if !defined(mad_f_mul)
cannam@90: # define mad_f_mul(x, y) \
cannam@90: ({ register mad_fixed64hi_t __hi; \
cannam@90: register mad_fixed64lo_t __lo; \
cannam@90: MAD_F_MLX(__hi, __lo, (x), (y)); \
cannam@90: mad_f_scale64(__hi, __lo); \
cannam@90: })
cannam@90: # endif
cannam@90:
cannam@90: # if !defined(MAD_F_MLA)
cannam@90: # define MAD_F_ML0(hi, lo, x, y) ((lo) = mad_f_mul((x), (y)))
cannam@90: # define MAD_F_MLA(hi, lo, x, y) ((lo) += mad_f_mul((x), (y)))
cannam@90: # define MAD_F_MLN(hi, lo) ((lo) = -(lo))
cannam@90: # define MAD_F_MLZ(hi, lo) ((void) (hi), (mad_fixed_t) (lo))
cannam@90: # endif
cannam@90:
cannam@90: # if !defined(MAD_F_ML0)
cannam@90: # define MAD_F_ML0(hi, lo, x, y) MAD_F_MLX((hi), (lo), (x), (y))
cannam@90: # endif
cannam@90:
cannam@90: # if !defined(MAD_F_MLN)
cannam@90: # define MAD_F_MLN(hi, lo) ((hi) = ((lo) = -(lo)) ? ~(hi) : -(hi))
cannam@90: # endif
cannam@90:
cannam@90: # if !defined(MAD_F_MLZ)
cannam@90: # define MAD_F_MLZ(hi, lo) mad_f_scale64((hi), (lo))
cannam@90: # endif
cannam@90:
cannam@90: # if !defined(mad_f_scale64)
cannam@90: # if defined(OPT_ACCURACY)
cannam@90: # define mad_f_scale64(hi, lo) \
cannam@90: ((((mad_fixed_t) \
cannam@90: (((hi) << (32 - (MAD_F_SCALEBITS - 1))) | \
cannam@90: ((lo) >> (MAD_F_SCALEBITS - 1)))) + 1) >> 1)
cannam@90: # else
cannam@90: # define mad_f_scale64(hi, lo) \
cannam@90: ((mad_fixed_t) \
cannam@90: (((hi) << (32 - MAD_F_SCALEBITS)) | \
cannam@90: ((lo) >> MAD_F_SCALEBITS)))
cannam@90: # endif
cannam@90: # define MAD_F_SCALEBITS MAD_F_FRACBITS
cannam@90: # endif
cannam@90:
cannam@90: /* C routines */
cannam@90:
cannam@90: mad_fixed_t mad_f_abs(mad_fixed_t);
cannam@90: mad_fixed_t mad_f_div(mad_fixed_t, mad_fixed_t);
cannam@90:
cannam@90: # endif
cannam@90:
cannam@90: /* Id: bit.h,v 1.12 2004/01/23 09:41:32 rob Exp */
cannam@90:
cannam@90: # ifndef LIBMAD_BIT_H
cannam@90: # define LIBMAD_BIT_H
cannam@90:
cannam@90: struct mad_bitptr {
cannam@90: unsigned char const *byte;
cannam@90: unsigned short cache;
cannam@90: unsigned short left;
cannam@90: };
cannam@90:
cannam@90: void mad_bit_init(struct mad_bitptr *, unsigned char const *);
cannam@90:
cannam@90: # define mad_bit_finish(bitptr) /* nothing */
cannam@90:
cannam@90: unsigned int mad_bit_length(struct mad_bitptr const *,
cannam@90: struct mad_bitptr const *);
cannam@90:
cannam@90: # define mad_bit_bitsleft(bitptr) ((bitptr)->left)
cannam@90: unsigned char const *mad_bit_nextbyte(struct mad_bitptr const *);
cannam@90:
cannam@90: void mad_bit_skip(struct mad_bitptr *, unsigned int);
cannam@90: unsigned long mad_bit_read(struct mad_bitptr *, unsigned int);
cannam@90: void mad_bit_write(struct mad_bitptr *, unsigned int, unsigned long);
cannam@90:
cannam@90: unsigned short mad_bit_crc(struct mad_bitptr, unsigned int, unsigned short);
cannam@90:
cannam@90: # endif
cannam@90:
cannam@90: /* Id: timer.h,v 1.16 2004/01/23 09:41:33 rob Exp */
cannam@90:
cannam@90: # ifndef LIBMAD_TIMER_H
cannam@90: # define LIBMAD_TIMER_H
cannam@90:
cannam@90: typedef struct {
cannam@90: signed long seconds; /* whole seconds */
cannam@90: unsigned long fraction; /* 1/MAD_TIMER_RESOLUTION seconds */
cannam@90: } mad_timer_t;
cannam@90:
cannam@90: extern mad_timer_t const mad_timer_zero;
cannam@90:
cannam@90: # define MAD_TIMER_RESOLUTION 352800000UL
cannam@90:
cannam@90: enum mad_units {
cannam@90: MAD_UNITS_HOURS = -2,
cannam@90: MAD_UNITS_MINUTES = -1,
cannam@90: MAD_UNITS_SECONDS = 0,
cannam@90:
cannam@90: /* metric units */
cannam@90:
cannam@90: MAD_UNITS_DECISECONDS = 10,
cannam@90: MAD_UNITS_CENTISECONDS = 100,
cannam@90: MAD_UNITS_MILLISECONDS = 1000,
cannam@90:
cannam@90: /* audio sample units */
cannam@90:
cannam@90: MAD_UNITS_8000_HZ = 8000,
cannam@90: MAD_UNITS_11025_HZ = 11025,
cannam@90: MAD_UNITS_12000_HZ = 12000,
cannam@90:
cannam@90: MAD_UNITS_16000_HZ = 16000,
cannam@90: MAD_UNITS_22050_HZ = 22050,
cannam@90: MAD_UNITS_24000_HZ = 24000,
cannam@90:
cannam@90: MAD_UNITS_32000_HZ = 32000,
cannam@90: MAD_UNITS_44100_HZ = 44100,
cannam@90: MAD_UNITS_48000_HZ = 48000,
cannam@90:
cannam@90: /* video frame/field units */
cannam@90:
cannam@90: MAD_UNITS_24_FPS = 24,
cannam@90: MAD_UNITS_25_FPS = 25,
cannam@90: MAD_UNITS_30_FPS = 30,
cannam@90: MAD_UNITS_48_FPS = 48,
cannam@90: MAD_UNITS_50_FPS = 50,
cannam@90: MAD_UNITS_60_FPS = 60,
cannam@90:
cannam@90: /* CD audio frames */
cannam@90:
cannam@90: MAD_UNITS_75_FPS = 75,
cannam@90:
cannam@90: /* video drop-frame units */
cannam@90:
cannam@90: MAD_UNITS_23_976_FPS = -24,
cannam@90: MAD_UNITS_24_975_FPS = -25,
cannam@90: MAD_UNITS_29_97_FPS = -30,
cannam@90: MAD_UNITS_47_952_FPS = -48,
cannam@90: MAD_UNITS_49_95_FPS = -50,
cannam@90: MAD_UNITS_59_94_FPS = -60
cannam@90: };
cannam@90:
cannam@90: # define mad_timer_reset(timer) ((void) (*(timer) = mad_timer_zero))
cannam@90:
cannam@90: int mad_timer_compare(mad_timer_t, mad_timer_t);
cannam@90:
cannam@90: # define mad_timer_sign(timer) mad_timer_compare((timer), mad_timer_zero)
cannam@90:
cannam@90: void mad_timer_negate(mad_timer_t *);
cannam@90: mad_timer_t mad_timer_abs(mad_timer_t);
cannam@90:
cannam@90: void mad_timer_set(mad_timer_t *, unsigned long, unsigned long, unsigned long);
cannam@90: void mad_timer_add(mad_timer_t *, mad_timer_t);
cannam@90: void mad_timer_multiply(mad_timer_t *, signed long);
cannam@90:
cannam@90: signed long mad_timer_count(mad_timer_t, enum mad_units);
cannam@90: unsigned long mad_timer_fraction(mad_timer_t, unsigned long);
cannam@90: void mad_timer_string(mad_timer_t, char *, char const *,
cannam@90: enum mad_units, enum mad_units, unsigned long);
cannam@90:
cannam@90: # endif
cannam@90:
cannam@90: /* Id: stream.h,v 1.20 2004/02/05 09:02:39 rob Exp */
cannam@90:
cannam@90: # ifndef LIBMAD_STREAM_H
cannam@90: # define LIBMAD_STREAM_H
cannam@90:
cannam@90:
cannam@90: # define MAD_BUFFER_GUARD 8
cannam@90: # define MAD_BUFFER_MDLEN (511 + 2048 + MAD_BUFFER_GUARD)
cannam@90:
cannam@90: enum mad_error {
cannam@90: MAD_ERROR_NONE = 0x0000, /* no error */
cannam@90:
cannam@90: MAD_ERROR_BUFLEN = 0x0001, /* input buffer too small (or EOF) */
cannam@90: MAD_ERROR_BUFPTR = 0x0002, /* invalid (null) buffer pointer */
cannam@90:
cannam@90: MAD_ERROR_NOMEM = 0x0031, /* not enough memory */
cannam@90:
cannam@90: MAD_ERROR_LOSTSYNC = 0x0101, /* lost synchronization */
cannam@90: MAD_ERROR_BADLAYER = 0x0102, /* reserved header layer value */
cannam@90: MAD_ERROR_BADBITRATE = 0x0103, /* forbidden bitrate value */
cannam@90: MAD_ERROR_BADSAMPLERATE = 0x0104, /* reserved sample frequency value */
cannam@90: MAD_ERROR_BADEMPHASIS = 0x0105, /* reserved emphasis value */
cannam@90:
cannam@90: MAD_ERROR_BADCRC = 0x0201, /* CRC check failed */
cannam@90: MAD_ERROR_BADBITALLOC = 0x0211, /* forbidden bit allocation value */
cannam@90: MAD_ERROR_BADSCALEFACTOR = 0x0221, /* bad scalefactor index */
cannam@90: MAD_ERROR_BADMODE = 0x0222, /* bad bitrate/mode combination */
cannam@90: MAD_ERROR_BADFRAMELEN = 0x0231, /* bad frame length */
cannam@90: MAD_ERROR_BADBIGVALUES = 0x0232, /* bad big_values count */
cannam@90: MAD_ERROR_BADBLOCKTYPE = 0x0233, /* reserved block_type */
cannam@90: MAD_ERROR_BADSCFSI = 0x0234, /* bad scalefactor selection info */
cannam@90: MAD_ERROR_BADDATAPTR = 0x0235, /* bad main_data_begin pointer */
cannam@90: MAD_ERROR_BADPART3LEN = 0x0236, /* bad audio data length */
cannam@90: MAD_ERROR_BADHUFFTABLE = 0x0237, /* bad Huffman table select */
cannam@90: MAD_ERROR_BADHUFFDATA = 0x0238, /* Huffman data overrun */
cannam@90: MAD_ERROR_BADSTEREO = 0x0239 /* incompatible block_type for JS */
cannam@90: };
cannam@90:
cannam@90: # define MAD_RECOVERABLE(error) ((error) & 0xff00)
cannam@90:
cannam@90: struct mad_stream {
cannam@90: unsigned char const *buffer; /* input bitstream buffer */
cannam@90: unsigned char const *bufend; /* end of buffer */
cannam@90: unsigned long skiplen; /* bytes to skip before next frame */
cannam@90:
cannam@90: int sync; /* stream sync found */
cannam@90: unsigned long freerate; /* free bitrate (fixed) */
cannam@90:
cannam@90: unsigned char const *this_frame; /* start of current frame */
cannam@90: unsigned char const *next_frame; /* start of next frame */
cannam@90: struct mad_bitptr ptr; /* current processing bit pointer */
cannam@90:
cannam@90: struct mad_bitptr anc_ptr; /* ancillary bits pointer */
cannam@90: unsigned int anc_bitlen; /* number of ancillary bits */
cannam@90:
cannam@90: unsigned char (*main_data)[MAD_BUFFER_MDLEN];
cannam@90: /* Layer III main_data() */
cannam@90: unsigned int md_len; /* bytes in main_data */
cannam@90:
cannam@90: int options; /* decoding options (see below) */
cannam@90: enum mad_error error; /* error code (see above) */
cannam@90: };
cannam@90:
cannam@90: enum {
cannam@90: MAD_OPTION_IGNORECRC = 0x0001, /* ignore CRC errors */
cannam@90: MAD_OPTION_HALFSAMPLERATE = 0x0002 /* generate PCM at 1/2 sample rate */
cannam@90: # if 0 /* not yet implemented */
cannam@90: MAD_OPTION_LEFTCHANNEL = 0x0010, /* decode left channel only */
cannam@90: MAD_OPTION_RIGHTCHANNEL = 0x0020, /* decode right channel only */
cannam@90: MAD_OPTION_SINGLECHANNEL = 0x0030 /* combine channels */
cannam@90: # endif
cannam@90: };
cannam@90:
cannam@90: void mad_stream_init(struct mad_stream *);
cannam@90: void mad_stream_finish(struct mad_stream *);
cannam@90:
cannam@90: # define mad_stream_options(stream, opts) \
cannam@90: ((void) ((stream)->options = (opts)))
cannam@90:
cannam@90: void mad_stream_buffer(struct mad_stream *,
cannam@90: unsigned char const *, unsigned long);
cannam@90: void mad_stream_skip(struct mad_stream *, unsigned long);
cannam@90:
cannam@90: int mad_stream_sync(struct mad_stream *);
cannam@90:
cannam@90: char const *mad_stream_errorstr(struct mad_stream const *);
cannam@90:
cannam@90: # endif
cannam@90:
cannam@90: /* Id: frame.h,v 1.20 2004/01/23 09:41:32 rob Exp */
cannam@90:
cannam@90: # ifndef LIBMAD_FRAME_H
cannam@90: # define LIBMAD_FRAME_H
cannam@90:
cannam@90:
cannam@90: enum mad_layer {
cannam@90: MAD_LAYER_I = 1, /* Layer I */
cannam@90: MAD_LAYER_II = 2, /* Layer II */
cannam@90: MAD_LAYER_III = 3 /* Layer III */
cannam@90: };
cannam@90:
cannam@90: enum mad_mode {
cannam@90: MAD_MODE_SINGLE_CHANNEL = 0, /* single channel */
cannam@90: MAD_MODE_DUAL_CHANNEL = 1, /* dual channel */
cannam@90: MAD_MODE_JOINT_STEREO = 2, /* joint (MS/intensity) stereo */
cannam@90: MAD_MODE_STEREO = 3 /* normal LR stereo */
cannam@90: };
cannam@90:
cannam@90: enum mad_emphasis {
cannam@90: MAD_EMPHASIS_NONE = 0, /* no emphasis */
cannam@90: MAD_EMPHASIS_50_15_US = 1, /* 50/15 microseconds emphasis */
cannam@90: MAD_EMPHASIS_CCITT_J_17 = 3, /* CCITT J.17 emphasis */
cannam@90: MAD_EMPHASIS_RESERVED = 2 /* unknown emphasis */
cannam@90: };
cannam@90:
cannam@90: struct mad_header {
cannam@90: enum mad_layer layer; /* audio layer (1, 2, or 3) */
cannam@90: enum mad_mode mode; /* channel mode (see above) */
cannam@90: int mode_extension; /* additional mode info */
cannam@90: enum mad_emphasis emphasis; /* de-emphasis to use (see above) */
cannam@90:
cannam@90: unsigned long bitrate; /* stream bitrate (bps) */
cannam@90: unsigned int samplerate; /* sampling frequency (Hz) */
cannam@90:
cannam@90: unsigned short crc_check; /* frame CRC accumulator */
cannam@90: unsigned short crc_target; /* final target CRC checksum */
cannam@90:
cannam@90: int flags; /* flags (see below) */
cannam@90: int private_bits; /* private bits (see below) */
cannam@90:
cannam@90: mad_timer_t duration; /* audio playing time of frame */
cannam@90: };
cannam@90:
cannam@90: struct mad_frame {
cannam@90: struct mad_header header; /* MPEG audio header */
cannam@90:
cannam@90: int options; /* decoding options (from stream) */
cannam@90:
cannam@90: mad_fixed_t sbsample[2][36][32]; /* synthesis subband filter samples */
cannam@90: mad_fixed_t (*overlap)[2][32][18]; /* Layer III block overlap data */
cannam@90: };
cannam@90:
cannam@90: # define MAD_NCHANNELS(header) ((header)->mode ? 2 : 1)
cannam@90: # define MAD_NSBSAMPLES(header) \
cannam@90: ((header)->layer == MAD_LAYER_I ? 12 : \
cannam@90: (((header)->layer == MAD_LAYER_III && \
cannam@90: ((header)->flags & MAD_FLAG_LSF_EXT)) ? 18 : 36))
cannam@90:
cannam@90: enum {
cannam@90: MAD_FLAG_NPRIVATE_III = 0x0007, /* number of Layer III private bits */
cannam@90: MAD_FLAG_INCOMPLETE = 0x0008, /* header but not data is decoded */
cannam@90:
cannam@90: MAD_FLAG_PROTECTION = 0x0010, /* frame has CRC protection */
cannam@90: MAD_FLAG_COPYRIGHT = 0x0020, /* frame is copyright */
cannam@90: MAD_FLAG_ORIGINAL = 0x0040, /* frame is original (else copy) */
cannam@90: MAD_FLAG_PADDING = 0x0080, /* frame has additional slot */
cannam@90:
cannam@90: MAD_FLAG_I_STEREO = 0x0100, /* uses intensity joint stereo */
cannam@90: MAD_FLAG_MS_STEREO = 0x0200, /* uses middle/side joint stereo */
cannam@90: MAD_FLAG_FREEFORMAT = 0x0400, /* uses free format bitrate */
cannam@90:
cannam@90: MAD_FLAG_LSF_EXT = 0x1000, /* lower sampling freq. extension */
cannam@90: MAD_FLAG_MC_EXT = 0x2000, /* multichannel audio extension */
cannam@90: MAD_FLAG_MPEG_2_5_EXT = 0x4000 /* MPEG 2.5 (unofficial) extension */
cannam@90: };
cannam@90:
cannam@90: enum {
cannam@90: MAD_PRIVATE_HEADER = 0x0100, /* header private bit */
cannam@90: MAD_PRIVATE_III = 0x001f /* Layer III private bits (up to 5) */
cannam@90: };
cannam@90:
cannam@90: void mad_header_init(struct mad_header *);
cannam@90:
cannam@90: # define mad_header_finish(header) /* nothing */
cannam@90:
cannam@90: int mad_header_decode(struct mad_header *, struct mad_stream *);
cannam@90:
cannam@90: void mad_frame_init(struct mad_frame *);
cannam@90: void mad_frame_finish(struct mad_frame *);
cannam@90:
cannam@90: int mad_frame_decode(struct mad_frame *, struct mad_stream *);
cannam@90:
cannam@90: void mad_frame_mute(struct mad_frame *);
cannam@90:
cannam@90: # endif
cannam@90:
cannam@90: /* Id: synth.h,v 1.15 2004/01/23 09:41:33 rob Exp */
cannam@90:
cannam@90: # ifndef LIBMAD_SYNTH_H
cannam@90: # define LIBMAD_SYNTH_H
cannam@90:
cannam@90:
cannam@90: struct mad_pcm {
cannam@90: unsigned int samplerate; /* sampling frequency (Hz) */
cannam@90: unsigned short channels; /* number of channels */
cannam@90: unsigned short length; /* number of samples per channel */
cannam@90: mad_fixed_t samples[2][1152]; /* PCM output samples [ch][sample] */
cannam@90: };
cannam@90:
cannam@90: struct mad_synth {
cannam@90: mad_fixed_t filter[2][2][2][16][8]; /* polyphase filterbank outputs */
cannam@90: /* [ch][eo][peo][s][v] */
cannam@90:
cannam@90: unsigned int phase; /* current processing phase */
cannam@90:
cannam@90: struct mad_pcm pcm; /* PCM output */
cannam@90: };
cannam@90:
cannam@90: /* single channel PCM selector */
cannam@90: enum {
cannam@90: MAD_PCM_CHANNEL_SINGLE = 0
cannam@90: };
cannam@90:
cannam@90: /* dual channel PCM selector */
cannam@90: enum {
cannam@90: MAD_PCM_CHANNEL_DUAL_1 = 0,
cannam@90: MAD_PCM_CHANNEL_DUAL_2 = 1
cannam@90: };
cannam@90:
cannam@90: /* stereo PCM selector */
cannam@90: enum {
cannam@90: MAD_PCM_CHANNEL_STEREO_LEFT = 0,
cannam@90: MAD_PCM_CHANNEL_STEREO_RIGHT = 1
cannam@90: };
cannam@90:
cannam@90: void mad_synth_init(struct mad_synth *);
cannam@90:
cannam@90: # define mad_synth_finish(synth) /* nothing */
cannam@90:
cannam@90: void mad_synth_mute(struct mad_synth *);
cannam@90:
cannam@90: void mad_synth_frame(struct mad_synth *, struct mad_frame const *);
cannam@90:
cannam@90: # endif
cannam@90:
cannam@90: /* Id: decoder.h,v 1.17 2004/01/23 09:41:32 rob Exp */
cannam@90:
cannam@90: # ifndef LIBMAD_DECODER_H
cannam@90: # define LIBMAD_DECODER_H
cannam@90:
cannam@90:
cannam@90: enum mad_decoder_mode {
cannam@90: MAD_DECODER_MODE_SYNC = 0,
cannam@90: MAD_DECODER_MODE_ASYNC
cannam@90: };
cannam@90:
cannam@90: enum mad_flow {
cannam@90: MAD_FLOW_CONTINUE = 0x0000, /* continue normally */
cannam@90: MAD_FLOW_STOP = 0x0010, /* stop decoding normally */
cannam@90: MAD_FLOW_BREAK = 0x0011, /* stop decoding and signal an error */
cannam@90: MAD_FLOW_IGNORE = 0x0020 /* ignore the current frame */
cannam@90: };
cannam@90:
cannam@90: struct mad_decoder {
cannam@90: enum mad_decoder_mode mode;
cannam@90:
cannam@90: int options;
cannam@90:
cannam@90: struct {
cannam@90: long pid;
cannam@90: int in;
cannam@90: int out;
cannam@90: } async;
cannam@90:
cannam@90: struct {
cannam@90: struct mad_stream stream;
cannam@90: struct mad_frame frame;
cannam@90: struct mad_synth synth;
cannam@90: } *sync;
cannam@90:
cannam@90: void *cb_data;
cannam@90:
cannam@90: enum mad_flow (*input_func)(void *, struct mad_stream *);
cannam@90: enum mad_flow (*header_func)(void *, struct mad_header const *);
cannam@90: enum mad_flow (*filter_func)(void *,
cannam@90: struct mad_stream const *, struct mad_frame *);
cannam@90: enum mad_flow (*output_func)(void *,
cannam@90: struct mad_header const *, struct mad_pcm *);
cannam@90: enum mad_flow (*error_func)(void *, struct mad_stream *, struct mad_frame *);
cannam@90: enum mad_flow (*message_func)(void *, void *, unsigned int *);
cannam@90: };
cannam@90:
cannam@90: void mad_decoder_init(struct mad_decoder *, void *,
cannam@90: enum mad_flow (*)(void *, struct mad_stream *),
cannam@90: enum mad_flow (*)(void *, struct mad_header const *),
cannam@90: enum mad_flow (*)(void *,
cannam@90: struct mad_stream const *,
cannam@90: struct mad_frame *),
cannam@90: enum mad_flow (*)(void *,
cannam@90: struct mad_header const *,
cannam@90: struct mad_pcm *),
cannam@90: enum mad_flow (*)(void *,
cannam@90: struct mad_stream *,
cannam@90: struct mad_frame *),
cannam@90: enum mad_flow (*)(void *, void *, unsigned int *));
cannam@90: int mad_decoder_finish(struct mad_decoder *);
cannam@90:
cannam@90: # define mad_decoder_options(decoder, opts) \
cannam@90: ((void) ((decoder)->options = (opts)))
cannam@90:
cannam@90: int mad_decoder_run(struct mad_decoder *, enum mad_decoder_mode);
cannam@90: int mad_decoder_message(struct mad_decoder *, void *, unsigned int *);
cannam@90:
cannam@90: # endif
cannam@90:
cannam@90: # ifdef __cplusplus
cannam@90: }
cannam@90: # endif