The R Project SVN – R – Diff – /trunk/src/main/plot3d.c

 	      int *npt)
+{
     double c;
     if (z1 > z2 ) {
-	if (z2 > high || z1 < low)
+	if (z2 > high || z1 < low) return;
-	    return;
 	if (z1 < high) {
 	    x[*npt] = x1;
 	    y[*npt] = y1;
 	    z[*npt] = z1;
 	    ++*npt;
+	} else if (z1 == R_PosInf) {
+	    x[*npt] = x2;
+	    y[*npt] = y1;
-	}
+	    z[*npt] = z2;
-	else {
+	    ++*npt;
+	} else { /* z1 >= high, z2 in range */
 	    c = (z1 - high) / (z1 - z2);
 	    x[*npt] = x1 + c * (x2 - x1);
 	    y[*npt] = y1;
 	    z[*npt] = z1 + c * (z2 - z1);
 	    ++*npt;
+	}
-	if (z2 > low) {
+	if (z2 == R_NegInf) {
+	    x[*npt] = x1;
+	    y[*npt] = y1;
-	}
+	    z[*npt] = z1;
-	else {
+	    ++*npt;
+	} else if (z2 <= low) { /* and z1 in range */
 	    c = (z2 -low) / (z2 - z1);
 	    x[*npt] = x2 - c * (x2 - x1);
 	    y[*npt] = y1;
 	    z[*npt] = z2 - c * (z2 - z1);
 	    ++*npt;
+	}
-    }
-    else if (z1 < z2) {
+    } else if (z1 < z2) {
-	if (z2 < low || z1 > high)
+	if (z2 < low || z1 > high) return;
-	    return;
 	if (z1 > low) {
 	    x[*npt] = x1;
 	    y[*npt] = y1;
 	    z[*npt] = z1;
 	    ++*npt;
+	} else if (z1 == R_NegInf) {
+	    x[*npt] = x2;
+	    y[*npt] = y1;
+	    z[*npt] = z2;;
-	}
+	    ++*npt;
-	else {
+	} else { /* and z2 in range */
 	    c = (z1 - low) / (z1 - z2);
 	    x[*npt] = x1 + c * (x2 - x1);
 	    y[*npt] = y1;
 	    z[*npt] = z1 + c * (z2 - z1);
 	    ++*npt;
 	    x[*npt] = x2;
 	    y[*npt] = y2;
 	    z[*npt] = z2;
 	    ++*npt;
 #endif
+	} else if (z2 == R_PosInf) {
+	    x[*npt] = x1;
+	    y[*npt] = y1;
-	}
+	    z[*npt] = z1;
-	else {
+	    ++*npt;
+	} else { /* z2 high, z1 in range */
 	    c = (z2 - high) / (z2 - z1);
 	    x[*npt] = x2 - c * (x2 - x1);
 	    y[*npt] = y1;
 	    z[*npt] = z2 - c * (z2 - z1);
 	    ++*npt;
+	}
-    }
-    else {
+    } else {
 	if(low <= z1 && z1 <= high) {
 	    x[*npt] = x1;
 	    y[*npt] = y1;
 	    z[*npt] = z1;
 	    ++*npt;
 				    z[i + (j - 1) * nx],
 				    z[i - 1 + j * nx],
 				    z[i + j * nx],
 				    px, py, pz, &npt);
 		if (npt > 2)
-		    GPolygon(npt, px, py, USER, col[(k-1)%ncol],
+		    GPolygon(npt, px, py, USER, col[(k-1) % ncol],
 			     R_TRANWHITE, dd);
+	    }
+	}
+    }
     GMode(0, dd);

The R Project SVN R

(root)/trunk/src/main/plot3d.c @ 63243 – Rev 56695 → 57771