113 void dct_III(FIXP_DBL *pDat, /*!< pointer to input/output */  in dct_III()  argument
139     cplxMultDiv2(&accu2, &accu1, pDat[L - i], pDat[i], sin_twiddle_L64[i*inc]);  in dct_III()
140     cplxMultDiv2(&accu4, &accu3, pDat[M+i], pDat[M-i], sin_twiddle_L64[(M-i)*inc]);  in dct_III()
155   xr     = fMultDiv2(pDat[M], sin_twiddle_L64[64/2].v.re );/* cos((PI/(2*L))*M); */  in dct_III()
156   tmp[0] = ((pDat[0]>>1) + xr)>>1;  in dct_III()
157   tmp[1] = ((pDat[0]>>1) - xr)>>1;  in dct_III()
159   cplxMultDiv2(&accu2, &accu1, pDat[L - (M/2)], pDat[M/2], sin_twiddle_L64[64/4]);  in dct_III()
175     *pDat++ = tmp1;  in dct_III()
176     *pDat++ = tmp3;  in dct_III()
177     *pDat++ = tmp2;  in dct_III()
178     *pDat++ = tmp4;  in dct_III()
186 void dct_II(FIXP_DBL *pDat, /*!< pointer to input/output */  in dct_II()  argument
206     FIXP_DBL *pdat  = &pDat[0];  in dct_II()
247       pDat[L - i] = accu4;  in dct_II()
248       pDat[i]     = accu3;  in dct_II()
251       pDat[M + i] = accu4;  in dct_II()
252       pDat[M - i] = accu3;  in dct_II()
257     pDat[L - (M/2)] = accu2;  in dct_II()
258     pDat[M/2]       = accu1;  in dct_II()
260     pDat[0] = (tmp[0]>>1)+(tmp[1]>>1);  in dct_II()
261     pDat[M] =  fMult(((tmp[0]>>1)-(tmp[1]>>1)), sin_twiddle_L64[64/2].v.re);/* cos((PI/(2*L))*M); */  in dct_II()
304 void dct_IV(FIXP_DBL *pDat,  in dct_IV()  argument
321     dct_IV_func1(M>>2, twiddle,  &pDat[0], &pDat[L-1]);  in dct_IV()
325     FIXP_DBL *RESTRICT pDat_0 = &pDat[0];  in dct_IV()
326     FIXP_DBL *RESTRICT pDat_1 = &pDat[L - 2];  in dct_IV()
355   fft(M, pDat, pDat_e);  in dct_IV()
360     dct_IV_func2(M>>2, sin_twiddle, &pDat[0], &pDat[L], sin_step);  in dct_IV()
364     FIXP_DBL *RESTRICT pDat_0 = &pDat[0];  in dct_IV()
365     FIXP_DBL *RESTRICT pDat_1 = &pDat[L - 2];  in dct_IV()
414 void dst_IV(FIXP_DBL *pDat,  in dst_IV()  argument
429     cplxMultDiv2(&tmp2, &tmp1, pDat[0], pDat[1], tab);  in dst_IV()
431     pDat[0] = tmp1;  in dst_IV()
432     pDat[1] = tmp2;  in dst_IV()
446     dst_IV_func1(M, twiddle, &pDat[0], &pDat[L]);  in dst_IV()
450     FIXP_DBL *RESTRICT pDat_0 = &pDat[0];  in dst_IV()
451     FIXP_DBL *RESTRICT pDat_1 = &pDat[L - 2];  in dst_IV()
481   fft(M, pDat, pDat_e);  in dst_IV()
485     dst_IV_func2(M>>2, sin_twiddle + sin_step, &pDat[0], &pDat[L - 1], sin_step);  in dst_IV()
494     pDat_0 = &pDat[0];  in dst_IV()
495     pDat_1 = &pDat[L - 2];  in dst_IV()