static void imdct36(int *out, int *buf, int *in, int *win)

{

    int i, j, t0, t1, t2, t3, s0, s1, s2, s3;

    int tmp[18], *tmp1, *in1;



    for(i=17;i>=1;i--)

        in[i] += in[i-1];

    for(i=17;i>=3;i-=2)

        in[i] += in[i-2];



    for(j=0;j<2;j++) {

        tmp1 = tmp + j;

        in1 = in + j;

#if 0

//more accurate but slower

        int64_t t0, t1, t2, t3;

        t2 = in1[2*4] + in1[2*8] - in1[2*2];

        

        t3 = (in1[2*0] + (int64_t)(in1[2*6]>>1))<<32;

        t1 = in1[2*0] - in1[2*6];

        tmp1[ 6] = t1 - (t2>>1);

        tmp1[16] = t1 + t2;



        t0 = MUL64(2*(in1[2*2] + in1[2*4]),    C2);

        t1 = MUL64(   in1[2*4] - in1[2*8] , -2*C8);

        t2 = MUL64(2*(in1[2*2] + in1[2*8]),   -C4);

        

        tmp1[10] = (t3 - t0 - t2) >> 32;

        tmp1[ 2] = (t3 + t0 + t1) >> 32;

        tmp1[14] = (t3 + t2 - t1) >> 32;

        

        tmp1[ 4] = MULH(2*(in1[2*5] + in1[2*7] - in1[2*1]), -C3);

        t2 = MUL64(2*(in1[2*1] + in1[2*5]),    C1);

        t3 = MUL64(   in1[2*5] - in1[2*7] , -2*C7);

        t0 = MUL64(2*in1[2*3], C3);



        t1 = MUL64(2*(in1[2*1] + in1[2*7]),   -C5);



        tmp1[ 0] = (t2 + t3 + t0) >> 32;

        tmp1[12] = (t2 + t1 - t0) >> 32;

        tmp1[ 8] = (t3 - t1 - t0) >> 32;

#else

        t2 = in1[2*4] + in1[2*8] - in1[2*2];

        

        t3 = in1[2*0] + (in1[2*6]>>1);

        t1 = in1[2*0] - in1[2*6];

        tmp1[ 6] = t1 - (t2>>1);

        tmp1[16] = t1 + t2;



        t0 = MULH(2*(in1[2*2] + in1[2*4]),    C2);

        t1 = MULH(   in1[2*4] - in1[2*8] , -2*C8);

        t2 = MULH(2*(in1[2*2] + in1[2*8]),   -C4);

        

        tmp1[10] = t3 - t0 - t2;

        tmp1[ 2] = t3 + t0 + t1;

        tmp1[14] = t3 + t2 - t1;

        

        tmp1[ 4] = MULH(2*(in1[2*5] + in1[2*7] - in1[2*1]), -C3);

        t2 = MULH(2*(in1[2*1] + in1[2*5]),    C1);

        t3 = MULH(   in1[2*5] - in1[2*7] , -2*C7);

        t0 = MULH(2*in1[2*3], C3);



        t1 = MULH(2*(in1[2*1] + in1[2*7]),   -C5);



        tmp1[ 0] = t2 + t3 + t0;

        tmp1[12] = t2 + t1 - t0;

        tmp1[ 8] = t3 - t1 - t0;

#endif

    }



    i = 0;

    for(j=0;j<4;j++) {

        t0 = tmp[i];

        t1 = tmp[i + 2];

        s0 = t1 + t0;

        s2 = t1 - t0;



        t2 = tmp[i + 1];

        t3 = tmp[i + 3];

        s1 = MULL(t3 + t2, icos36[j]);

        s3 = MULL(t3 - t2, icos36[8 - j]);

        

        t0 = (s0 + s1) << 5;

        t1 = (s0 - s1) << 5;

        out[(9 + j)*SBLIMIT] =  MULH(t1, win[9 + j]) + buf[9 + j];

        out[(8 - j)*SBLIMIT] =  MULH(t1, win[8 - j]) + buf[8 - j];

        buf[9 + j] = MULH(t0, win[18 + 9 + j]);

        buf[8 - j] = MULH(t0, win[18 + 8 - j]);

        

        t0 = (s2 + s3) << 5;

        t1 = (s2 - s3) << 5;

        out[(9 + 8 - j)*SBLIMIT] =  MULH(t1, win[9 + 8 - j]) + buf[9 + 8 - j];

        out[(        j)*SBLIMIT] =  MULH(t1, win[        j]) + buf[        j];

        buf[9 + 8 - j] = MULH(t0, win[18 + 9 + 8 - j]);

        buf[      + j] = MULH(t0, win[18         + j]);

        i += 4;

    }



    s0 = tmp[16];

    s1 = MULL(tmp[17], icos36[4]);

    t0 = (s0 + s1) << 5;

    t1 = (s0 - s1) << 5;

    out[(9 + 4)*SBLIMIT] =  MULH(t1, win[9 + 4]) + buf[9 + 4];

    out[(8 - 4)*SBLIMIT] =  MULH(t1, win[8 - 4]) + buf[8 - 4];

    buf[9 + 4] = MULH(t0, win[18 + 9 + 4]);

    buf[8 - 4] = MULH(t0, win[18 + 8 - 4]);

}

