for(v=0;v<Y;v++) {
//dprintf(fourier->outfd,"[fourier] (u=%d,v=%d)\n",v,x);
for(v=0;v<Y;v++) {
//dprintf(fourier->outfd,"[fourier] (u=%d,v=%d)\n",v,x);
- fourier->ftdata[off_f+x].r+=(cos(arg)*data[off_r+x].r-sin(arg)*data[off_r+x].i);
- fourier->ftdata[off_f+x].r+=(sin(arg)*data[off_r+x].r+cos(arg)*data[off_r+x].i);
+ fourier->ftdata[off_f].r+=(cos(arg)*data[off_r].r-sin(arg)*data[off_r].i);
+ fourier->ftdata[off_f].r+=(sin(arg)*data[off_r].r+cos(arg)*data[off_r].i);