From 6e50a75e0929c975aefd5e74e08fbce65df8ebf9 Mon Sep 17 00:00:00 2001
From: Lee Johnston <lee.johnston@ni.com>
Date: Tue, 21 Jan 2020 13:21:04 -0600
Subject: [PATCH 01/14] update covert to handle sequences and make the
 interface optional

---
 unyt/__init__.py                 |  2 +-
 unyt/mpl_interface.py            | 39 ++++++++++++++++++++++++----
 unyt/tests/test_mpl_interface.py | 44 +++++++++++++++++++++++++++++++-
 3 files changed, 78 insertions(+), 7 deletions(-)

diff --git a/unyt/__init__.py b/unyt/__init__.py
index f8989964..e3720101 100644
--- a/unyt/__init__.py
+++ b/unyt/__init__.py
@@ -96,7 +96,7 @@
 from unyt.unit_systems import UnitSystem  # NOQA: F401
 from unyt.testing import assert_allclose_units  # NOQA: F401
 from unyt.dimensions import accepts, returns  # NOQA: F401
-from unyt import mpl_interface  # NOQA: F401
+from unyt.mpl_interface import MplUnitsCM  # NOQA: F401
 
 
 # function to only import quantities into this namespace
diff --git a/unyt/mpl_interface.py b/unyt/mpl_interface.py
index 8d055910..1c2ad466 100644
--- a/unyt/mpl_interface.py
+++ b/unyt/mpl_interface.py
@@ -18,6 +18,7 @@
         ConversionInterface,
         AxisInfo,
         registry,
+        ConversionError,
     )
 except ImportError:
     pass
@@ -81,7 +82,7 @@ def convert(value, unit, axis):
             Parameters
             ----------
 
-            value : unyt_array
+            value : unyt_array, unyt_quantity, or sequence there of
             unit : Unit, string or tuple
                 This parameter comes from unyt_arrayConverter.default_units() or from
                 user code such as Axes.plot(), Axis.set_units(), etc. In user code, it
@@ -98,11 +99,39 @@ def convert(value, unit, axis):
             Raises
             ------
 
-            ConversionError if unit does not have the same dimensions as value
+            ConversionError if unit does not have the same dimensions as value or
+            if we don't know how to convert value.
             """
             if isinstance(unit, str) or isinstance(unit, Unit):
                 unit = (unit,)
-            return value.to(*unit)
+            if isinstance(value, (unyt_array, unyt_quantity)):
+                converted_value = value.to(*unit)
+            elif isinstance(value, (list, tuple)):
+                value_type = type(value)
+                converted_value = []
+                for obj in value:
+                    converted_value.append(obj.to(*unit))
+                converted_value = value_type(converted_value)
+            else:
+                raise ConversionError(f"unable to convert {value}")
+            return converted_value
+
+    class MplUnitsCM:
+        """Context manager for experimenting with Unyt in Matplotlib"""
+
+        def __enter__(self):
+            registry[unyt_array] = unyt_arrayConverter()
+            registry[unyt_quantity] = unyt_arrayConverter()
+
+        def __exit__(self, exc_type, exc_val, exc_tb):
+            registry.pop(unyt_array)
+            registry.pop(unyt_quantity)
+
+        def enable(self):
+            self.__enter__()
+
+        def disable(self):
+            self.__exit__()
 
-    registry[unyt_array] = unyt_arrayConverter()
-    registry[unyt_quantity] = unyt_arrayConverter()
+        def __repr__(self):
+            return "<MplUnitsCM>"
diff --git a/unyt/tests/test_mpl_interface.py b/unyt/tests/test_mpl_interface.py
index 94cb6cb2..c975fe91 100644
--- a/unyt/tests/test_mpl_interface.py
+++ b/unyt/tests/test_mpl_interface.py
@@ -2,8 +2,9 @@
 import numpy as np
 import pytest
 from unyt._on_demand_imports import _matplotlib, NotAModule
-from unyt import s, K
+from unyt import s, K, unyt_array, unyt_quantity
 from unyt.exceptions import UnitConversionError
+from unyt.mpl_interface import MplUnitsCM
 
 
 check_matplotlib = pytest.mark.skipif(
@@ -13,9 +14,12 @@
 
 @pytest.fixture
 def ax(scope="module"):
+    mpl_units = MplUnitsCM()
+    mpl_units.enable()
     fig, ax = _matplotlib.pyplot.subplots()
     yield ax
     _matplotlib.pyplot.close()
+    mpl_units.disable()
 
 
 @check_matplotlib
@@ -101,3 +105,41 @@ def test_list_label(ax):
     assert ax.xaxis.get_label().get_text() == expected_xlabel
     expected_ylabel = ""
     assert ax.yaxis.get_label().get_text() == expected_ylabel
+
+
+@check_matplotlib
+def test_errorbar(ax):
+    x = unyt_array([8, 9, 10], "cm")
+    y = unyt_array([8, 9, 10], "kg")
+    y_scatter = [
+        unyt_array([0.1, 0.2, 0.3], "kg"),
+        unyt_array([0.1, 0.2, 0.3], "kg"),
+    ]
+    x_lims = (unyt_quantity(5, "cm"), unyt_quantity(12, "cm"))
+    y_lims = (unyt_quantity(5, "kg"), unyt_quantity(12, "kg"))
+
+    ax.errorbar(x, y, yerr=y_scatter)
+    x_lims = (unyt_quantity(5, "cm"), unyt_quantity(12, "cm"))
+    y_lims = (unyt_quantity(5, "kg"), unyt_quantity(12, "kg"))
+    ax.set_xlim(*x_lims)
+    ax.set_ylim(*y_lims)
+
+
+@check_matplotlib
+def test_hist2d(ax):
+    x = np.random.normal(size=50000) * s
+    y = 3 * x + np.random.normal(size=50000) * s
+    ax.hist2d(x, y, bins=(50, 50))
+
+
+@check_matplotlib
+def test_imshow(ax):
+    data = np.reshape(np.random.normal(size=10000), (100, 100))
+    ax.imshow(data, vmin=data.min(), vmax=data.max())
+
+
+@check_matplotlib
+def test_hist(ax):
+    data = np.random.normal(size=10000) * s
+    bin_edges = np.linspace(data.min(), data.max(), 50)
+    ax.hist(data, bins=bin_edges)

From 12472f8ccdfd4a251a6bdf451674a8f292bbaf18 Mon Sep 17 00:00:00 2001
From: Lee Johnston <lee.johnston@ni.com>
Date: Tue, 21 Jan 2020 13:52:26 -0600
Subject: [PATCH 02/14] update documentation and add context manager test

---
 docs/usage.rst                   | 17 ++++++++++-------
 unyt/mpl_interface.py            |  2 +-
 unyt/tests/test_mpl_interface.py | 14 +++++++++++++-
 3 files changed, 24 insertions(+), 9 deletions(-)

diff --git a/docs/usage.rst b/docs/usage.rst
index b9ae7b6c..9efe05bf 100644
--- a/docs/usage.rst
+++ b/docs/usage.rst
@@ -1184,25 +1184,28 @@ calculation involves many operations on arrays with only a few elements.
 
 Plotting with Matplotlib
 ++++++++++++++++++++++++
+.. note::
+  - This is an experimental feature. Please report issues.
+  - The context manager ``MplUnitsCM`` will temporarily enable this feature
 
 Matplotlib is Unyt aware. With no additional effort, Matplotlib will label the x and y
 axes with the units.
 
   >>> import matplotlib.pyplot as plt
-  >>> from unyt import s, K
+  >>> from unyt import s, K, MplUnitsCM
   >>> x = [0.0, 0.01, 0.02]*s
   >>> y = [298.15, 308.15, 318.15]*K
-  >>> plt.plot(x, y)
-  [<matplotlib.lines.Line2D object at ...>]
-  >>> plt.show()
+  >>> with MplUnitsCM():
+  ... plt.plot(x, y)
+  ... plt.show()
 
 .. image:: _static/mpl_fig1.png
 
 You can change the plotted units without affecting the original data.
 
-  >>> plt.plot(x, y, xunits="ms", yunits=("J", "thermal"))
-  [<matplotlib.lines.Line2D object at ...>]
-  >>> plt.show()
+  >>> with MplUnitsCM():
+  ... plt.plot(x, y, xunits="ms", yunits=("J", "thermal"))
+  ... plt.show()
 
 .. image:: _static/mpl_fig2.png
 
diff --git a/unyt/mpl_interface.py b/unyt/mpl_interface.py
index 1c2ad466..7e50c65a 100644
--- a/unyt/mpl_interface.py
+++ b/unyt/mpl_interface.py
@@ -131,7 +131,7 @@ def enable(self):
             self.__enter__()
 
         def disable(self):
-            self.__exit__()
+            self.__exit__(None, None, None)
 
         def __repr__(self):
             return "<MplUnitsCM>"
diff --git a/unyt/tests/test_mpl_interface.py b/unyt/tests/test_mpl_interface.py
index c975fe91..fa40235a 100644
--- a/unyt/tests/test_mpl_interface.py
+++ b/unyt/tests/test_mpl_interface.py
@@ -4,7 +4,7 @@
 from unyt._on_demand_imports import _matplotlib, NotAModule
 from unyt import s, K, unyt_array, unyt_quantity
 from unyt.exceptions import UnitConversionError
-from unyt.mpl_interface import MplUnitsCM
+from unyt.mpl_interface import MplUnitsCM, unyt_arrayConverter
 
 
 check_matplotlib = pytest.mark.skipif(
@@ -143,3 +143,15 @@ def test_hist(ax):
     data = np.random.normal(size=10000) * s
     bin_edges = np.linspace(data.min(), data.max(), 50)
     ax.hist(data, bins=bin_edges)
+
+
+@check_matplotlib
+def test_MplUnitsCM():
+    mplunits = MplUnitsCM()
+    with pytest.raises(KeyError):
+        _matplotlib.units.registry[unyt_array]
+    mplunits.enable()
+    assert isinstance(_matplotlib.units.registry[unyt_array], unyt_arrayConverter)
+    mplunits.disable()
+    assert unyt_array not in _matplotlib.units.registry.keys()
+    assert unyt_quantity not in _matplotlib.units.registry.keys()

From 031fa5c30fe82e650b827888286dc9a53129a0ff Mon Sep 17 00:00:00 2001
From: Lee Johnston <lee.johnston@ni.com>
Date: Tue, 21 Jan 2020 14:04:57 -0600
Subject: [PATCH 03/14] fix usage doctest

---
 docs/usage.rst | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/docs/usage.rst b/docs/usage.rst
index 9efe05bf..995c50ce 100644
--- a/docs/usage.rst
+++ b/docs/usage.rst
@@ -1196,16 +1196,16 @@ axes with the units.
   >>> x = [0.0, 0.01, 0.02]*s
   >>> y = [298.15, 308.15, 318.15]*K
   >>> with MplUnitsCM():
-  ... plt.plot(x, y)
-  ... plt.show()
+      plt.plot(x, y)
+      plt.show()
 
 .. image:: _static/mpl_fig1.png
 
 You can change the plotted units without affecting the original data.
 
   >>> with MplUnitsCM():
-  ... plt.plot(x, y, xunits="ms", yunits=("J", "thermal"))
-  ... plt.show()
+      plt.plot(x, y, xunits="ms", yunits=("J", "thermal"))
+      plt.show()
 
 .. image:: _static/mpl_fig2.png
 

From 972cb4bba6e09eea2527846a4ca6ebe30297c3ff Mon Sep 17 00:00:00 2001
From: Lee Johnston <lee.johnston@ni.com>
Date: Tue, 21 Jan 2020 14:33:00 -0600
Subject: [PATCH 04/14] fix doctest in usage

---
 docs/usage.rst | 10 ++++++----
 1 file changed, 6 insertions(+), 4 deletions(-)

diff --git a/docs/usage.rst b/docs/usage.rst
index 995c50ce..f328ff7c 100644
--- a/docs/usage.rst
+++ b/docs/usage.rst
@@ -1196,16 +1196,18 @@ axes with the units.
   >>> x = [0.0, 0.01, 0.02]*s
   >>> y = [298.15, 308.15, 318.15]*K
   >>> with MplUnitsCM():
-      plt.plot(x, y)
-      plt.show()
+  ...   plt.plot(x, y)
+  ...   plt.show()
+  [<matplotlib.lines.Line2D object at ...>]
 
 .. image:: _static/mpl_fig1.png
 
 You can change the plotted units without affecting the original data.
 
   >>> with MplUnitsCM():
-      plt.plot(x, y, xunits="ms", yunits=("J", "thermal"))
-      plt.show()
+  ...   plt.plot(x, y, xunits="ms", yunits=("J", "thermal"))
+  ...   plt.show()
+  [<matplotlib.lines.Line2D object at ...>]
 
 .. image:: _static/mpl_fig2.png
 

From c56bb7b42cabd7940f321275278a550392ada292 Mon Sep 17 00:00:00 2001
From: Lee Johnston <lee.johnston@ni.com>
Date: Tue, 21 Jan 2020 17:15:26 -0600
Subject: [PATCH 05/14] make context manager callable, update docs and tests

---
 docs/usage.rst                   |  8 ++---
 unyt/__init__.py                 |  4 ++-
 unyt/mpl_interface.py            | 59 +++++++++++++++++++++++++++++---
 unyt/tests/test_mpl_interface.py | 42 ++++++++++++++++++-----
 4 files changed, 95 insertions(+), 18 deletions(-)

diff --git a/docs/usage.rst b/docs/usage.rst
index f328ff7c..faa877de 100644
--- a/docs/usage.rst
+++ b/docs/usage.rst
@@ -1186,16 +1186,16 @@ Plotting with Matplotlib
 ++++++++++++++++++++++++
 .. note::
   - This is an experimental feature. Please report issues.
-  - The context manager ``MplUnitsCM`` will temporarily enable this feature
+  - The context manager ``matplotlib_support`` will temporarily enable this feature
 
 Matplotlib is Unyt aware. With no additional effort, Matplotlib will label the x and y
 axes with the units.
 
   >>> import matplotlib.pyplot as plt
-  >>> from unyt import s, K, MplUnitsCM
+  >>> from unyt import s, K, matplotlib_support
   >>> x = [0.0, 0.01, 0.02]*s
   >>> y = [298.15, 308.15, 318.15]*K
-  >>> with MplUnitsCM():
+  >>> with matplotlib_support:
   ...   plt.plot(x, y)
   ...   plt.show()
   [<matplotlib.lines.Line2D object at ...>]
@@ -1204,7 +1204,7 @@ axes with the units.
 
 You can change the plotted units without affecting the original data.
 
-  >>> with MplUnitsCM():
+  >>> with matplotlib_support:
   ...   plt.plot(x, y, xunits="ms", yunits=("J", "thermal"))
   ...   plt.show()
   [<matplotlib.lines.Line2D object at ...>]
diff --git a/unyt/__init__.py b/unyt/__init__.py
index e3720101..dcc63318 100644
--- a/unyt/__init__.py
+++ b/unyt/__init__.py
@@ -96,7 +96,9 @@
 from unyt.unit_systems import UnitSystem  # NOQA: F401
 from unyt.testing import assert_allclose_units  # NOQA: F401
 from unyt.dimensions import accepts, returns  # NOQA: F401
-from unyt.mpl_interface import MplUnitsCM  # NOQA: F401
+from unyt.mpl_interface import matplotlib_support  # NOQA: F401
+
+matplotlib_support = matplotlib_support()
 
 
 # function to only import quantities into this namespace
diff --git a/unyt/mpl_interface.py b/unyt/mpl_interface.py
index 7e50c65a..2629f769 100644
--- a/unyt/mpl_interface.py
+++ b/unyt/mpl_interface.py
@@ -28,6 +28,8 @@
     class unyt_arrayConverter(ConversionInterface):
         """Matplotlib interface for unyt_array"""
 
+        _labelstyle = "()"
+
         @staticmethod
         def axisinfo(unit, axis):
             """Set the axis label based on unit
