mergeron 2025.739290.6__py3-none-any.whl → 2025.739290.7__py3-none-any.whl

mergeron/__init__.py

@@ -12,7 +12,7 @@ from ruamel import yaml
 
 _PKG_NAME: str = Path(__file__).parent.stem
 
-VERSION = "2025.739290.6"
+VERSION = "2025.739290.7"
 
 __version__ = VERSION
 

mergeron/core/__init__.py

@@ -9,7 +9,14 @@ import mpmath  # type: ignore
 import numpy as np
 from attrs import cmp_using, field, frozen
 
-from .. import VERSION, ArrayBIGINT, this_yaml, yaml_rt_mapper  # noqa: TID252
+from .. import (  # noqa: TID252
+    VERSION,
+    ArrayBIGINT,
+    ArrayDouble,
+    this_yaml,
+    yamelize_attrs,
+    yaml_rt_mapper,
+)
 
 __version__ = VERSION
 
@@ -17,6 +24,17 @@ type MPFloat = mpmath.ctx_mp_python.mpf
 type MPMatrix = mpmath.matrix  # type: ignore
 
 
+@frozen
+class GuidelinesBoundary:
+    """Output of a Guidelines boundary function."""
+
+    coordinates: ArrayDouble
+    """Market-share pairs as Cartesian coordinates of points on the boundary."""
+
+    area: float
+    """Area under the boundary."""
+
+
 @frozen
 class INVTableData:
     industry_group: str
@@ -29,6 +47,7 @@ type INVData = MappingProxyType[
 ]
 type INVData_in = Mapping[str, Mapping[str, Mapping[str, INVTableData]]]
 
+yamelize_attrs(INVTableData)
 
 (_, _) = (
     this_yaml.representer.add_representer(
@@ -59,6 +78,16 @@ type INVData_in = Mapping[str, Mapping[str, Mapping[str, INVTableData]]]
     ),
 )
 
+_, _ = (
+    this_yaml.representer.add_representer(
+        MappingProxyType,
+        lambda _r, _d: _r.represent_mapping("!mappingproxy", dict(_d.items())),
+    ),
+    this_yaml.constructor.add_constructor(
+        "!mappingproxy", lambda _c, _n: MappingProxyType(dict(**yaml_rt_mapper(_c, _n)))
+    ),
+)
+
 
 def _dict_from_mapping(_p: Mapping[Any, Any], /) -> dict[Any, Any]:
     retval: dict[Any, Any] = {}
@@ -76,34 +105,3 @@ def _mappingproxy_from_mapping(_p: Mapping[Any, Any], /) -> MappingProxyType[Any
             else {_k: _v}
         )
     return MappingProxyType(retval)
-
-
-_, _ = (
-    this_yaml.representer.add_representer(
-        MappingProxyType,
-        lambda _r, _d: _r.represent_mapping("!mappingproxy", dict(_d.items())),
-    ),
-    this_yaml.constructor.add_constructor(
-        "!mappingproxy", lambda _c, _n: MappingProxyType(yaml_rt_mapper(_c, _n))
-    ),
-)
-
-
-for _typ in (INVTableData,):
-    _, _ = (
-        this_yaml.representer.add_representer(
-            _typ,
-            lambda _r, _d: _r.represent_mapping(
-                f"!{_d.__class__.__name__}",
-                {
-                    _a.name: getattr(_d, _a.name)
-                    for _a in _d.__attrs_attrs__
-                    if _a.name not in {"coordinates", "area"}
-                },
-            ),
-        ),
-        this_yaml.constructor.add_constructor(
-            f"!{_typ.__name__}",
-            lambda _c, _n: globals()[_n.tag.lstrip("!")](**yaml_rt_mapper(_c, _n)),
-        ),
-    )

mergeron/core/ftc_merger_investigations_data.py

@@ -12,7 +12,7 @@ from __future__ import annotations
 
 import re
 import shutil
-from collections.abc import Sequence
+from collections.abc import Mapping, Sequence
 from operator import itemgetter
 from pathlib import Path
 from types import MappingProxyType
@@ -94,6 +94,10 @@ CNT_FCOUNT_DICT = {
 }
 
 