static void imdct36(int *out, int *buf, int *in, int *win)

{

    int i, j, t0, t1, t2, t3, s0, s1, s2, s3;

    int tmp[18], *tmp1, *in1;



    for(i=17;i>=1;i--)

        in[i] += in[i-1];

    for(i=17;i>=3;i-=2)

        in[i] += in[i-2];



    for(j=0;j<2;j++) {

        tmp1 = tmp + j;

        in1 = in + j;

#if 0



        int64_t t0, t1, t2, t3;

        t2 = in1[2*4] + in1[2*8] - in1[2*2];

        

        t3 = (in1[2*0] + (int64_t)(in1[2*6]>>1))<<32;

        t1 = in1[2*0] - in1[2*6];

        tmp1[ 6] = t1 - (t2>>1);

        tmp1[16] = t1 + t2;



        t0 = MUL64(2*(in1[2*2] + in1[2*4]),    C2);

        t1 = MUL64(   in1[2*4] - in1[2*8] , -2*C8);

        t2 = MUL64(2*(in1[2*2] + in1[2*8]),   -C4);

        

        tmp1[10] = (t3 - t0 - t2) >> 32;

        tmp1[ 2] = (t3 + t0 + t1) >> 32;

        tmp1[14] = (t3 + t2 - t1) >> 32;

        

        tmp1[ 4] = MULH(2*(in1[2*5] + in1[2*7] - in1[2*1]), -C3);

        t2 = MUL64(2*(in1[2*1] + in1[2*5]),    C1);

        t3 = MUL64(   in1[2*5] - in1[2*7] , -2*C7);

        t0 = MUL64(2*in1[2*3], C3);



        t1 = MUL64(2*(in1[2*1] + in1[2*7]),   -C5);



        tmp1[ 0] = (t2 + t3 + t0) >> 32;

        tmp1[12] = (t2 + t1 - t0) >> 32;

        tmp1[ 8] = (t3 - t1 - t0) >> 32;

#else

        t2 = in1[2*4] + in1[2*8] - in1[2*2];

        

        t3 = in1[2*0] + (in1[2*6]>>1);

        t1 = in1[2*0] - in1[2*6];

        tmp1[ 6] = t1 - (t2>>1);

        tmp1[16] = t1 + t2;



        t0 = MULH(2*(in1[2*2] + in1[2*4]),    C2);

        t1 = MULH(   in1[2*4] - in1[2*8] , -2*C8);

        t2 = MULH(2*(in1[2*2] + in1[2*8]),   -C4);

        

        tmp1[10] = t3 - t0 - t2;

        tmp1[ 2] = t3 + t0 + t1;

        tmp1[14] = t3 + t2 - t1;

        

        tmp1[ 4] = MULH(2*(in1[2*5] + in1[2*7] - in1[2*1]), -C3);

        t2 = MULH(2*(in1[2*1] + in1[2*5]),    C1);

        t3 = MULH(   in1[2*5] - in1[2*7] , -2*C7);

        t0 = MULH(2*in1[2*3], C3);



        t1 = MULH(2*(in1[2*1] + in1[2*7]),   -C5);



        tmp1[ 0] = t2 + t3 + t0;

        tmp1[12] = t2 + t1 - t0;

        tmp1[ 8] = t3 - t1 - t0;

#endif

    }



    i = 0;

    for(j=0;j<4;j++) {

        t0 = tmp[i];

        t1 = tmp[i + 2];

        s0 = t1 + t0;

        s2 = t1 - t0;



        t2 = tmp[i + 1];

        t3 = tmp[i + 3];

        s1 = MULL(t3 + t2, icos36[j]);

        s3 = MULL(t3 - t2, icos36[8 - j]);

        

        t0 = (s0 + s1) << 5;

        t1 = (s0 - s1) << 5;

        out[(9 + j)*SBLIMIT] =  MULH(t1, win[9 + j]) + buf[9 + j];

        out[(8 - j)*SBLIMIT] =  MULH(t1, win[8 - j]) + buf[8 - j];

        buf[9 + j] = MULH(t0, win[18 + 9 + j]);

        buf[8 - j] = MULH(t0, win[18 + 8 - j]);

        

        t0 = (s2 + s3) << 5;

        t1 = (s2 - s3) << 5;

        out[(9 + 8 - j)*SBLIMIT] =  MULH(t1, win[9 + 8 - j]) + buf[9 + 8 - j];

        out[(        j)*SBLIMIT] =  MULH(t1, win[        j]) + buf[        j];

        buf[9 + 8 - j] = MULH(t0, win[18 + 9 + 8 - j]);

        buf[      + j] = MULH(t0, win[18         + j]);

        i += 4;

    }



    s0 = tmp[16];

    s1 = MULL(tmp[17], icos36[4]);

    t0 = (s0 + s1) << 5;

    t1 = (s0 - s1) << 5;

    out[(9 + 4)*SBLIMIT] =  MULH(t1, win[9 + 4]) + buf[9 + 4];

    out[(8 - 4)*SBLIMIT] =  MULH(t1, win[8 - 4]) + buf[8 - 4];

    buf[9 + 4] = MULH(t0, win[18 + 9 + 4]);

    buf[8 - 4] = MULH(t0, win[18 + 8 - 4]);

}

        t0 = (s0 + s1) << 5;         t1 = (s0 - s1) << 5;         t0 = (s2 + s3) << 5;         t1 = (s2 - s3) << 5;     t0 = (s0 + s1) << 5;     t1 = (s0 - s1) << 5;
