seal 0.1.0

data/mpg123/src/libmpg123/synths.h

@@ -0,0 +1,52 @@
+#ifndef MPG123_SYNTH_H
+#define MPG123_SYNTH_H
+
+/* This is included inside frame.h, which is included in mpg123lib_intern.h,
+   at the appropriate place.
+   Explicit header inclusions here would cause circular dependencies. */
+
+/* The handle needs these types for selecting the decoding routine at runtime.
+   Not just for optimization, mainly for XtoY, mono/stereo. */
+typedef int (*func_synth)(real *,int, mpg123_handle *,int );
+typedef int (*func_synth_mono)(real *, mpg123_handle *);
+typedef int (*func_synth_stereo)(real *, real *, mpg123_handle *);
+enum synth_channel  { c_plain=0, c_stereo, c_m2s, c_mono, c_limit };
+enum synth_resample
+{
+	 r_none=-1
+	,r_1to1=0
+#	ifndef NO_DOWNSAMPLE
+	,r_2to1
+	,r_4to1
+#	endif
+#	ifndef NO_NTOM
+	,r_ntom
+#	endif
+	,r_limit
+};
+enum synth_format
+{
+	 f_none=-1
+#	ifndef NO_16BIT
+	,f_16
+#	endif
+#	ifndef NO_8BIT
+	,f_8
+#	endif
+#	ifndef NO_REAL
+	,f_real
+#	endif
+#	ifndef NO_32BIT
+	,f_32
+#	endif
+	,f_limit
+};
+struct synth_s
+{
+	func_synth              plain[r_limit][f_limit];
+	func_synth_stereo      stereo[r_limit][f_limit];
+	func_synth_mono   mono2stereo[r_limit][f_limit];
+	func_synth_mono          mono[r_limit][f_limit];
+};
+
+#endif

data/mpg123/src/libmpg123/tabinit.c

