1 /*
2  * libmad - MPEG audio decoder library
3  * Copyright (C) 2000-2004 Underbit Technologies, Inc.
4  *
5  * This program is free software; you can redistribute it and/or modify
6  * it under the terms of the GNU General Public License as published by
7  * the Free Software Foundation; either version 2 of the License, or
8  * (at your option) any later version.
9  *
10  * This program is distributed in the hope that it will be useful,
11  * but WITHOUT ANY WARRANTY; without even the implied warranty of
12  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
13  * GNU General Public License for more details.
14  *
15  * You should have received a copy of the GNU General Public License
16  * along with this program; if not, write to the Free Software
17  * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
18  *
19  * If you would like to negotiate alternate licensing terms, you may do
20  * so by contacting: Underbit Technologies, Inc. <info@underbit.com>
21  */
22 
23 #ifdef __cplusplus
24 extern "C" {
25 #endif
26 
27 /* configuration for windows */
28 #if defined(_WIN64)
29 # define FPM_64BIT
30 #elif defined(_M_IX86) || defined(__i386__) || defined(__386__)
31 # define FPM_INTEL
32 #else
33 # define FPM_DEFAULT
34 #endif
35 
36 /* Id: version.h,v 1.26 2004/01/23 09:41:33 rob Exp */
37 
38 # ifndef LIBMAD_VERSION_H
39 # define LIBMAD_VERSION_H
40 
41 # define MAD_VERSION_MAJOR	0
42 # define MAD_VERSION_MINOR	15
43 # define MAD_VERSION_PATCH	1
44 # define MAD_VERSION_EXTRA	" (beta)"
45 
46 # define MAD_VERSION_STRINGIZE(str)	#str
47 # define MAD_VERSION_STRING(num)	MAD_VERSION_STRINGIZE(num)
48 
49 # define MAD_VERSION		MAD_VERSION_STRING(MAD_VERSION_MAJOR) "."  \
50 				MAD_VERSION_STRING(MAD_VERSION_MINOR) "."  \
51 				MAD_VERSION_STRING(MAD_VERSION_PATCH)  \
52 				MAD_VERSION_EXTRA
53 
54 # define MAD_PUBLISHYEAR	"2000-2004"
55 # define MAD_AUTHOR		"Underbit Technologies, Inc."
56 # define MAD_EMAIL		"info@underbit.com"
57 
58 extern char const mad_version[];
59 extern char const mad_copyright[];
60 extern char const mad_author[];
61 extern char const mad_build[];
62 
63 # endif
64 
65 /* Id: fixed.h,v 1.38 2004/02/17 02:02:03 rob Exp */
66 
67 # ifndef LIBMAD_FIXED_H
68 # define LIBMAD_FIXED_H
69 
70 #  include <limits.h>
71 
72 # if INT_MAX >= 2147483647
73 typedef   signed int mad_fixed_t;
74 
75 typedef   signed int mad_fixed64hi_t;
76 typedef unsigned int mad_fixed64lo_t;
77 # else
78 typedef   signed long mad_fixed_t;
79 
80 typedef   signed long mad_fixed64hi_t;
81 typedef unsigned long mad_fixed64lo_t;
82 # endif
83 /* compile-time assert: */
84 typedef int _mad_check_fixed_t[2*(sizeof(mad_fixed_t)>=4) - 1];
85 
86 # if defined(_MSC_VER)
87 #  define mad_fixed64_t  signed __int64
88 # elif 1 || defined(__GNUC__)
89 #  define mad_fixed64_t  signed long long
90 # endif
91 
92 # if defined(FPM_FLOAT)
93 typedef double mad_sample_t;
94 # else
95 typedef mad_fixed_t mad_sample_t;
96 # endif
97 
98 /*
99  * Fixed-point format: 0xABBBBBBB
100  * A == whole part      (sign + 3 bits)
101  * B == fractional part (28 bits)
102  *
103  * Values are signed two's complement, so the effective range is:
104  * 0x80000000 to 0x7fffffff
105  *       -8.0 to +7.9999999962747097015380859375
106  *
107  * The smallest representable value is:
108  * 0x00000001 == 0.0000000037252902984619140625 (i.e. about 3.725e-9)
109  *
110  * 28 bits of fractional accuracy represent about
111  * 8.6 digits of decimal accuracy.
112  *
113  * Fixed-point numbers can be added or subtracted as normal
114  * integers, but multiplication requires shifting the 64-bit result
115  * from 56 fractional bits back to 28 (and rounding.)
116  *
117  * Changing the definition of MAD_F_FRACBITS is only partially
118  * supported, and must be done with care.
119  */
120 
121 # define MAD_F_FRACBITS		28
122 
123 # if MAD_F_FRACBITS == 28
124 #  define MAD_F(x)		((mad_fixed_t) (x##L))
125 # else
126 #  if MAD_F_FRACBITS < 28
127 #   warning "MAD_F_FRACBITS < 28"
128 #   define MAD_F(x)		((mad_fixed_t)  \
129 				 (((x##L) +  \
130 				   (1L << (28 - MAD_F_FRACBITS - 1))) >>  \
131 				  (28 - MAD_F_FRACBITS)))
132 #  elif MAD_F_FRACBITS > 28
133 #   error "MAD_F_FRACBITS > 28 not currently supported"
134 #   define MAD_F(x)		((mad_fixed_t)  \
135 				 ((x##L) << (MAD_F_FRACBITS - 28)))
136 #  endif
137 # endif
138 
139 # define MAD_F_MIN		((mad_fixed_t) -0x80000000L)
140 # define MAD_F_MAX		((mad_fixed_t) +0x7fffffffL)
141 
142 # define MAD_F_ONE		MAD_F(0x10000000)
143 
144 # define mad_f_tofixed(x)	((mad_fixed_t)  \
145 				 ((x) * (double) (1L << MAD_F_FRACBITS) + 0.5))
146 # define mad_f_todouble(x)	((double)  \
147 				 ((x) / (double) (1L << MAD_F_FRACBITS)))
148 
149 # define mad_f_intpart(x)	((x) >> MAD_F_FRACBITS)
150 # define mad_f_fracpart(x)	((x) & ((1L << MAD_F_FRACBITS) - 1))
151 				/* (x should be positive) */
152 
153 # define mad_f_fromint(x)	((x) << MAD_F_FRACBITS)
154 
155 # define mad_f_add(x, y)	((x) + (y))
156 # define mad_f_sub(x, y)	((x) - (y))
157 
158 # if defined(FPM_FLOAT)
159 #  error "FPM_FLOAT not yet supported"
160 
161 #  undef MAD_F
162 #  define MAD_F(x)		mad_f_todouble(x)
163 
164 #  define mad_f_mul(x, y)	((x) * (y))
165 #  define mad_f_scale64
166 
167 #  undef ASO_ZEROCHECK
168 
169 # elif defined(FPM_64BIT)
170 
171 /*
172  * This version should be the most accurate if 64-bit types are supported by
173  * the compiler, although it may not be the most efficient.
174  */
175 #  if defined(OPT_ACCURACY)
176 #   define mad_f_mul(x, y)  \
177     ((mad_fixed_t)  \
178      ((((mad_fixed64_t) (x) * (y)) +  \
179        (1L << (MAD_F_SCALEBITS - 1))) >> MAD_F_SCALEBITS))
180 #  else
181 #   define mad_f_mul(x, y)  \
182     ((mad_fixed_t) (((mad_fixed64_t) (x) * (y)) >> MAD_F_SCALEBITS))
183 #  endif
184 
185 #  define MAD_F_SCALEBITS  MAD_F_FRACBITS
186 
187 /* --- Intel --------------------------------------------------------------- */
188 
189 # elif defined(FPM_INTEL)
190 
191 #  if defined(_MSC_VER)
192 #   pragma warning(push)
193 #   pragma warning(disable: 4035)  /* no return value */
194 static __forceinline
mad_f_mul_inline(mad_fixed_t x,mad_fixed_t y)195 mad_fixed_t mad_f_mul_inline(mad_fixed_t x, mad_fixed_t y)
196 {
197   enum {
198     fracbits = MAD_F_FRACBITS
199   };
200 
201   __asm {
202     mov eax, x
203     imul y
204     shrd eax, edx, fracbits
205   }
206 
207   /* implicit return of eax */
208 }
209 #   pragma warning(pop)
210 
211 #   define mad_f_mul		mad_f_mul_inline
212 #   define mad_f_scale64
213 
214 #  elif defined(__WATCOMC__) && defined(__386__)
215 mad_fixed_t mad_f_mul_inl(mad_fixed_t,mad_fixed_t);
216 /* 28 == MAD_F_FRACBITS */
217 #pragma aux mad_f_mul_inl =		\
218 	"imul ebx",			\
219 	"shrd eax,edx,28"		\
220 	parm	[eax] [ebx]		\
221 	value	[eax]			\
222 	modify exact [eax edx]
223 
224 #   define mad_f_mul		mad_f_mul_inl
225 #   define mad_f_scale64
226 
227 #  else
228 /*
229  * This Intel version is fast and accurate; the disposition of the least
230  * significant bit depends on OPT_ACCURACY via mad_f_scale64().
231  */
232 #   define MAD_F_MLX(hi, lo, x, y)  \
233     asm ("imull %3"  \
234 	 : "=a" (lo), "=d" (hi)  \
235 	 : "%a" (x), "rm" (y)  \
236 	 : "cc")
237 
238 #   if defined(OPT_ACCURACY)
239 /*
240  * This gives best accuracy but is not very fast.
241  */
242 #    define MAD_F_MLA(hi, lo, x, y)  \
243     ({ mad_fixed64hi_t __hi;  \
244        mad_fixed64lo_t __lo;  \
245        MAD_F_MLX(__hi, __lo, (x), (y));  \
246        asm ("addl %2,%0\n\t"  \
247 	    "adcl %3,%1"  \
248 	    : "=rm" (lo), "=rm" (hi)  \
249 	    : "r" (__lo), "r" (__hi), "0" (lo), "1" (hi)  \
250 	    : "cc");  \
251     })
252 #   endif  /* OPT_ACCURACY */
253 
254 #   if defined(OPT_ACCURACY)
255 /*
256  * Surprisingly, this is faster than SHRD followed by ADC.
257  */
258 #    define mad_f_scale64(hi, lo)  \
259     ({ mad_fixed64hi_t __hi_;  \
260        mad_fixed64lo_t __lo_;  \
261        mad_fixed_t __result;  \
262        asm ("addl %4,%2\n\t"  \
263 	    "adcl %5,%3"  \
264 	    : "=rm" (__lo_), "=rm" (__hi_)  \
265 	    : "0" (lo), "1" (hi),  \
266 	      "ir" (1L << (MAD_F_SCALEBITS - 1)), "ir" (0)  \
267 	    : "cc");  \
268        asm ("shrdl %3,%2,%1"  \
269 	    : "=rm" (__result)  \
270 	    : "0" (__lo_), "r" (__hi_), "I" (MAD_F_SCALEBITS)  \
271 	    : "cc");  \
272        __result;  \
273     })
274 #   elif defined(OPT_INTEL)
275 /*
276  * Alternate Intel scaling that may or may not perform better.
277  */
278 #    define mad_f_scale64(hi, lo)  \
279     ({ mad_fixed_t __result;  \
280        asm ("shrl %3,%1\n\t"  \
281 	    "shll %4,%2\n\t"  \
282 	    "orl %2,%1"  \
283 	    : "=rm" (__result)  \
284 	    : "0" (lo), "r" (hi),  \
285 	      "I" (MAD_F_SCALEBITS), "I" (32 - MAD_F_SCALEBITS)  \
286 	    : "cc");  \
287        __result;  \
288     })
289 #   else
290 #    define mad_f_scale64(hi, lo)  \
291     ({ mad_fixed_t __result;  \
292        asm ("shrdl %3,%2,%1"  \
293 	    : "=rm" (__result)  \
294 	    : "0" (lo), "r" (hi), "I" (MAD_F_SCALEBITS)  \
295 	    : "cc");  \
296        __result;  \
297     })
298 #   endif  /* OPT_ACCURACY */
299 
300 #   define MAD_F_SCALEBITS  MAD_F_FRACBITS
301 #  endif
302 
303 /* --- ARM ----------------------------------------------------------------- */
304 
305 # elif defined(FPM_ARM)
306 
307 /*
308  * This ARM V4 version is as accurate as FPM_64BIT but much faster. The
309  * least significant bit is properly rounded at no CPU cycle cost!
310  */
311 # if 1
312 /*
313  * This is faster than the default implementation via MAD_F_MLX() and
314  * mad_f_scale64().
315  */
316 #  define mad_f_mul(x, y)  \
317     ({ mad_fixed64hi_t __hi;  \
318        mad_fixed64lo_t __lo;  \
319        mad_fixed_t __result;  \
320        asm ("smull	%0, %1, %3, %4\n\t"  \
321 	    "movs	%0, %0, lsr %5\n\t"  \
322 	    "adc	%2, %0, %1, lsl %6"  \
323 	    : "=&r" (__lo), "=&r" (__hi), "=r" (__result)  \
324 	    : "%r" (x), "r" (y),  \
325 	      "M" (MAD_F_SCALEBITS), "M" (32 - MAD_F_SCALEBITS)  \
326 	    : "cc");  \
327        __result;  \
328     })
329 # endif
330 
331 #  define MAD_F_MLX(hi, lo, x, y)  \
332     asm ("smull	%0, %1, %2, %3"  \
333 	 : "=&r" (lo), "=&r" (hi)  \
334 	 : "%r" (x), "r" (y))
335 
336 #  define MAD_F_MLA(hi, lo, x, y)  \
337     asm ("smlal	%0, %1, %2, %3"  \
338 	 : "+r" (lo), "+r" (hi)  \
339 	 : "%r" (x), "r" (y))
340 
341 #ifdef __thumb__
342 /* In Thumb-2, the RSB-immediate instruction is only allowed with a zero
343 	operand.  If needed this code can also support Thumb-1
344 	(simply append "s" to the end of the second two instructions). */
345 #  define MAD_F_MLN(hi, lo)  \
346     asm ("rsbs        %0, %0, #0\n\t"  \
347 	 "sbc   %1, %1, %1\n\t"  \
348 	 "sub   %1, %1, %2"  \
349 	 : "+&r" (lo), "=&r" (hi)  \
350 	 : "r" (hi)  \
351 	 : "cc")
352 #else /* ! __thumb__ */
353 #  define MAD_F_MLN(hi, lo)  \
354     asm ("rsbs	%0, %2, #0\n\t"  \
355 	 "rsc	%1, %3, #0"  \
356 	 : "=&r" (lo), "=r" (hi)  \
357 	 : "0" (lo), "1" (hi)  \
358 	 : "cc")
359 #endif /* __thumb__ */
360 
361 #  define mad_f_scale64(hi, lo)  \
362     ({ mad_fixed_t __result;  \
363        asm ("movs	%0, %1, lsr %3\n\t"  \
364 	    "adc	%0, %0, %2, lsl %4"  \
365 	    : "=&r" (__result)  \
366 	    : "r" (lo), "r" (hi),  \
367 	      "M" (MAD_F_SCALEBITS), "M" (32 - MAD_F_SCALEBITS)  \
368 	    : "cc");  \
369        __result;  \
370     })
371 
372 #  define MAD_F_SCALEBITS  MAD_F_FRACBITS
373 
374 /* --- MIPS ---------------------------------------------------------------- */
375 
376 # elif defined(FPM_MIPS)
377 
378 #if defined (__GNUC__) && (__GNUC__ > 4 || (__GNUC__ == 4 && __GNUC_MINOR__ >= 4))
379   typedef unsigned int u64_di_t __attribute__ ((mode (DI)));
380 # define MAD_F_MLX(hi, lo, x, y) \
381    do { \
382       u64_di_t __ll = (u64_di_t) (x) * (y); \
383       hi = __ll >> 32; \
384       lo = __ll; \
385    } while (0)
386 #else
387 /*
388  * This MIPS version is fast and accurate; the disposition of the least
389  * significant bit depends on OPT_ACCURACY via mad_f_scale64().
390  */
391 #  define MAD_F_MLX(hi, lo, x, y)  \
392     asm ("mult	%2,%3"  \
393 	 : "=l" (lo), "=h" (hi)  \
394 	 : "%r" (x), "r" (y))
395 
396 # if defined(HAVE_MADD_ASM)
397 #  define MAD_F_MLA(hi, lo, x, y)  \
398     asm ("madd	%2,%3"  \
399 	 : "+l" (lo), "+h" (hi)  \
400 	 : "%r" (x), "r" (y))
401 # elif defined(HAVE_MADD16_ASM)
402 /*
403  * This loses significant accuracy due to the 16-bit integer limit in the
404  * multiply/accumulate instruction.
405  */
406 #  define MAD_F_ML0(hi, lo, x, y)  \
407     asm ("mult	%2,%3"  \
408 	 : "=l" (lo), "=h" (hi)  \
409 	 : "%r" ((x) >> 12), "r" ((y) >> 16))
410 #  define MAD_F_MLA(hi, lo, x, y)  \
411     asm ("madd16	%2,%3"  \
412 	 : "+l" (lo), "+h" (hi)  \
413 	 : "%r" ((x) >> 12), "r" ((y) >> 16))
414 #  define MAD_F_MLZ(hi, lo)  ((mad_fixed_t) (lo))
415 # endif
416 #endif /* MIPS / gcc-4.4. */
417 
418 # if defined(OPT_SPEED)
419 #  define mad_f_scale64(hi, lo)  \
420     ((mad_fixed_t) ((hi) << (32 - MAD_F_SCALEBITS)))
421 #  define MAD_F_SCALEBITS  MAD_F_FRACBITS
422 # endif
423 
424 /* --- SPARC --------------------------------------------------------------- */
425 
426 # elif defined(FPM_SPARC)
427 
428 /*
429  * This SPARC V8 version is fast and accurate; the disposition of the least
430  * significant bit depends on OPT_ACCURACY via mad_f_scale64().
431  */
432 #  define MAD_F_MLX(hi, lo, x, y)  \
433     asm ("smul %2, %3, %0\n\t"  \
434 	 "rd %%y, %1"  \
435 	 : "=r" (lo), "=r" (hi)  \
436 	 : "%r" (x), "rI" (y))
437 
438 /* --- PowerPC ------------------------------------------------------------- */
439 
440 # elif defined(FPM_PPC)
441 
442 /*
443  * This PowerPC version is fast and accurate; the disposition of the least
444  * significant bit depends on OPT_ACCURACY via mad_f_scale64().
445  */
446 #  define MAD_F_MLX(hi, lo, x, y)  \
447     do {  \
448       asm ("mullw %0,%1,%2"  \
449 	   : "=r" (lo)  \
450 	   : "%r" (x), "r" (y));  \
451       asm ("mulhw %0,%1,%2"  \
452 	   : "=r" (hi)  \
453 	   : "%r" (x), "r" (y));  \
454     }  \
455     while (0)
456 
457 #  if defined(OPT_ACCURACY)
458 /*
459  * This gives best accuracy but is not very fast.
460  */
461 #   define MAD_F_MLA(hi, lo, x, y)  \
462     ({ mad_fixed64hi_t __hi;  \
463        mad_fixed64lo_t __lo;  \
464        MAD_F_MLX(__hi, __lo, (x), (y));  \
465        asm ("addc %0,%2,%3\n\t"  \
466 	    "adde %1,%4,%5"  \
467 	    : "=r" (lo), "=r" (hi)  \
468 	    : "0" (lo), "r" (__lo), \
469 	      "1" (hi), "r" (__hi) \
470 	    : "xer");  \
471     })
472 #  endif
473 
474 #  if defined(OPT_ACCURACY)
475 /*
476  * This is slower than the truncating version below it.
477  */
478 #   define mad_f_scale64(hi, lo)  \
479     ({ mad_fixed_t __result, __round;  \
480        asm ("rotrwi %0,%1,%2"  \
481 	    : "=r" (__result)  \
482 	    : "r" (lo), "i" (MAD_F_SCALEBITS));  \
483        asm ("extrwi %0,%1,1,0"  \
484 	    : "=r" (__round)  \
485 	    : "r" (__result));  \
486        asm ("insrwi %0,%1,%2,0"  \
487 	    : "+r" (__result)  \
488 	    : "r" (hi), "i" (MAD_F_SCALEBITS));  \
489        asm ("add %0,%1,%2"  \
490 	    : "=r" (__result)  \
491 	    : "%r" (__result), "r" (__round));  \
492        __result;  \
493     })
494 #  else
495 #   define mad_f_scale64(hi, lo)  \
496     ({ mad_fixed_t __result;  \
497        asm ("rotrwi %0,%1,%2"  \
498 	    : "=r" (__result)  \
499 	    : "r" (lo), "i" (MAD_F_SCALEBITS));  \
500        asm ("insrwi %0,%1,%2,0"  \
501 	    : "+r" (__result)  \
502 	    : "r" (hi), "i" (MAD_F_SCALEBITS));  \
503        __result;  \
504     })
505 #  endif
506 
507 #  define MAD_F_SCALEBITS  MAD_F_FRACBITS
508 
509 /* --- Default ------------------------------------------------------------- */
510 
511 # elif defined(FPM_DEFAULT)
512 
513 /*
514  * This version is the most portable but it loses significant accuracy.
515  * Furthermore, accuracy is biased against the second argument, so care
516  * should be taken when ordering operands.
517  *
518  * The scale factors are constant as this is not used with SSO.
519  *
520  * Pre-rounding is required to stay within the limits of compliance.
521  */
522 #  if defined(OPT_SPEED)
523 #   define mad_f_mul(x, y)	(((x) >> 12) * ((y) >> 16))
524 #  else
525 #   define mad_f_mul(x, y)	((((x) + (1L << 11)) >> 12) *  \
526 				 (((y) + (1L << 15)) >> 16))
527 #  endif
528 
529 /* ------------------------------------------------------------------------- */
530 
531 # else
532 #  error "no FPM selected"
533 # endif
534 
535 /* default implementations */
536 
537 # if !defined(mad_f_mul)
538 #  define mad_f_mul(x, y)  \
539     ({ register mad_fixed64hi_t __hi;  \
540        register mad_fixed64lo_t __lo;  \
541        MAD_F_MLX(__hi, __lo, (x), (y));  \
542        mad_f_scale64(__hi, __lo);  \
543     })
544 # endif
545 
546 # if !defined(MAD_F_MLA)
547 #  define MAD_F_ML0(hi, lo, x, y)	((lo)  = mad_f_mul((x), (y)))
548 #  define MAD_F_MLA(hi, lo, x, y)	((lo) += mad_f_mul((x), (y)))
549 #  define MAD_F_MLN(hi, lo)		((lo)  = -(lo))
550 #  define MAD_F_MLZ(hi, lo)		((void) (hi), (mad_fixed_t) (lo))
551 # endif
552 
553 # if !defined(MAD_F_ML0)
554 #  define MAD_F_ML0(hi, lo, x, y)	MAD_F_MLX((hi), (lo), (x), (y))
555 # endif
556 
557 # if !defined(MAD_F_MLN)
558 #  define MAD_F_MLN(hi, lo)		((hi) = ((lo) = -(lo)) ? ~(hi) : -(hi))
559 # endif
560 
561 # if !defined(MAD_F_MLZ)
562 #  define MAD_F_MLZ(hi, lo)		mad_f_scale64((hi), (lo))
563 # endif
564 
565 # if !defined(mad_f_scale64)
566 #  if defined(OPT_ACCURACY)
567 #   define mad_f_scale64(hi, lo)  \
568     ((((mad_fixed_t)  \
569        (((hi) << (32 - (MAD_F_SCALEBITS - 1))) |  \
570 	((lo) >> (MAD_F_SCALEBITS - 1)))) + 1) >> 1)
571 #  else
572 #   define mad_f_scale64(hi, lo)  \
573     ((mad_fixed_t)  \
574      (((hi) << (32 - MAD_F_SCALEBITS)) |  \
575       ((lo) >> MAD_F_SCALEBITS)))
576 #  endif
577 #  define MAD_F_SCALEBITS  MAD_F_FRACBITS
578 # endif
579 
580 /* C routines */
581 
582 mad_fixed_t mad_f_abs(mad_fixed_t);
583 mad_fixed_t mad_f_div(mad_fixed_t, mad_fixed_t);
584 
585 # endif
586 
587 /* Id: bit.h,v 1.12 2004/01/23 09:41:32 rob Exp */
588 
589 # ifndef LIBMAD_BIT_H
590 # define LIBMAD_BIT_H
591 
592 struct mad_bitptr {
593   unsigned char const *byte;
594   unsigned short cache;
595   unsigned short left;
596 };
597 
598 void mad_bit_init(struct mad_bitptr *, unsigned char const *);
599 
600 # define mad_bit_finish(bitptr)		/* nothing */
601 
602 unsigned int mad_bit_length(struct mad_bitptr const *,
603 			    struct mad_bitptr const *);
604 
605 # define mad_bit_bitsleft(bitptr)  ((bitptr)->left)
606 unsigned char const *mad_bit_nextbyte(struct mad_bitptr const *);
607 
608 void mad_bit_skip(struct mad_bitptr *, unsigned int);
609 unsigned long mad_bit_read(struct mad_bitptr *, unsigned int);
610 void mad_bit_write(struct mad_bitptr *, unsigned int, unsigned long);
611 
612 unsigned short mad_bit_crc(struct mad_bitptr, unsigned int, unsigned short);
613 
614 # endif
615 
616 /* Id: timer.h,v 1.16 2004/01/23 09:41:33 rob Exp */
617 
618 # ifndef LIBMAD_TIMER_H
619 # define LIBMAD_TIMER_H
620 
621 typedef struct {
622   signed long seconds;		/* whole seconds */
623   unsigned long fraction;	/* 1/MAD_TIMER_RESOLUTION seconds */
624 } mad_timer_t;
625 
626 extern mad_timer_t const mad_timer_zero;
627 
628 # define MAD_TIMER_RESOLUTION	352800000UL
629 
630 enum mad_units {
631   MAD_UNITS_HOURS	 =    -2,
632   MAD_UNITS_MINUTES	 =    -1,
633   MAD_UNITS_SECONDS	 =     0,
634 
635   /* metric units */
636 
637   MAD_UNITS_DECISECONDS	 =    10,
638   MAD_UNITS_CENTISECONDS =   100,
639   MAD_UNITS_MILLISECONDS =  1000,
640 
641   /* audio sample units */
642 
643   MAD_UNITS_8000_HZ	 =  8000,
644   MAD_UNITS_11025_HZ	 = 11025,
645   MAD_UNITS_12000_HZ	 = 12000,
646 
647   MAD_UNITS_16000_HZ	 = 16000,
648   MAD_UNITS_22050_HZ	 = 22050,
649   MAD_UNITS_24000_HZ	 = 24000,
650 
651   MAD_UNITS_32000_HZ	 = 32000,
652   MAD_UNITS_44100_HZ	 = 44100,
653   MAD_UNITS_48000_HZ	 = 48000,
654 
655   /* video frame/field units */
656 
657   MAD_UNITS_24_FPS	 =    24,
658   MAD_UNITS_25_FPS	 =    25,
659   MAD_UNITS_30_FPS	 =    30,
660   MAD_UNITS_48_FPS	 =    48,
661   MAD_UNITS_50_FPS	 =    50,
662   MAD_UNITS_60_FPS	 =    60,
663 
664   /* CD audio frames */
665 
666   MAD_UNITS_75_FPS	 =    75,
667 
668   /* video drop-frame units */
669 
670   MAD_UNITS_23_976_FPS	 =   -24,
671   MAD_UNITS_24_975_FPS	 =   -25,
672   MAD_UNITS_29_97_FPS	 =   -30,
673   MAD_UNITS_47_952_FPS	 =   -48,
674   MAD_UNITS_49_95_FPS	 =   -50,
675   MAD_UNITS_59_94_FPS	 =   -60
676 };
677 
678 # define mad_timer_reset(timer)	((void) (*(timer) = mad_timer_zero))
679 
680 int mad_timer_compare(mad_timer_t, mad_timer_t);
681 
682 # define mad_timer_sign(timer)	mad_timer_compare((timer), mad_timer_zero)
683 
684 void mad_timer_negate(mad_timer_t *);
685 mad_timer_t mad_timer_abs(mad_timer_t);
686 
687 void mad_timer_set(mad_timer_t *, unsigned long, unsigned long, unsigned long);
688 void mad_timer_add(mad_timer_t *, mad_timer_t);
689 void mad_timer_multiply(mad_timer_t *, signed long);
690 
691 signed long mad_timer_count(mad_timer_t, enum mad_units);
692 unsigned long mad_timer_fraction(mad_timer_t, unsigned long);
693 void mad_timer_string(mad_timer_t, char *, char const *,
694 		      enum mad_units, enum mad_units, unsigned long);
695 
696 # endif
697 
698 /* Id: stream.h,v 1.20 2004/02/05 09:02:39 rob Exp */
699 
700 # ifndef LIBMAD_STREAM_H
701 # define LIBMAD_STREAM_H
702 
703 
704 # define MAD_BUFFER_GUARD	8
705 # define MAD_BUFFER_MDLEN	(511 + 2048 + MAD_BUFFER_GUARD)
706 
707 enum mad_error {
708   MAD_ERROR_NONE	   = 0x0000,	/* no error */
709 
710   MAD_ERROR_BUFLEN	   = 0x0001,	/* input buffer too small (or EOF) */
711   MAD_ERROR_BUFPTR	   = 0x0002,	/* invalid (null) buffer pointer */
712 
713   MAD_ERROR_NOMEM	   = 0x0031,	/* not enough memory */
714 
715   MAD_ERROR_LOSTSYNC	   = 0x0101,	/* lost synchronization */
716   MAD_ERROR_BADLAYER	   = 0x0102,	/* reserved header layer value */
717   MAD_ERROR_BADBITRATE	   = 0x0103,	/* forbidden bitrate value */
718   MAD_ERROR_BADSAMPLERATE  = 0x0104,	/* reserved sample frequency value */
719   MAD_ERROR_BADEMPHASIS	   = 0x0105,	/* reserved emphasis value */
720 
721   MAD_ERROR_BADCRC	   = 0x0201,	/* CRC check failed */
722   MAD_ERROR_BADBITALLOC	   = 0x0211,	/* forbidden bit allocation value */
723   MAD_ERROR_BADSCALEFACTOR = 0x0221,	/* bad scalefactor index */
724   MAD_ERROR_BADMODE        = 0x0222,	/* bad bitrate/mode combination */
725   MAD_ERROR_BADFRAMELEN	   = 0x0231,	/* bad frame length */
726   MAD_ERROR_BADBIGVALUES   = 0x0232,	/* bad big_values count */
727   MAD_ERROR_BADBLOCKTYPE   = 0x0233,	/* reserved block_type */
728   MAD_ERROR_BADSCFSI	   = 0x0234,	/* bad scalefactor selection info */
729   MAD_ERROR_BADDATAPTR	   = 0x0235,	/* bad main_data_begin pointer */
730   MAD_ERROR_BADPART3LEN	   = 0x0236,	/* bad audio data length */
731   MAD_ERROR_BADHUFFTABLE   = 0x0237,	/* bad Huffman table select */
732   MAD_ERROR_BADHUFFDATA	   = 0x0238,	/* Huffman data overrun */
733   MAD_ERROR_BADSTEREO	   = 0x0239	/* incompatible block_type for JS */
734 };
735 
736 # define MAD_RECOVERABLE(error)	((error) & 0xff00)
737 
738 struct mad_stream {
739   unsigned char const *buffer;		/* input bitstream buffer */
740   unsigned char const *bufend;		/* end of buffer */
741   unsigned long skiplen;		/* bytes to skip before next frame */
742 
743   int sync;				/* stream sync found */
744   unsigned long freerate;		/* free bitrate (fixed) */
745 
746   unsigned char const *this_frame;	/* start of current frame */
747   unsigned char const *next_frame;	/* start of next frame */
748   struct mad_bitptr ptr;		/* current processing bit pointer */
749 
750   struct mad_bitptr anc_ptr;		/* ancillary bits pointer */
751   unsigned int anc_bitlen;		/* number of ancillary bits */
752 
753   unsigned char (*main_data)[MAD_BUFFER_MDLEN];
754 					/* Layer III main_data() */
755   unsigned int md_len;			/* bytes in main_data */
756 
757   int options;				/* decoding options (see below) */
758   enum mad_error error;			/* error code (see above) */
759 };
760 
761 enum {
762   MAD_OPTION_IGNORECRC      = 0x0001,	/* ignore CRC errors */
763   MAD_OPTION_HALFSAMPLERATE = 0x0002	/* generate PCM at 1/2 sample rate */
764 # if 0  /* not yet implemented */
765   MAD_OPTION_LEFTCHANNEL    = 0x0010,	/* decode left channel only */
766   MAD_OPTION_RIGHTCHANNEL   = 0x0020,	/* decode right channel only */
767   MAD_OPTION_SINGLECHANNEL  = 0x0030	/* combine channels */
768 # endif
769 };
770 
771 void mad_stream_init(struct mad_stream *);
772 void mad_stream_finish(struct mad_stream *);
773 
774 # define mad_stream_options(stream, opts)  \
775     ((void) ((stream)->options = (opts)))
776 
777 void mad_stream_buffer(struct mad_stream *,
778 		       unsigned char const *, unsigned long);
779 void mad_stream_skip(struct mad_stream *, unsigned long);
780 
781 int mad_stream_sync(struct mad_stream *);
782 
783 char const *mad_stream_errorstr(struct mad_stream const *);
784 
785 # endif
786 
787 /* Id: frame.h,v 1.20 2004/01/23 09:41:32 rob Exp */
788 
789 # ifndef LIBMAD_FRAME_H
790 # define LIBMAD_FRAME_H
791 
792 
793 enum mad_layer {
794   MAD_LAYER_I   = 1,			/* Layer I */
795   MAD_LAYER_II  = 2,			/* Layer II */
796   MAD_LAYER_III = 3			/* Layer III */
797 };
798 
799 enum mad_mode {
800   MAD_MODE_SINGLE_CHANNEL = 0,		/* single channel */
801   MAD_MODE_DUAL_CHANNEL	  = 1,		/* dual channel */
802   MAD_MODE_JOINT_STEREO	  = 2,		/* joint (MS/intensity) stereo */
803   MAD_MODE_STEREO	  = 3		/* normal LR stereo */
804 };
805 
806 enum mad_emphasis {
807   MAD_EMPHASIS_NONE	  = 0,		/* no emphasis */
808   MAD_EMPHASIS_50_15_US	  = 1,		/* 50/15 microseconds emphasis */
809   MAD_EMPHASIS_CCITT_J_17 = 3,		/* CCITT J.17 emphasis */
810   MAD_EMPHASIS_RESERVED   = 2		/* unknown emphasis */
811 };
812 
813 struct mad_header {
814   enum mad_layer layer;			/* audio layer (1, 2, or 3) */
815   enum mad_mode mode;			/* channel mode (see above) */
816   int mode_extension;			/* additional mode info */
817   enum mad_emphasis emphasis;		/* de-emphasis to use (see above) */
818 
819   unsigned long bitrate;		/* stream bitrate (bps) */
820   unsigned int samplerate;		/* sampling frequency (Hz) */
821 
822   unsigned short crc_check;		/* frame CRC accumulator */
823   unsigned short crc_target;		/* final target CRC checksum */
824 
825   int flags;				/* flags (see below) */
826   int private_bits;			/* private bits (see below) */
827 
828   mad_timer_t duration;			/* audio playing time of frame */
829 };
830 
831 struct mad_frame {
832   struct mad_header header;		/* MPEG audio header */
833 
834   int options;				/* decoding options (from stream) */
835 
836   mad_fixed_t sbsample[2][36][32];	/* synthesis subband filter samples */
837   mad_fixed_t (*overlap)[2][32][18];	/* Layer III block overlap data */
838 };
839 
840 # define MAD_NCHANNELS(header)		((header)->mode ? 2 : 1)
841 # define MAD_NSBSAMPLES(header)  \
842   ((header)->layer == MAD_LAYER_I ? 12 :  \
843    (((header)->layer == MAD_LAYER_III &&  \
844      ((header)->flags & MAD_FLAG_LSF_EXT)) ? 18 : 36))
845 
846 enum {
847   MAD_FLAG_NPRIVATE_III	= 0x0007,	/* number of Layer III private bits */
848   MAD_FLAG_INCOMPLETE	= 0x0008,	/* header but not data is decoded */
849 
850   MAD_FLAG_PROTECTION	= 0x0010,	/* frame has CRC protection */
851   MAD_FLAG_COPYRIGHT	= 0x0020,	/* frame is copyright */
852   MAD_FLAG_ORIGINAL	= 0x0040,	/* frame is original (else copy) */
853   MAD_FLAG_PADDING	= 0x0080,	/* frame has additional slot */
854 
855   MAD_FLAG_I_STEREO	= 0x0100,	/* uses intensity joint stereo */
856   MAD_FLAG_MS_STEREO	= 0x0200,	/* uses middle/side joint stereo */
857   MAD_FLAG_FREEFORMAT	= 0x0400,	/* uses free format bitrate */
858 
859   MAD_FLAG_LSF_EXT	= 0x1000,	/* lower sampling freq. extension */
860   MAD_FLAG_MC_EXT	= 0x2000,	/* multichannel audio extension */
861   MAD_FLAG_MPEG_2_5_EXT	= 0x4000	/* MPEG 2.5 (unofficial) extension */
862 };
863 
864 enum {
865   MAD_PRIVATE_HEADER	= 0x0100,	/* header private bit */
866   MAD_PRIVATE_III	= 0x001f	/* Layer III private bits (up to 5) */
867 };
868 
869 void mad_header_init(struct mad_header *);
870 
871 # define mad_header_finish(header)  /* nothing */
872 
873 int mad_header_decode(struct mad_header *, struct mad_stream *);
874 
875 void mad_frame_init(struct mad_frame *);
876 void mad_frame_finish(struct mad_frame *);
877 
878 int mad_frame_decode(struct mad_frame *, struct mad_stream *);
879 
880 void mad_frame_mute(struct mad_frame *);
881 
882 # endif
883 
884 /* Id: synth.h,v 1.15 2004/01/23 09:41:33 rob Exp */
885 
886 # ifndef LIBMAD_SYNTH_H
887 # define LIBMAD_SYNTH_H
888 
889 
890 struct mad_pcm {
891   unsigned int samplerate;		/* sampling frequency (Hz) */
892   unsigned short channels;		/* number of channels */
893   unsigned short length;		/* number of samples per channel */
894   mad_fixed_t samples[2][1152];		/* PCM output samples [ch][sample] */
895 };
896 
897 struct mad_synth {
898   mad_fixed_t filter[2][2][2][16][8];	/* polyphase filterbank outputs */
899   					/* [ch][eo][peo][s][v] */
900 
901   unsigned int phase;			/* current processing phase */
902 
903   struct mad_pcm pcm;			/* PCM output */
904 };
905 
906 /* single channel PCM selector */
907 enum {
908   MAD_PCM_CHANNEL_SINGLE = 0
909 };
910 
911 /* dual channel PCM selector */
912 enum {
913   MAD_PCM_CHANNEL_DUAL_1 = 0,
914   MAD_PCM_CHANNEL_DUAL_2 = 1
915 };
916 
917 /* stereo PCM selector */
918 enum {
919   MAD_PCM_CHANNEL_STEREO_LEFT  = 0,
920   MAD_PCM_CHANNEL_STEREO_RIGHT = 1
921 };
922 
923 void mad_synth_init(struct mad_synth *);
924 
925 # define mad_synth_finish(synth)  /* nothing */
926 
927 void mad_synth_mute(struct mad_synth *);
928 
929 void mad_synth_frame(struct mad_synth *, struct mad_frame const *);
930 
931 # endif
932 
933 /* Id: decoder.h,v 1.17 2004/01/23 09:41:32 rob Exp */
934 
935 # ifndef LIBMAD_DECODER_H
936 # define LIBMAD_DECODER_H
937 
938 
939 enum mad_decoder_mode {
940   MAD_DECODER_MODE_SYNC  = 0,
941   MAD_DECODER_MODE_ASYNC
942 };
943 
944 enum mad_flow {
945   MAD_FLOW_CONTINUE = 0x0000,	/* continue normally */
946   MAD_FLOW_STOP     = 0x0010,	/* stop decoding normally */
947   MAD_FLOW_BREAK    = 0x0011,	/* stop decoding and signal an error */
948   MAD_FLOW_IGNORE   = 0x0020	/* ignore the current frame */
949 };
950 
951 struct mad_decoder {
952   enum mad_decoder_mode mode;
953 
954   int options;
955 
956   struct {
957     long pid;
958     int in;
959     int out;
960   } async;
961 
962   struct {
963     struct mad_stream stream;
964     struct mad_frame frame;
965     struct mad_synth synth;
966   } *sync;
967 
968   void *cb_data;
969 
970   enum mad_flow (*input_func)(void *, struct mad_stream *);
971   enum mad_flow (*header_func)(void *, struct mad_header const *);
972   enum mad_flow (*filter_func)(void *,
973 			       struct mad_stream const *, struct mad_frame *);
974   enum mad_flow (*output_func)(void *,
975 			       struct mad_header const *, struct mad_pcm *);
976   enum mad_flow (*error_func)(void *, struct mad_stream *, struct mad_frame *);
977   enum mad_flow (*message_func)(void *, void *, unsigned int *);
978 };
979 
980 void mad_decoder_init(struct mad_decoder *, void *,
981 		      enum mad_flow (*)(void *, struct mad_stream *),
982 		      enum mad_flow (*)(void *, struct mad_header const *),
983 		      enum mad_flow (*)(void *,
984 					struct mad_stream const *,
985 					struct mad_frame *),
986 		      enum mad_flow (*)(void *,
987 					struct mad_header const *,
988 					struct mad_pcm *),
989 		      enum mad_flow (*)(void *,
990 					struct mad_stream *,
991 					struct mad_frame *),
992 		      enum mad_flow (*)(void *, void *, unsigned int *));
993 int mad_decoder_finish(struct mad_decoder *);
994 
995 # define mad_decoder_options(decoder, opts)  \
996     ((void) ((decoder)->options = (opts)))
997 
998 int mad_decoder_run(struct mad_decoder *, enum mad_decoder_mode);
999 int mad_decoder_message(struct mad_decoder *, void *, unsigned int *);
1000 
1001 # endif
1002 
1003 #ifdef __cplusplus
1004 }
1005 #endif
1006