avoid a possible sign-extension for offsets/sizes in SGI
[imager.git] / gaussian.im
index 8489273..7e93769 100644 (file)
@@ -1,9 +1,10 @@
+#define IMAGER_NO_CONTEXT
 #include "imager.h"
 #include <math.h>
 
 static double
 gauss(int x, double std) {
-  return 1.0/(sqrt(2.0*PI)*std)*exp(-(double)(x)*(float)(x)/(2*std*std));
+  return 1.0/(sqrt(2.0*PI)*std)*exp(-(double)(x)*(double)(x)/(2*std*std));
 }
 
 /* Counters are as follows
@@ -18,14 +19,16 @@ gauss(int x, double std) {
 
 int
 i_gaussian(i_img *im, double stddev) {
-  int i,l,c,ch;
+  int i, c, ch;
+  i_img_dim x, y;
   double pc;
   double *coeff;
   double res[MAXCHANNELS];
   i_img *timg;
   int radius, diameter;
+  dIMCTXim(im);
 
-  mm_log((1,"i_gaussian(im %p, stdev %.2f)\n",im,stddev));
+  im_log((aIMCTX, 1,"i_gaussian(im %p, stdev %.2f)\n",im,stddev));
   i_clear_error();
 
   if (stddev <= 0) {
@@ -59,13 +62,13 @@ i_gaussian(i_img *im, double stddev) {
 
 #code im->bits <= 8
   IM_COLOR rcolor;
-  for(l=0;l<im->ysize;l++) {
-    for(i=0;i<im->xsize;i++) {
+  for(y = 0; y < im->ysize; y++) {
+    for(x = 0; x < im->xsize; x++) {
       pc=0.0;
       for(ch=0;ch<im->channels;ch++) 
        res[ch]=0; 
       for(c = 0;c < diameter; c++)
-       if (IM_GPIX(im,i+c-radius,l,&rcolor)!=-1) {
+       if (IM_GPIX(im,x+c-radius,y,&rcolor)!=-1) {
          for(ch=0;ch<im->channels;ch++) 
            res[ch]+= rcolor.channel[ch] * coeff[c];
          pc+=coeff[c];
@@ -74,17 +77,17 @@ i_gaussian(i_img *im, double stddev) {
        double value = res[ch] / pc;
        rcolor.channel[ch] = value > IM_SAMPLE_MAX ? IM_SAMPLE_MAX : IM_ROUND(value);
       }
-      IM_PPIX(timg,i,l,&rcolor);
+      IM_PPIX(timg, x, y, &rcolor);
     }
   }
   
-  for(l=0;l<im->xsize;l++) {
-    for(i=0;i<im->ysize;i++) {
+  for(x = 0;x < im->xsize; x++) {
+    for(y = 0; y < im->ysize; y++) {
       pc=0.0;
       for(ch=0; ch<im->channels; ch++)
        res[ch]=0; 
       for(c=0; c < diameter; c++)
-       if (IM_GPIX(timg,l,i+c-radius,&rcolor)!=-1) {
+       if (IM_GPIX(timg, x, y+c-radius, &rcolor)!=-1) {
          for(ch=0;ch<im->channels;ch++) 
            res[ch]+= rcolor.channel[ch] * coeff[c];
          pc+=coeff[c];
@@ -93,7 +96,7 @@ i_gaussian(i_img *im, double stddev) {
        double value = res[ch]/pc;
        rcolor.channel[ch] = value > IM_SAMPLE_MAX ? IM_SAMPLE_MAX : IM_ROUND(value);
       }
-      IM_PPIX(im,l,i,&rcolor);
+      IM_PPIX(im, x, y, &rcolor);
     }
   }
 #/code