matplotlib
diff --git a/‎lib/matplotlib/axes/_base.py
Lines changed: 33 additions & 1 deletion b/‎lib/matplotlib/axes/_base.py
Lines changed: 33 additions & 1 deletion
diff --git a/‎lib/matplotlib/tests/test_cycles.py
Lines changed: 29 additions & 0 deletions b/‎lib/matplotlib/tests/test_cycles.py
Lines changed: 29 additions & 0 deletions
@@ -208,6 +208,38 @@ def _process_plot_format(fmt, *, ambiguous_fmt_datakey=False):
     return linestyle, marker, color
 
 
+def to_shareable_cycler(cycler):
+    """
+    Return a copy of *cycler*, that can be shared across multiple Axes.
+
+    Plotting on any Axes will advance the joint property iterators of all Axes.
+    """
+    from cycler import Cycler
+    new = Cycler(cycler)
+    # This doesn't carry forward any previous shares, but that can easily be
+    # changed by checking first whether the dict exists.
+    new._cycling_iterators_cache = {}
+    return new
+
+
+def _get_cycling_iterator(cycler, name):
+    """
+    Return an `itertools.cycle()` over this *cycler*, intended for plotting
+    command *name*.
+
+    If this cycler is shareable (per `to_shareable_cycler`), then multiple
+    invocations with the same *cycler* and *name* will return the same
+    `itertools.cycle` instance (*name* only serves as a key for this purpose).
+    If this cycler is not shareable, then a new `itertools.cycle` instance is
+    returned every time.
+    """
+    if not hasattr(cycler, "_cycling_iterators_cache"):
+        return itertools.cycle(cycler)
+    if name not in cycler._cycling_iterators_cache:
+        cycler._cycling_iterators_cache[name] = itertools.cycle(cycler)
+    return cycler._cycling_iterators_cache[name]
+
+
 class _process_plot_var_args:
     """
     Process variable length arguments to `~.Axes.plot`, to support ::
@@ -235,7 +267,7 @@ def __setstate__(self, state):
     def set_prop_cycle(self, cycler):
         if cycler is None:
             cycler = mpl.rcParams['axes.prop_cycle']
-        self.prop_cycler = itertools.cycle(cycler)
+        self.prop_cycler = _get_cycling_iterator(cycler, self.command)
         self._prop_keys = cycler.keys  # This should make a copy
 
     def __call__(self, *args, data=None, **kwargs):
 
@@ -158,3 +158,32 @@ def test_invalid_input_forms():
         ax.set_prop_cycle(cycler(foobar=[1, 2]))
     with pytest.raises(ValueError):
         ax.set_prop_cycle(cycler(color='rgb', c='cmy'))
+
+
+@mpl.style.context("default")
+def test_cycler_sharing():
+    from matplotlib.axes._base import to_shareable_cycler
+
+    shared = to_shareable_cycler(plt.rcParams["axes.prop_cycle"])
+    fig, axs = plt.subplots(2)
+    axs[0].set_prop_cycle(shared)
+    axs[1].set_prop_cycle(shared)
+
+    axs[0].plot([0, 1])
+    axs[1].plot([0, 1])
+    # Check that fill uses a separate iterator.
+    axs[1].fill([2, 3, 4], [0, 1, 0])
+    axs[0].fill([2, 3, 4], [0, 1, 0])
+    # Check that the second axes goes back to the first color, even though it
+    # didn't use it on the first cycle.
+    for i in range(10):
+        axs[1].plot([5+i, 5+i], [0, 1])
+
+    assert mpl.colors.same_color([l.get_color() for l in axs[0].lines],
+                                 ["C0"])
+    assert mpl.colors.same_color([l.get_color() for l in axs[1].lines],
+                                 [f"C{i}" for i in "12345678901"])
+    assert mpl.colors.same_color([p.get_facecolor() for p in axs[0].patches],
+                                 ["C1"])
+    assert mpl.colors.same_color([p.get_facecolor() for p in axs[1].patches],
+                                 ["C0"])