Check for isfinite in get_path_extents.

efiring · efiring · commit 6f0fa53b2d1b · 2008-08-06T18:06:21.000Z
This is part of an experiment to make nan and inf be ignored
equally throughout plotting; it may turn out to be unnecessary,
in which case it should be reverted for the sake of efficiency.

svn path=/trunk/matplotlib/; revision=5981
diff --git a/examples/pylab_examples/nan_test.py b/examples/pylab_examples/nan_test.py
@@ -3,10 +3,12 @@
 Example: simple line plots with NaNs inserted.
 """
 from pylab import *
+import numpy as np
 
 t = arange(0.0, 1.0+0.01, 0.01)
 s = cos(2*2*pi*t)
 t[41:60] = NaN
+#t[50:60] = np.inf
 
 subplot(2,1,1)
 plot(t, s, '-', lw=2)
@@ -20,6 +22,7 @@
 t[0] = NaN
 t[-1] = NaN
 plot(t, s, '-', lw=2)
+title('Also with NaN in first and last point')
 
 xlabel('time (s)')
 ylabel('more nans')
diff --git a/src/_path.cpp b/src/_path.cpp
@@ -280,6 +280,8 @@ void get_path_extents(PathIterator& path, const agg::trans_affine& trans,
     {
         if ((code & agg::path_cmd_end_poly) == agg::path_cmd_end_poly)
             continue;
+        if (MPL_notisfinite64(x) || MPL_notisfinite64(y))
+            continue;
         if (x < *x0) *x0 = x;
         if (y < *y0) *y0 = y;
         if (x > *x1) *x1 = x;

Original file line number	Diff line number	Diff line change
`@@ -280,6 +280,8 @@ void get_path_extents(PathIterator& path, const agg::trans_affine& trans,`
`280`	`280`	`{`
`281`	`281`	`if ((code & agg::path_cmd_end_poly) == agg::path_cmd_end_poly)`
`282`	`282`	`continue;`
	`283`	`+ if (MPL_notisfinite64(x) \|\| MPL_notisfinite64(y))`
	`284`	`+ continue;`
`283`	`285`	`if (x < x0) x0 = x;`
`284`	`286`	`if (y < y0) y0 = y;`
`285`	`287`	`if (x > x1) x1 = x;`