@@ -0,0 +1,294 @@
+/*
+	tabinit.c: initialize tables...
+
+	copyright ?-2008 by the mpg123 project - free software under the terms of the LGPL 2.1
+	see COPYING and AUTHORS files in distribution or http://mpg123.org
+	initially written by Michael Hipp
+*/
+
+#include "mpg123lib_intern.h"
+#include "debug.h"
+
+/* That altivec alignment part here should not hurt generic code, I hope */
+#ifdef OPT_ALTIVEC
+static ALIGNED(16) real cos64[16];
+static ALIGNED(16) real cos32[8];
+static ALIGNED(16) real cos16[4];
+static ALIGNED(16) real cos8[2];
+static ALIGNED(16) real cos4[1];
+#elif defined(REAL_IS_FIXED) && defined(PRECALC_TABLES)
+static real cos64[16] = 
+{
+	8398725,8480395,8647771,8909416,9279544,9780026,10443886,11321405,
+	12491246,14081950,16316987,19619946,24900150,34523836,57170182,170959967
+};
+static real cos32[8] =
+{
+	8429197,8766072,9511743,10851869,13223040,17795219,28897867,85583072
+};
+static real cos16[4] =
+{
+	8552951,10088893,15099095,42998586
+};
+static real cos8[2] =
+{
+	9079764,21920489
+};
+static real cos4[1] =
+{
+	11863283
+};
+#else
+static real cos64[16],cos32[8],cos16[4],cos8[2],cos4[1];
+#endif
+
+real *pnts[] = { cos64,cos32,cos16,cos8,cos4 };
+
+
+static long intwinbase[] = {
+     0,    -1,    -1,    -1,    -1,    -1,    -1,    -2,    -2,    -2,
+    -2,    -3,    -3,    -4,    -4,    -5,    -5,    -6,    -7,    -7,
+    -8,    -9,   -10,   -11,   -13,   -14,   -16,   -17,   -19,   -21,
+   -24,   -26,   -29,   -31,   -35,   -38,   -41,   -45,   -49,   -53,
+   -58,   -63,   -68,   -73,   -79,   -85,   -91,   -97,  -104,  -111,
+  -117,  -125,  -132,  -139,  -147,  -154,  -161,  -169,  -176,  -183,
+  -190,  -196,  -202,  -208,  -213,  -218,  -222,  -225,  -227,  -228,
+  -228,  -227,  -224,  -221,  -215,  -208,  -200,  -189,  -177,  -163,
+  -146,  -127,  -106,   -83,   -57,   -29,     2,    36,    72,   111,
+   153,   197,   244,   294,   347,   401,   459,   519,   581,   645,
+   711,   779,   848,   919,   991,  1064,  1137,  1210,  1283,  1356,
+  1428,  1498,  1567,  1634,  1698,  1759,  1817,  1870,  1919,  1962,
+  2001,  2032,  2057,  2075,  2085,  2087,  2080,  2063,  2037,  2000,
+  1952,  1893,  1822,  1739,  1644,  1535,  1414,  1280,  1131,   970,
+   794,   605,   402,   185,   -45,  -288,  -545,  -814, -1095, -1388,
+ -1692, -2006, -2330, -2663, -3004, -3351, -3705, -4063, -4425, -4788,
+ -5153, -5517, -5879, -6237, -6589, -6935, -7271, -7597, -7910, -8209,
+ -8491, -8755, -8998, -9219, -9416, -9585, -9727, -9838, -9916, -9959,
+ -9966, -9935, -9863, -9750, -9592, -9389, -9139, -8840, -8492, -8092,
+ -7640, -7134, -6574, -5959, -5288, -4561, -3776, -2935, -2037, -1082,
+   -70,   998,  2122,  3300,  4533,  5818,  7154,  8540,  9975, 11455,
+ 12980, 14548, 16155, 17799, 19478, 21189, 22929, 24694, 26482, 28289,
+ 30112, 31947, 33791, 35640, 37489, 39336, 41176, 43006, 44821, 46617,
+ 48390, 50137, 51853, 53534, 55178, 56778, 58333, 59838, 61289, 62684,
+ 64019, 65290, 66494, 67629, 68692, 69679, 70590, 71420, 72169, 72835,
+ 73415, 73908, 74313, 74630, 74856, 74992, 75038 };
+
+void prepare_decode_tables()
+{
+#if !defined(REAL_IS_FIXED) || !defined(PRECALC_TABLES)
+  int i,k,kr,divv;
+  real *costab;
+
+  for(i=0;i<5;i++)
+  {
+    kr=0x10>>i; divv=0x40>>i;
+    costab = pnts[i];
+    for(k=0;k<kr;k++)
+      costab[k] = DOUBLE_TO_REAL(1.0 / (2.0 * cos(M_PI * ((double) k * 2.0 + 1.0) / (double) divv)));
+  }
+#endif
+}
+
+#ifdef OPT_MMXORSSE
+#if !defined(OPT_X86_64) && !defined(OPT_NEON)
+void make_decode_tables_mmx_asm(long scaleval, float* decwin_mmx, float *decwins);
+void make_decode_tables_mmx(mpg123_handle *fr)
+{
+	debug("MMX decode tables");
+	/* Take care: The scale should be like before, when we didn't have float output all around. */
+	make_decode_tables_mmx_asm((long)((fr->lastscale < 0 ? fr->p.outscale : fr->lastscale)*SHORT_SCALE), fr->decwin_mmx, fr->decwins);
+	debug("MMX decode tables done");
+}
+#else
+
+/* This mimics round() as defined in C99. We stay C89. */
+static int rounded(double f)
+{
+	return (int)(f>0 ? floor(f+0.5) : ceil(f-0.5));
+}
+
+/* x86-64 doesn't use asm version */
+void make_decode_tables_mmx(mpg123_handle *fr)
+{
+	int i,j,val;
+	int idx = 0;
+	short *ptr = (short *)fr->decwins;
+	/* Scale is always based on 1.0 . */
+	double scaleval = -0.5*(fr->lastscale < 0 ? fr->p.outscale : fr->lastscale);
+	debug1("MMX decode tables with scaleval %g", scaleval);
+	for(i=0,j=0;i<256;i++,j++,idx+=32)
+	{
+		if(idx < 512+16)
+		fr->decwin_mmx[idx+16] = fr->decwin_mmx[idx] = DOUBLE_TO_REAL((double) intwinbase[j] * scaleval);
+		
+		if(i % 32 == 31)
+		idx -= 1023;
+		if(i % 64 == 63)
+		scaleval = - scaleval;
+	}
+	
+	for( /* i=256 */ ;i<512;i++,j--,idx+=32)
+	{
+		if(idx < 512+16)
+		fr->decwin_mmx[idx+16] = fr->decwin_mmx[idx] = DOUBLE_TO_REAL((double) intwinbase[j] * scaleval);
+		
+		if(i % 32 == 31)
+		idx -= 1023;
+		if(i % 64 == 63)
+		scaleval = - scaleval;
+	}
+	
+	for(i=0; i<512; i++) {
+		if(i&1) val = rounded(fr->decwin_mmx[i]*0.5);
+		else val = rounded(fr->decwin_mmx[i]*-0.5);
+		if(val > 32767) val = 32767;
+		else if(val < -32768) val = -32768;
+		ptr[i] = val;
+	}
+	for(i=512; i<512+32; i++) {
+		if(i&1) val = rounded(fr->decwin_mmx[i]*0.5);
+		else val = 0;
+		if(val > 32767) val = 32767;
+		else if(val < -32768) val = -32768;
+		ptr[i] = val;
+	}
+	for(i=0; i<512; i++) {
+		val = rounded(fr->decwin_mmx[511-i]*-0.5);
+		if(val > 32767) val = 32767;
+		else if(val < -32768) val = -32768;
+		ptr[512+32+i] = val;
+	}
+	debug("decode tables done");
+}
+#endif
+#endif
+
+void make_decode_tables(mpg123_handle *fr)
+{
+	int i,j;
+	int idx = 0;
+	/* Scale is always based on 1.0 . */
+	double scaleval = -0.5*(fr->lastscale < 0 ? fr->p.outscale : fr->lastscale);
+	debug1("decode tables with scaleval %g", scaleval);
+#ifdef REAL_IS_FIXED
+	long scaleval_long = DOUBLE_TO_REAL_15(scaleval);
+#endif
+	for(i=0,j=0;i<256;i++,j++,idx+=32)
+	{
+		if(idx < 512+16)
+#ifdef REAL_IS_FIXED
+		fr->decwin[idx+16] = fr->decwin[idx] = REAL_SCALE_WINDOW(intwinbase[j] * scaleval_long);
+#else
+		fr->decwin[idx+16] = fr->decwin[idx] = DOUBLE_TO_REAL((double) intwinbase[j] * scaleval);
+#endif
+
+		if(i % 32 == 31)
+		idx -= 1023;
+		if(i % 64 == 63)
+#ifdef REAL_IS_FIXED
+		scaleval_long = - scaleval_long;
+#else
+		scaleval = - scaleval;
+#endif
+	}
+
+	for( /* i=256 */ ;i<512;i++,j--,idx+=32)
+	{
+		if(idx < 512+16)
+#ifdef REAL_IS_FIXED
+		fr->decwin[idx+16] = fr->decwin[idx] = REAL_SCALE_WINDOW(intwinbase[j] * scaleval_long);
+#else
+		fr->decwin[idx+16] = fr->decwin[idx] = DOUBLE_TO_REAL((double) intwinbase[j] * scaleval);
+#endif
+
+		if(i % 32 == 31)
+		idx -= 1023;
+		if(i % 64 == 63)
+#ifdef REAL_IS_FIXED
+		scaleval_long = - scaleval_long;
+#else
+		scaleval = - scaleval;
+#endif
+	}
+#if defined(OPT_X86_64) || defined(OPT_ALTIVEC) || defined(OPT_SSE) || defined(OPT_ARM) || defined(OPT_NEON)
+	if(fr->cpu_opts.type == x86_64 || fr->cpu_opts.type == altivec || fr->cpu_opts.type == sse || fr->cpu_opts.type == arm || fr->cpu_opts.type == neon)
+	{ /* for float SSE / AltiVec / ARM decoder */
+		for(i=512; i<512+32; i++)
+		{
+			fr->decwin[i] = (i&1) ? fr->decwin[i] : 0;
+		}
+		for(i=0; i<512; i++)
+		{
+			fr->decwin[512+32+i] = -fr->decwin[511-i];
+		}
+#ifdef OPT_NEON
+		if(fr->cpu_opts.type == neon)
+		{
+			for(i=0; i<512; i+=2)
+			{
+				fr->decwin[i] = -fr->decwin[i];
+			}
+		}
+#endif
+	}
+#endif
+	debug("decode tables done");
+}
+
+#ifndef NO_8BIT
+int make_conv16to8_table(mpg123_handle *fr)
+{
+  int i;
+	int mode = fr->af.encoding;
+
+  /*
+   * ????: 8.0 is right but on SB cards '2.0' is a better value ???
+   */
+  const double mul = 8.0;
+
+  if(!fr->conv16to8_buf){
+    fr->conv16to8_buf = (unsigned char *) malloc(8192);
+    if(!fr->conv16to8_buf) {
+      fr->err = MPG123_ERR_16TO8TABLE;
+      if(NOQUIET) error("Can't allocate 16 to 8 converter table!");
+      return -1;
+    }
+    fr->conv16to8 = fr->conv16to8_buf + 4096;
+  }
+
+  if(fr->af.encoding == MPG123_ENC_ULAW_8){
+    double m=127.0 / log(256.0);
+    int c1;
+
+    for(i=-4096;i<4096;i++) {
+/* dunno whether this is a valid transformation rule ?!?!? */
+      if(i < 0)
+        c1 = 127 - (int) (log( 1.0 - 255.0 * (double) i*mul / 32768.0 ) * m);
+      else
+        c1 = 255 - (int) (log( 1.0 + 255.0 * (double) i*mul / 32768.0 ) * m);
+      if((c1 < 0 || c1 > 255) && NOQUIET) error2("Converror %d %d",i,c1);
+
+      if(c1 == 0)
+        c1 = 2;
+      fr->conv16to8[i] = (unsigned char) c1;
+    }
+  }
+  else if(mode == MPG123_ENC_SIGNED_8) {
+    for(i=-4096;i<4096;i++) {
+      fr->conv16to8[i] = i>>5;
+    }
+  }
+  else if(mode == MPG123_ENC_UNSIGNED_8) {
+    for(i=-4096;i<4096;i++) {
+      fr->conv16to8[i] = (i>>5)+128;
+    }
+  }
+  else {
+    for(i=-4096;i<4096;i++) {
+      fr->conv16to8[i] = 0;
+    }
+  }
+	return 0;
+}
+#endif
+

