comparison Plugins/Input/mpg123/dct64_i386.c @ 970:0f294f2b0a9b trunk

[svn] - integer-mode decoding accuracy improvements

969:971bd2b03ff5

 */

#include "mpg123.h"

static void
mpg123_dct64_1(real * out0, real * out1, real * b1, real * b2, real * samples)
{

    {
        register real *costab = mpg123_pnts[0];

        b1[0x00] = samples[0x00] + samples[0x1F];
        b1[0x1F] = (samples[0x00] - samples[0x1F]) * costab[0x0];

        b1[0x01] = samples[0x01] + samples[0x1E];

        b1[0x0F] = samples[0x0F] + samples[0x10];
        b1[0x10] = (samples[0x0F] - samples[0x10]) * costab[0xF];
    }

    {
        register real *costab = mpg123_pnts[1];

        b2[0x00] = b1[0x00] + b1[0x0F];
        b2[0x0F] = (b1[0x00] - b1[0x0F]) * costab[0];
        b2[0x01] = b1[0x01] + b1[0x0E];
        b2[0x0E] = (b1[0x01] - b1[0x0E]) * costab[1];

        b2[0x17] = b1[0x17] + b1[0x18];
        b2[0x18] = (b1[0x18] - b1[0x17]) * costab[7];
    }

    {
        register real *costab = mpg123_pnts[2];

        b1[0x00] = b2[0x00] + b2[0x07];
        b1[0x07] = (b2[0x00] - b2[0x07]) * costab[0];
        b1[0x01] = b2[0x01] + b2[0x06];
        b1[0x06] = (b2[0x01] - b2[0x06]) * costab[1];

        b1[0x1B] = b2[0x1B] + b2[0x1C];
        b1[0x1C] = (b2[0x1C] - b2[0x1B]) * costab[3];
    }

    {
        register real const cos0 = mpg123_pnts[3][0];
        register real const cos1 = mpg123_pnts[3][1];

        b2[0x00] = b1[0x00] + b1[0x03];
        b2[0x03] = (b1[0x00] - b1[0x03]) * cos0;
        b2[0x01] = b1[0x01] + b1[0x02];
        b2[0x02] = (b1[0x01] - b1[0x02]) * cos1;

        b2[0x1D] = b1[0x1D] + b1[0x1E];
        b2[0x1E] = (b1[0x1E] - b1[0x1D]) * cos1;
    }

    {
        register real const cos0 = mpg123_pnts[4][0];

        b1[0x00] = b2[0x00] + b2[0x01];
        b1[0x01] = (b2[0x00] - b2[0x01]) * cos0;
        b1[0x02] = b2[0x02] + b2[0x03];
        b1[0x03] = (b2[0x03] - b2[0x02]) * cos0;

/*
 * the call via mpg123_dct64 is a trick to force GCC to use
 * (new) registers for the b1,b2 pointer to the bufs[xx] field
 */
void
mpg123_dct64(real * a, real * b, real * c)
{
    real bufs[0x40];

    mpg123_dct64_1(a, b, bufs, bufs + 0x20, c);
}

970:0f294f2b0a9b

 */

#include "mpg123.h"

static void
mpg123_dct64_1(mpgdec_real * out0, mpgdec_real * out1, mpgdec_real * b1, mpgdec_real * b2, mpgdec_real * samples)
{

    {
        register mpgdec_real *costab = mpg123_pnts[0];

        b1[0x00] = samples[0x00] + samples[0x1F];
        b1[0x1F] = (samples[0x00] - samples[0x1F]) * costab[0x0];

        b1[0x01] = samples[0x01] + samples[0x1E];

        b1[0x0F] = samples[0x0F] + samples[0x10];
        b1[0x10] = (samples[0x0F] - samples[0x10]) * costab[0xF];
    }

    {
        register mpgdec_real *costab = mpg123_pnts[1];

        b2[0x00] = b1[0x00] + b1[0x0F];
        b2[0x0F] = (b1[0x00] - b1[0x0F]) * costab[0];
        b2[0x01] = b1[0x01] + b1[0x0E];
        b2[0x0E] = (b1[0x01] - b1[0x0E]) * costab[1];

        b2[0x17] = b1[0x17] + b1[0x18];
        b2[0x18] = (b1[0x18] - b1[0x17]) * costab[7];
    }

    {
        register mpgdec_real *costab = mpg123_pnts[2];

        b1[0x00] = b2[0x00] + b2[0x07];
        b1[0x07] = (b2[0x00] - b2[0x07]) * costab[0];
        b1[0x01] = b2[0x01] + b2[0x06];
        b1[0x06] = (b2[0x01] - b2[0x06]) * costab[1];

        b1[0x1B] = b2[0x1B] + b2[0x1C];
        b1[0x1C] = (b2[0x1C] - b2[0x1B]) * costab[3];
    }

    {
        register mpgdec_real const cos0 = mpg123_pnts[3][0];
        register mpgdec_real const cos1 = mpg123_pnts[3][1];

        b2[0x00] = b1[0x00] + b1[0x03];
        b2[0x03] = (b1[0x00] - b1[0x03]) * cos0;
        b2[0x01] = b1[0x01] + b1[0x02];
        b2[0x02] = (b1[0x01] - b1[0x02]) * cos1;

        b2[0x1D] = b1[0x1D] + b1[0x1E];
        b2[0x1E] = (b1[0x1E] - b1[0x1D]) * cos1;
    }

    {
        register mpgdec_real const cos0 = mpg123_pnts[4][0];

        b1[0x00] = b2[0x00] + b2[0x01];
        b1[0x01] = (b2[0x00] - b2[0x01]) * cos0;
        b1[0x02] = b2[0x02] + b2[0x03];
        b1[0x03] = (b2[0x03] - b2[0x02]) * cos0;

/*
 * the call via mpg123_dct64 is a trick to force GCC to use
 * (new) registers for the b1,b2 pointer to the bufs[xx] field
 */
void
mpg123_dct64(mpgdec_real * a, mpgdec_real * b, mpgdec_real * c)
{
    mpgdec_real bufs[0x40];

    mpg123_dct64_1(a, b, bufs, bufs + 0x20, c);
}