void pvmp3_split(int32 *vect)
{
int32 i;
const int32 *pt_cosTerms = &CosTable_dct32[15];
int32 *pt_vect = vect;
int32 *pt_vect_2 = pt_vect - 1;
for (i = 3; i != 0; i--)
{
int32 tmp2 = *(pt_vect);
int32 tmp1 = *(pt_vect_2);
int32 cosx = *(pt_cosTerms--);
*(pt_vect_2--) = (tmp1 + tmp2);
*(pt_vect++) = fxp_mul32_Q27((tmp1 - tmp2), cosx);
tmp2 = *(pt_vect);
tmp1 = *(pt_vect_2);
cosx = *(pt_cosTerms--);
*(pt_vect_2--) = (tmp1 + tmp2);
*(pt_vect++) = fxp_mul32_Q27((tmp1 - tmp2), cosx);
}
for (i = 5; i != 0; i--)
{
int32 tmp2 = *(pt_vect);
int32 tmp1 = *(pt_vect_2);
int32 cosx = *(pt_cosTerms--);
*(pt_vect_2--) = (tmp1 + tmp2);
*(pt_vect++) = fxp_mul32_Q32((tmp1 - tmp2) << 1, cosx);
tmp2 = *(pt_vect);
tmp1 = *(pt_vect_2);
cosx = *(pt_cosTerms--);
*(pt_vect_2--) = (tmp1 + tmp2);
*(pt_vect++) = fxp_mul32_Q32((tmp1 - tmp2) << 1, cosx);
}
}
void pvmp3_mdct_18(int32 vec[], int32 *history, const int32 *window)
{
int32 i;
int32 tmp;
int32 tmp1;
int32 tmp2;
int32 tmp3;
int32 tmp4;
const int32 *pt_cos_split = cosTerms_dct18;
const int32 *pt_cos = cosTerms_1_ov_cos_phi;
const int32 *pt_cos_x = &cosTerms_1_ov_cos_phi[17];
int32 *pt_vec = vec;
int32 *pt_vec_o = &vec[17];
for (i = 9; i != 0; i--)
{
tmp = *(pt_vec);
tmp1 = *(pt_vec_o);
tmp = fxp_mul32_Q32(tmp << 1, *(pt_cos++));
tmp1 = fxp_mul32_Q27(tmp1, *(pt_cos_x--));
*(pt_vec++) = tmp + tmp1 ;
*(pt_vec_o--) = fxp_mul32_Q28((tmp - tmp1), *(pt_cos_split++));
}
pvmp3_dct_9(vec); // Even terms
pvmp3_dct_9(&vec[9]); // Odd terms
tmp3 = vec[16]; //
vec[16] = vec[ 8];
tmp4 = vec[14]; //
vec[14] = vec[ 7];
tmp = vec[12];
vec[12] = vec[ 6];
tmp2 = vec[10]; // vec[10]
vec[10] = vec[ 5];
vec[ 8] = vec[ 4];
vec[ 6] = vec[ 3];
vec[ 4] = vec[ 2];
vec[ 2] = vec[ 1];
vec[ 1] = vec[ 9] - tmp2; // vec[9] + vec[10]
vec[ 3] = vec[11] - tmp2;
vec[ 5] = vec[11] - tmp;
vec[ 7] = vec[13] - tmp;
vec[ 9] = vec[13] - tmp4;
vec[11] = vec[15] - tmp4;
vec[13] = vec[15] - tmp3;
vec[15] = vec[17] - tmp3;
/* overlap and add */
tmp2 = vec[0];
tmp3 = vec[9];
for (i = 0; i < 6; i++)
{
tmp = history[ i];
tmp4 = vec[i+10];
vec[i+10] = tmp3 + tmp4;
tmp1 = vec[i+1];
vec[ i] = fxp_mac32_Q32(tmp, (vec[i+10]), window[ i]);
tmp3 = tmp4;
history[i ] = -(tmp2 + tmp1);
tmp2 = tmp1;
}
tmp = history[ 6];
tmp4 = vec[16];
vec[16] = tmp3 + tmp4;
tmp1 = vec[7];
vec[ 6] = fxp_mac32_Q32(tmp, vec[16] << 1, window[ i]);
tmp = history[ 7];
history[6] = -(tmp2 + tmp1);
history[7] = -(tmp1 + vec[8]);
tmp1 = history[ 8];
tmp4 = vec[17] + tmp4;
vec[ 7] = fxp_mac32_Q32(tmp, tmp4 << 1, window[ 7]);
history[8] = -(vec[8] + vec[9]);
vec[ 8] = fxp_mac32_Q32(tmp1, vec[17] << 1, window[ 8]);
tmp = history[9];
tmp1 = history[17];
tmp2 = history[16];
vec[ 9] = fxp_mac32_Q32(tmp, vec[17] << 1, window[ 9]);
vec[17] = fxp_mac32_Q32(tmp1, vec[10] << 1, window[17]);
vec[10] = -vec[ 16];
vec[16] = fxp_mac32_Q32(tmp2, vec[11] << 1, window[16]);
tmp1 = history[15];
tmp2 = history[14];
vec[11] = -vec[ 15];
vec[15] = fxp_mac32_Q32(tmp1, vec[12] << 1, window[15]);
vec[12] = -vec[ 14];
vec[14] = fxp_mac32_Q32(tmp2, vec[13] << 1, window[14]);
tmp = history[13];
tmp1 = history[12];
tmp2 = history[11];
tmp3 = history[10];
vec[13] = fxp_mac32_Q32(tmp, vec[12] << 1, window[13]);
vec[12] = fxp_mac32_Q32(tmp1, vec[11] << 1, window[12]);
vec[11] = fxp_mac32_Q32(tmp2, vec[10] << 1, window[11]);
vec[10] = fxp_mac32_Q32(tmp3, tmp4 << 1, window[10]);
/* next iteration overlap */
tmp1 = history[ 8];
tmp3 = history[ 7];
tmp2 = history[ 1];
tmp = history[ 0];
tmp1 <<= 1;
tmp3 <<= 1;
history[ 0] = fxp_mul32_Q32(tmp1, window[18]);
history[17] = fxp_mul32_Q32(tmp1, window[35]);
history[ 1] = fxp_mul32_Q32(tmp3, window[19]);
history[16] = fxp_mul32_Q32(tmp3, window[34]);
tmp2 <<= 1;
tmp <<= 1;
history[ 7] = fxp_mul32_Q32(tmp2, window[25]);
history[10] = fxp_mul32_Q32(tmp2, window[28]);
history[ 8] = fxp_mul32_Q32(tmp, window[26]);
history[ 9] = fxp_mul32_Q32(tmp, window[27]);
tmp1 = history[ 6];
tmp3 = history[ 5];
tmp4 = history[ 4];
tmp2 = history[ 3];
tmp = history[ 2];
tmp1 <<= 1;
tmp3 <<= 1;
tmp4 <<= 1;
history[ 2] = fxp_mul32_Q32(tmp1, window[20]);
history[15] = fxp_mul32_Q32(tmp1, window[33]);
history[ 3] = fxp_mul32_Q32(tmp3, window[21]);
history[14] = fxp_mul32_Q32(tmp3, window[32]);
history[ 4] = fxp_mul32_Q32(tmp4, window[22]);
history[13] = fxp_mul32_Q32(tmp4, window[31]);
tmp2 <<= 1;
tmp <<= 1;
history[ 5] = fxp_mul32_Q32(tmp2, window[23]);
history[12] = fxp_mul32_Q32(tmp2, window[30]);
history[ 6] = fxp_mul32_Q32(tmp, window[24]);
history[11] = fxp_mul32_Q32(tmp, window[29]);
}