static void FUNC_0(int *VAR_0, int *VAR_1, int *VAR_2, int *VAR_3)

{

    int VAR_4, VAR_5, VAR_6, VAR_7, VAR_8, VAR_9, VAR_10, VAR_11, VAR_12, VAR_13;

    int VAR_14[18], *VAR_15, *VAR_16;



    for(VAR_4=17;VAR_4>=1;VAR_4--)

        VAR_2[VAR_4] += VAR_2[VAR_4-1];

    for(VAR_4=17;VAR_4>=3;VAR_4-=2)

        VAR_2[VAR_4] += VAR_2[VAR_4-2];



    for(VAR_5=0;VAR_5<2;VAR_5++) {

        VAR_15 = VAR_14 + VAR_5;

        VAR_16 = VAR_2 + VAR_5;

#if 0



        int64_t VAR_6, VAR_7, VAR_8, VAR_9;

        VAR_8 = VAR_16[2*4] + VAR_16[2*8] - VAR_16[2*2];

        

        VAR_9 = (VAR_16[2*0] + (int64_t)(VAR_16[2*6]>>1))<<32;

        VAR_7 = VAR_16[2*0] - VAR_16[2*6];

        VAR_15[ 6] = VAR_7 - (VAR_8>>1);

        VAR_15[16] = VAR_7 + VAR_8;



        VAR_6 = MUL64(2*(VAR_16[2*2] + VAR_16[2*4]),    C2);

        VAR_7 = MUL64(   VAR_16[2*4] - VAR_16[2*8] , -2*C8);

        VAR_8 = MUL64(2*(VAR_16[2*2] + VAR_16[2*8]),   -C4);

        

        VAR_15[10] = (VAR_9 - VAR_6 - VAR_8) >> 32;

        VAR_15[ 2] = (VAR_9 + VAR_6 + VAR_7) >> 32;

        VAR_15[14] = (VAR_9 + VAR_8 - VAR_7) >> 32;

        

        VAR_15[ 4] = MULH(2*(VAR_16[2*5] + VAR_16[2*7] - VAR_16[2*1]), -C3);

        VAR_8 = MUL64(2*(VAR_16[2*1] + VAR_16[2*5]),    C1);

        VAR_9 = MUL64(   VAR_16[2*5] - VAR_16[2*7] , -2*C7);

        VAR_6 = MUL64(2*VAR_16[2*3], C3);



        VAR_7 = MUL64(2*(VAR_16[2*1] + VAR_16[2*7]),   -C5);



        VAR_15[ 0] = (VAR_8 + VAR_9 + VAR_6) >> 32;

        VAR_15[12] = (VAR_8 + VAR_7 - VAR_6) >> 32;

        VAR_15[ 8] = (VAR_9 - VAR_7 - VAR_6) >> 32;

#else

        VAR_8 = VAR_16[2*4] + VAR_16[2*8] - VAR_16[2*2];

        

        VAR_9 = VAR_16[2*0] + (VAR_16[2*6]>>1);

        VAR_7 = VAR_16[2*0] - VAR_16[2*6];

        VAR_15[ 6] = VAR_7 - (VAR_8>>1);

        VAR_15[16] = VAR_7 + VAR_8;



        VAR_6 = MULH(2*(VAR_16[2*2] + VAR_16[2*4]),    C2);

        VAR_7 = MULH(   VAR_16[2*4] - VAR_16[2*8] , -2*C8);

        VAR_8 = MULH(2*(VAR_16[2*2] + VAR_16[2*8]),   -C4);

        

        VAR_15[10] = VAR_9 - VAR_6 - VAR_8;

        VAR_15[ 2] = VAR_9 + VAR_6 + VAR_7;

        VAR_15[14] = VAR_9 + VAR_8 - VAR_7;

        

        VAR_15[ 4] = MULH(2*(VAR_16[2*5] + VAR_16[2*7] - VAR_16[2*1]), -C3);

        VAR_8 = MULH(2*(VAR_16[2*1] + VAR_16[2*5]),    C1);

        VAR_9 = MULH(   VAR_16[2*5] - VAR_16[2*7] , -2*C7);

        VAR_6 = MULH(2*VAR_16[2*3], C3);



        VAR_7 = MULH(2*(VAR_16[2*1] + VAR_16[2*7]),   -C5);



        VAR_15[ 0] = VAR_8 + VAR_9 + VAR_6;

        VAR_15[12] = VAR_8 + VAR_7 - VAR_6;

        VAR_15[ 8] = VAR_9 - VAR_7 - VAR_6;

#endif

    }



    VAR_4 = 0;

    for(VAR_5=0;VAR_5<4;VAR_5++) {

        VAR_6 = VAR_14[VAR_4];

        VAR_7 = VAR_14[VAR_4 + 2];

        VAR_10 = VAR_7 + VAR_6;

        VAR_12 = VAR_7 - VAR_6;



        VAR_8 = VAR_14[VAR_4 + 1];

        VAR_9 = VAR_14[VAR_4 + 3];

        VAR_11 = MULL(VAR_9 + VAR_8, icos36[VAR_5]);

        VAR_13 = MULL(VAR_9 - VAR_8, icos36[8 - VAR_5]);

        

        VAR_6 = (VAR_10 + VAR_11) << 5;

        VAR_7 = (VAR_10 - VAR_11) << 5;

        VAR_0[(9 + VAR_5)*SBLIMIT] =  MULH(VAR_7, VAR_3[9 + VAR_5]) + VAR_1[9 + VAR_5];

        VAR_0[(8 - VAR_5)*SBLIMIT] =  MULH(VAR_7, VAR_3[8 - VAR_5]) + VAR_1[8 - VAR_5];

        VAR_1[9 + VAR_5] = MULH(VAR_6, VAR_3[18 + 9 + VAR_5]);

        VAR_1[8 - VAR_5] = MULH(VAR_6, VAR_3[18 + 8 - VAR_5]);

        

        VAR_6 = (VAR_12 + VAR_13) << 5;

        VAR_7 = (VAR_12 - VAR_13) << 5;

        VAR_0[(9 + 8 - VAR_5)*SBLIMIT] =  MULH(VAR_7, VAR_3[9 + 8 - VAR_5]) + VAR_1[9 + 8 - VAR_5];

        VAR_0[(        VAR_5)*SBLIMIT] =  MULH(VAR_7, VAR_3[        VAR_5]) + VAR_1[        VAR_5];

        VAR_1[9 + 8 - VAR_5] = MULH(VAR_6, VAR_3[18 + 9 + 8 - VAR_5]);

        VAR_1[      + VAR_5] = MULH(VAR_6, VAR_3[18         + VAR_5]);

        VAR_4 += 4;

    }



    VAR_10 = VAR_14[16];

    VAR_11 = MULL(VAR_14[17], icos36[4]);

    VAR_6 = (VAR_10 + VAR_11) << 5;

    VAR_7 = (VAR_10 - VAR_11) << 5;

    VAR_0[(9 + 4)*SBLIMIT] =  MULH(VAR_7, VAR_3[9 + 4]) + VAR_1[9 + 4];

    VAR_0[(8 - 4)*SBLIMIT] =  MULH(VAR_7, VAR_3[8 - 4]) + VAR_1[8 - 4];

    VAR_1[9 + 4] = MULH(VAR_6, VAR_3[18 + 9 + 4]);

    VAR_1[8 - 4] = MULH(VAR_6, VAR_3[18 + 8 - 4]);

}

