Я ищу простой алгоритм для быстрого DCT и IDCT матрицы

Вот мой расчет 1D FDCT и IFDCT с помощью FFT с одинаковой длиной:

 //--------------------------------------------------------------------------- void DFCTrr(double *dst,double *src,double *tmp,int n) { // exact normalized DCT II by N DFFT int i,j; double nn=n,a,da=(M_PI*(nn-0.5))/nn,a0,b0,a1,b1,m; for (j= 0,i=n-1;i>=0;i-=2,j++) dst[j]=src[i]; for (j=n-1,i=n-2;i>=0;i-=2,j--) dst[j]=src[i]; DFFTcr(tmp,dst,n); m=2.0*sqrt(2.0); for (a=0.0,j=0,i=0;i=0;i-=2,j++) dst[i]=src[j]-m; for (j=n-1,i=n-2;i>=0;i-=2,j--) dst[i]=src[j]-m; } //--------------------------------------------------------------------------- 

• dst – целевой вектор [n]
• src – исходный вектор [n]
• tmp – это временный вектор [2n]

Эти массивы не должны пересекаться !!! Это взято из classа трансформации мины, поэтому я, надеюсь, не забыл что-то копировать.

• XXXrr означает, что пункт назначения является реальным, а источник также является реальным доменом
• XXXrc означает, что назначение является реальным, а источник – сложным доменом
• XXXcr означает, что назначение является сложным, а источник – это реальный домен

Все данные представляют собой double массивы, для первого числа сложных доменов – Реальная и вторая мнимая часть, поэтому массив имеет размер 2N . Обе функции используют FFT и iFFT, если вам нужен код и для них комментировать меня. Просто, чтобы быть уверенным, я добавил не быструю их реализацию ниже. Гораздо проще скопировать это, потому что быстрые используют слишком большую иерархию classов преобразования

медленная реализация DFT, iDFT для тестирования:

 //--------------------------------------------------------------------------- void transform::DFTcr(double *dst,double *src,int n) { int i,j; double a,b,a0,_n,q,qq,dq; dq=+2.0*M_PI/double(n); _n=2.0/double(n); for (q=0.0,j=0;j