@@ -55,7 +57,18 @@ def axisinfo(unit, axis):
                 label = ""
             else:
                 unit_str = unit_obj.latex_representation()
-                label = "$\\left(" + unit_str + "\\right)$"
+                if unyt_arrayConverter._labelstyle == "[]":
+                    label = "$\\left[" + unit_str + "\\right]$"
+                elif unyt_arrayConverter._labelstyle == "/":
+                    axsym = axis.axis_name
+                    if "/" in unit_str:
+                        label = (
+                            "$q_{" + axsym + "}\\;/\\;\\left(" + unit_str + "\\right)$"
+                        )
+                    else:
+                        label = "$q_{" + axsym + "}\\;/\\;" + unit_str + "$"
+                else:
+                    label = "$\\left(" + unit_str + "\\right)$"
             return AxisInfo(label=label)
 
         @staticmethod
@@ -113,11 +126,49 @@ def convert(value, unit, axis):
                     converted_value.append(obj.to(*unit))
                 converted_value = value_type(converted_value)
             else:
-                raise ConversionError(f"unable to convert {value}")
+                raise ConversionError("unable to convert {%s}".format(value))
             return converted_value
 
-    class MplUnitsCM:
-        """Context manager for experimenting with Unyt in Matplotlib"""
+    class matplotlib_support:
+        """Context manager for experimenting with Unyt in Matplotlib
+
+        Parameters
+        ----------
+
+        label_style : string, one from the set {'()', '[]', '/'}
+            The axis label style.
+            '()' -> '(unit)'
+            '[]' -> '[unit]'
+            '/' -> 'q / unit'
+            SI standard where label is a mathematical expression.
+            'q' is a generic quantity symbol and for a value on the axis, x,
+            the equation is q = x * unit.
+        """
+
+        def __init__(self, label_style="()"):
+            self._labelstyle = label_style
+            unyt_arrayConverter._labelstyle = label_style
+
+        def __call__(self):
+            self.__enter__()
+
+        @property
+        def label_style(self):
+            """label_style : string, one from the set {'()', '[]', '/'}
+                The axis label style.
+                '()' -> '(unit)'
+                '[]' -> '[unit]'
+                '/' -> 'q / unit'
+                SI standard where label is a mathematical expression.
+                'q' is a generic quantity symbol and for a value on the axis, x,
+                the equation is q = x * unit.
+            """
+            return self._labelstyle
+
+        @label_style.setter
+        def label_style(self, label_style="()"):
+            self._labelstyle = label_style
+            unyt_arrayConverter._labelstyle = label_style
 
         def __enter__(self):
             registry[unyt_array] = unyt_arrayConverter()
diff --git a/unyt/tests/test_mpl_interface.py b/unyt/tests/test_mpl_interface.py
index fa40235a..39f182d7 100644
--- a/unyt/tests/test_mpl_interface.py
+++ b/unyt/tests/test_mpl_interface.py
@@ -2,9 +2,9 @@
 import numpy as np
 import pytest
 from unyt._on_demand_imports import _matplotlib, NotAModule
-from unyt import s, K, unyt_array, unyt_quantity
+from unyt import s, K, unyt_array, unyt_quantity, matplotlib_support
 from unyt.exceptions import UnitConversionError