data/mpg123/src/libmpg123/tabinit_mmx.S

@@ -0,0 +1,210 @@
+/*
+	tabinit_mmx: make_decode_tables_mmx
+
+	copyright ?-2006 by the mpg123 project - free software under the terms of the LGPL 2.1
+	see COPYING and AUTHORS files in distribution or http://mpg123.org
+	initially written by the mysterious higway (apparently)
+*/
+
+#include "mangle.h"
+
+.data
+	ALIGN32
+.globl ASM_NAME(costab_mmxsse)
+ASM_NAME(costab_mmxsse):
+	.long 1056974725
+	.long 1057056395
+	.long 1057223771
+	.long 1057485416
+	.long 1057855544
+	.long 1058356026
+	.long 1059019886
+	.long 1059897405
+	.long 1061067246
+	.long 1062657950
+	.long 1064892987
+	.long 1066774581
+	.long 1069414683
+	.long 1073984175
+	.long 1079645762
+	.long 1092815430
+	.long 1057005197
+	.long 1057342072
+	.long 1058087743
+	.long 1059427869
+	.long 1061799040
+	.long 1065862217
+	.long 1071413542
+	.long 1084439708
+	.long 1057128951
+	.long 1058664893
+	.long 1063675095
+	.long 1076102863
+	.long 1057655764
+	.long 1067924853
+	.long 1060439283
+	ALIGN32
+intwinbase:
+	.short      0,    -1,    -1,    -1,    -1,    -1,    -1,    -2
+	.short     -2,    -2,    -2,    -3,    -3,    -4,    -4,    -5
+	.short     -5,    -6,    -7,    -7,    -8,    -9,   -10,   -11
+	.short    -13,   -14,   -16,   -17,   -19,   -21,   -24,   -26
+	.short    -29,   -31,   -35,   -38,   -41,   -45,   -49,   -53
+	.short    -58,   -63,   -68,   -73,   -79,   -85,   -91,   -97
+	.short   -104,  -111,  -117,  -125,  -132,  -139,  -147,  -154
+	.short   -161,  -169,  -176,  -183,  -190,  -196,  -202,  -208
+	.short   -213,  -218,  -222,  -225,  -227,  -228,  -228,  -227
+	.short   -224,  -221,  -215,  -208,  -200,  -189,  -177,  -163
+	.short   -146,  -127,  -106,   -83,   -57,   -29,     2,    36
+	.short     72,   111,   153,   197,   244,   294,   347,   401
+	.short    459,   519,   581,   645,   711,   779,   848,   919
+	.short    991,  1064,  1137,  1210,  1283,  1356,  1428,  1498
+	.short   1567,  1634,  1698,  1759,  1817,  1870,  1919,  1962
+	.short   2001,  2032,  2057,  2075,  2085,  2087,  2080,  2063
+	.short   2037,  2000,  1952,  1893,  1822,  1739,  1644,  1535
+	.short   1414,  1280,  1131,   970,   794,   605,   402,   185
+	.short    -45,  -288,  -545,  -814, -1095, -1388, -1692, -2006
+	.short  -2330, -2663, -3004, -3351, -3705, -4063, -4425, -4788
+	.short  -5153, -5517, -5879, -6237, -6589, -6935, -7271, -7597
+	.short  -7910, -8209, -8491, -8755, -8998, -9219, -9416, -9585
+	.short  -9727, -9838, -9916, -9959, -9966, -9935, -9863, -9750
+	.short  -9592, -9389, -9139, -8840, -8492, -8092, -7640, -7134
+	.short  -6574, -5959, -5288, -4561, -3776, -2935, -2037, -1082
+	.short    -70,   998,  2122,  3300,  4533,  5818,  7154,  8540
+	.short   9975, 11455, 12980, 14548, 16155, 17799, 19478, 21189
+	.short  22929, 24694, 26482, 28289, 30112, 31947,-26209,-24360
+	.short -22511,-20664,-18824,-16994,-15179,-13383,-11610, -9863
+	.short  -8147, -6466, -4822, -3222, -1667,  -162,  1289,  2684
+	.short   4019,  5290,  6494,  7629,  8692,  9679, 10590, 11420
+	.short  12169, 12835, 13415, 13908, 14313, 14630, 14856, 14992
+	.short  15038
+
+intwindiv:
+	.long 0x47800000			# 65536.0
+.text
+	ALIGN32
+/* void make_decode_tables_mmx_asm(long scaleval, float* decwin_mmx, float *decwins); */
+.globl ASM_NAME(make_decode_tables_mmx_asm)
+ASM_NAME(make_decode_tables_mmx_asm):
+	pushl %edi
+	pushl %esi
+	pushl %ebx
+/* stack: 0=ebx, 4=esi, 8=edi, 12=back, 16=scaleval, 20=decwin_mmx, 24=decwins */
+
+	xorl %ecx,%ecx
+	xorl %ebx,%ebx
+	movl $32,%esi
+	movl $intwinbase,%edi
+	negl 16(%esp)	/* scaleval */
+	pushl $2	/* intwinbase step */
+/* stack: 20=scaleval 24=decwin_mmx, 28=decwins */
+.L00:
+	cmpl $528,%ecx
+	jnc .L02
+	movswl (%edi),%eax
+	cmpl $intwinbase+444,%edi
+	jc .L01
+	addl $60000,%eax
+.L01:
+	pushl %eax
+/* stack: 24=scaleval 28=decwin_mmx, 32=decwins */
+	fildl (%esp)
+	fdivs intwindiv
+	fimull 24(%esp) /* scaleval */
+/* eax used to be popped the line before... I'll just use it here a bit */
+	movl 28(%esp),%eax /* decwin_mmx */
+	fsts    (%eax,%ecx,4)
+	fstps 64(%eax,%ecx,4)
+	popl %eax
+/* stack: 20=scaleval 24=decwin_mmx, 28=decwins */
+.L02:
+	leal -1(%esi),%edx
+	andl %ebx,%edx
+	cmpl $31,%edx
+	jnz .L03
+	addl $-1023,%ecx
+	testl %esi,%ebx
+	jz  .L03
+	negl 20(%esp)
+.L03:
+	addl %esi,%ecx
+	addl (%esp),%edi
+	incl %ebx
+	cmpl $intwinbase,%edi
+	jz .L04
+	cmpl $256,%ebx
+	jnz .L00
+	negl (%esp)
+	jmp .L00
+.L04:
+	popl %eax
+
+	xorl %ecx,%ecx
+	xorl %ebx,%ebx
+	pushl $2 /* paired with popl above */
+.L05:
+	cmpl $528,%ecx
+	jnc .L11
+	movswl (%edi),%eax
+	cmpl $intwinbase+444,%edi
+	jc .L06
+	addl $60000,%eax
+.L06:
+	cltd
+	imull 20(%esp)
+	shrdl $17,%edx,%eax
+	cmpl $32767,%eax
+	movl $1055,%edx
+	jle .L07
+	movl $32767,%eax
+	jmp .L08
+.L07:
+	cmpl $-32767,%eax
+	jge .L08
+	movl $-32767,%eax
+.L08:
+/* going to use ebx for decwins, watch the jumps */
+	pushl %ebx 
+/* stack: 24=scaleval 28=decwin_mmx, 32=decwins */
+	movl 32(%esp),%ebx
+	cmpl $512,%ecx
+	jnc .L09
+	subl %ecx,%edx
+	movw %ax,(%ebx,%edx,2) /* decwins */
+	movw %ax,-32(%ebx,%edx,2)
+.L09:
+	testl $1,%ecx
+	jnz .L10
+	negl %eax
+.L10:
+	movw %ax,(%ebx,%ecx,2)
+	movw %ax,32(%ebx,%ecx,2)
+	popl %ebx /* that has to match the pushl before */
+.L11:
+	leal -1(%esi),%edx
+	andl %ebx,%edx
+	cmpl $31,%edx
+	jnz .L12
+	addl $-1023,%ecx
+	testl %esi,%ebx
+	jz  .L12
+	negl 20(%esp)
+.L12:
+	addl %esi,%ecx
+	addl (%esp),%edi
+	incl %ebx
+	cmpl $intwinbase,%edi
+	jz .L13
+	cmpl $256,%ebx
+	jnz .L05
+	negl (%esp)
+	jmp .L05
+.L13:
+	popl %eax
+	
+	popl %ebx
+	popl %esi
+	popl %edi
+	ret
+
+NONEXEC_STACK