gegl r2557 - branches/branch2_zhangjb/operations/frequency/tools
- From: zhangjb svn gnome org
- To: svn-commits-list gnome org
- Subject: gegl r2557 - branches/branch2_zhangjb/operations/frequency/tools
- Date: Sun, 3 Aug 2008 02:05:35 +0000 (UTC)
Author: zhangjb
Date: Sun Aug 3 02:05:35 2008
New Revision: 2557
URL: http://svn.gnome.org/viewvc/gegl?rev=2557&view=rev
Log:
Modified:
branches/branch2_zhangjb/operations/frequency/tools/filters.c
Modified: branches/branch2_zhangjb/operations/frequency/tools/filters.c
==============================================================================
--- branches/branch2_zhangjb/operations/frequency/tools/filters.c (original)
+++ branches/branch2_zhangjb/operations/frequency/tools/filters.c Sun Aug 3 02:05:35 2008
@@ -35,33 +35,31 @@
gint index,h_index;
gint max_x = FFT_HALF(width);
-#include <stdio.h>
-
- for(x=0;x<max_x;x++)
- {
- for(y=0;y<height/2;y++)
- {
- index = y*max_x+x;
- h_index = (height/2-y-1)*max_x+width/2-x-1;
-
- Yr= Xr[index]*Hr[h_index] - Xi[index]*Hi[h_index];
- Yi= Xi[index]*Hr[h_index] + Xr[index]*Hi[h_index];
- Xr[index] = Yr;
- Xi[index] = Yi;
- }
-
- for(y=height/2;y<height;y++)
- {
- index = (y*max_x)+x;
- h_index = (3*height/2-y-1)*max_x+width/2-x-1;
-
- Yr= Xr[index]*Hr[h_index] - Xi[index]*Hi[h_index];
- Yi= Xi[index]*Hr[h_index] + Xr[index]*Hi[h_index];
- Xr[index] = Yr;
- Xi[index] = Yi;
- }
- }
- return TRUE;
+ for(x=0;x<max_x;x++)
+ {
+ for(y=0;y<height/2;y++)
+ {
+ index = y*max_x+x;
+ h_index = (height/2-y-1)*max_x+width/2-x-1;
+
+ Yr= Xr[index]*Hr[h_index] - Xi[index]*Hi[h_index];
+ Yi= Xi[index]*Hr[h_index] + Xr[index]*Hi[h_index];
+ Xr[index] = Yr;
+ Xi[index] = Yi;
+ }
+
+ for(y=height/2;y<height;y++)
+ {
+ index = (y*max_x)+x;
+ h_index = (3*height/2-y-1)*max_x+width/2-x-1;
+
+ Yr= Xr[index]*Hr[h_index] - Xi[index]*Hi[h_index];
+ Yi= Xi[index]*Hr[h_index] + Xr[index]*Hi[h_index];
+ Xr[index] = Yr;
+ Xi[index] = Yi;
+ }
+ }
+ return TRUE;
}
gboolean
@@ -77,8 +75,8 @@
{
index = ELEM_ID_HALF_MATRIX(x, y, width);
Hi[index] = 0;
- Hr[index] = exp( -((gdouble)(x+1-width/2)*(x+1-width/2)
- +(y+1-height/2)*(y+1-height/2))/(2*cutoff*cutoff) );
+ Hr[index] = exp( -((gdouble)(x+1-width/2)*(x+1-width/2)
+ +(y+1-height/2)*(y+1-height/2))/(2*cutoff*cutoff) );
}
}
[
Date Prev][
Date Next] [
Thread Prev][
Thread Next]
[
Thread Index]
[
Date Index]
[
Author Index]