Index: trunk/psLib/src/astro/psCoord.c
===================================================================
--- trunk/psLib/src/astro/psCoord.c	(revision 6230)
+++ trunk/psLib/src/astro/psCoord.c	(revision 6251)
@@ -10,6 +10,6 @@
 *  @author GLG, MHPCC
 *
-*  @version $Revision: 1.102 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2006-01-28 01:31:44 $
+*  @version $Revision: 1.103 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2006-01-30 20:28:33 $
 *
 *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -1085,6 +1085,4 @@
     PS_ASSERT_PTR_NON_NULL(input, NULL);
     PS_ASSERT_PTR_NON_NULL(inToOut, NULL);
-    PS_ASSERT_PTR_NON_NULL(inToOut->x, NULL);
-    PS_ASSERT_PTR_NON_NULL(inToOut->y, NULL);
     if (out == NULL) {
         //XXX: Should the length (nalloc) be 1 and append be used everytime a pixel is added?
@@ -1104,8 +1102,8 @@
         deriv = psPlaneTransformDeriv(deriv, inToOut, coord);
         fxnVal = psPlaneTransformApply(fxnVal, inToOut, coord);
-        if (fabs(fxnVal->x - coord->x) < fabs(deriv->x) &&
-                fabs(fxnVal->y - coord->y) < fabs(deriv->y)) {
-            int x = (int)(round(fabs(deriv->x)));
-            int y = (int)(round(fabs(deriv->y)));
+        if (fabs(fxnVal->x - coord->x) <= fabs(deriv->x) &&
+                fabs(fxnVal->y - coord->y) <= fabs(deriv->y)) {
+            int x = (int)(ceil(fabs(deriv->x)));
+            int y = (int)(ceil(fabs(deriv->y)));
             for (int j = -x; j <= x; j++) {
                 for (int k = -y; k <= y; k++) {