+def reverse_map(_dict: Mapping[Any, Any]) -> Mapping[Any, Any]:
+    return {_v: _k for _k, _v in _dict.items()}
+
+
 def construct_data(
     _archive_path: Path = INVDATA_ARCHIVE_PATH,
     *,

mergeron/core/guidelines_boundaries.py

@@ -146,6 +146,15 @@ class GuidelinesThresholds:
 class ConcentrationBoundary:
     """Concentration parameters, boundary coordinates, and area under concentration boundary."""
 
+    threshold: float = field(kw_only=False, default=0.01)
+
+    @threshold.validator
+    def _tv(
+        _instance: ConcentrationBoundary, _attribute: Attribute[float], _value: float, /
+    ) -> None:
+        if not 0 <= _value <= 1:
+            raise ValueError("Concentration threshold must lie between 0 and 1.")
+
     measure_name: Literal[
         "ΔHHI",
         "Combined share",
@@ -165,17 +174,8 @@ class ConcentrationBoundary:
         }:
             raise ValueError(f"Invalid name for a concentration measure, {_value!r}.")
 
-    threshold: float = field(kw_only=False, default=0.01)
-
-    @threshold.validator
-    def _tv(
-        _instance: ConcentrationBoundary, _attribute: Attribute[float], _value: float, /
-    ) -> None:
-        if not 0 <= _value <= 1:
-            raise ValueError("Concentration threshold must lie between 0 and 1.")
-
     precision: int = field(
-        kw_only=False, default=5, validator=validators.instance_of(int)
+        kw_only=True, default=5, validator=validators.instance_of(int)
     )
 
     area: float = field(init=False, kw_only=True)
@@ -284,7 +284,7 @@ class DiversionRatioBoundary:
 
     agg_method: UPPAggrSelector = field(
         kw_only=True,
-        default=UPPAggrSelector.MAX,
+        default=UPPAggrSelector.MIN,
         validator=validators.instance_of(UPPAggrSelector),
     )
     """

mergeron/core/guidelines_boundary_functions.py

@@ -1,5 +1,6 @@
 import decimal
-from typing import Any, Literal, TypedDict
+from collections.abc import Callable
+from typing import Literal, TypedDict
 
 import matplotlib as mpl
 import matplotlib.axes as mpa
@@ -7,11 +8,10 @@ import matplotlib.patches as mpp
 import matplotlib.pyplot as plt
 import matplotlib.ticker as mpt
 import numpy as np
-from attrs import frozen
 from mpmath import mp, mpf  # type: ignore
 
 from .. import DEFAULT_REC_RATIO, VERSION, ArrayBIGINT, ArrayDouble  # noqa: TID252
-from . import MPFloat
+from . import GuidelinesBoundary, MPFloat
 
 __version__ = VERSION
 
@@ -28,17 +28,6 @@ class ShareRatioBoundaryKeywords(TypedDict, total=False):
     weighting: Literal["own-share", "cross-product-share", None]
 
 
-@frozen
-class GuidelinesBoundary:
-    """Output of a Guidelines boundary function."""
-
-    coordinates: ArrayDouble
-    """Market-share pairs as Cartesian coordinates of points on the boundary."""
-
-    area: float
-    """Area under the boundary."""
-
-
 def dh_area(_delta_bound: float | MPFloat = 0.01, /, *, dps: int = 9) -> float:
     R"""
     Area under the ΔHHI boundary.
@@ -56,7 +45,7 @@ def dh_area(_delta_bound: float | MPFloat = 0.01, /, *, dps: int = 9) -> float:
     .. math::
 
         2 s1 s_2 &= ΔHHI\\
-        _s_1 + s_2 &= 1
+        s_1 + s_2 &= 1
 
     Parameters
     ----------
@@ -113,7 +102,7 @@ def hhi_delta_boundary(
         np.column_stack((_s_1, _delta_bound / (2 * _s_1))).astype(float),
         np.array([(mpf("0.0"), mpf("1.0"))], float),
     ))
-    bdry = np.vstack((half_bdry[::-1], half_bdry[1:, ::-1]))
+    bdry = np.vstack((half_bdry[::-1], half_bdry[1:, ::-1]), dtype=float)
 
     return GuidelinesBoundary(bdry, dh_area(_delta_bound, dps=dps))
 
