求救 缩减多重for循环

程序需要优化部分如下：
for ( k = 0; k < M ; k ++)
{
for ( l = 0; l < N ; l ++)
{
for ( i = 0; i < M ; i ++)
{
for ( j = 0; j < N ; j ++)
{
gray1 = fre[j * step2 + i] ;
gray2 = fure[j * step2 + i] ;
co = cos ( dPi * ( k * i / M + l * j / N));
si = sin ( dPi * ( k * i / M + l * j / N));
re = gray1 * co - gray2 * si;
ure= gray1 * si + gray2 * co;
renum += re;
urenum += ure;
}
}
printf("%d",k);
c[ l * step3 + k] = (uchar)(sqrt(renum * renum + urenum
* urenum));
renum = 0;
urenum = 0;
}
}
在取值的时候M,N都要取到300-400之间，结果运行的时候程序需要半个小时以上才能完成，求救高人帮忙，看程序是否可以优化，那怕优化后10分钟运行出来也可以接受啊～ 小弟在此谢谢了
补充一下 这个是对输入图像进行傅立叶变换，其中输入就是图像的X,Y轴的像素点，有一种快速傅立叶变换办法，可是时间仓促来不及学了，能不能在这个算法的基础上改变循环？

C程序for多重循环 vb 多重循环 求救急 C语言for多重循环语句 c中的for多重循环问题,高手指教.谢谢 菜鸟,求救js for循环问题 我不明白pascal里面for的多重循环,请给我讲解一下 VB的多重循环 关于for循环，while循环 C语言 for的多重循环是怎么一回事啊 能给小弟详细讲解以下吗? 多重循环结构的运算？