/*

 * tabinit.c

 *

 * Copyright (C) 1999-2010 The L.A.M.E. project

 *

 * Initially written by Michael Hipp, see also AUTHORS and README.

 *  

 * This library is free software; you can redistribute it and/or

 * modify it under the terms of the GNU Library General Public

 * License as published by the Free Software Foundation; either

 * version 2 of the License, or (at your option) any later version.

 *

 * This library 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

 * Library General Public License for more details.

 *

 * You should have received a copy of the GNU Library General Public

 * License along with this library; if not, write to the

 * Free Software Foundation, Inc., 59 Temple Place - Suite 330,

 * Boston, MA 02111-1307, USA.

 */

/* $Id: tabinit.c,v 1.17 2017/09/06 15:07:30 robert Exp $ */

#ifdef HAVE_CONFIG_H

# include <config.h>

#endif

#include <stdlib.h>

#include "tabinit.h"

#include "mpg123.h"

#ifdef WITH_DMALLOC

#include <dmalloc.h>

#endif

static int table_init_called = 0;

real    decwin[512 + 32];

static real cos64[16], cos32[8], cos16[4], cos8[2], cos4[1];

real   *pnts[] = { cos64, cos32, cos16, cos8, cos4 };

/* *INDENT-OFF* */

static const double dewin[512] = {

   0.000000000,-0.000015259,-0.000015259,-0.000015259,

  -0.000015259,-0.000015259,-0.000015259,-0.000030518,

  -0.000030518,-0.000030518,-0.000030518,-0.000045776,

  -0.000045776,-0.000061035,-0.000061035,-0.000076294,

  -0.000076294,-0.000091553,-0.000106812,-0.000106812,

  -0.000122070,-0.000137329,-0.000152588,-0.000167847,

  -0.000198364,-0.000213623,-0.000244141,-0.000259399,

  -0.000289917,-0.000320435,-0.000366211,-0.000396729,

  -0.000442505,-0.000473022,-0.000534058,-0.000579834,

  -0.000625610,-0.000686646,-0.000747681,-0.000808716,

  -0.000885010,-0.000961304,-0.001037598,-0.001113892,

  -0.001205444,-0.001296997,-0.001388550,-0.001480103,

  -0.001586914,-0.001693726,-0.001785278,-0.001907349,

  -0.002014160,-0.002120972,-0.002243042,-0.002349854,

  -0.002456665,-0.002578735,-0.002685547,-0.002792358,

  -0.002899170,-0.002990723,-0.003082275,-0.003173828,

  -0.003250122,-0.003326416,-0.003387451,-0.003433228,

  -0.003463745,-0.003479004,-0.003479004,-0.003463745,

  -0.003417969,-0.003372192,-0.003280640,-0.003173828,

  -0.003051758,-0.002883911,-0.002700806,-0.002487183,

  -0.002227783,-0.001937866,-0.001617432,-0.001266479,

  -0.000869751,-0.000442505, 0.000030518, 0.000549316,

   0.001098633, 0.001693726, 0.002334595, 0.003005981,

   0.003723145, 0.004486084, 0.005294800, 0.006118774,

   0.007003784, 0.007919312, 0.008865356, 0.009841919,

   0.010848999, 0.011886597, 0.012939453, 0.014022827,

   0.015121460, 0.016235352, 0.017349243, 0.018463135,

   0.019577026, 0.020690918, 0.021789551, 0.022857666,

   0.023910522, 0.024932861, 0.025909424, 0.026840210,

   0.027725220, 0.028533936, 0.029281616, 0.029937744,

   0.030532837, 0.031005859, 0.031387329, 0.031661987,

   0.031814575, 0.031845093, 0.031738281, 0.031478882,

   0.031082153, 0.030517578, 0.029785156, 0.028884888,

   0.027801514, 0.026535034, 0.025085449, 0.023422241,

   0.021575928, 0.019531250, 0.017257690, 0.014801025,

   0.012115479, 0.009231567, 0.006134033, 0.002822876,

  -0.000686646,-0.004394531,-0.008316040,-0.012420654,

  -0.016708374,-0.021179199,-0.025817871,-0.030609131,

  -0.035552979,-0.040634155,-0.045837402,-0.051132202,

  -0.056533813,-0.061996460,-0.067520142,-0.073059082,

  -0.078628540,-0.084182739,-0.089706421,-0.095169067,

  -0.100540161,-0.105819702,-0.110946655,-0.115921021,

  -0.120697021,-0.125259399,-0.129562378,-0.133590698,

  -0.137298584,-0.140670776,-0.143676758,-0.146255493,

  -0.148422241,-0.150115967,-0.151306152,-0.151962280,

  -0.152069092,-0.151596069,-0.150497437,-0.148773193,

  -0.146362305,-0.143264771,-0.139450073,-0.134887695,

  -0.129577637,-0.123474121,-0.116577148,-0.108856201,

  -0.100311279,-0.090927124,-0.080688477,-0.069595337,

  -0.057617187,-0.044784546,-0.031082153,-0.016510010,

  -0.001068115, 0.015228271, 0.032379150, 0.050354004,

   0.069168091, 0.088775635, 0.109161377, 0.130310059,

   0.152206421, 0.174789429, 0.198059082, 0.221984863,

   0.246505737, 0.271591187, 0.297210693, 0.323318481,

   0.349868774, 0.376800537, 0.404083252, 0.431655884,

   0.459472656, 0.487472534, 0.515609741, 0.543823242,

   0.572036743, 0.600219727, 0.628295898, 0.656219482,

   0.683914185, 0.711318970, 0.738372803, 0.765029907,

   0.791213989, 0.816864014, 0.841949463, 0.866363525,

   0.890090942, 0.913055420, 0.935195923, 0.956481934,

   0.976852417, 0.996246338, 1.014617920, 1.031936646,

   1.048156738, 1.063217163, 1.077117920, 1.089782715,

   1.101211548, 1.111373901, 1.120223999, 1.127746582,

   1.133926392, 1.138763428, 1.142211914, 1.144287109,

   1.144989014

};

/* *INDENT-ON* */

void

make_decode_tables(long scaleval)

{

    int     i, j, k, kr, divv;

    real   *table, *costab;

    if (table_init_called)

        return;

    table_init_called = 1;

    for (i = 0; i < 5; i++) {

        kr = 0x10 >> i;

        divv = 0x40 >> i;

        costab = pnts[i];

        for (k = 0; k < kr; k++)

            costab[k] = (real) (1.0 / (2.0 * cos(M_PI * ((double) k * 2.0 + 1.0) / (double) divv)));

    }

    table = decwin;

    scaleval = -scaleval;

    for (i = 0, j = 0; i < 256; i++, j++, table += 32) {

        if (table < decwin + 512 + 16)

            table[16] = table[0] = (real) (dewin[j] * scaleval);

        if (i % 32 == 31)

            table -= 1023;

        if (i % 64 == 63)

            scaleval = -scaleval;

    }

    for ( /* i=256 */ ; i < 512; i++, j--, table += 32) {

        if (table < decwin + 512 + 16)

            table[16] = table[0] = (real) (dewin[j] * scaleval);

        if (i % 32 == 31)

            table -= 1023;

        if (i % 64 == 63)

            scaleval = -scaleval;

    }

}