@@ -141,12 +130,12 @@ def hhi_pre_contrib_boundary(
 
     step_size = mp.power(10, -dps)
     # Range-limit is 0 less a step, which is -1 * step-size
-    _s_1 = np.array(mp.arange(_s_mid, -step_size, -step_size))
-    s_2 = np.sqrt(_hhi_bound - _s_1**2)
-    half_bdry = np.column_stack((_s_1, s_2)).astype(float)
+    s_1 = np.array(mp.arange(_s_mid, -step_size, -step_size))
+    s_2 = np.sqrt(_hhi_bound - s_1**2)
+    half_bdry = np.column_stack((s_1, s_2)).astype(float)
 
     return GuidelinesBoundary(
-        np.vstack((half_bdry[::-1], half_bdry[1:, ::-1])),
+        np.vstack((half_bdry[::-1], half_bdry[1:, ::-1]), dtype=float),
         round(float(mp.pi * _hhi_bound / 4), dps),
     )
 
@@ -178,7 +167,8 @@ def combined_share_boundary(
 
     _s1 = np.array([0, _s_mid, _s_intcpt], float)
     return GuidelinesBoundary(
-        np.column_stack((_s1, _s1[::-1])), round(float(_s_intcpt * _s_mid), dps)
+        np.array(list(zip(_s1, _s1[::-1])), float),
+        round(float(_s_intcpt * _s_mid), dps),
     )
 
 
@@ -211,6 +201,7 @@ def hhi_post_contrib_boundary(
     )
 
 
+# hand-rolled root finding
 def shrratio_boundary_wtd_avg(  # noqa: PLR0914
     _delta_star: float = 0.075,
     _r_val: float = DEFAULT_REC_RATIO,
@@ -252,58 +243,57 @@ def shrratio_boundary_wtd_avg(  # noqa: PLR0914
     is derived and plotted from y-intercept to the ray of symmetry as follows::
 
         from sympy import plot as symplot, solve, symbols
-        _s_1, s_2 = symbols("_s_1 s_2", positive=True)
+        s_1, s_2 = symbols("s_1 s_2", positive=True)
 
         g_val, r_val, m_val = 0.06, 0.80, 0.30
         delta_star = g_val / (r_val * m_val)
 
         # recapture_form == "inside-out"
         oswag = solve(
-            _s_1 * s_2 / (1 - _s_1)
-            + s_2 * _s_1 / (1 - (r_val * s_2 + (1 - r_val) * _s_1))
-            - (_s_1 + s_2) * delta_star,
+            s_1 * s_2 / (1 - s_1)
+            + s_2 * s_1 / (1 - (r_val * s_2 + (1 - r_val) * s_1))
+            - (s_1 + s_2) * delta_star,
             s_2
         )[0]
         symplot(
             oswag,
-            (_s_1, 0., d_hat / (1 + d_hat)),
+            (s_1, 0., d_hat / (1 + d_hat)),
             ylabel=s_2
         )
 
         cpswag = solve(
-            s_2 * s_2 / (1 - _s_1)
-            + _s_1 * _s_1 / (1 - (r_val * s_2 + (1 - r_val) * _s_1))
-            - (_s_1 + s_2) * delta_star,
+            s_2 * s_2 / (1 - s_1)
+            + s_1 * s_1 / (1 - (r_val * s_2 + (1 - r_val) * s_1))
+            - (s_1 + s_2) * delta_star,
             s_2
         )[1]
         symplot(
             cpwag,
-            (_s_1, 0., d_hat / (1 + d_hat)),
-            ylabel=s_2
+            (s_1, 0.0, d_hat / (1 + d_hat)), ylabel=s_2
         )
 
         # recapture_form == "proportional"
         oswag = solve(
-            _s_1 * s_2 / (1 - _s_1)
-            + s_2 * _s_1 / (1 - s_2)
-            - (_s_1 + s_2) * delta_star,
+            s_1 * s_2 / (1 - s_1)
+            + s_2 * s_1 / (1 - s_2)
+            - (s_1 + s_2) * delta_star,
              s_2
         )[0]
         symplot(
             oswag,
-            (_s_1, 0., d_hat / (1 + d_hat)),
+            (s_1, 0., d_hat / (1 + d_hat)),
             ylabel=s_2
         )
 
         cpswag = solve(
-            s_2 * s_2 / (1 - _s_1)
-            + _s_1 * _s_1 / (1 - s_2)
-            - (_s_1 + s_2) * delta_star,
+            s_2 * s_2 / (1 - s_1)
+            + s_1 * s_1 / (1 - s_2)
+            - (s_1 + s_2) * delta_star,
              s_2
         )[1]
         symplot(
             cpswag,
-            (_s_1, 0.0, d_hat / (1 + d_hat)),
+            (s_1, 0.0, d_hat / (1 + d_hat)),
             ylabel=s_2
         )
 
@@ -321,26 +311,22 @@ def shrratio_boundary_wtd_avg(  # noqa: PLR0914
     # parameters for iteration
     _step_size = mp.power(10, -dps)
     theta_ = _step_size * (10 if weighting == "cross-product-share" else 1)
-    for _s_1 in mp.arange(_s_mid - _step_size, 0, -_step_size):
+    for s_1 in mp.arange(_s_mid - _step_size, 0, -_step_size):
         # The wtd. avg. GUPPI is not always convex to the origin, so we
         #   increment s_2 after each iteration in which our algorithm
         #   finds (s1, s2) on the boundary
         s_2 = s_2_pre * (1 + theta_)
 
-        if (_s_1 + s_2) > mpf("0.99875"):
-            # Loss of accuracy at 3-9s and up
-            break
-
         while True:
-            de_1 = s_2 / (1 - _s_1)
+            de_1 = s_2 / (1 - s_1)
             de_2 = (
-                _s_1 / (1 - lerp(_s_1, s_2, _r_val))
+                s_1 / (1 - lerp(s_1, s_2, _r_val))
                 if recapture_form == "inside-out"
-                else _s_1 / (1 - s_2)
+                else s_1 / (1 - s_2)
             )
 
             r_ = (
-                mp.fdiv(_s_1 if weighting == "cross-product-share" else s_2, _s_1 + s_2)
+                mp.fdiv(s_1 if weighting == "cross-product-share" else s_2, s_1 + s_2)
                 if weighting
                 else 0.5
             )
@@ -365,7 +351,7 @@ def shrratio_boundary_wtd_avg(  # noqa: PLR0914
                 break
 
         # Build-up boundary points
-        bdry.append((_s_1, s_2))
+        bdry.append((s_1, s_2))
 
         # Build up area terms
         s_2_oddsum += s_2 if s_2_oddval else 0
@@ -374,7 +360,11 @@ def shrratio_boundary_wtd_avg(  # noqa: PLR0914
 
         # Hold share points
         s_2_pre = s_2
-        s_1_pre = _s_1
+        s_1_pre = s_1
+
+        if (s_1 + s_2) > mpf("0.99875"):
+            # Loss of accuracy at 3-9s and up
+            break
 
     if s_2_oddval:
         s_2_evnsum -= s_2_pre
@@ -412,7 +402,7 @@ def shrratio_boundary_wtd_avg(  # noqa: PLR0914
 
     # Points defining boundary to point-of-symmetry
     return GuidelinesBoundary(
-        np.vstack((bdry_array[::-1], bdry_array[1:, ::-1])),
+        np.vstack((bdry_array[::-1], bdry_array[1:, ::-1]), dtype=float),
         round(float(bdry_area_total), dps),
     )
 
@@ -436,29 +426,29 @@ def shrratio_boundary_xact_avg(  # noqa: PLR0914
 
         from sympy import latex, plot as symplot, solve, symbols
 
-        _s_1, s_2 = symbols("_s_1 s_2")
+        s_1, s_2 = symbols("s_1 s_2")
 
         g_val, r_val, m_val = 0.06, 0.80, 0.30
         d_hat = g_val / (r_val * m_val)
 
         # recapture_form = "inside-out"
         sag = solve(
-            (s_2 / (1 - _s_1))
-            + (_s_1 / (1 - (r_val * s_2 + (1 - r_val) * _s_1)))
+            (s_2 / (1 - s_1))
+            + (s_1 / (1 - (r_val * s_2 + (1 - r_val) * s_1)))
             - 2 * d_hat,
             s_2
         )[0]
         symplot(
             sag,
-            (_s_1, 0., d_hat / (1 + d_hat)),
+            (s_1, 0., d_hat / (1 + d_hat)),
             ylabel=s_2
         )
 
         # recapture_form = "proportional"
-        sag = solve((s_2/(1 - _s_1)) + (_s_1/(1 - s_2)) - 2 * d_hat, s_2)[0]
+        sag = solve((s_2/(1 - s_1)) + (s_1/(1 - s_2)) - 2 * d_hat, s_2)[0]
         symplot(
             sag,
-            (_s_1, 0., d_hat / (1 + d_hat)),
+            (s_1, 0., d_hat / (1 + d_hat)),
             ylabel=s_2
         )
 
@@ -614,7 +604,7 @@ def shrratio_boundary_min(
     _s_mid = _delta_star / (1 + _delta_star)
 
     if recapture_form == "inside-out":
-        # ## Plot envelope of GUPPI boundaries with rk_ = r_bar if sk_ = min(_s_1, s_2)
+        # ## Plot envelope of GUPPI boundaries with rk_ = r_bar if sk_ = min(s_1, s_2)
         # ## See (si_, sj_) in equation~(44), or thereabouts, in paper
         smin_nr = _delta_star * (1 - _r_val)
         smax_nr = 1 - _delta_star * _r_val
@@ -626,7 +616,7 @@ def shrratio_boundary_min(
         s_1, bdry_area = np.array((0, _s_mid, _s_intcpt), float), _s_mid
 
     return GuidelinesBoundary(
-        np.column_stack((s_1, s_1[::-1])), round(float(bdry_area), dps)
+        np.array(list(zip(s_1, s_1[::-1])), float), round(float(bdry_area), dps)
     )
 
 
@@ -659,7 +649,7 @@ def shrratio_boundary_max(
     _s1_pts = (0, _s_mid, _s_intcpt)
 
     return GuidelinesBoundary(
-        np.column_stack((np.array(_s1_pts, float), np.array(_s1_pts[::-1], float))),
+        np.array(list(zip(_s1_pts, _s1_pts[::-1])), float),
         round(float(_s_intcpt * _s_mid), dps),  # simplified calculation
     )
 
@@ -797,51 +787,57 @@ def round_cust(
     return float(f_ * (n_ / f_).quantize(e_, rounding=rounding_mode))
 
 
-def boundary_plot(*, mktshares_plot_flag: bool = True) -> tuple[Any, ...]:
+def boundary_plot(
+    *,
+    mktshare_plot_flag: bool = True,
+    mktshare_axes_flag: bool = True,
+    backend: Literal["pgf"] | str | None = "pgf",
+) -> tuple[mpl.pyplot, mpl.pyplot.Figure, mpl.axes.Axes, Callable[..., mpl.axes.Axes]]:
     """Setup basic figure and axes for plots of safe harbor boundaries.
 
     See, https://matplotlib.org/stable/tutorials/text/pgf.html
     """
 
-    mpl.use("pgf")
-
-    plt.rcParams.update({
-        "pgf.rcfonts": False,
-        "pgf.texsystem": "lualatex",
-        "pgf.preamble": "\n".join([
-            R"\pdfvariable minorversion=7",
-            R"\usepackage{fontspec}",
-            R"\usepackage{luacode}",
-            R"\begin{luacode}",
-            R"local function embedfull(tfmdata)",
-            R'    tfmdata.embedding = "full"',
-            R"end",
-            R"",
-            R"luatexbase.add_to_callback("
-            R'  "luaotfload.patch_font", embedfull, "embedfull"'
-            R")",
-            R"\end{luacode}",
-            R"\usepackage{mathtools}",
-            R"\usepackage{unicode-math}",
-            R"\setmathfont[math-style=ISO]{STIX Two Math}",
-            R"\setmainfont{STIX Two Text}",
-            r"\setsansfont{Fira Sans Light}",
-            R"\setmonofont[Scale=MatchLowercase,]{Fira Mono}",
-            R"\defaultfontfeatures[\rmfamily]{",
-            R"  Ligatures={TeX, Common},",
-            R"  Numbers={Proportional, Lining},",
-            R"  }",
-            R"\defaultfontfeatures[\sffamily]{",
-            R"  Ligatures={TeX, Common},",
-            R"  Numbers={Monospaced, Lining},",
-            R"  LetterSpace=0.50,",
-            R"  }",
-            R"\usepackage[",
-            R"  activate={true, nocompatibility},",
-            R"  tracking=true,",
-            R"  ]{microtype}",
-        ]),
-    })
+    if backend == "pgf":
+        mpl.use("pgf")
+
+        plt.rcParams.update({
+            "pgf.rcfonts": False,
+            "pgf.texsystem": "lualatex",
+            "pgf.preamble": "\n".join([
+                R"\pdfvariable minorversion=7",
+                R"\usepackage{fontspec}",
+                R"\usepackage{luacode}",
+                R"\begin{luacode}",
+                R"local function embedfull(tfmdata)",
+                R'    tfmdata.embedding = "full"',
+                R"end",
+                R"",
+                R"luatexbase.add_to_callback("
+                R'  "luaotfload.patch_font", embedfull, "embedfull"'
+                R")",
+                R"\end{luacode}",
+                R"\usepackage{mathtools}",
+                R"\usepackage{unicode-math}",
+                R"\setmathfont[math-style=ISO]{STIX Two Math}",
+                R"\setmainfont{STIX Two Text}",
+                r"\setsansfont{Fira Sans Light}",
+                R"\setmonofont[Scale=MatchLowercase,]{Fira Mono}",
+                R"\defaultfontfeatures[\rmfamily]{",
+                R"  Ligatures={TeX, Common},",
+                R"  Numbers={Proportional, Lining},",
+                R"  }",
+                R"\defaultfontfeatures[\sffamily]{",
+                R"  Ligatures={TeX, Common},",
+                R"  Numbers={Monospaced, Lining},",
+                R"  LetterSpace=0.50,",
+                R"  }",
+                R"\usepackage[",
+                R"  activate={true, nocompatibility},",
+                R"  tracking=true,",
+                R"  ]{microtype}",
+            ]),
+        })
 
     # Initialize a canvas with a single figure (set of axes)
     fig_ = plt.figure(figsize=(5, 5), dpi=600)
@@ -851,8 +847,8 @@ def boundary_plot(*, mktshares_plot_flag: bool = True) -> tuple[Any, ...]:
         ax1_: mpa.Axes,
         /,
         *,
-        mktshares_plot_flag: bool = False,
-        mktshares_axlbls_flag: bool = False,
+        mktshare_plot_flag: bool = False,
+        mktshare_axes_flag: bool = False,
     ) -> mpa.Axes:
         # Set the width of axis grid lines, and tick marks:
         # both axes, both major and minor ticks
@@ -883,16 +879,7 @@ def boundary_plot(*, mktshares_plot_flag: bool = True) -> tuple[Any, ...]:
             ax1_.yaxis.get_majorticklabels(), horizontalalignment="right", fontsize=6
         )
 
-        if mktshares_plot_flag:
-            # Axis labels
-            if mktshares_axlbls_flag:
-                # x-axis
-                ax1_.set_xlabel("Firm 1 Market Share, $_s_1$", fontsize=10)
-                ax1_.xaxis.set_label_coords(0.75, -0.1)
-                # y-axis
-                ax1_.set_ylabel("Firm 2 Market Share, $s_2$", fontsize=10)
-                ax1_.yaxis.set_label_coords(-0.1, 0.75)
-
+        if mktshare_plot_flag:
             # Plot the ray of symmetry
             ax1_.plot(
                 [0, 1], [0, 1], linewidth=0.5, linestyle=":", color="grey", zorder=1
@@ -944,8 +931,21 @@ def boundary_plot(*, mktshares_plot_flag: bool = True) -> tuple[Any, ...]:
             for axl_ in ax1_.get_xticklabels(), ax1_.get_yticklabels():
                 plt.setp(axl_[::2], visible=False)
 
+            # Axis labels
+            if mktshare_axes_flag:
+                # x-axis
+                ax1_.set_xlabel("Firm 1 Market Share, $s_1$", fontsize=10)
+                ax1_.xaxis.set_label_coords(0.75, -0.1)
+                # y-axis
+                ax1_.set_ylabel("Firm 2 Market Share, $s_2$", fontsize=10)
+                ax1_.yaxis.set_label_coords(-0.1, 0.75)
+
         return ax1_
 
-    ax_out = _set_axis_def(ax_out, mktshares_plot_flag=mktshares_plot_flag)
+    ax_out = _set_axis_def(
+        ax_out,
+        mktshare_plot_flag=mktshare_plot_flag,
+        mktshare_axes_flag=mktshare_axes_flag,
+    )
 
     return plt, fig_, ax_out, _set_axis_def