static void FUNC_0(int *VAR_0, int *VAR_1, int *VAR_2, int *VAR_3)
{ int VAR_4, VAR_5, VAR_6, VAR_7, VAR_8, VAR_9, VAR_10, VAR_11, VAR_12, VAR_13; int VAR_14[18], *VAR_15, *VAR_16; for(VAR_4=17;VAR_4>=1;VAR_4--) VAR_2[VAR_4] += VAR_2[VAR_4-1]; for(VAR_4=17;VAR_4>=3;VAR_4-=2) VAR_2[VAR_4] += VAR_2[VAR_4-2]; for(VAR_5=0;VAR_5<2;VAR_5++) { VAR_15 = VAR_14 + VAR_5; VAR_16 = VAR_2 + VAR_5; #if 0
int64_t VAR_6, VAR_7, VAR_8, VAR_9; VAR_8 = VAR_16[2*4] + VAR_16[2*8] - VAR_16[2*2]; VAR_9 = (VAR_16[2*0] + (int64_t)(VAR_16[2*6]>>1))<<32; VAR_7 = VAR_16[2*0] - VAR_16[2*6]; VAR_15[ 6] = VAR_7 - (VAR_8>>1); VAR_15[16] = VAR_7 + VAR_8; VAR_6 = MUL64(2*(VAR_16[2*2] + VAR_16[2*4]),    C2); VAR_7 = MUL64(   VAR_16[2*4] - VAR_16[2*8] , -2*C8); VAR_8 = MUL64(2*(VAR_16[2*2] + VAR_16[2*8]),   -C4); VAR_15[10] = (VAR_9 - VAR_6 - VAR_8) >> 32; VAR_15[ 2] = (VAR_9 + VAR_6 + VAR_7) >> 32; VAR_15[14] = (VAR_9 + VAR_8 - VAR_7) >> 32; VAR_15[ 4] = MULH(2*(VAR_16[2*5] + VAR_16[2*7] - VAR_16[2*1]), -C3); VAR_8 = MUL64(2*(VAR_16[2*1] + VAR_16[2*5]),    C1); VAR_9 = MUL64(   VAR_16[2*5] - VAR_16[2*7] , -2*C7); VAR_6 = MUL64(2*VAR_16[2*3], C3); VAR_7 = MUL64(2*(VAR_16[2*1] + VAR_16[2*7]),   -C5); VAR_15[ 0] = (VAR_8 + VAR_9 + VAR_6) >> 32; VAR_15[12] = (VAR_8 + VAR_7 - VAR_6) >> 32; VAR_15[ 8] = (VAR_9 - VAR_7 - VAR_6) >> 32; #else
VAR_8 = VAR_16[2*4] + VAR_16[2*8] - VAR_16[2*2]; VAR_9 = VAR_16[2*0] + (VAR_16[2*6]>>1); VAR_7 = VAR_16[2*0] - VAR_16[2*6]; VAR_15[ 6] = VAR_7 - (VAR_8>>1); VAR_15[16] = VAR_7 + VAR_8; VAR_6 = MULH(2*(VAR_16[2*2] + VAR_16[2*4]),    C2); VAR_7 = MULH(   VAR_16[2*4] - VAR_16[2*8] , -2*C8); VAR_8 = MULH(2*(VAR_16[2*2] + VAR_16[2*8]),   -C4); VAR_15[10] = VAR_9 - VAR_6 - VAR_8; VAR_15[ 2] = VAR_9 + VAR_6 + VAR_7; VAR_15[14] = VAR_9 + VAR_8 - VAR_7; VAR_15[ 4] = MULH(2*(VAR_16[2*5] + VAR_16[2*7] - VAR_16[2*1]), -C3); VAR_8 = MULH(2*(VAR_16[2*1] + VAR_16[2*5]),    C1); VAR_9 = MULH(   VAR_16[2*5] - VAR_16[2*7] , -2*C7); VAR_6 = MULH(2*VAR_16[2*3], C3); VAR_7 = MULH(2*(VAR_16[2*1] + VAR_16[2*7]),   -C5); VAR_15[ 0] = VAR_8 + VAR_9 + VAR_6; VAR_15[12] = VAR_8 + VAR_7 - VAR_6; VAR_15[ 8] = VAR_9 - VAR_7 - VAR_6; #endif
} VAR_4 = 0; for(VAR_5=0;VAR_5<4;VAR_5++) { VAR_6 = VAR_14[VAR_4]; VAR_7 = VAR_14[VAR_4 + 2]; VAR_10 = VAR_7 + VAR_6; VAR_12 = VAR_7 - VAR_6; VAR_8 = VAR_14[VAR_4 + 1]; VAR_9 = VAR_14[VAR_4 + 3]; VAR_11 = MULL(VAR_9 + VAR_8, icos36[VAR_5]); VAR_13 = MULL(VAR_9 - VAR_8, icos36[8 - VAR_5]); VAR_6 = (VAR_10 + VAR_11) << 5; VAR_7 = (VAR_10 - VAR_11) << 5; VAR_0[(9 + VAR_5)*SBLIMIT] =  MULH(VAR_7, VAR_3[9 + VAR_5]) + VAR_1[9 + VAR_5]; VAR_0[(8 - VAR_5)*SBLIMIT] =  MULH(VAR_7, VAR_3[8 - VAR_5]) + VAR_1[8 - VAR_5]; VAR_1[9 + VAR_5] = MULH(VAR_6, VAR_3[18 + 9 + VAR_5]); VAR_1[8 - VAR_5] = MULH(VAR_6, VAR_3[18 + 8 - VAR_5]); VAR_6 = (VAR_12 + VAR_13) << 5; VAR_7 = (VAR_12 - VAR_13) << 5; VAR_0[(9 + 8 - VAR_5)*SBLIMIT] =  MULH(VAR_7, VAR_3[9 + 8 - VAR_5]) + VAR_1[9 + 8 - VAR_5]; VAR_0[(        VAR_5)*SBLIMIT] =  MULH(VAR_7, VAR_3[        VAR_5]) + VAR_1[        VAR_5]; VAR_1[9 + 8 - VAR_5] = MULH(VAR_6, VAR_3[18 + 9 + 8 - VAR_5]); VAR_1[      + VAR_5] = MULH(VAR_6, VAR_3[18         + VAR_5]); VAR_4 += 4; } VAR_10 = VAR_14[16]; VAR_11 = MULL(VAR_14[17], icos36[4]); VAR_6 = (VAR_10 + VAR_11) << 5; VAR_7 = (VAR_10 - VAR_11) << 5; VAR_0[(9 + 4)*SBLIMIT] =  MULH(VAR_7, VAR_3[9 + 4]) + VAR_1[9 + 4]; VAR_0[(8 - 4)*SBLIMIT] =  MULH(VAR_7, VAR_3[8 - 4]) + VAR_1[8 - 4]; VAR_1[9 + 4] = MULH(VAR_6, VAR_3[18 + 9 + 4]); VAR_1[8 - 4] = MULH(VAR_6, VAR_3[18 + 8 - 4]); }