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

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

mergeron/core/guidelines_boundary_functions_extra.py

@@ -17,6 +17,7 @@ from scipy.spatial.distance import minkowski as distance_function  # type: ignor
 from sympy import lambdify, simplify, solve, symbols  # type: ignore
 
 from .. import DEFAULT_REC_RATIO, VERSION, ArrayDouble  # noqa: TID252
+from . import GuidelinesBoundary, MPFloat
 from . import guidelines_boundary_functions as gbf
 
 __version__ = VERSION
@@ -529,3 +530,209 @@ def shrratio_boundary_xact_avg_mp(  # noqa: PLR0914
     ) - mp.power(_s_mid, 2)
 
     return gbf.GuidelinesBoundary(bdry, float(mp.nstr(bdry_area_simpson, dps)))
+
+
+# shrratio_boundary_wtd_avg_autoroot
+# this function is about half as fast as the manual one! ... and a touch less precise
+def _shrratio_boundary_wtd_avg_autoroot(  # noqa: PLR0914
+    _delta_star: float = 0.075,
+    _r_val: float = DEFAULT_REC_RATIO,
+    /,
+    *,
+    agg_method: Literal[
+        "arithmetic mean", "geometric mean", "distance"
+    ] = "arithmetic mean",
+    weighting: Literal["own-share", "cross-product-share", None] = "own-share",
+    recapture_form: Literal["inside-out", "proportional"] = "inside-out",
+    dps: int = 5,
+) -> GuidelinesBoundary:
+    """
+    Share combinations on the share-weighted average diversion ratio boundary.
+
+    Parameters
+    ----------
+    _delta_star
+        Share ratio (:math:`\\overline{d} / \\overline{r}`)
+    _r_val
+        recapture ratio
+    agg_method
+        Whether "arithmetic mean", "geometric mean", or "distance".
+    weighting
+        Whether "own-share" or "cross-product-share"  (or None for simple, unweighted average).
+    recapture_form
+        Whether recapture-ratio is MNL-consistent ("inside-out") or has fixed
+        value for both merging firms ("proportional").
+    dps
+        Number of decimal places for rounding returned shares and area.
+
+    Returns
+    -------
+        Array of share-pairs, area under boundary.
+
+    Notes
+    -----
+    An analytical expression for the share-weighted arithmetic mean boundary
+    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)
+
+        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_2
+        )[0]
+        symplot(
+            oswag,
+            (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
+        )[1]
+        symplot(
+            cpwag,
+            (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_2
+        )[0]
+        symplot(
+            oswag,
+            (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
+        )[1]
+        symplot(
+            cpswag,
+            (s_1, 0.0, d_hat / (1 + d_hat)),
+            ylabel=s_2
+        )
+
+
+    """
+
+    _delta_star, _r_val = (mpf(f"{_v}") for _v in (_delta_star, _r_val))
+    _s_mid = mp.fdiv(_delta_star, 1 + _delta_star)
+
+    # initial conditions
+    bdry = [(_s_mid, _s_mid)]
+    s_1_pre, s_2_pre = _s_mid, _s_mid
+    s_2_oddval, s_2_oddsum, s_2_evnsum = True, 0.0, 0.0
+
+    # 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):
+
+        def delta_test(x: MPFloat) -> MPFloat:
+            _de_1 = x / (1 - s_1)
+            _de_2 = (
+                s_1 / (1 - gbf.lerp(s_1, x, _r_val))
+                if recapture_form == "inside-out"
+                else s_1 / (1 - x)
+            )
+            _w = (
+                mp.fdiv(s_1 if weighting == "cross-product-share" else x, s_1 + x)
+                if weighting
+                else 0.5
+            )
+
+            match agg_method:
+                case "geometric mean":
+                    delta_test = mp.expm1(
+                        gbf.lerp(mp.log1p(_de_1), mp.log1p(_de_2), _w)
+                    )
+                case "distance":
+                    delta_test = mp.sqrt(gbf.lerp(_de_1**2, _de_2**2, _w))
+                case _:
+                    delta_test = gbf.lerp(_de_1, _de_2, _w)
+
+            return _delta_star - delta_test
+
+        try:
+            s_2 = mp.findroot(
+                delta_test,
+                x0=(s_2_pre * (1 - theta_), s_2_pre * (1 + theta_)),
+                tol=mp.sqrt(_step_size),
+                solver="ridder",
+            )
+        except (mp.ComplexResult, ValueError, ZeroDivisionError) as _e:
+            print(s_1, s_2_pre)
+            raise _e
+
+        # Build-up boundary points
+        bdry.append((s_1, s_2))
+
+        # Build up area terms
+        s_2_oddsum += s_2 if s_2_oddval else 0
+        s_2_evnsum += s_2 if not s_2_oddval else 0
+        s_2_oddval = not s_2_oddval
+
+        # Hold share points
+        s_2_pre = s_2
+        s_1_pre = s_1
+
+        if (s_1_pre + s_2_pre) > mpf("0.99875"):
+            # Loss of accuracy at 3-9s and up
+            break
+
+    if s_2_oddval:
+        s_2_evnsum -= s_2_pre
+    else:
+        s_2_oddsum -= s_1_pre
+
+    _s_intcpt = gbf._shrratio_boundary_intcpt(
+        s_2_pre,
+        _delta_star,
+        _r_val,
+        recapture_form=recapture_form,
+        agg_method=agg_method,
+        weighting=weighting,
+    )
+
+    if weighting == "own-share":
+        gbd_prtlarea = (
+            _step_size * (4 * s_2_oddsum + 2 * s_2_evnsum + _s_mid + s_2_pre) / 3
+        )
+        # Area under boundary
+        bdry_area_total = float(
+            2 * (s_1_pre + gbd_prtlarea)
+            - (mp.power(_s_mid, "2") + mp.power(s_1_pre, "2"))
+        )
+
+    else:
+        gbd_prtlarea = (
+            _step_size * (4 * s_2_oddsum + 2 * s_2_evnsum + _s_mid + _s_intcpt) / 3
+        )
+        # Area under boundary
+        bdry_area_total = float(2 * gbd_prtlarea - mp.power(_s_mid, "2"))
+
+    bdry.append((mpf("0.0"), _s_intcpt))
+    bdry_array = np.array(bdry, float)
+
+    # Points defining boundary to point-of-symmetry
+    return GuidelinesBoundary(
+        np.vstack((bdry_array[::-1], bdry_array[1:, ::-1]), dtype=float),
+        round(float(bdry_area_total), dps),
+    )