Prefer geomspace() to logspace().

anntzer · anntzer · commit 784c759328ec · 2020-02-12T23:59:43.000+01:00
geomspace() ("geometric progression from `a` to `b`") is, I think,
easier to read than logspace ("geometric progression from `10**a` to
`10**b`") -- see eg changes in ticker.py, colorbar.py.
diff --git a/examples/text_labels_and_annotations/custom_legends.py b/examples/text_labels_and_annotations/custom_legends.py
@@ -27,8 +27,7 @@
 np.random.seed(19680801)
 
 N = 10
-data = np.transpose([np.logspace(0, 1, 100) + np.random.randn(100) + ii
-                     for ii in range(N)])
+data = (np.geomspace(1, 10, 100) + np.random.randn(N, 100)).T
 cmap = plt.cm.coolwarm
 rcParams['axes.prop_cycle'] = cycler(color=cmap(np.linspace(0, 1, N)))
 
diff --git a/lib/matplotlib/ticker.py b/lib/matplotlib/ticker.py
@@ -968,7 +968,7 @@ def set_locs(self, locs=None):
             # Add labels between bases at log-spaced coefficients;
             # include base powers in case the locations include
             # "major" and "minor" points, as in colorbar.
-            c = np.logspace(0, 1, int(b)//2 + 1, base=b)
+            c = np.geomspace(1, b, int(b)//2 + 1)
             self._sublabels = set(np.round(c))
             # For base 10, this yields (1, 2, 3, 4, 6, 10).
         else:
diff --git a/lib/mpl_toolkits/axes_grid1/colorbar.py b/lib/mpl_toolkits/axes_grid1/colorbar.py
@@ -634,7 +634,7 @@ def _uniform_y(self, N):
         """
         vmin, vmax = self._get_colorbar_limits()
         if isinstance(self.norm, colors.LogNorm):
-            y = np.logspace(np.log10(vmin), np.log10(vmax), N)
+            y = np.geomspace(vmin, vmax, N)
         else:
             y = np.linspace(vmin, vmax, N)
         return y
diff --git a/tutorials/introductory/usage.py b/tutorials/introductory/usage.py
@@ -670,7 +670,7 @@ def my_plotter(ax, data1, data2, param_dict):
 #   # Setup, and create the data to plot
 #   y = np.random.rand(100000)
 #   y[50000:] *= 2
-#   y[np.logspace(1, np.log10(50000), 400).astype(int)] = -1
+#   y[np.geomspace(10, 50000, 400).astype(int)] = -1
 #   mpl.rcParams['path.simplify'] = True
 #
 #   mpl.rcParams['path.simplify_threshold'] = 0.0
@@ -749,7 +749,7 @@ def my_plotter(ax, data1, data2, param_dict):
 #   # Setup, and create the data to plot
 #   y = np.random.rand(100000)
 #   y[50000:] *= 2
-#   y[np.logspace(1, np.log10(50000), 400).astype(int)] = -1
+#   y[np.geomspace(10, 50000, 400).astype(int)] = -1
 #   mpl.rcParams['path.simplify'] = True
 #
 #   mpl.rcParams['agg.path.chunksize'] = 0
