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]); }