[Kst] extragear/graphics/kst/src/libkstmath

Duncan Hanson duncan.hanson at gmail.com
Tue Jun 13 23:40:40 CEST 2006


SVN commit 551186 by dhanson:

BUG:124158 sqrt(sW/len) for all apodizations.

 M  +9 -1      kstpsdgenerator.cpp  


--- trunk/extragear/graphics/kst/src/libkstmath/kstpsdgenerator.cpp #551185:551186
@@ -126,6 +126,7 @@
         w[i] = 1.0 - fabs(x) / a;
         sW += w[i] * w[i];
       }
+      sW = sqrt(sW / double(len));
       for (int i = 0; i < len; ++i) {
         w[i] /= sW;
       }
@@ -137,6 +138,7 @@
         w[i] = 0.42 + 0.5 * cos(M_PI * x / a) + 0.08 * cos(2 * M_PI * x/a);
         sW += w[i] * w[i];
       }
+      sW = sqrt(sW / double(len));
       for (int i = 0; i < len; ++i) {
         w[i] /= sW;
       }
@@ -148,6 +150,7 @@
         w[i] = pow(1.0 - (x * x) / (a * a), 2.0);
         sW += w[i] * w[i];
       }
+      sW = sqrt(sW / double(len));
       for (int i = 0; i < len; ++i) {
         w[i] /= sW;
       }
@@ -159,6 +162,7 @@
         w[i] = cos(M_PI * x / (2.0 * a));
         sW += w[i] * w[i];
       }
+      sW = sqrt(sW / double(len));
       for (int i = 0; i < len; ++i) {
         w[i] /= sW;
       }
@@ -169,6 +173,7 @@
         x = i - a;
         w[i] = exp(-1.0 * x * x/(2.0 * gaussianSigma * gaussianSigma));
       }
+      sW = sqrt(sW / double(len));
       for (int i = 0; i < len; ++i) {
         w[i] /= sW;
       }
@@ -179,7 +184,8 @@
         x = i - a;
         w[i] = 0.53836 + 0.46164 * cos(M_PI * x / a);
         sW += w[i] * w[i];
-      }  
+      }
+      sW = sqrt(sW / double(len));
       for (int i = 0; i < len; ++i) {
         w[i] /= sW;
       }
@@ -191,6 +197,7 @@
         w[i] = pow(cos(M_PI * x/(2.0 * a)), 2.0);
         sW += w[i] * w[i];
       }
+      sW = sqrt(sW / double(len));
       for (int i = 0; i < len; ++i) {
         w[i] /= sW;
       }
@@ -202,6 +209,7 @@
         w[i] = 1.0 - x * x / (a * a);
         sW += w[i] * w[i];
       }
+      sW = sqrt(sW / double(len));
       for (int i = 0; i < len; ++i) {
         w[i] /= sW;
       }


More information about the Kst mailing list