-from unyt.mpl_interface import MplUnitsCM, unyt_arrayConverter
+from unyt.mpl_interface import unyt_arrayConverter
 
 
 check_matplotlib = pytest.mark.skipif(
@@ -14,12 +14,11 @@
 
 @pytest.fixture
 def ax(scope="module"):
-    mpl_units = MplUnitsCM()
-    mpl_units.enable()
+    matplotlib_support.enable()
     fig, ax = _matplotlib.pyplot.subplots()
     yield ax
     _matplotlib.pyplot.close()
-    mpl_units.disable()
+    matplotlib_support.disable()
 
 
 @check_matplotlib
@@ -146,12 +145,37 @@ def test_hist(ax):
 
 
 @check_matplotlib
-def test_MplUnitsCM():
-    mplunits = MplUnitsCM()
+def test_matplotlib_support():
     with pytest.raises(KeyError):
         _matplotlib.units.registry[unyt_array]
-    mplunits.enable()
+    matplotlib_support.enable()
     assert isinstance(_matplotlib.units.registry[unyt_array], unyt_arrayConverter)
-    mplunits.disable()
+    matplotlib_support.disable()
     assert unyt_array not in _matplotlib.units.registry.keys()
     assert unyt_quantity not in _matplotlib.units.registry.keys()
+    # test as a callable
+    matplotlib_support()
+    assert isinstance(_matplotlib.units.registry[unyt_array], unyt_arrayConverter)
+
+
+@check_matplotlib
+def test_labelstyle():
+    x = [0, 1, 2] * s
+    y = [3, 4, 5] * K
+    matplotlib_support.label_style = "[]"
+    matplotlib_support.enable()
+    assert unyt_arrayConverter._labelstyle == "[]"
+    fig, ax = _matplotlib.pyplot.subplots()
+    ax.plot(x, y)
+    expected_xlabel = "$\\left[\\rm{s}\\right]$"
+    assert ax.xaxis.get_label().get_text() == expected_xlabel
+    expected_ylabel = "$\\left[\\rm{K}\\right]$"
+    assert ax.yaxis.get_label().get_text() == expected_ylabel
+    matplotlib_support.label_style = "/"
+    ax.clear()
+    ax.plot(x, y)
+    expected_xlabel = "$q_{x}\\;/\\;\\rm{s}$"
+    assert ax.xaxis.get_label().get_text() == expected_xlabel
+    expected_ylabel = "$q_{y}\\;/\\;\\rm{K}$"
+    assert ax.yaxis.get_label().get_text() == expected_ylabel
+    _matplotlib.pyplot.close()

From 20287125503d9a5c48ee6dcf9d5f491ebd8585c9 Mon Sep 17 00:00:00 2001
From: Lee Johnston <lee.johnston@ni.com>
Date: Tue, 21 Jan 2020 17:58:58 -0600
Subject: [PATCH 06/14] fix import issue when matplotlib isn't available

---
 unyt/__init__.py      | 8 ++++++--
 unyt/mpl_interface.py | 3 ---
 2 files changed, 6 insertions(+), 5 deletions(-)

diff --git a/unyt/__init__.py b/unyt/__init__.py
index dcc63318..a83476ab 100644
--- a/unyt/__init__.py
+++ b/unyt/__init__.py
@@ -96,9 +96,13 @@
 from unyt.unit_systems import UnitSystem  # NOQA: F401
 from unyt.testing import assert_allclose_units  # NOQA: F401
 from unyt.dimensions import accepts, returns  # NOQA: F401
-from unyt.mpl_interface import matplotlib_support  # NOQA: F401
 
-matplotlib_support = matplotlib_support()
+try:
+    from unyt.mpl_interface import matplotlib_support  # NOQA: F401
+except ImportError:
+    pass
+else:
+    matplotlib_support = matplotlib_support()
 
 
 # function to only import quantities into this namespace
diff --git a/unyt/mpl_interface.py b/unyt/mpl_interface.py
index 2629f769..3c8b37d9 100644
--- a/unyt/mpl_interface.py
+++ b/unyt/mpl_interface.py
@@ -183,6 +183,3 @@ def enable(self):
 
         def disable(self):
             self.__exit__(None, None, None)
-
-        def __repr__(self):
-            return "<MplUnitsCM>"

From 15609539bf135e6b7da96559f4507cc00d481c53 Mon Sep 17 00:00:00 2001
From: Lee Johnston <lee.johnston@ni.com>
Date: Tue, 21 Jan 2020 18:42:53 -0600
Subject: [PATCH 07/14] fix import issue and add some more documentation

---
 docs/usage.rst                   | 30 ++++++++++++++++++++++++++++++
 unyt/tests/test_mpl_interface.py |  8 ++++++--
 2 files changed, 36 insertions(+), 2 deletions(-)

diff --git a/docs/usage.rst b/docs/usage.rst
index faa877de..f57e5706 100644
--- a/docs/usage.rst
+++ b/docs/usage.rst
@@ -1211,6 +1211,36 @@ You can change the plotted units without affecting the original data.
 
 .. image:: _static/mpl_fig2.png
 
+There are three ways to use the context manager.
+
+1. As a conventional context manager in a ``with`` statement as shown above
+
+2. As a feature toggle in an interactive session
+
+  >>> import matplotlib.pyplot as plt
+  >>> from unyt import s, K, matplotlib_support
+  >>> matplotlib_support.enable()
+  >>> plt.plot([0, 1, 2]*s, [3, 4, 5]*K)
+  [<matplotlib.lines.Line2D object at ...>]
+  >>> plt.show()
+  >>> matplotlib_support.disable()
+
+3. As an enable for a complete session
+
+  >>> import unyt
+  >>> unyt.matplotlib_support()
+  >>> import matplotlib.pyplot as plt
+
+It is possible to set the label style.
+
+  >>> import matplotlib.pyplot as plt
+  >>> from unyt import s, K, matplotlib_support
+  >>> matplotlib_support.label_style = "[]"
+  >>> with matplotlib_support:
+  ...   plt.plot([0, 1, 2]*s, [3, 4, 5]*K)
+  ...   plt.show()
+  [<matplotlib.lines.Line2D object at ...>]
+
 .. note::
   
   - This feature works in Matplotlib versions 2.2.4 and above
diff --git a/unyt/tests/test_mpl_interface.py b/unyt/tests/test_mpl_interface.py
index 39f182d7..7aee8a22 100644
--- a/unyt/tests/test_mpl_interface.py
+++ b/unyt/tests/test_mpl_interface.py
@@ -2,10 +2,14 @@
 import numpy as np
 import pytest
 from unyt._on_demand_imports import _matplotlib, NotAModule
-from unyt import s, K, unyt_array, unyt_quantity, matplotlib_support
+from unyt import s, K, unyt_array, unyt_quantity
 from unyt.exceptions import UnitConversionError
-from unyt.mpl_interface import unyt_arrayConverter
 
+try:
+    from unyt import matplotlib_support
+    from unyt.mpl_interface import unyt_arrayConverter
+except ImportError:
+    pass
 
 check_matplotlib = pytest.mark.skipif(
     isinstance(_matplotlib.pyplot, NotAModule), reason="matplotlib not installed"

From e24853cede8cf4d9a20e90d31d44778f9b6dfcf1 Mon Sep 17 00:00:00 2001
From: Lee Johnston <lee.johnston@ni.com>
Date: Tue, 21 Jan 2020 19:24:37 -0600
Subject: [PATCH 08/14] fix matplotlib 2.2.4 issues

---
 unyt/_on_demand_imports.py       | 11 +++++++++++
 unyt/mpl_interface.py            |  4 ++--
 unyt/tests/test_mpl_interface.py |  3 ++-
 3 files changed, 15 insertions(+), 3 deletions(-)

diff --git a/unyt/_on_demand_imports.py b/unyt/_on_demand_imports.py
index 0b09d8b3..fe6dd854 100644
--- a/unyt/_on_demand_imports.py
+++ b/unyt/_on_demand_imports.py
@@ -131,6 +131,7 @@ class matplotlib_imports(object):
     _name = "matplotlib"
     _pyplot = None
     _units = None
+    _use = None
 
     @property
     def __version__(self):
@@ -163,5 +164,15 @@ def units(self):
             self._units = units
         return self._units
 
+    @property
+    def use(self):
+        if self._use is None:
+            try:
+                from matplotlib import use
+            except ImportError:
+                use = NotAModule(self._name)
+            self._use = use
+        return self._use
+
 
 _matplotlib = matplotlib_imports()
diff --git a/unyt/mpl_interface.py b/unyt/mpl_interface.py
index 3c8b37d9..53b3d016 100644
--- a/unyt/mpl_interface.py
+++ b/unyt/mpl_interface.py
@@ -18,12 +18,12 @@
         ConversionInterface,
         AxisInfo,
         registry,
-        ConversionError,
     )
 except ImportError:
     pass
 else:
     from unyt import unyt_array, unyt_quantity, Unit
+    from unyt.exceptions import UnitConversionError
 
     class unyt_arrayConverter(ConversionInterface):
         """Matplotlib interface for unyt_array"""
@@ -126,7 +126,7 @@ def convert(value, unit, axis):
                     converted_value.append(obj.to(*unit))
                 converted_value = value_type(converted_value)
             else:
-                raise ConversionError("unable to convert {%s}".format(value))
+                raise UnitConversionError("unable to convert {%s}".format(value))
             return converted_value
 
     class matplotlib_support:
diff --git a/unyt/tests/test_mpl_interface.py b/unyt/tests/test_mpl_interface.py
index 7aee8a22..367ab67d 100644
--- a/unyt/tests/test_mpl_interface.py
+++ b/unyt/tests/test_mpl_interface.py
@@ -17,7 +17,8 @@
 
 
 @pytest.fixture
-def ax(scope="module"):
+def ax():
+    _matplotlib.use("agg")
     matplotlib_support.enable()
     fig, ax = _matplotlib.pyplot.subplots()
     yield ax

From f1b5f3c6f07ec4d595ecc0c546e493a860f8610c Mon Sep 17 00:00:00 2001
From: Lee Johnston <lee.johnston@ni.com>
Date: Tue, 21 Jan 2020 19:33:28 -0600
Subject: [PATCH 09/14] fix label_style test issue

---
 unyt/tests/test_mpl_interface.py | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/unyt/tests/test_mpl_interface.py b/unyt/tests/test_mpl_interface.py
index 367ab67d..6c55b6a2 100644
--- a/unyt/tests/test_mpl_interface.py
+++ b/unyt/tests/test_mpl_interface.py
@@ -30,6 +30,7 @@ def ax():
 def test_label(ax):
     x = [0, 1, 2] * s
     y = [3, 4, 5] * K
+    matplotlib_support.label_style = "()"
     ax.plot(x, y)
     expected_xlabel = "$\\left(\\rm{s}\\right)$"
     assert ax.xaxis.get_label().get_text() == expected_xlabel
@@ -93,6 +94,7 @@ def test_conversionerror(ax):
 def test_ndarray_label(ax):
     x = [0, 1, 2] * s
     y = np.arange(3, 6)
+    matplotlib_support.label_style = "()"
     ax.plot(x, y)
     expected_xlabel = "$\\left(\\rm{s}\\right)$"
     assert ax.xaxis.get_label().get_text() == expected_xlabel
@@ -104,6 +106,7 @@ def test_ndarray_label(ax):
 def test_list_label(ax):
     x = [0, 1, 2] * s
     y = [3, 4, 5]
+    matplotlib_support.label_style = "()"
     ax.plot(x, y)
     expected_xlabel = "$\\left(\\rm{s}\\right)$"
     assert ax.xaxis.get_label().get_text() == expected_xlabel

From 63b9902b4c571a4655f87a56f1438aa96f42d0ec Mon Sep 17 00:00:00 2001
From: Lee Johnston <lee.johnston@ni.com>
Date: Tue, 21 Jan 2020 20:24:27 -0600
Subject: [PATCH 10/14] increase code coverage

---
 unyt/mpl_interface.py            | 9 +++------
 unyt/tests/test_mpl_interface.py | 8 +++++++-
 2 files changed, 10 insertions(+), 7 deletions(-)

diff --git a/unyt/mpl_interface.py b/unyt/mpl_interface.py
index 53b3d016..d579cdaf 100644
--- a/unyt/mpl_interface.py
+++ b/unyt/mpl_interface.py
@@ -23,7 +23,6 @@
     pass
 else:
     from unyt import unyt_array, unyt_quantity, Unit
-    from unyt.exceptions import UnitConversionError
 
     class unyt_arrayConverter(ConversionInterface):
         """Matplotlib interface for unyt_array"""
@@ -112,21 +111,19 @@ def convert(value, unit, axis):
             Raises
             ------
 
-            ConversionError if unit does not have the same dimensions as value or
-            if we don't know how to convert value.
+            ConversionError if unit does not have the same dimensions as value
             """
+            converted_value = value
             if isinstance(unit, str) or isinstance(unit, Unit):
                 unit = (unit,)
             if isinstance(value, (unyt_array, unyt_quantity)):
                 converted_value = value.to(*unit)
-            elif isinstance(value, (list, tuple)):
+            else:
                 value_type = type(value)
                 converted_value = []
                 for obj in value:
                     converted_value.append(obj.to(*unit))
                 converted_value = value_type(converted_value)
-            else:
-                raise UnitConversionError("unable to convert {%s}".format(value))
             return converted_value
 
     class matplotlib_support:
diff --git a/unyt/tests/test_mpl_interface.py b/unyt/tests/test_mpl_interface.py
index 6c55b6a2..ace8b403 100644
--- a/unyt/tests/test_mpl_interface.py
+++ b/unyt/tests/test_mpl_interface.py
@@ -2,7 +2,7 @@
 import numpy as np
 import pytest
 from unyt._on_demand_imports import _matplotlib, NotAModule
-from unyt import s, K, unyt_array, unyt_quantity
+from unyt import m, s, K, unyt_array, unyt_quantity
 from unyt.exceptions import UnitConversionError
 
 try:
@@ -171,6 +171,7 @@ def test_labelstyle():
     x = [0, 1, 2] * s
     y = [3, 4, 5] * K
     matplotlib_support.label_style = "[]"
+    assert matplotlib_support.label_style == "[]"
     matplotlib_support.enable()
     assert unyt_arrayConverter._labelstyle == "[]"
     fig, ax = _matplotlib.pyplot.subplots()
@@ -186,4 +187,9 @@ def test_labelstyle():
     assert ax.xaxis.get_label().get_text() == expected_xlabel
     expected_ylabel = "$q_{y}\\;/\\;\\rm{K}$"
     assert ax.yaxis.get_label().get_text() == expected_ylabel
+    x = [0, 1, 2] * m / s
+    ax.clear()
+    ax.plot(x, y)
+    expected_xlabel = "$q_{x}\\;/\\;\\left(\\rm{m} / \\rm{s}\\right)$"
+    assert ax.xaxis.get_label().get_text() == expected_xlabel
     _matplotlib.pyplot.close()

From 3341aa34e14504f845af8f9ff1b60cf4c1f2854b Mon Sep 17 00:00:00 2001
From: Nathan Goldbaum <nathan12343@gmail.com>
Date: Tue, 21 Jan 2020 20:15:11 -0700
Subject: [PATCH 11/14] tweak docs and docstring formatting

---
 docs/_static/mpl_fig3.png        | Bin 0 -> 22200 bytes
 docs/usage.rst                   |  42 +++++++++++++++----------------
 unyt/mpl_interface.py            |  34 ++++++++++++-------------
 unyt/tests/test_mpl_interface.py |   2 ++
 4 files changed, 39 insertions(+), 39 deletions(-)
 create mode 100644 docs/_static/mpl_fig3.png

diff --git a/docs/_static/mpl_fig3.png b/docs/_static/mpl_fig3.png
new file mode 100644
index 0000000000000000000000000000000000000000..54ff844f7a41a2f821b4414f6fafaff2e224052a
GIT binary patch
literal 22200
zcmeFZcT`o)mo0h#Ns>e)M-KujNEXQ$P(Y%flCywh$skctqJW4>21NuBNs@C;l9Zft
zkensMTl@HR_pkfjzPHDCW4yl}V|@DIK4(|$s@gTznscrSR8f*4IZbyOg+h_Y$x5rD
zP}n{w6xKKa9{kJIU&GVz*C~gaa`y?~$AjQe0Q^q)L{`%Qg(5RV{=@nrm1++EDC~G!
z!%@xF)X~Mz-UMZB==j9S*3rt`h{f5&-of0~Mu3}-o1cs2v7_S?Q68TExXx{BZ^pwG
z{P8Xd#e$NPzH#3*ZfV3_=dIJ|$$H<DJE!U>IPTqB=i(S4?rKw1YrWfdE6BOjsjqj$
zVn9ptX{w1<RD8E!u6@RcSTa{`fp1EIb7$+RQ<pmM1#u}~V!g#O%^YJsdao<gf!#;E
z=_Sdu=_=UgF*oVC{iS$RCo#I(h{o?(0Q?k@kmcea|08>(igk{hoZL>}vJV=5Z0_Ky
zNJ~rCwS7SaA)lkAybf2%uYJP42v;Xf@1t(R<>6PSKOo<Zz4!lr`TvUniyXkxiLRX(
z|NZOJDLg!*9F5jcR+x^+DqiBdckjYPRMpl>WW7v7>An0U6Lyx8(!`k<W~~_8pNu)>
zLHa|}M&6gF#~%<76QA<hBfEI<V$-xBeBb`kw}+hf?%nH|vvWztcBMQP`9pI0G-g5A
z@vR*Gxw6CH5sa(?zW-ONAo8M*vgF<cs^HHJ84h#p)E@0W-f_+k7YlfKd2Rjv<*O`p
zl-DrjXUOX?t58#0yS}qyHeU5?X&^V1o9q^{7%vq`W1afln?c_g%9vi^x!HDYPPf#&
zi(OS!wP-Q3w5hfA(eBy=F|%YGqnPvktmtdDOhnY&(>?0?CN1x;oH})i*LBU{M@!4%
z)DyU!w1P1oaV@V{UA_3;bi_Bk@|bVmzR7DznYPmW`0-<RcPcEsP}hAqzpk-yb0Iy;
zhz1XaQhe*@E%C#(TGTmeY6?2KAF;L-&HXvXJo9rJ84sWP`d;GZ?j6IXA}3EX@iXjh
z9$W0$8+w_ES{cao4+tO_%(raZU7u=hX%Sl;Y>T}n@MOG-ZU~Qzsdct37EAcaIBIWi
z?`rQUeEhCbm(_)<b4IxX*;*M_ZD5jA5BN8g1RhoQF^FrLI={thW;<=LI$G`qqoA2<
zN+SQBYoh8I&emdI*qb+R95OpPIyeLc^BM|%{#4o6+)QUanJd`L5Q*&Wz7O3@w}t6e
zchG?QZZH$qj&EoE5!ml(TWaMoTrm@C3+-lci)I{GC8p(1-cK|(Hl80Ty=Fa_H!<Qe
zHj`Pp^0knmD@D=x$J<LiI~fnt@4_Y^dqqU|VRbEMVXonv9m-~bykR&zhdtFIN+$_*
zoXV~k^E1a-=FiHAdwXMTO${u)FF6?>KKNv3XD>9K6kib+FKdX|-QC?h*3sA3_we-0
z5~_d|cj(VdnlSl{GMeEEcVlCVFYPs#nLnc$i_;Q5lW5lQllxj)KPo&AEZ!D+N4|Yq
z=SM<sVWnMY!zd>whmw|&p`fCAXdKdGflPmoB}{qpiU<Q<mep`eRA0r^V>^*wZ8*E*
z%AUJfmP20-d(KA|VehWhk~y?oDg2qDNFT~1@q<EH{8?%$$E8b`Pzq7J&zhQKRn^q6
zqQu?!`}Dmfc*@qWZ~2q5Nl3&qd+k2!({m?51=rX3l1Vnji8urV2A&rcjiD7VJwr}j
zy*5#+D5eB!$VFAd!uUt&(4PGcA<C~_W^HTL{bzq4+d`ncL!fu9Oen*7Hnx{FH8-{L
z9^V+N@Z4J|U48s32^Zz&=H}{pxM*aWDCsFHFF&>3$ZWR1HD6a(m;H3xqQqrYuiSK@
zJMHbecXk~TJ1Y0@OT)L3#Q`tawwDJ?_KIoC@l7V#j^rwf+P)TF%|d76R@)@Uo~(qT
zx*e1#r)SvffPki%95G_Oy+%Dftx9@-tYz&}*{&Yw_}Yb4@vobO!8(||zR#sBIZmoF
zqQ20b2#Y4VKWLz~^jDLNQ*y=pSP{%zcd_kG!Xx{qR2=fwyYBA%?yfF(wyjK@xe7mH
zAjGL*g;yinkc`JW2fx-Q^kayoICpr|&+ECLimzgJnKTOAUVA2u@aFe8sVJI**oz58
zpXrY!6WBQ@{$6*|Uy-bS>+<~~S2{~oV}#kAEdvdY<bj;&6E&oV;aX}YE^9R<BU}ST
zWIo&zH@KK0ua{DVC@Qg^C}F_OqW)J>u6Q;&f|)cPZ~Tv?_U7)-%w0Fj=|w9HyiI@g
z%*pk+v}%+3TY)Xavph=dJN~uCJ5h}C?YA&D$kg~eQte>rvOcJOfc2Z4wl&F)FwdOI
zz@ESO45p98Xc?tMLsczyP5)q2(I;i?l#y;)eV!1H_;?I(>J;P`?}>l67wg@){QR&i
z>e}Vd=@moM#y9fGm_dsCmh|Jsf|*p#;gBz??)Fu~dv~$#Dbh#ND{u$@y}_HZ8UY2Z
z@Ed=^jZnw^JTWqOE)9(>wbPI6i~l}Dl4P4m!>3gfHKX(-Iab>QX3EB@%#Y9n=iesg
z9_X3R%`1CN+Kw?|e=v(Zqo*Xtk`!3k+&pc8!Tu%@Mq$l{x<Y~7^(ac8gZ)CYA0dqE
zf4ZCcdeV^Cpef5EYFnU7eMaS>VXC<g8P6z{%!$&KwotZwY0M~neWyMcy}hIJ(8HU?
z$1*AQ`x7f{!)h;fj)=j(T|g~K=>32O`X8t~ug8w?>U-9VW3<-UR9O8ttQZ}0e$lAD
zI&Sqjem2B2T8cfPg0Btk?}WL3T}XV|7o%DK;<CbJ0+i?AHRI!cS1mk?>7mz{{uQN-
z;Ib)ri1|u1%CBcvO=B@eJKn+<m*6s$i4aXehI!HycEAy4(#<-V;<xIek7<9mJkKG#
zkZ)`1ESScDiHB^GyC2q1F;8cUE?z_(8O3aNM!CjZ`*T$!x|;(t)7RFg?r}zR^7ngx
zCsrmSA|qL7Iyo$mV4~YAVMIO<j@QmhJUaTOQH}XT_Vzp9H3KQClnZOI1&VE<?ChB5
zG_Sh*gVW9G0cqkRW}I=q_A72)jtq)89tukArWh<X|6m__PQd&&p8x5N@fkxy3ww^0
zl_4co{l8y-3{R1vi0?Mr6@R~<QCetyZpM!ALIEkCSd1{?-##SoTCs&yM7G5Yd!THV
z=zeT|nk2B9wzlWDW5+xuiFEG-n*FoMlj$wZGKnh+^qnuovMht^lftvQr|rE*+*XE*
zHMpplJl3i=Gv~)DrKo~u)C)h}@_(3BJQXVDwe4zeZ|}+Q9N-{>pV;?W#c0PDwpd~N
zeXfjHp~<6z{k1XAPY;f#<8&5IkN0M1YnC`(%E`&0)+w~;xePyFSE;F}&Qeo1LV#oW
zO<qaKbZ=v(XUEfXn=K0FNh_%R{@U7)<!>kR!l-N81jZ{(8qW^}H{NYJL}P0J6acom
z9c}kT3z)Uz_6b>QZMuhac6FI|fBBVl`a@P0kD;Mqe!>HQ0ZA({$Do5Qei1F5(6C5+
z0h^T{%jMdOUlVoSw!HG;P>O5%`SBJl0J4dc%fNS&SBzrseSCcAhOXHR%L10ncc)S8
zjJPnreEE`kIIPXI)*`>RC=fQiYyP$4Jv=*hL!GO_S~^OOn<&X$HRj{RL7~gu_HsIt
zrwyJy6{V%6wHhhe{+@Vnv@@1#+(`QViuOZecYAw&);q6Q05Cic{_uJ3JNgmR)<9ga
zzO^MjR{^LKebus;lF=6i*ZLw?{`$`!%T~iq*go1*XjfBdGP;A3x=XroFY0-TYnRSY
z+wZO!KHXlZFL7G@I(E=k>Fqr^H5D!8eSAC7d$oM$yW6Je_uO2bBAb)rqdeQOgzd#X
z0fQ<Zhs>M-tA4A3g@cFb@4kGw7!({#c=~k7t50VkjP%g{L%nhRDLd{%`U*WJaW~6r
z4ae+v)sMI%dMwz4h5Iu)#*|JLYYMD~guM=h0|NtlhL)Vg^gXvJgQ>36nC8as1nFvO
z($3G%i*I+UQ6Gsc59WVz+nik*Dcvc<{n^*2&Wvo_Pt9Ct$~MBv%ksRgk>U&%5*V!S
zN{jjzo4JLj&u!s8sq(>^o}TXRcb<!P=`40IWAfbkwfSDk=N+dy#7jj-e};JW9YPsI
zZW<b%pO~0XkJp45FNnLCvp(nYeb;nV4t3|k9e#HW7i?HM>El}`+Yn55mb>rJ4Pj$r
zKkogex&4Qbg8d9Ft(=90MTiqjxG}_)f|kAC8vb;4s+zXN)Hj8*5)u&^1HLN6i@qae
z6ot@oyhMTpx-Xw+zADEK$QXE*{&+M&{OULQPwuy)8Rh9vsXRWBXdWIivalHD$D{d{
z)MFl-B#GX~ds)5iPoKVtk8jhqEx&#59s~*G#B_oI1G#3Fb)0B)O0F0&l~sRji~dqI
z?+4?@##CYQ7uW@VE<NtNE_~#FJuZ4w`Z5}wpHZdNKE0(-G30&Php#-zb)(m2_aN@%
zin6_%<FO5kI+C=}J<oDZ9Pe0tbWh{q4OV95ggs9(KMv^_=pPGaD~otAz=#Yzgh~E}
z*CW?hps^3S!40rq;@w4d^g8LrppK2IAcX=ZagwUt(Tr$22ld0Us~hi1FzKPu8}B9Y
z_O~v%!k#H_8Y=W!83%l0jagem*=(*O-%n0>Fvuug=}awOIY@SWr-_B?^M%d=x2oF-
zmsG_tsI~AbyvwY}`!;*D+ileFuhc7DHVdz^BJWhFdjGxNJO{YBp(j*t#!0xfIaz!!
z=YmV0{)adFsmO+|+IaaZL$cQ&ik3<Vg6lf>DWHKfJ$d29)pgbwOz8gjFR;0?Lf_Us
z-M;JUrZXLlc9Ow#G1TL|!>cq5Y@++eVyW*m@?~3o8OV{wv~Qv9_Bq~6@+h#<M)K74
z4x#g03UVyp5zv(O`iSX4Flii@v8oncu}lW)uw#4K?XikXnZ(r%%DKPS<s}Eb-p;4n
zU(A2-OR&nH+9r6)C0+|AVOky2LZ&LTpy_=QPB}J*^+a(f1!;k+zhn9ldJY3x_LP5^
zcub%Lu}%4G4TH{AQC_#}#CmjlK@^zqr}$<%f9Tro0%KfSc2`*r8Ba)Hs=fkIu5qBw
zTO@R7G7(j~S5Lt_``vzs9JQZk+i=LYo=4o-q^{%p_j<psN)3X_4rS-ah!mK9hRa=Z
ztm>*Yauy7i|JwyZ;#z|lGMT8oU!*@c8S{v@Ose}&y&1eFAQsbP^0yZO{MwXZS)<kX
z(NZ`+c~V4>ZnzoiM??hbU?_!fqTTtZcY(HM*(Z-^6PTz{hWo=0GHt{B?^I0IU`C?7
znwPk<qQkg%%sUyPA&u9}JYn+V*Q{e>7!{|#Vy5%Z|INnTl#Cw}Q%lCE!%Q_Iy`n1#
zDOi#aDyA`FfbGPzv7=t$Zp*mA1Z|l(mxAI--H7Lq`7|d<5c04Dm9IQbO7_kudv8!s
zI^%~c9%@FtLFUy|h9b=gxx?h2ppJ(d{Raac1+2_I^re0q!(TLF1Yx-7*pau!&T@#1
znh(Kd<@oxcuA02_x<#8M1sYs*dv|Lw5wA@%FY)ol@{L(om+zc?zoZ;<-Icq*{WWml
zm-Ip$(u<<_BmqkCGXB_YDM}P?vOA4~OvBlwRe<O4eaxz{F&X;UHJtdq`q!krhqy$P
zoddN^WkCoG)k9>c>aGYwa^qQ^IM=&Uv%>p$=o@)(16!vAfB({B^mcpg5-QU2UjZCR
z9dk6euG&C)aJ84luevo_fRK!=acnH{84mt;vyOz4P#!dzUHHMnYVXSYcg2azbt}9@
zY~0T)$HvrWl!xETQ&weY7k+n~S2Z#=<|z@ncrn1ZG4#P0?1%TK)?sJQc`|WR&vW}h
zfz`mx>(^1f_+%+%Wo4O@z$Oq@NIeOPrSfQN(B)X$-e7T(CZh@s;rOY%x<j!u!;hPt
zm-lvks&R9FF^l4?{W{M~o(9*8moFcT-Mn=xomqg`cW-lU6OwcRtA1W+eQ#0!<j1d1
zGMeR2ophP}0E;R7YYl!;yVD#Dnr>vw<XMteEd=(3+bdd@k0}o)w4AF7q~PFp{H=mV
z%HS%uC@(K>-jzZEX)!Uq(5vC$;e|$02I1>4G$Ee2n0(w;$ad6Y)N{~E=-Wffkl26#
z;Y12Hq2ugB_k<G<0IE9|q&IftJm=3cJ#rQ#PaQwYj4h|E%)!d)3uC0$xvvy2x*l11
z__DuW15nlNcyDgLE0wIcxH!+EC+P4{zheGqEbH{~=CKsB7?P#~%s9>9L(=AC(^V0n
zpwuxb9Y`W&OJ+*8*#2AZo@wJ*yUR-YC^BZr*6HR5%|e@zEXgK2Nl1M!u(3TDTZNnJ
zmN**DcYbb7kT?j+RfrV|xulkIf2>d|L;CjoZCEh0{fu~!fcd|q|1<udaRsf^I8QYG
z-y-?M{JdMWr>nbr#Ov2pqvh^c1VltQA0VG!TU%Rj!o|e}(j?JEs-mJ|cYj~6Z31Ez
zMq!&PW{?lY?oeYt(aw6<T9PkzEJ-aX6;Q;P>P<H8gN4)7+<cysGkw1|U7e>Pj9K~N
z!`|C1%!jHp4)*q6Tl5?q1tEa5>d*f8^q-vmY9<3@Ae#qII|yyGcCW&%JDVCKA{*-X
zZ12bYNNmCtKBBCus=AJ(#>U#P0NY}&rjJS*)|^7KvdSqc;*O`KCQJ@mW_^Yv{v0Lc
z`mjTX-S400S1Jx{d=^&8;{?n|`S|#NNf0n@Acm=+k@RpqVuE*7C;W1$@8050w(FE0
zF^;G|?srI$m$e0j#9{)v9QLtYFKujWtgY?#)wgz~DsMJQ9fuztdfdEu^Y=h5N34)_
zl4W0Jk7bHd{9Oo$c7`3uZr{GWwLBCFA)e000Wr*tIDjm_+ij{bW^vrAA^=%81Z16h
zl&4iSfVOHSo%jl@U3GOe!bF`*&Z0U)o<?@k?-qAy)zoeuld#G%{_@qHL<8`0OV6$1
z*(1nCS8=!0cdl>6XdcN~nP}^h-bxP&k@1Q|2P|Uph}0d2Uz#h*Ei9~A1U_Q>c_xF3
z@$1{L(M{jQ<b!J3JH}XKS;I=}19!2H%3N*kb2X3D@*GQ(i|}A}!n>3YNE38(BTU8^
zP%N1bcwrAZ3*ZabsFR3nbHeepsqMU%uJ|s?1Ys4W+X{xUBRT7QWyQlfw!cki=1$Rs
zl|?GbxlvXhkBUo)hQf}kcNmDaZ^mqNm-@dhiElP5m*wxBK?P>Su5A=aqhoht;}nu7
zFb~niD6ZL%Oi8=6J8-(HysvO-;zs-$<gDt++?W>9jw_BXc0HU?u}Dusa@Jw0##PaL
zL9zF;z8GjrxtNu=+1)7hrhtvi$8u1}#PBZUtosphXeUKv6c<J+C_W=eW&9qS1{d$V
z3uBZx3mzo)w-;U_&)A(##vis;$-FeK@DfQ{UCkB5kBqc{I;hJ)Iv`svtZg_K80NgT
zLMDTfwdi|=YdmF6n;f1c>xlqoxM*SOhZj=i@zn`apP>PO_MNdazN!B8u$9}$VY=|j
zkmlR4^NVQal@Dc|H-fM0Q@sjFRX&&n%(AmYFiVx=ht-2tjdI@y{|opAU^-cn7^)fl
z_I-5)g^8=i_{SUZMXURNTpv$!PM$-CIT^u(pj9$~#ns~euM976_i}jIo<Rl=$SZvu
zg1JQyg1p|@zLmC8r>kHH$6pRjn$wLN+w+s}DI=Y*?{PlVX-&=}CZ$TKZmF$d^%Yxh
z^ja3&ptMuPT(5_;^`t_qZYj(=sb&SDC@CVfv%d$}>)=W4u`nOG4cN<X{Xi&jb{dB_
zTmE84Sm<XtRv7F!42_~wu<Rckm526)bv3zcK9O@?cXxGh&?pCalGDhOOuwP^tJbfO
zXd3SDO6s&Z)6<=J`%)kfDv%oo=-?sCw)nzrv&*j2<yO5%8mOw*D~7F%wDQ@Ng|KX=
z?U5*>>9)yq{Ooa_f)}@Fa9ApS*x2yPlM)1d(cGA@xV>e9*Y#-eNW9RQ`T?nLK^tAp
zV;%Get1o2uX!RS&oVCy2*TjW)r=40X4Cjfw&Z&;W3bAKi5VRGEBAwWKmw~Tr=dH!Q
z*K6-lF77({-X>Xx;xnuY983Vv!$9iGroS>(?^FsD-C#asMn%#kFYu#bLM+0O2|+#i
z-rbNm);zf=!GT2+ex`Y8VuWz5C+%kP$8}z@s+kxu7jCNHm-x&cYh{*t-bYUR+lw-?
zvIOVPpC@Jze!W_GQqd42wo-%>T(3IyJ}yqDZ~eA(T<5Ik%i+9n%ZgElvw%ReAwe@R
z&-ag|=`Kd@OSv!hW^n12yzTGrC(0ZuvSS6ZbZfo~UrI_U!SjzM9toYajSY8rczDm^
zn}`U_3J<q!PksF{(H3s><f|o8-^j}T+8Pq~Cp2fezVih)a&+<`nGIUiUnd3i{{4G?
zuR~X)r?aim4%PX#l}FAY2exJFc+gi#dLb)=niqsUL!O?VkdLKDl%dh>E3R@PWGQAj
z=KAMRh1(YO-?U4oR;SEiNv2tH3t%S*+o55*Hvw%S@@N>S0u@b7IiRyGt;V45SF{U~
z)l!v)UB*1mQBmnVkrIj}N$=5`Vn<H~RSY=1K_ZqE&zFK5>(li^lfc+1%_!F?7S&N!
zZEZR7^YiyZ!c8f0@=J!94UJB^<E8*4MqL0oKuE6bm}H7_V*Bg!xBAQ61QHVFys#7$
z6-^-DRDXK|!1(?EW2=6pe?9$<5zS8#uGhDx&d&e*o&)h=GLt4|%eOm^2s!Fmdf&7Q
zGxy8Aj}CS?Bqb9d64Kw=14+eXw9L(FxM=<PD&ugax76}bq1n;?HZcwFc$}q`5WP@r
z#(o8)>X$DEO->p;f4)UYh@KqY2ySs>k%(9&i^De1S?z5tE9UeXEm`1Lskg{K%fPVI
zJ7=3(iCttpBoxXlrDSAev@KLR(yp%WS)X(}ID5cquH0znaI-_eWkshqOW*q!S!kYR
zUnme4xh|`tn^f0<&-aA2uhZF}+ZesE;;i)U0kGY{z}JrV77;)~5+D~>fg+{h^@)p%
zqZ?W&+f)*DS!w+gCi!%9<C^WLc!Jlyd2gn!e3o80_rE6g`ERt5Q5$3|TidhIJF;g+
zji8@v0F7MGh6AJ-+SGm@n}&u)c4uY8(u#<H-~&+MJwqG{QP&?mx{#KZh7vzm%59Ak
ziO9<n94xSw_wtfNh0tG1v72gm3(L;3m+MX_gPev2&B32#W8l>tfB$?|?Tg>HiWL;l
z9iAB6j>&-;MM<id-wg6Dr*Ol29I#1UB_AFc>F?iIDfzSPQS(`#wP>zd+z)Y57GFPw
z=wAs536IRoHhZ;g?*T(EEq%%tkHqND(g59hkgRX^aU(X9m#r9VMv&MVjW%icZT)<4
zOWJE<iy7O_^eIJJfN*TkW&F+FrRfk`Qq-Bg%TbYPXU#{xApjfa-#$h%jqC$e6KZC#
zB}Tx`h=R>9$<0dE%%e<amQ(+Basq-tf+$O8hN<wK@0TuRAsOzROrxD~p9(S^5eVC8
zXIOAqtixl2!qys6T%$#HNI*Y42!P$r-V523xfjNp`Pnm?3>rHs4>f<dXW2H^`if30
zr-f$)=tLo1Aler`xENNRcq=U<==Gh^GFNT!%EjT%B@Of>8rhU}WFf^JU_NOY(%wMr
zDh{0iLi7yNk}bTac5Ij?jK%_0*(>6!+YirjQAPDv;*gQPs+MQfCofjOTsLr<?NB&j
z@fV<Uj=|mK^|ty4W2OJ5xD;Z5vCU!jZgc*%s0;a)4bO8R!`%>oP^X9=*={-V0b~Qs
zX_vS!)*C&;w+WbdhQ;Ul6jW7qG&iCV)DdPJPVfpywax(dO!cgSi-tlukrd&c(IhXB
zXQa54`+#q%mhsz<U0NS6zBae@0#~h5`!R6}(sXArP2>oSVYqFVBw{imaltPXJGx^`
zWPn}C?-5jqYK82!kZHD?#F6XYya3SUoChKVL>tBsb@WMJy9yvBsDeB(`{n!xe4Y<T
znXRQxUk8%i6`G5|+l@k+!#q*M(r$5-oH;E6mpB9qFVp@VY{NwnW<<5)-@sPcD+!Wm
z>1-sKP6pOL_fud%tP0F~noRG;T_$&_SY?L|ESi3{x3Dlelfz-&|8g%2a(}`RiujGo
z0`X`96%$azCl#4jsDhDCR3h`b5M@Zib@GQLE%8YD^9ZaNso+Mov(!BZpiib?KB2N6
zeV572RxI&E19jcLkT-#m@*k|Ni$_MKJG`~53)1!Jxuo1(*QaZ{ZMuj=2Y8N7LBgXn
z6H@=`@3y3COEK~HMTiFlE$b%j0v*GMk?@<~1M;NlN|6uNfFkwC*#xe*(Jw+Y2_J*B
z0h^e|WRqeLuRr4C7kP+bfR#r(L*7~@L0CebijY|L@BTB?UL`|F$bX=AIu3b?%$jCZ
zdbLs53QLtyYX>i6XNgype6sZPZr8c>403kMFnc-<bQG^Q(1Vhll|g_;Z!eS{y+j13
zDK??QG6`1;4VODRRsin-!sS|=KDGX2mEX1joh97={LIuX9SmYYtMULq-qv3-srDNY
zCJrzkh^n8Wn5Ut83@K)~ylT1tWR4zM1y*iymN+=uH!q<x=K%e@y*Pqxcc`#=OUyV;
zNfz(0GQ)zabLt-C>}5V)^KYS>8}&L^K7~6^(utlFtr&5C(>$*k$CdvJvL}(+QUP?%
zJD|FxxnpIlmmdQKAL;lYD=+Vy51ah^pJi^FEyp4NHuXk1YvaCTdQV3Fhcf=$_ssy=
zT6Iw6Cg)0vJO`+h-Vgp}plXtd*#CMM7|`-HADdvUQ=Uh6ZT-!gv!IH5)E|tMIC08H
za<d5u{dvA9WIZT-_wH$z5hv>C=x9XepQ!Z*0zGu$>gH^l@kGtbvi*g0tAXz?eSA>x
zeFZnSJ<2(XaKBXGMT@Q#WUmF6)8HlRbj9>gWPQ5T0+c2Nd9e{yE$^$r!RX3LDQL+A
z!imbK-~@Pe!x!(KhDK9+ycm!B*yqomFR-(3{NT*uHEU=1he{Xe3#AqC`~IC@@_2V@
z!5Y-xE9>j)NFQkYZotM)Cb4yfOzMyqRbOBK!q0DGHrAHkd8uFi$D@2LLQ>L{q9Tnc
zt8l!>Q<+l=ZZ2Y}KQ0rjtzU3@#8udE8&U=7%s1_0`5w4yw%!NyJEQKaSrunfpkshR
z4wH3+$LnF<I=?FZMXb2pXER#1r@r_<F%)#JLAW7aYRcilHTEOz`{|O~r`FeTjv%H>
zuKv;!mJvdCH7F?w4c3a?=?LZDp!JxTm|WoCP*YK<iq!X#b9LP{S-l3j_p4XL1vVp3
z<rX(LHwj5dGIFj63cfu&JVa#uWe>}CwyRg)0r2W7Ig;D)$@~RPFJVVIA3{0nC@(c?
zco9I7MXS179Jv!Ai8t^h_1T8Ch!7tw;=nEw#x&EyQ+WYBx#8_D>9NZj!J%wAQsT6`
zx7Xu5+`?}?D43!Ut?D@j(Fzz|jL_t1iu7KW$3MOsx+JP5_aTeABQ$6xg(>nxK0P2H
zH$8$UGe*?ua;j3i5eNx+)<f?B^E<!xyz|tu87i=5kn)NHL4Cp6ximhjVn671FqL8M
zz}nWgy;wI01%TSPT0nX3+!=5d9PLUK9vx6NKYxx53!)>ozORoO;#y8oQgS;zNP=fG
z2Ab*rNrHJ&FLah@vI6fJ&>f<6#Rrt0FLq-?L$MGi3Ak-c(+#O-XgM_>Fz_0nL_|cc
z7L2&f#XGJHN5d|%r9MN(a_Lea3B6GF?$;=%#oq9@Z%<wS4?3h|O*7*mP;51xAWAh4
zNY`yMN;>%+s6Tq+R7Bu&Qo$Y5=;zqbnEr4KeFZ&v-y$D{hmYS^l=$Y2ea`_gF>$i{
z&dSvS+*7B1e$&idSf*fCVBzO~0~8dBT{)q3yHCG^h^MRy6fxkl=KPtx#`F91aK0Qf
zRf{;xo}r+)1w9ckZDj$L^!rJ2d_Ym=e^RwwLB{kgc?klDbaH!-1(psZ0M;ujZ)Rp@
zCMMRi#(EF7=5J|f(S=DJ3j<fe!Org2ldj%(@ePFGpJBPXuN3>kO5Fdm+!AxuQo%d%
zV?eir_Ph9QWtWZF6kpVZJe!)9ON4o+e>JH~58ul0Br|Y!78)sa$+}m^ec{@*NLWx4
zGczULQfgzY7bc&QQoiTAbR{A#F&L|8r*PTT^$7=tnYW`B_3oS8L9t)BprWDS+PoO?
z=1o&eOL_(~58D17h_X|&FUMYQI`EwjDcfsdERS+lqkbnvZYP%Syw}D{Y}y)i<I|l;
zp_^%6zI@4GzJj*DG4sxO)?<V9WM_n@((yROngB9e0<u5XAOHyvvB*Fj|0%3A>C4oO
zER(!+iO}1-d^bpTWY9b5ZJ>Gdp<`=MJ+}A4Sb`MUQb1(peqaE7dOmg$1gC?8??7MO
z<y|yO$@=l}Ta8jKidm12qQh3KfLH8&K)}4t4h0PwJ!e-?)Bde;?oGuuFRHJ=@&}08
zFYxBoj?riYAxT<Crx0B{LKhb(jI}BfK5<?aJY}+A$Ty_p{2;m!q*=8Vt9jKUCai6q
zdiw(pQA1r<>MaK(L=;9uZ96x!%U`2%L?d?uLQy{pZ#h6Gl!LDyE!&WXJ0LVnWC3R-
zkQn75o*X~cva8<Gi*2MHGp_xFd73T*2(4p-v|LP^3=!TDf&p(2@dlI+T?GLlAfO9D
z`bag4^hs~2rit2bpHMnuIk<?uAh!XDv~^A+eYjLC6I?8HCE|#qWo;{at!Z3oxGO6j
zf)R8il5*=jH>5cndoOP4eaH&0B>MRg$UN8*LxM<970`n8GS%!#)}z4#5sx<x#6Hh-
zwSd_dzKM;<?&!%Fq?c(X%$JP7#Nn3O4zfGY3-QIkT45p{v^tXd$GtJ6DJ!05>G0S`
zsQriMkGbd~Wj*dwAW>D3IMPkgdb7*X`^V0`vs5VFBJGQtck0P8)INE*AJWgbyQZV&
z5LRsf`3BVf?$}4&sMp<hAl41nLR=C`pn{C<?xQ*~2jt$>eo)N41L{GtOu7NmgrG14
zhejn`CZHbJi_tqX;Tt@*5_XW;V)?0z+@S=S@09(C@k=ZE;mtW99+V;6LguT|lze2Q
zoMfQ~3KLIg%JzFC^>#h3Is4xzoy9^PqA!32(}F9spt^cm-WmZaFy3HQt?d%b9#-Io
zcy7PBfL3MRvf%t1-vCjKxQjZ5d{b?ffvC<vVyZM5<oif0Q-9Lma=2Y;NAvTe2_j=A
zt7&{dR@`tXthnytDYluzy?CyQ5p9C1F9{hCB3+}~`^_dRT}S+Mn=)K<dwCQ?g+EFo
z@wQclh^xIRAp;0b^kgqgRk9Kg+C)N26Ei%23}}8{5`2wBOMi1{>EQ&Vg%WH_5wdSh
zk6NCop@3VdrzmuUzIK7GhQ>zv*itd<`5f^z0VZO5h`WKZ4i+sC=A{2pUfn@J{8INW
zn>IIiq8v%cj36vbip57laOvX+{{bxS`ZsLy0U;r5AclPSzyfBx-KE@)_R<RPlVc68
z_xdA#q>NhSv!P-uyFF9Q&CT0(!f29qa3a0(TrOSM*wODRFziX({{j(_@Fg&$^?!Z9
zIolpDub@EeXSfDgBOeGJ`~m{=y%|GUy;)lMLGZJ109J&$`g*>nTfgc46U$N}Mi7p-
z5FC$-OdKYw+@2h`mowbWIt-Z@Fl^lSP`1T0@9B0DA<cfzrC9^ahf6zlf1Y_4-A>nk
zbx$yB)R92O-Kw1%(N`&nx^9(jfT*|sKwDEiI5ZDk06Ia7SAF>Kj;d-Wn}pfFJ1AnN
zz93-x-9#Rk@S+Jn`Y0gSfs+?o_4*HbqWrt>89O`B(I$>5qaQJHV~@8l4sP_nZr^=6
zuLi3%v&!VOfszs_2pw-DBfE<1j3sx+aY~$*pZ=QjJUQA~>d)~7!(UH-LR?&H06A+w
zXs9Wo(<e&CpTFgo4aNs?i4M=R@_(t9e0y<priHbQWZ1Xi(EPF>7Ac+tdd9|ao02y`
zPu=D6Go&%~F9L0Va(ZUwo{mnd`)YY3y=_@F1VL`U&bs~|y%ZvSlvV#>qB!(DYj+;9
zY6q#I<_%ji>@&<gFM#8DObPayfFi2&p1V(;WNRl@>dv125Ee$pDC&s%03sZCMc&~O
zQDr;McWz3i0w4m>l<;&REqJ9;x71nCb~FyG3kij&{r&v~U-abrf5o>sz_&*TzUfW|
zpOTiImMT)fI#ee0|J2Y^<e*t>&n6TbG6a-b(t()$Gy$xk7ygS~!RV8fh$Un3KVCv?
z1(#z4H)et9D*EVv$lub_mtl(Igskaq`4NAFP~`2~x3Ay7Z)gai*DQ6ug8KK^p_{Hl
zyTFPTa+;GKe6hRo^59;fVwDNKz{T|n@wYfF+!uBHT{T|qySQ8=9j6IgfJowD(}`y!
z{`IWCf;z^91}$0?OJj+Hp#<(d4qjVu$5PPKw}6)x$*r&eqHAkE&0e$sAoT!6(*uOH
z%E3hNYz_<zlx_EF!`y+LqWZ%P-zW+9f;(mCs%H@<2zSpMcC;b4BeiYY0mw9oq+znY
z<NQFg+yR<JP@38FsCDas$lC|~=2YXg!N!xt6~>5u9?(q}U#W|*EZyZBcg2*w-pcRy
z)}nNVGA|J|u%DruHUkEOJk1(*lbU~p*e53>d<sneM&hGY@nORn9zGcQtI}uWfl3b$
z+u)$<KMUF(NC2?!Po}irzt6W>t1~*}e9C80RB-AuqSAwqA)WjY@X_v8>4tn%)cfku
zsn0;+ci;fRjyWL!Tc-#tPZ4=7uj`I~ZR;T~we~yTpDdseD2@+0w{agoA}2TJMS36&
z^h~2vlc7`ZJyhmu|8r{F;7^*yo;#ld0!jfM?Qo=+^7%on`3GRvJ}5z-oe`A))*9S^
zju+Z+muFJ;#2{>q?~$tr+cQ0kJ`Wyic!<vHnk}5s1(*GmzlUM*esGXaU}lt&^$7#Q
zO1jPphKWR_J=_Smq;VCDIxjAiEa0*o(qM!Hbco@C?i@5>*>xjf&ynLFd0vqF*^6F8
zo(C2u32EY7V%npHgCWgWZB*5V%#Y$=r5wLA3Jw5rax+NiJ5i5jfxQPKv4|%XvocxA
zSa-D4mU%i%dO%zR%q0;h%Ai;ofc-hmo+i&`FrXgj*I$<r?HoN>2ffhgp8tNDNv}Ab
z1ju{dJQ?9MR`)a0smL8Lgcs4E^Y3#zsV4S$MGVY5WTvCT@7Bw;En%F{k>FIUdkN}1
zI9@t76IVf<53pgX0=uYVFfv)v+=!4746ZJ|nUQ6W;)7GD+DTT|{T+$NOX(*eAYg3#
zZaqJVU9iM6{cKC&^fw_y9RO^!S)@ePZGa9r1Q{Hr|HU^o{%&1j{u~v;6#$<H1Ff_m
z9T<Xu+ckG(Y+n8x%3Ob!PvK$5g_Pvs>7N7w11R>L!Ml`dX7|0@ud<K>E}Px&fU*Ht
zj(9iW83O+|`1<>&Dg+tiX%a{SN&n(N%Bc*T4Ad@}_{+VMn?Hs@jH`+)Rpy<E=Q`bj
zB)ikm<PtPg%;wHqupB?Yb%px~Jy0@I_T7#drCZ|*bnY|fD|Ddjfpm#bToDYMk}Wu#
z`~y$d1W|?hLotp1&h;JIm>y_D57>8&%&^GHOXIr`26Y?l^oYzRGz1eD7l8{h{ONsj
zm$l(R-LgMDUZ|(ic)5>tu&rZoG59?kc>05^goyn<@lICL1%dyfe}d{4`G3Io<F;#K
zbw)YeMPT9Y(J-dL17+cYxHyBiw>LNpni3@37xq=^6f@e|=2CF=Bcj!%Unf#$jgCfe
z*e<7RaDd6{N#MVMeG3ThyA%G*2t?8dS<!$=V9rxXMdb{p@Bx6~OT*NE$9TBuGKUt0
z14W#BoX}G4J>&<VR`iA4W8=Ein^J#7DUu$Jl7Hp%zkt3BI0)MJ7=r6>YV;Q?e}_5@
z?_&=p??XFq)#Fv^pX{rH^8)M@SAlTJ7{NJH_*UK3PmYtnNm7-#iuvI7X45#BsEX@W
zb<yZfpW@-@CY!>DgajoJxB;r?bk()A?t=0i5EQiE(^G?kPx?ABu>(pL(hL8${)(nz
zj~ItnABig24`lfi>D+Rc$jXth_`xxz5~PmuI_L!q5MEv6;Yo$82t29yj?07Fah76E
z3w8R(8`R(mfpU!5y}86x<wVKp+1XAArJxo9ivVbf^|GHoe=aos8$Ho8dPB}v?cc4?
zM;q^b^{j7dfD1^P!W+rj0}g}#Ix8cSl9Ha6m5C!B4gi1GsgP^C8!bE?+Sn($CMIlP
zxj@|e;M~iF$`q)T5QQp&yRNRT+a*Ibr6f=sAhufO(0<f!C1fR}z3mC8fB*hDAc#)k
z;_4pV2$C0j^8FR~q1L8y^fQUa`YI*VY?cBFQ;`HvnBspBBLqwK0j==j#VIJZAW}IR
zt$P0aIUQ~`cr%dlij~H{Bt^OivH1td*8OP2gO;A%g?P|l<M&Q|kG$;UA{>$Bf)h`K
zEp$Rk&y<cR!H?5GaA+typf}JE_p@@eGJt2BdA`bdZciEVTyS5Rg1N(l(m6BG=0_^o
zmf9Zdys5Ki)lsR6u??VkJsv5!%Erc){P}Z^-P9e>bBL;55K;uZdUfIIRZ5@Kgr`B4
zS>^A+IRteJbq*a8jbJM-I#?-vrZ0K^5O_v@;N3ucLZZ2bSG1s;@Otb%=?YcTX#wm*
z;2e|3J~>$xU90k;DMPXja*nyPxNBevMk+U8)sizXC|Fr>f%Rghf!_8M3&x4*9q$OM
zqmB7wR|eJwAtB*MB<32-OKUj51(N80^=BL|5PiDcK67L(<La?Z+Wv6Re6&Au&v<g4
zz6%plqa*7~Bu`Dg+G1DTaY#uDie{eO7%gB$ZJCd~+Cx<SFZv81hrMhc;|51Y9il!M
zP7cv7iRcCS<1#VyhX+<J`utV;0Xc5ts!XH*dw)9kO*sMo?6E!<TxJ3(q(@T*YuCDw
z3e*;cHu}XyA9z`ufKrcZ43o(KM&00~nIF9GAVjd9Yz|?>e$aA2>zI}T&T??at#$t(
z?EZZbpeYdbq0)S)8v#IcKJ5O}V{hN8Z$P3h*(`}yMG~*QQS<xXGG*ZJVD~<*jAN~b
za)j!SH3F_~8m@xH2$1(82uH_2dEL$($zW1^#W&}bE-#USV`Fop=+s$l|0rZvN>+l1
zMVnJnyBW$Aa@nsj``pmI<z{oAp*o}tVF=NxF85LO)+UVHkc{I&VFlu|+lnoDB9n;m
z*@3H~7R(zqjl9I~Jso_&t(N@5@MnoRD^-X*OThz7G9J)P1l{07q<+^K=iU=X2A|lc
zFZf&@aos>v_<(=|5C;l#mVaIX2h21S2ubQOwH6g~7Vg+&AUInDyodQ3NCSOfZHLMa
zLl+w9Yh(t$icSH48jCzd-MPy?ogXOQ54t{{aTZz;#j47%BNrD)LK5XfC}}F<iO2n1
zJ@~!yHGZeSrR+hlz8UmpjL>YxUIdO_x^6iQb=TMUzX17P12isX^RBz}1wLV)RELrI
zW6KelpQd@qaD)&ONhHU#2Xz`)4FaXV`p7P+4uNrHlEsF<jxmS>L@)bHZa^ZAOaEtW
zgW8>!_&3cUhze%|2}uMIP9#7L<iE&=28daUeQa3@w<=2S9^&;r^c2TT!z55Q^(LPA
z2d5C3q2~K(x;@v8)>U_2@$-ndhej(Y<ErQ=y&I{PUv%00(cnpg=ST>7a?z>5{Kr;E
zEJ7X$gar4D#=0g_-_z@xibXS4dKFCgL5YP3AAybpRz;^l%9+XRN2z9<q?(+uKF@6E
zQgvGnIc9sZ6M=`=y1Y!9pkEsNS$`S-gZra1!K!D?hs3LpOdg#H7MZxg|HecN*Adz#
zuJ?aqA||_rFlYfS5L>6J-J4<A&P~)fUcP~vOp3|A$oA|@N=PdSEXB_2uvXQYImK(i
zWb}yn65jYq#V>G`z6y)<N5WcZI#O6?bcr(z#tZHPZ~MOz{Mng(`OA)YiOvFvb;9u)
z*|^iWG{3(9Hy~`#1gB4<f#L0tJXk7Q)c>E@GW^kVDthK-v;Ja5tINc9tsZPyQPiYH
z9Haqn5YSEje<coVhl_wkJ_qGvzKMy<&@dOs&}d*wo0`0j(z2(gAHC<%*Yh0f>nqd?
zl{ucN0=4(E$!quS(^V%rV+sCvPIW_FY+l!3Htpm<K_Sk)r4yK<7`uU3==yYBbR9EB
zY&b4lcmV`D9~kb&BRPG*<91m2XGMbR>(ddOT^AJo>z0e?YqwAFxK5<J7fcC@VfjA7
z)Vy)?bMG&UgWr-*e`L(7p&JK3nERPpSj!LLPLm9mp4$-8x}$+jIC}E}`HYc`A>QlY
z*44vbrUk5hHw<jXIyPHBZECM-mGM++X7=Uxwt{3)D2D%8l`vYtZ)3w=SX8vYU2JxN
zS!%52E6&_KO>Hs`u}ea2r|YEksXbG{AI~1zXv->sv-@Fx|LM|DA=$fk?_RgeSZk$z
z{d(HT$?1u`eJ7pz#dA;W>~3@<NUV3ODNigHj@<!sLd=Mgq9O^D3<!3-yK8a~#Ou)G
zv>k9fcl$!z6%^{#%Ib?N;jH8AxxXAtVTH_AKc=f%TC!@IvA?gbUOPI<1W^(Mdm6F5
zx{#QQi;JH*dp3zSmfz$wP|H}hwzfb(7;Z205Rj3{0A0>?CHDCBL{t<FSW8@FC1w{a
z%yHk>nT^}$oIXpK{b=$!SG*v>!?I~0uDVTL&7OT?ZGGkR>C-BDdZAz&nL;!rAjW$v
zC}`;Eg8|d8e^g;%Vfipn1n_7G%80P@5)G7Xc0W>3Geb#8NVM&r`_4^jdTz39n&Ey=
zZ|#6yP+P7quZ@Z=X4p*qB7H)B6MwjWxgKIVEYznvAvM702NV?vFAW#x4mb(f3<tns
zXApD#RzFLI^6>B&Us`&VuAU_{gwoR1-V8x6t2}ULzUR+vFjjmMRTYI*zIMFCQY5_b
zh6i34Ng73*UEX}H3@4!#L$Li+K^Wx=5{3Ame&Y<AtnBPFE-S;&O2sNIzOf#uCe=*x
z@W0m9!XhHfcI-6Za<gFJx^iJ(DFlq~R<yTqilE-;g}?tN{Yo#DwEF>JVg5HiUbnEg
zj4}WR&2MndLZpF0E%s$Kg5+?9nK|tE=-^XP5e>?~)|SlO-TnUk`%tJu*CD=13|8US
z;P;^9d}siUBJja>X~(BX6Fl>jqBLrgvRBbf^$YgBIyw8w$K(2@?f3KG3d%G*M`B=Z
zehIW6MA>M1#1FN}ED#2L{P+>&^Sn+(Q4)osJVwq1c#f^1aX6=VYN*gw1||^kQ)Mhm
z-@a`Kr8D4=t^?~&j#&p2lzdj#*5ZIP)&Tw^*T6!Y`@O@BvCSf~KgtMEG}oGFW6kiZ
z2NEi^Ve`9XAZUxtCtVhpmlIt0R^x2qGt0^cHJvw~_ua%d{<VvhJ<g_&M)RK$_o8p6
zz0U7(lb|l~%JnTNo}@n&bT*+v{rkfN;-h0}+wWkj`68c#l`NC;HM91s@4)VY2RX_h
z$=22?hVm4LSXzPrl&2@9WWy)kIOH+}PYFCZQlk9EbwPu#o|>p}e+@IWSjITA%jg$X
zX6AkC?R}9d+F5rJnF{*PieP^#d%7TY6}oJJOMksnbFW<V+pss&C2Fq|*Z9{TGy8=F
zP!c7Zn49n!;K+(Lm)!GXTc^EP=1w-k#il_id5%Y7fcY#T&mP3Iv_wKa-Rw6BmiD=s
z7uC2wW$IORHB&K5hf04iBlAqX1Mj4$5sqxe<@B%h8=iA^OR9t@mf;SPoQ_l1O@|aa
zqG7AM=MK+86<>PHb0{K7QTQYMe12Y)wl_l%b;8;tGHPhn2E5lX$q$J<wwAA!#XiN|
zws>}2UNYr7cc&&ROkpeyZt{y2nh6<B`%=Al;N>^@!x~kUcb)8T+gwVx|L=%gf``#+
zpFJV+R%*H)_#LnJA(h5~el6deowD8r_{#XtTcHT6RwF&6?(5XM%XmN2LP?e{XD=6r
z(=6B|=)ZwD<k63Q)LEN!i?zW~E|1s)dXzy!fTRZX1iqMMg*YIb0(sb)>B>QRrQv*}
zW>Lo1XNuQ|4#ym&S{ZsGBRybZ6=>+M3p~i8>bDoC)>fLnNnr3gSVGdn>jGY}=x@xV
zEIhJ8wW*|ir<141eeJwViZbXm_9Q>j_sV{!wICt@gRyWLKG8a_Z!_gx*;*6-L(4Hz
zYcf9dXl*N?VqXb4QlR7vv|?AO7t+?~RqY=nPJl&h=TL>;5mI-Z+WqgO=)NNPKzs3W
zrT=F_P2q#ZN))OCcc(T+iqu!}+TYQTu|g>hH@AFtwYH^NdDv4uI(V8-;E?ai&s$^W
zFy7%zIaq@CoVG2j6w05+r)FWH4)=v8eVf96uUe(fwL&js$NEt*A7V{**5-DJADw%P
zg_><C&wjd~qM^NxtS)l6frt#xdxO1-#Cu*@q@+R*`_~juHAf@c@lfkePJ2ekhaGu4
za%8xtA_VqOK2;U`9&5w0uA)wTnBD%-V`#yL{^V=Tn_c2V9Su^H25&DtXDact?jJWX
zODp`m5zgc3&_~pU_0Kj{9*p_77SIc|6HEW*RmW_w(s(cq;a@Y;8}-TbYv5hMLa`of
zSG{i{|2-6enXxY>R;ZAlS(Oir>^$(ehgGprdW@^SGZCJ@+TQ&AfD?S6Rm!X4We<Mt
zy?YMK5~VJ+?zzu-?E-=;x?sN0tTlK)W$*X(u<j>Phd(Jmp?FIwDJIdr<J`GDf1ktZ
z)nv@XpXDL)GiOjTGBSXtH@PFHpd9W>CSq`e69Y~O2??$3>@<QmBmX_o`E1iRx}wt3
z+TLC&Fqt;BwJo>u4_;Yy#zhyz->4bDmc3B*CJ=}G#GW@+3VL8c{pm3Z40pkR6>(xN
zq*6y)6u@?00{Whjnp@{F|M9_Ur3@JOsJJwIT3dk@k_O=!8#MEnC_pg-DAogx7$3?$
z4b9B_v$J^tTwNY$TEBU#8DOqmSvzAy8Mpj-cd1gw*e#0Zn`<RA%wK+lG{O0;cxiWG
z5fS8!gZk&VM08hTCESHBb8@c3`4Q#~39gfbkfIv`19Rpk9J7IwRe3B5ArdJH_y_;=
zbWZ2h(GbwX^Dh0TV8dC#@Z{mYz?#o+D)r>-sp|S)DYv{yVQm)Aw@wx2@RoGRZf|eP
zDkw+;hcIpL7Qtt3Zcac*_=}E+hzL=~Q7lqY%s>-p@$p9!gX=#gDy6%-`+C}a3klI$
z9=6L0C*#C_;yGVPV{;MFzUZLF%iTlvUQ$EY>zmu#{dBlxoSd#MEiXgV_rkE&?<%df
zjt;T9`ZC!qXD6rIc6K?YZnYW?e2P%EdBf3h>(kp`?96dRLBpM=x16)@49Bs&q8%Te
zmx_SbrPcpe9M;zM_U(rcsbHXU9{)nOCQ30=^XM&gGCa)CdI^D6+Z0Fw@f!}kiaWKr
zpY|u$PF@*r%D`+fwm)u1-TkZzcH=WBUW@LFhfu^55FA_!MjQJ2lH%fz$;r}|-;$G0
z+1S{`%#awbtgLA9-MDe1ZPvirn#9s_Q6}^A)$c2n-k<4V>h78G#t8V3sipUPNlP=g
zkJjOCgsP$Z24G#j7MuMX8ls2WXz|_sd=|=1+eX!8p$m^LznUZ-|N0>P3=>nR*U|os
z&#LDoC7DjY#7119K2W#x5hwvTy#WgazUFgNjbV)t+46#90Un^q)I_fuaAllAy_UV8
z0TKsTsjQ(`lVEgo)O9hVaC~9{j@BV4?~#L+x+z~JVR2o#QV*^v;5AM`5DM0$<dhV&
zpx`-_ii!$~DEK6<l|jIi7%U{~U=T8b2#-ea@fmPRPt47w25ROQ;UbF=oJSwNeIo`I
z<2G2JtgLpG5Gd!jzH=br<L8seiaeHw>!Hx3I|?@MfG?h~{loA?(+U9vAT4d`N5TTX
zhDRIXU`D|7|8n`&WcoPqgWAY17vJ5LFuMtNu<m&>EHuY+%NS>u;QnPwR8?1(T1PY-
zY!RLbzdquRuITzn<VC&4s*>JQ<g`e|z`(!~E-UId1Vay&S1>O4$$2TCUlj04(PRDR
zPfMxur{K#F>lzpuCb-nD6OXgPi6GZOt`plfgS{nYI&{Qmn^sCKEBv27^TUc0_ACt&
ztdVUwoEAa{<jGCBK>aCwI3}mk1*47R8ycF?v+thvnjI>X-}&cHndSOx&!EFFD0<kC
z#G&~GtT2PsN=$zpenkt^;=%D30ZpyhSW9T|0^NezXhUOT-|=ySIFOF&cR{!^Y<hh@
zJ16H^Wu+7;v*fpY9fD_2vr{!P!dOsH059wfSZ-$Tf1-644q#pL*$fMBEe%kBqfc5^
zR`$~_wY~@W*+A(<K9e7B+qpm8F6k};^Ua;Tq5l5I^)#CdN>RKY?|GgLJO$>Lo;sRN
zh6G7ZIxwgE`upQwI7<ktLR(knql&?FQ+RekK?8)^Imq|<N>=*f<@)Hzn8X9%z?ZF+
z5&F>3P{><k6%~7eK$)}tB*MV|NY?EAFD~6hwTtHzEA9jb2j@=`n~*wXWThR+qr(lL
z#yN)OR!k-26yjiLd-CLoPjE01jLVamAMYOa2^1KA<TJ$S$E`TrkT3f)XbDX2%{zDS
zuyJt6si^!qJ3HB@Do>8xrSIG^s(y|Or}4ndoPw{vj)?daIFN0K0~V#{a54{zgao6r
zvonOsw;-WXfA{&Lr6;_5P#_@ews8UOGd@2b*c`zz1@W<#omz^*jeGaPXTx8={z*9_
z2oL#eOsf32`Y=nsGBbZTYHkfQM6LXp>FKJ@&a*(ecLxc)6#)aEw1EK&@aVlMyI_%o
zDs4&%4kc;Wv7p-C?<Enn#Z}p5H$_KwT}yG#eL<kM?Vw^FK^*v49f;qrqodO(uK~xH
zA#w-@u!yoAJ)GCl_vKX2d_1!{IhmwsD(mpY($mh)PRnVe$PNc6(vUwlZbetZ%o(mv
z)+7A=lda$Px_f#IH>R7veEVkQ%xEk|MaOs}-Ug04GJ?;<%)zD)M;+Y)$C)6l6gbWi
zT{I&H7Th`<IduB`3wW7o+S(M$fC!D50S8{q;Goup3l|)3;$G>>ukatK^wt*^77l>*
z^|)cSBasIR3`Qr43nDYGdhO4r-eixsS}LIyTMK6m5#i(GLv`avT7f(fE7WzU5zME_
zHL5?&hdg2q@YTTpmTj9dNAXT61>OB`0V4u`|9;%I2?91b4NZS|no~eQx68(~9M<*W
za`(aknp^R;e_HibRo}oyI0>iW^t_^7I3zqTZEl{1VnK3Ct5-RUL?I*{DCh6cZDedA
z0*K)QM0@e#1<REy8O=^G+j)~%EL>a)6p4w6VNjNkK?9rQ`1rVK;eh6#+F0$^uU}Qv
z)$utvI5O;B-o&D&riPM2qhCq4iyY>Z{AWiWE@$K<C1LC7>8)*T)j{2f-9!y;8$UcW
zsN;t1aU$8dxm|SE-imZTuskr3Icz5m(aP|<``IFod3JUd1|LUKQqmPf<9jM9A5w{j
zI4)i!0LTH2pdQYV5S&U8@-KB>z6mulP`R4dK;m=^CrDv^__ervqrs6&@QxKqYycE=
zKcEYWN=ht3N=r(FN~k~IfQ1pfPU5}ry{)}7VV|kMeg;ti9xvCyg0418^z;K*F^1|4
zVW)+2(2INE9Q?Rfc-!3qfVG|@(c?MPs7!2jNcn4JG(_<kgC&iiDsOfTj>Cd7$?BSy
zMEMOQ=A=$N!;)X_U$%pjph(?D=wrZM0r~`Xu4!B6E09%;q4q&mPOcgPws06;kVydR
z=J~n`>|9(Rk`aURr?41=pP*oaj^Z`MhWZthV{4_Z-y|H)A#6V_`^odC%A9h>Lu0UT
zxNbCa03`KU#qq(cl2>ZK+`23zL<RAJ%Efb9wr{}85jwkie0&UN<^g|z^FiezuX*PM
z&*ZNZ#X6+LI+<>gaI_XNqv$y?mz6q5Q@5ft8;9T!gu|mF@U1+~urp4)<!}~~ZCG{;
zqK>kQOvP(!Puj+Cu(8iDFc9<b@W4S4)I?NVKV`mJ&n_+!L1jtXCZBNwLgdF@vuzkV
zId+jS?Olgr3STCPrw|w3v>aBMvxPSe#_ZS7QNWg10@l)__V#wjou%dE0()KM<Om^y
zkI`vtYs;wL^^<ZMFf}t*_9Xd-u0R|G@H~(CdK(+@;d@1JrUvWWw{)4boExB*@fy_;
zC?2Qq`vXe^!JAcAwREDxm-N{`Hy}3d9~^XUvaC^opuq^d$^E4+wx3$8$rgWo!Uu=0
z)j5a$w@|QbX@BjV$hn*25z3O)!^IBg-`&-Ib@j!b_V-7G75fXEapG=g0d}nUdZR*L
zHUy9ua75hw3<V|Upcp5U7Q!AlYNM!^HVMCZ#y-OpCH!>j3;XDQiF#~ubMsOq$9G(8
zpc0uBs<T{|avHUcM8NL>FA?eR$cWelTiv|E2K=7BK2|X??TE@tt>@PcWCS6^9Dhlz
z_jr5mw*#!TLCJpdEW%r_U%&Q=)O!kA385DAL|0r;&}rbNXzLk-*f=@mBTCa+C(_hX
zpFkudbi|>MxXJON3OrG-A@GE~wQ6Xnw8!byyvF0VTSwgI99~C7UHq1jQNJ?s|5S4B
zF-@Lf6fe<@sw@;7mk`p?iXCZJ3NlG3u#$2aVu6f{6l{ZI$S4;#Wa5R1j1_1L6_{bn
zZf1%@W!Pc{h-)idEOlrVbx77&u+ZH^3hgG6GI87a%>G)kB}@88lQ#X@=Kaq5JI{Nb
zM`tp%1ejR;dOA`}G@P$BQusGc#OCP{^H8~=>|_>;MI0Vg&3(7o$Iovw)(`_jOtwa1
zhD<kd!gfSd?W-$~icr<OVzus^((S)BGvA_A`epUzIBbSOL_s7}W@j3iCIyw4yEp-$
zY&wruh^myJH8>^<9G6NMhJ%L}dLy~q0+5;lXMI(RQ!$Y^e*z3_3Kv~EDUbX)rXr%o
zA&jhDcBK5!6y1yuEiNfhz>s((=#kZ?M@40e40U2MZ4_Z48UkE>azo96>cRcXI=8v-
z4|=#*8;aeD0bBvlDHvylZ25T;LsvN?ZoV)-UtVdwC|PKR;;TM^w2H%CT#p-=13BA-
z=#+F|U?7am-s=7xJ?XCi3X=Ei(Tt6c`y;^a=kM=thXo?G@<j1`x$DXXjtqLr7?z|y
z8V0J<j$M=TPu)IrIX5>~Qq=lXSXyBr)pTHKD8Rz&dUKd5ZB%xhd)D*jS3v0&l$gp*
z&S_nq7d1|p+LJ|W7{4MET+Ytxo3yKk@B3;r8e(4(s9xPB`g(5<D^2oI_+tz!TSrn-
zvE}f|G>G^pn$(jaQ^~k(q+MNI2UdRSdR()s^5k|UJln*mR;yL;R%gA@nVnc(do`|^
zPp(?|ZMR1PaSJLh?k=sixH>|-c|m0=-jR_&Ufuhr_n~qw1wTT-?Gu^2%E9T6>)sZp
ze;@R1+e@}0J;iWjNvYH6+S7Wv!Sh6qmHpvlpLnOmFjrth<D&^GB~>bwsOkD5Ttq>X
zzI0HWc-*+GPcdJ=E>)}5?gyvaR8wbYHw&|eJZ|6X7t9Jv1)+TM9m&}Mg0shFT-1NE
zV<?#CKhW}DU%qh<2clqRJlZ=uC(Bk_I3wrenHd?24PJ5?mkLPTQp+){Zb9I020~m`
zArb}Sxi*AF?;d}v$GKJ}TT%lD->t4L&pt0;d?;Yd|JW*%aRY9_48#`$XJO2n)oM0W
ztg*eFlCxvS!&o+4qA?0C!8?q}<ilA<S6J#{a)Owh>uq&U*is5nH$7*(;1Nnclz`zh
zi#9_WX4-vYHg6!!IGz2!oW~{nIUXB*;zJ@LD90G-4WQ$a@z3Na00*K8aB_bEe4hdi
z<jxwb_bO0PaLnGd#-<mH*k@+v=UYH%M6u`Yn}jWiPoj6sTzmQRB;?8wbc7<pvwq<m
z#8teU&%OkPC_(;?f)R}-({nI(mkvLPMy~Qp=Uw7W_@lZ#;U4Uufdfbcp2}e&0N$AX
zX*4Y^PDIq;_5~8Q$UBhxUKlzYmQFnH^Tt_we0@^!@@uQEM2po571T{ycu(gtFANVm
z8;D=nmgJg9myfkh>kn?kwz`IdS&QCTy@Mt<W&*(`k?M82Sck*GjeniqeH3NDXI8@H
j8}Wa9*Wb46_@AEUcg*Y$ni?#aGE3T)oROsA?`rr99N!{Z

literal 0
HcmV?d00001

diff --git a/docs/usage.rst b/docs/usage.rst
index f57e5706..63ebc032 100644
--- a/docs/usage.rst
+++ b/docs/usage.rst
@@ -1186,10 +1186,12 @@ Plotting with Matplotlib
 ++++++++++++++++++++++++
 .. note::
   - This is an experimental feature. Please report issues.
-  - The context manager ``matplotlib_support`` will temporarily enable this feature
+  - This feature works in Matplotlib versions 2.2.4 and above
+  - Matplotlib is not a dependency of Unyt
 
-Matplotlib is Unyt aware. With no additional effort, Matplotlib will label the x and y
-axes with the units.
+Matplotlib is Unyt aware. After enabling support in :mod:`unyt` using the
+:class:`unyt.matplotlib_support <unyt.mpl_interface.matplotlib_support>` context
+manager, Matplotlib will label the x and y axes with the units.
 
   >>> import matplotlib.pyplot as plt
   >>> from unyt import s, K, matplotlib_support
@@ -1211,11 +1213,24 @@ You can change the plotted units without affecting the original data.
 
 .. image:: _static/mpl_fig2.png
 
-There are three ways to use the context manager.
+It is also possible to set the label style, the choices ``"()"``, ``"[]"`` and
+``"/"`` are supported.
+
+  >>> import matplotlib.pyplot as plt
+  >>> from unyt import s, K, matplotlib_support
+  >>> matplotlib_support.label_style = "[]"
+  >>> with matplotlib_support:
+  ...   plt.plot([0, 1, 2]*s, [3, 4, 5]*K)
+  ...   plt.show()
+  [<matplotlib.lines.Line2D object at ...>]
+
+.. image:: _static/mpl_fig3.png
+
+There are three ways to use the context manager:
 
 1. As a conventional context manager in a ``with`` statement as shown above
 
-2. As a feature toggle in an interactive session
+2. As a feature toggle in an interactive session:
 
   >>> import matplotlib.pyplot as plt
   >>> from unyt import s, K, matplotlib_support
@@ -1225,23 +1240,8 @@ There are three ways to use the context manager.
   >>> plt.show()
   >>> matplotlib_support.disable()
 
-3. As an enable for a complete session
+3. As an enable for a complete session:
 
   >>> import unyt
   >>> unyt.matplotlib_support()
   >>> import matplotlib.pyplot as plt
-
-It is possible to set the label style.
-
-  >>> import matplotlib.pyplot as plt
-  >>> from unyt import s, K, matplotlib_support
-  >>> matplotlib_support.label_style = "[]"
-  >>> with matplotlib_support:
-  ...   plt.plot([0, 1, 2]*s, [3, 4, 5]*K)
-  ...   plt.show()
-  [<matplotlib.lines.Line2D object at ...>]
-
-.. note::
-  
-  - This feature works in Matplotlib versions 2.2.4 and above
-  - Matplotlib is not a dependency of Unyt
diff --git a/unyt/mpl_interface.py b/unyt/mpl_interface.py
index d579cdaf..e82d6db8 100644
--- a/unyt/mpl_interface.py
+++ b/unyt/mpl_interface.py
@@ -111,7 +111,8 @@ def convert(value, unit, axis):
             Raises
             ------
 
-            ConversionError if unit does not have the same dimensions as value
+            UnitConversionError if unit does not have the same dimensions as value or
+            if we don't know how to convert value.
             """
             converted_value = value
             if isinstance(unit, str) or isinstance(unit, Unit):
@@ -127,19 +128,18 @@ def convert(value, unit, axis):
             return converted_value
 
     class matplotlib_support:
-        """Context manager for experimenting with Unyt in Matplotlib
+        """Context manager for setting up integration with Unyt in Matplotlib
 
         Parameters
         ----------
 
-        label_style : string, one from the set {'()', '[]', '/'}
-            The axis label style.
-            '()' -> '(unit)'
-            '[]' -> '[unit]'
-            '/' -> 'q / unit'
-            SI standard where label is a mathematical expression.
-            'q' is a generic quantity symbol and for a value on the axis, x,
-            the equation is q = x * unit.
+        label_style : str
+          One of the following set, ``{'()', '[]', '/'}``. These choices correspond to the following
+          unit labels:
+
+            * ``'()'`` -> ``'(unit)'``
+            * ``'[]'`` -> ``'[unit]'``
+            * ``'/'`` -> ``'q_x / unit'``
         """
 
         def __init__(self, label_style="()"):
@@ -151,14 +151,12 @@ def __call__(self):
 
         @property
         def label_style(self):
-            """label_style : string, one from the set {'()', '[]', '/'}
-                The axis label style.
-                '()' -> '(unit)'
-                '[]' -> '[unit]'
-                '/' -> 'q / unit'
-                SI standard where label is a mathematical expression.
-                'q' is a generic quantity symbol and for a value on the axis, x,
-                the equation is q = x * unit.
+            """str: One of the following set, ``{'()', '[]', '/'}``. These choices correspond to the following
+               unit labels:
+
+                 * ``'()'`` -> ``'(unit)'``
+                 * ``'[]'`` -> ``'[unit]'``
+                 * ``'/'`` -> ``'q_x / unit'``
             """
             return self._labelstyle
 
diff --git a/unyt/tests/test_mpl_interface.py b/unyt/tests/test_mpl_interface.py
index ace8b403..72b93c47 100644
--- a/unyt/tests/test_mpl_interface.py
+++ b/unyt/tests/test_mpl_interface.py
@@ -36,6 +36,7 @@ def test_label(ax):
     assert ax.xaxis.get_label().get_text() == expected_xlabel
     expected_ylabel = "$\\left(\\rm{K}\\right)$"
     assert ax.yaxis.get_label().get_text() == expected_ylabel
+    _matplotlib.pyplot.close()
 
 
 @check_matplotlib
@@ -193,3 +194,4 @@ def test_labelstyle():
     expected_xlabel = "$q_{x}\\;/\\;\\left(\\rm{m} / \\rm{s}\\right)$"
     assert ax.xaxis.get_label().get_text() == expected_xlabel
     _matplotlib.pyplot.close()
+    matplotlib_support.disable()

From 2a4f8748787f6fb71a1aca62623c37637b28ba72 Mon Sep 17 00:00:00 2001
From: Nathan Goldbaum <nathan12343@gmail.com>
Date: Tue, 21 Jan 2020 21:09:22 -0700
Subject: [PATCH 12/14] fix formatting

---
 unyt/mpl_interface.py            | 14 +++++---------
 unyt/tests/test_mpl_interface.py |  5 +----
 2 files changed, 6 insertions(+), 13 deletions(-)

diff --git a/unyt/mpl_interface.py b/unyt/mpl_interface.py
index e82d6db8..16c2ce49 100644
--- a/unyt/mpl_interface.py
+++ b/unyt/mpl_interface.py
@@ -14,11 +14,7 @@
 
 
 try:
-    from matplotlib.units import (
-        ConversionInterface,
-        AxisInfo,
-        registry,
-    )
+    from matplotlib.units import ConversionInterface, AxisInfo, registry
 except ImportError:
     pass
 else:
@@ -134,8 +130,8 @@ class matplotlib_support:
         ----------
 
         label_style : str
-          One of the following set, ``{'()', '[]', '/'}``. These choices correspond to the following
-          unit labels:
+          One of the following set, ``{'()', '[]', '/'}``. These choices
+          correspond to the following unit labels:
 
             * ``'()'`` -> ``'(unit)'``
             * ``'[]'`` -> ``'[unit]'``
@@ -151,8 +147,8 @@ def __call__(self):
 
         @property
         def label_style(self):
-            """str: One of the following set, ``{'()', '[]', '/'}``. These choices correspond to the following
-               unit labels:
+            """str: One of the following set, ``{'()', '[]', '/'}``.
+               These choices correspond to the following unit labels:
 
                  * ``'()'`` -> ``'(unit)'``
                  * ``'[]'`` -> ``'[unit]'``
diff --git a/unyt/tests/test_mpl_interface.py b/unyt/tests/test_mpl_interface.py
index 72b93c47..198c8a45 100644
--- a/unyt/tests/test_mpl_interface.py
+++ b/unyt/tests/test_mpl_interface.py
@@ -119,10 +119,7 @@ def test_list_label(ax):
 def test_errorbar(ax):
     x = unyt_array([8, 9, 10], "cm")
     y = unyt_array([8, 9, 10], "kg")
-    y_scatter = [
-        unyt_array([0.1, 0.2, 0.3], "kg"),
-        unyt_array([0.1, 0.2, 0.3], "kg"),
-    ]
+    y_scatter = [unyt_array([0.1, 0.2, 0.3], "kg"), unyt_array([0.1, 0.2, 0.3], "kg")]
     x_lims = (unyt_quantity(5, "cm"), unyt_quantity(12, "cm"))
     y_lims = (unyt_quantity(5, "kg"), unyt_quantity(12, "kg"))
 

From a2eb1217714540edaa2ac9b6ed7b976094d734f6 Mon Sep 17 00:00:00 2001
From: Nathan Goldbaum <nathan12343@gmail.com>
Date: Tue, 21 Jan 2020 21:19:40 -0700
Subject: [PATCH 13/14] add matplotlib to doc requirements

---
 docs/requirements.txt | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/docs/requirements.txt b/docs/requirements.txt
index 677daf5c..df8aee05 100644
--- a/docs/requirements.txt
+++ b/docs/requirements.txt
@@ -2,4 +2,5 @@ numpy
 sympy
 astropy
 pint
-sphinx==1.7.9
\ No newline at end of file
+matplotlib
+sphinx
\ No newline at end of file

From 131a2c8adab5ad6e5494ae3f67e68a9173c6e005 Mon Sep 17 00:00:00 2001
From: Nathan Goldbaum <nathan12343@gmail.com>
Date: Tue, 21 Jan 2020 21:22:05 -0700
Subject: [PATCH 14/14] slim down doc requirements

---
 docs/requirements.txt | 2 --
 1 file changed, 2 deletions(-)

diff --git a/docs/requirements.txt b/docs/requirements.txt
index df8aee05..5b14b3f3 100644
--- a/docs/requirements.txt
+++ b/docs/requirements.txt
@@ -1,6 +1,4 @@
 numpy
 sympy
-astropy
-pint
 matplotlib
 sphinx
\ No newline at end of file