Apr 16, 2012 · Mar 22, 2012
diff --git a/lib/matplotlib/streamplot.py b/lib/matplotlib/streamplot.py
 import numpy as np
 import matplotlib
 import matplotlib.patches as patches
 import matplotlib.cbook as cbook


 __all__ = ['streamplot']
        return -dxi, -dyi

    def integrate(x0, y0):
        """Return x, y coordinates of trajectory based on starting point.
        """Return x, y grid-coordinates of trajectory based on starting point.

        Integrate both forward and backward in time from starting point in
        grid coordinates.

        Integrate both forward and backward in time from starting point.
        Integration is terminated when a trajectory reaches a domain boundary
        or when it crosses into an already occupied cell in the StreamMask. The
        resulting trajectory is None if it is shorter than `minlength`.

 def _euler_step(xf_traj, yf_traj, dmap, f):
    """Simple Euler integration step."""
 nx, ny = dmap.grid.shape
 ny, nx = dmap.grid.shape
    xi = xf_traj[-1]
    yi = yf_traj[-1]
    cx, cy = f(xi, yi) # ds.cx is in data coordinates, ds in axis coord.
    cx, cy = f(xi, yi)
    if cx > 0:
        dsx = (nx - 1 - xi) / cx
    else:
Original file line number	Diff line number	Diff line change
Expand Up		@@ -5,7 +5,6 @@
		import numpy as np
		import matplotlib
		import matplotlib.patches as patches
		import matplotlib.cbook as cbook


		__all__ = ['streamplot']
Expand DownExpand Up		@@ -332,9 +331,11 @@ def backward_time(xi, yi):
		return -dxi, -dyi

		def integrate(x0, y0):
		"""Return x, y coordinates of trajectory based on starting point.
		"""Return x, y grid-coordinates of trajectory based on starting point.

		Integrate both forward and backward in time from starting point in
		grid coordinates.

		Integrate both forward and backward in time from starting point.
		Integration is terminated when a trajectory reaches a domain boundary
		or when it crosses into an already occupied cell in the StreamMask. The
		resulting trajectory is None if it is shorter than `minlength`.
Expand DownExpand Up		@@ -446,10 +447,10 @@ def _integrate_rk12(x0, y0, dmap, f):

		def _euler_step(xf_traj, yf_traj, dmap, f):
		"""Simple Euler integration step."""
		nx, ny = dmap.grid.shape
		ny, nx = dmap.grid.shape
		xi = xf_traj[-1]
		yi = yf_traj[-1]
		cx, cy = f(xi, yi) # ds.cx is in data coordinates, ds in axis coord.
		cx, cy = f(xi, yi)
		if cx > 0:
		dsx = (nx - 1 - xi) / cx
		else:
Expand Down