Index: trunk/psLib/src/imageops/psImageConvolve.c
===================================================================
--- trunk/psLib/src/imageops/psImageConvolve.c	(revision 15437)
+++ trunk/psLib/src/imageops/psImageConvolve.c	(revision 15958)
@@ -7,6 +7,6 @@
 /// @author Eugene Magnier, IfA
 ///
-/// @version $Revision: 1.60 $ $Name: not supported by cvs2svn $
-/// @date $Date: 2007-11-02 03:23:46 $
+/// @version $Revision: 1.61 $ $Name: not supported by cvs2svn $
+/// @date $Date: 2007-12-31 03:02:08 $
 ///
 /// Copyright 2004-2007 Institute for Astronomy, University of Hawaii
@@ -826,9 +826,15 @@
                 g += *sg;
             }
-            *vo = (g > 0) ? s / g : 0.0;
+            *vo = (g > 0.25) ? s / g : NAN;
         }
         memcpy(image->data.F32[j], calculation->data.F32, Nx*sizeof(psF32));
     }
     psFree(calculation);
+
+    // XXX test
+    if (0) {
+	psFree(gaussnorm);
+	return true;
+    }
 
     /** Smooth in Y direction  **/
@@ -855,6 +861,6 @@
             double g = gauss[n - j];
             for (int i = 0; i < Nx; i++, vi++, vo++, vm++, vs++) {
-                if (*vm & maskVal)
-                    continue;
+                if (*vm & maskVal) continue;
+                if (!isfinite(*vi)) continue;
                 *vo += *vi * g;
                 *vs += g;
@@ -865,5 +871,5 @@
         psF32 *vs = outsum->data.F32;
         for (int i = 0; i < Nx; i++, vo++, vs++) {
-            *vo = (*vs > 0) ? *vo / *vs : 0.0;
+            *vo = (*vs > 0.25) ? *vo / *vs : NAN;
         }