1 /*
2 * tabinit.c
3 *
4 * Copyright (C) 1999-2010 The L.A.M.E. project
5 *
6 * Initially written by Michael Hipp, see also AUTHORS and README.
7 *
8 * This library is free software; you can redistribute it and/or
9 * modify it under the terms of the GNU Library General Public
10 * License as published by the Free Software Foundation; either
11 * version 2 of the License, or (at your option) any later version.
12 *
13 * This library is distributed in the hope that it will be useful,
14 * but WITHOUT ANY WARRANTY; without even the implied warranty of
15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
16 * Library General Public License for more details.
17 *
18 * You should have received a copy of the GNU Library General Public
19 * License along with this library; if not, write to the
20 * Free Software Foundation, Inc., 59 Temple Place - Suite 330,
21 * Boston, MA 02111-1307, USA.
22 */
23 /* $Id: tabinit.c,v 1.17 2017/09/06 15:07:30 robert Exp $ */
24
25 #ifdef HAVE_CONFIG_H
26 # include <config.h>
27 #endif
28
29 #include <stdlib.h>
30 #include "tabinit.h"
31 #include "mpg123.h"
32
33 #ifdef WITH_DMALLOC
34 #include <dmalloc.h>
35 #endif
36
37 static int table_init_called = 0;
38
39 real decwin[512 + 32];
40 static real cos64[16], cos32[8], cos16[4], cos8[2], cos4[1];
41 real *pnts[] = { cos64, cos32, cos16, cos8, cos4 };
42
43 /* *INDENT-OFF* */
44
45 static const double dewin[512] = {
46 0.000000000,-0.000015259,-0.000015259,-0.000015259,
47 -0.000015259,-0.000015259,-0.000015259,-0.000030518,
48 -0.000030518,-0.000030518,-0.000030518,-0.000045776,
49 -0.000045776,-0.000061035,-0.000061035,-0.000076294,
50 -0.000076294,-0.000091553,-0.000106812,-0.000106812,
51 -0.000122070,-0.000137329,-0.000152588,-0.000167847,
52 -0.000198364,-0.000213623,-0.000244141,-0.000259399,
53 -0.000289917,-0.000320435,-0.000366211,-0.000396729,
54 -0.000442505,-0.000473022,-0.000534058,-0.000579834,
55 -0.000625610,-0.000686646,-0.000747681,-0.000808716,
56 -0.000885010,-0.000961304,-0.001037598,-0.001113892,
57 -0.001205444,-0.001296997,-0.001388550,-0.001480103,
58 -0.001586914,-0.001693726,-0.001785278,-0.001907349,
59 -0.002014160,-0.002120972,-0.002243042,-0.002349854,
60 -0.002456665,-0.002578735,-0.002685547,-0.002792358,
61 -0.002899170,-0.002990723,-0.003082275,-0.003173828,
62 -0.003250122,-0.003326416,-0.003387451,-0.003433228,
63 -0.003463745,-0.003479004,-0.003479004,-0.003463745,
64 -0.003417969,-0.003372192,-0.003280640,-0.003173828,
65 -0.003051758,-0.002883911,-0.002700806,-0.002487183,
66 -0.002227783,-0.001937866,-0.001617432,-0.001266479,
67 -0.000869751,-0.000442505, 0.000030518, 0.000549316,
68 0.001098633, 0.001693726, 0.002334595, 0.003005981,
69 0.003723145, 0.004486084, 0.005294800, 0.006118774,
70 0.007003784, 0.007919312, 0.008865356, 0.009841919,
71 0.010848999, 0.011886597, 0.012939453, 0.014022827,
72 0.015121460, 0.016235352, 0.017349243, 0.018463135,
73 0.019577026, 0.020690918, 0.021789551, 0.022857666,
74 0.023910522, 0.024932861, 0.025909424, 0.026840210,
75 0.027725220, 0.028533936, 0.029281616, 0.029937744,
76 0.030532837, 0.031005859, 0.031387329, 0.031661987,
77 0.031814575, 0.031845093, 0.031738281, 0.031478882,
78 0.031082153, 0.030517578, 0.029785156, 0.028884888,
79 0.027801514, 0.026535034, 0.025085449, 0.023422241,
80 0.021575928, 0.019531250, 0.017257690, 0.014801025,
81 0.012115479, 0.009231567, 0.006134033, 0.002822876,
82 -0.000686646,-0.004394531,-0.008316040,-0.012420654,
83 -0.016708374,-0.021179199,-0.025817871,-0.030609131,
84 -0.035552979,-0.040634155,-0.045837402,-0.051132202,
85 -0.056533813,-0.061996460,-0.067520142,-0.073059082,
86 -0.078628540,-0.084182739,-0.089706421,-0.095169067,
87 -0.100540161,-0.105819702,-0.110946655,-0.115921021,
88 -0.120697021,-0.125259399,-0.129562378,-0.133590698,
89 -0.137298584,-0.140670776,-0.143676758,-0.146255493,
90 -0.148422241,-0.150115967,-0.151306152,-0.151962280,
91 -0.152069092,-0.151596069,-0.150497437,-0.148773193,
92 -0.146362305,-0.143264771,-0.139450073,-0.134887695,
93 -0.129577637,-0.123474121,-0.116577148,-0.108856201,
94 -0.100311279,-0.090927124,-0.080688477,-0.069595337,
95 -0.057617187,-0.044784546,-0.031082153,-0.016510010,
96 -0.001068115, 0.015228271, 0.032379150, 0.050354004,
97 0.069168091, 0.088775635, 0.109161377, 0.130310059,
98 0.152206421, 0.174789429, 0.198059082, 0.221984863,
99 0.246505737, 0.271591187, 0.297210693, 0.323318481,
100 0.349868774, 0.376800537, 0.404083252, 0.431655884,
101 0.459472656, 0.487472534, 0.515609741, 0.543823242,
102 0.572036743, 0.600219727, 0.628295898, 0.656219482,
103 0.683914185, 0.711318970, 0.738372803, 0.765029907,
104 0.791213989, 0.816864014, 0.841949463, 0.866363525,
105 0.890090942, 0.913055420, 0.935195923, 0.956481934,
106 0.976852417, 0.996246338, 1.014617920, 1.031936646,
107 1.048156738, 1.063217163, 1.077117920, 1.089782715,
108 1.101211548, 1.111373901, 1.120223999, 1.127746582,
109 1.133926392, 1.138763428, 1.142211914, 1.144287109,
110 1.144989014
111 };
112 /* *INDENT-ON* */
113
114 void
make_decode_tables(long scaleval)115 make_decode_tables(long scaleval)
116 {
117 int i, j, k, kr, divv;
118 real *table, *costab;
119
120 if (table_init_called)
121 return;
122
123 table_init_called = 1;
124
125 for (i = 0; i < 5; i++) {
126 kr = 0x10 >> i;
127 divv = 0x40 >> i;
128 costab = pnts[i];
129 for (k = 0; k < kr; k++)
130 costab[k] = (real) (1.0 / (2.0 * cos(M_PI * ((double) k * 2.0 + 1.0) / (double) divv)));
131 }
132
133 table = decwin;
134 scaleval = -scaleval;
135 for (i = 0, j = 0; i < 256; i++, j++, table += 32) {
136 if (table < decwin + 512 + 16)
137 table[16] = table[0] = (real) (dewin[j] * scaleval);
138 if (i % 32 == 31)
139 table -= 1023;
140 if (i % 64 == 63)
141 scaleval = -scaleval;
142 }
143
144 for ( /* i=256 */ ; i < 512; i++, j--, table += 32) {
145 if (table < decwin + 512 + 16)
146 table[16] = table[0] = (real) (dewin[j] * scaleval);
147 if (i % 32 == 31)
148 table -= 1023;
149 if (i % 64 == 63)
150 scaleval = -scaleval;
151 }
152 }
153