freealg 0.6.3__py3-none-any.whl → 0.7.1__py3-none-any.whl

freealg/distributions/_marchenko_pastur.py

@@ -13,8 +13,9 @@
 
 import numpy
 from scipy.interpolate import interp1d
-from .._plot_util import plot_density, plot_hilbert, plot_stieltjes, \
-    plot_stieltjes_on_disk, plot_samples
+from .._freeform._plot_util import plot_density, plot_hilbert, \
+    plot_stieltjes, plot_stieltjes_on_disk, plot_samples
+from ..visualization import glue_branches
 
 try:
     from scipy.integrate import cumtrapz
@@ -94,15 +95,20 @@ class MarchenkoPastur(object):
     # init
     # ====
 
-    def __init__(self, lam):
+    def __init__(self, lam, sigma=1.0):
         """
         Initialization.
         """
 
         self.lam = lam
-        self.lam_p = (1 + numpy.sqrt(self.lam))**2
-        self.lam_m = (1 - numpy.sqrt(self.lam))**2
-        self.support = (self.lam_m, self.lam_p)
+        self.sigma = sigma
+
+        # self.lam_p = (1 + numpy.sqrt(self.lam))**2
+        # self.lam_m = (1 - numpy.sqrt(self.lam))**2
+        self.lam_p = sigma**2 * (1.0 + numpy.sqrt(lam))**2
+        self.lam_m = sigma**2 * (1.0 - numpy.sqrt(lam))**2
+
+        self.supp = (self.lam_m, self.lam_p)
 
     # =======
     # density
@@ -117,7 +123,7 @@ class MarchenkoPastur(object):
 
         x : numpy.array, default=None
             The locations where density is evaluated at. If `None`, an interval
-            slightly larger than the support interval of the spectral density
+            slightly larger than the supp interval of the spectral density
             is used.
 
         rho : numpy.array, default=None
@@ -132,7 +138,7 @@ class MarchenkoPastur(object):
 
         save : bool, default=False
             If not `False`, the plot is saved. If a string is given, it is
-            assumed to the save filename (with the file extension). This option
+        assumed to the save filename (with the file extension). This option
             is relevant only if ``plot=True``.
 
         eig : numpy.array, default=None
@@ -168,11 +174,21 @@ class MarchenkoPastur(object):
             x_max = numpy.ceil(center + radius * scale)
             x = numpy.linspace(x_min, x_max, 500)
 
+        # Unpack parameters
+        lam = self.lam
+        lam_p = self.lam_p
+        lam_m = self.lam_m
+        sigma = self.sigma
+
         rho = numpy.zeros_like(x)
-        mask = numpy.logical_and(x >= self.lam_m, x <= self.lam_p)
+        # mask = numpy.logical_and(x >= self.lam_m, x <= self.lam_p)
+        mask = (x > lam_m) & (x < lam_p)
+
+        # rho[mask] = (1.0 / (2.0 * numpy.pi * x[mask] * self.lam)) * \
+        #     numpy.sqrt((self.lam_p - x[mask]) * (x[mask] - self.lam_m))
 
-        rho[mask] = (1.0 / (2.0 * numpy.pi * x[mask] * self.lam)) * \
-            numpy.sqrt((self.lam_p - x[mask]) * (x[mask] - self.lam_m))
+        rho[mask] = numpy.sqrt((lam_p - x[mask]) * (x[mask] - lam_m)) / \
+            (lam * x[mask] * 2 * numpy.pi * sigma**2)
 
         if plot:
             if eig is not None:
@@ -183,6 +199,17 @@ class MarchenkoPastur(object):
 
         return rho
 
+    # =======
+    # support
+    # =======
+
+    def support(self):
+        """
+        supp
+        """
+
+        return [self.supp]
+
     # =======
     # hilbert
     # =======
@@ -196,7 +223,7 @@ class MarchenkoPastur(object):
 
         x : numpy.array, default=None
             The locations where Hilbert transform is evaluated at. If `None`,
-            an interval slightly larger than the support interval of the
+            an interval slightly larger than the supp interval of the
             spectral density is used.
 
         plot : bool, default=False
@@ -257,7 +284,7 @@ class MarchenkoPastur(object):
         hilb = -hilb
 
         if plot:
-            plot_hilbert(x, hilb, support=self.support, latex=latex, save=save)
+            plot_hilbert(x, hilb, support=self.supp, latex=latex, save=save)
 
         return hilb
 
@@ -265,59 +292,126 @@ class MarchenkoPastur(object):
     # m mp numeric vectorized
     # =======================
 
-    def _m_mp_numeric_vectorized(self, z, alt_branch=False, tol=1e-8):
+    # def _m_mp_numeric_vectorized(self, z, alt_branch=False, tol=1e-8):
+    #     """
+    #     Stieltjes transform (principal or secondary branch)
+    #     for Marchenko-Pastur distribution on upper half-plane.
+    #     """
+    #
+    #     sigma = 1.0
+    #     m = numpy.empty_like(z, dtype=complex)
+    #
+    #     # When z is too small, do not use quadratic form.
+    #     mask = numpy.abs(z) < tol
+    #     m[mask] = 1 / (sigma**2 * (1 - self.lam))
+    #
+    #     # Use quadratic form
+    #     not_mask = ~mask
+    #     if numpy.any(not_mask):
+    #
+    #         sign = -1 if alt_branch else 1
+    #         A = self.lam * sigma**2 * z[not_mask]
+    #         B = z[not_mask] - sigma**2 * (1 - self.lam)
+    #         D = B**2 - 4 * A
+    #         sqrtD = numpy.sqrt(D)
+    #         m1 = (-B + sqrtD) / (2 * A)
+    #         m2 = (-B - sqrtD) / (2 * A)
+    #
+    #         # pick correct branch only for non-masked entries
+    #         upper = z[not_mask].imag >= 0
+    #         branch = numpy.empty_like(m1)
+    #         branch[upper] = numpy.where(sign*m1[upper].imag > 0, m1[upper],
+    #                                     m2[upper])
+    #         branch[~upper] = numpy.where(sign*m1[~upper].imag < 0,
+    #                                      m1[~upper], m2[~upper])
+    #         m[not_mask] = branch
+    #
+    #     return m
+
+    # =============
+    # sqrt pos imag
+    # =============
+
+    def _sqrt_pos_imag(self, z):
         """
-        Stieltjes transform (principal or secondary branch)
-        for Marchenko-Pastur distribution on upper half-plane.
+        Square root on a branch cut with always positive imaginary part.
         """
 
-        sigma = 1.0
-        m = numpy.empty_like(z, dtype=complex)
+        sq = numpy.sqrt(z)
+        sq = numpy.where(sq.imag < 0, -sq, sq)
 
-        # When z is too small, do not use quadratic form.
-        mask = numpy.abs(z) < tol
-        m[mask] = 1 / (sigma**2 * (1 - self.lam))
-
-        # Use quadratic form
-        not_mask = ~mask
-        if numpy.any(not_mask):
-
-            sign = -1 if alt_branch else 1
-            A = self.lam * sigma**2 * z[not_mask]
-            B = z[not_mask] - sigma**2 * (1 - self.lam)
-            D = B**2 - 4 * A
-            sqrtD = numpy.sqrt(D)
-            m1 = (-B + sqrtD) / (2 * A)
-            m2 = (-B - sqrtD) / (2 * A)
-
-            # pick correct branch only for non-masked entries
-            upper = z[not_mask].imag >= 0
-            branch = numpy.empty_like(m1)
-            branch[upper] = numpy.where(sign*m1[upper].imag > 0, m1[upper],
-                                        m2[upper])
-            branch[~upper] = numpy.where(sign*m1[~upper].imag < 0, m1[~upper],
-                                         m2[~upper])
-            m[not_mask] = branch
-
-        return m
+        return sq
 
     # ============
     # m mp reflect
     # ============
 
-    def _m_mp_reflect(self, z, alt_branch=False):
+    # def _m_mp_reflect(self, z, alt_branch=False):
+    #     """
+    #     Analytic continuation using Schwarz reflection.
+    #     """
+    #
+    #     mask_p = z.imag >= 0.0
+    #     mask_n = z.imag < 0.0
+    #
+    #     m = numpy.zeros_like(z)
+    #
+    #     f = self._m_mp_numeric_vectorized
+    #     m[mask_p] = f(z[mask_p], alt_branch=False)
+    #     m[mask_n] = f(z[mask_n], alt_branch=alt_branch)
+    #
+    #     return m
+
+    # ================
+    # stieltjes branch
+    # ================
+
+    def _stieltjes_branch(self, z, alt_branch=False, tol=1e-8):
         """
-        Analytic continuation using Schwarz reflection.
         """
 
-        mask_p = z.imag >= 0.0
-        mask_n = z.imag < 0.0
+        # Unpack parameters
+        lam = self.lam
+        sigma = self.sigma
+
+        z = numpy.asarray(z, dtype=complex)
+        m = numpy.empty_like(z, dtype=complex)
+
+        def _eval_upper(zu):
+            mu = numpy.empty_like(zu, dtype=complex)
+
+            mask = numpy.abs(zu) < tol
+            if numpy.any(mask):
+                if alt_branch:
+                    mu[mask] = numpy.inf + 0.0j
+                else:
+                    mu[mask] = 1.0 / (sigma**2 * (1.0 - lam))
 
-        m = numpy.zeros_like(z)
+            not_mask = ~mask
+            if numpy.any(not_mask):
+                sign = -1 if alt_branch else 1
 
-        f = self._m_mp_numeric_vectorized
-        m[mask_p] = f(z[mask_p], alt_branch=False)
-        m[mask_n] = f(z[mask_n], alt_branch=alt_branch)
+                A = lam * sigma**2 * zu[not_mask]
+                B = zu[not_mask] - sigma**2 * (1.0 - lam)
+                D = B**2 - 4.0 * A
+
+                sqrtD = self._sqrt_pos_imag(D)
+
+                r1 = (-B + sqrtD) / (2.0 * A)
+                r2 = (-B - sqrtD) / (2.0 * A)
+
+                mu[not_mask] = numpy.where(sign * r1.imag > 0.0, r1, r2)
+
+            return mu
+
+        mask_p = numpy.imag(z) >= 0.0
+        if numpy.any(mask_p):
+            m[mask_p] = _eval_upper(z[mask_p])
+
+        mask_n = ~mask_p
+        if numpy.any(mask_n):
+            z_ref = numpy.conjugate(z[mask_n])
+            m[mask_n] = numpy.conjugate(_eval_upper(z_ref))
 
         return m
 
@@ -325,8 +419,8 @@ class MarchenkoPastur(object):
     # stieltjes
     # =========
 
-    def stieltjes(self, x=None, y=None, plot=False, on_disk=False, latex=False,
-                  save=False):
+    def stieltjes(self, z=None, x=None, y=None, alt_branch='both', plot=False,
+                  on_disk=False, latex=False, save=False):
         """
         Stieltjes transform of distribution.
 
@@ -335,13 +429,17 @@ class MarchenkoPastur(object):
 
         x : numpy.array, default=None
             The x axis of the grid where the Stieltjes transform is evaluated.
-            If `None`, an interval slightly larger than the support interval of
+            If `None`, an interval slightly larger than the supp interval of
             the spectral density is used.
 
         y : numpy.array, default=None
             The y axis of the grid where the Stieltjes transform is evaluated.
             If `None`, a grid on the interval ``[-1, 1]`` is used.
 
+        alt_branch : {``True``, ``False``, ``'both'``} default=``'both'``
+            If `True`, returns non-physical branch. If `False`, returns
+            physical branch. If ``'both'``, returns both.
+
         plot : bool, default=False
             If `True`, Stieltjes transform is plotted.
 
@@ -407,38 +505,57 @@ class MarchenkoPastur(object):
             # Cayley transform mapping zeta on D to z on H
             z_H = 1j * (1 + zeta) / (1 - zeta)
 
-            m1_D = self._m_mp_reflect(z_H, alt_branch=False)
-            m2_D = self._m_mp_reflect(z_H, alt_branch=True)
+            # m1_D = self._m_mp_reflect(z_H, alt_branch=False)
+            # m2_D = self._m_mp_reflect(z_H, alt_branch=True)
+            m1_D = self._stieltjes_branch(z_H, alt_branch=False)
+            m2_D = self._stieltjes_branch(z_H, alt_branch=True)
+            m12_D = glue_branches(z_H, m1_D, m2_D)
 
-            plot_stieltjes_on_disk(r, t, m1_D, m2_D, support=self.support,
+            plot_stieltjes_on_disk(r, t, m1_D, m12_D, support=self.supp,
                                    latex=latex, save=save)
 
-            return m1_D, m2_D
-
-        # Create x if not given
-        if x is None:
-            radius = 0.5 * (self.lam_p - self.lam_m)
-            center = 0.5 * (self.lam_p + self.lam_m)
-            scale = 2.0
-            x_min = numpy.floor(2.0 * (center - 2.0 * radius * scale)) / 2.0
-            x_max = numpy.ceil(2.0 * (center + 2.0 * radius * scale)) / 2.0
-            x = numpy.linspace(x_min, x_max, 500)
-
-        # Create y if not given
-        if y is None:
-            y = numpy.linspace(-1, 1, 400)
-
-        x_grid, y_grid = numpy.meshgrid(x, y)
-        z = x_grid + 1j * y_grid              # shape (Ny, Nx)
-
-        m1 = self._m_mp_reflect(z, alt_branch=False)
-        m2 = self._m_mp_reflect(z, alt_branch=True)
+            if alt_branch == 'both':
+                return m1_D, m2_D
+            elif alt_branch is True:
+                return m2_D
+            else:
+                return m1_D
+
+        if z is None:
+            # Create x if not given
+            if x is None:
+                radius = 0.5 * (self.lam_p - self.lam_m)
+                center = 0.5 * (self.lam_p + self.lam_m)
+                scale = 2.0
+                x_min = numpy.floor(
+                    2.0 * (center - 2.0 * radius * scale)) / 2.0
+                x_max = numpy.ceil(
+                    2.0 * (center + 2.0 * radius * scale)) / 2.0
+                x = numpy.linspace(x_min, x_max, 500)
+
+            # Create y if not given
+            if y is None:
+                y = numpy.linspace(-1, 1, 400)
+
+            x_grid, y_grid = numpy.meshgrid(x, y)
+            z = x_grid + 1j * y_grid              # shape (Ny, Nx)
+
+        # m1 = self._m_mp_reflect(z, alt_branch=False)
+        # m2 = self._m_mp_reflect(z, alt_branch=True)
+        m1 = self._stieltjes_branch(z, alt_branch=False)
+        m2 = self._stieltjes_branch(z, alt_branch=True)
 
         if plot:
-            plot_stieltjes(x, y, m1, m2, support=self.support, latex=latex,
+            m12 = glue_branches(z, m1, m2)
+            plot_stieltjes(x, y, m1, m12, support=self.supp, latex=latex,
                            save=save)
 
-        return m1, m2
+        if alt_branch == 'both':
+            return m1, m2
+        elif alt_branch is True:
+            return m2
+        else:
+            return m1
 
     # ======
     # sample
@@ -456,11 +573,11 @@ class MarchenkoPastur(object):
             Size of sample.
 
         x_min : float, default=None
-            Minimum of sample values. If `None`, the left edge of the support
+            Minimum of sample values. If `None`, the left edge of the supp
             is used.
 
         x_max : float, default=None
-            Maximum of sample values. If `None`, the right edge of the support
+            Maximum of sample values. If `None`, the right edge of the supp
             is used.
 
         method : {``'mc'``, ``'qmc'``}, default= ``'qmc'``

freealg/distributions/_meixner.py

@@ -13,8 +13,8 @@
 
 import numpy
 from scipy.interpolate import interp1d
-from .._plot_util import plot_density, plot_hilbert, plot_stieltjes, \
-    plot_stieltjes_on_disk, plot_samples
+from .._freeform._plot_util import plot_density, plot_hilbert, \
+    plot_stieltjes, plot_stieltjes_on_disk, plot_samples
 
 try:
     from scipy.integrate import cumtrapz

freealg/distributions/_wachter.py

@@ -13,8 +13,8 @@
 
 import numpy
 from scipy.interpolate import interp1d
-from .._plot_util import plot_density, plot_hilbert, plot_stieltjes, \
-    plot_stieltjes_on_disk, plot_samples
+from .._freeform._plot_util import plot_density, plot_hilbert, \
+    plot_stieltjes, plot_stieltjes_on_disk, plot_samples
 
 try:
     from scipy.integrate import cumtrapz

freealg/distributions/_wigner.py

@@ -13,8 +13,8 @@
 
 import numpy
 from scipy.interpolate import interp1d
-from .._plot_util import plot_density, plot_hilbert, plot_stieltjes, \
-    plot_stieltjes_on_disk, plot_samples
+from .._freeform._plot_util import plot_density, plot_hilbert, \
+    plot_stieltjes, plot_stieltjes_on_disk, plot_samples
 
 try:
     from scipy.integrate import cumtrapz

freealg/visualization/__init__.py

@@ -0,0 +1,12 @@
+# SPDX-FileCopyrightText: Copyright 2026, Siavash Ameli <sameli@berkeley.edu>
+# SPDX-License-Identifier: BSD-3-Clause
+# SPDX-FileType: SOURCE
+#
+# This program is free software: you can redistribute it and/or modify it
+# under the terms of the license found in the LICENSE.txt file in the root
+# directory of this source tree.
+
+from ._rgb_hsv import rgb_hsv
+from ._glue_util import glue_branches
+
+__all__ = ['rgb_hsv', 'glue_branches']

freealg/visualization/_glue_util.py

@@ -0,0 +1,32 @@
+# SPDX-FileCopyrightText: Copyright 2026, Siavash Ameli <sameli@berkeley.edu>
+# SPDX-License-Identifier: BSD-3-Clause
+# SPDX-FileType: SOURCE
+#
+# This program is free software: you can redistribute it and/or modify it
+# under the terms of the license found in the LICENSE.txt file in the root
+# directory of this source tree.
+
+
+# =======
+# Imports
+# =======
+
+import numpy
+
+__all__ = ['glue_branches']
+
+
+# =============
+# glue branches
+# =============
+
+def glue_branches(z, m1, m2):
+    """
+    m12 is the mixing of m1 and m2 where it contains m1 on C^+ and m2 on C^-.
+    """
+
+    m12 = numpy.array(m2, copy=True)
+    mask_p = numpy.imag(z) >= 0.0
+    m12[mask_p] = m1[mask_p]
+
+    return m12

freealg/visualization/_rgb_hsv.py

@@ -0,0 +1,125 @@
+# SPDX-FileCopyrightText: Copyright 2026, Siavash Ameli <sameli@berkeley.edu>
+# SPDX-License-Identifier: BSD-3-Clause
+# SPDX-FileType: SOURCE
+#
+# This program is free software: you can redistribute it and/or modify it under
+# the terms of the license found in the LICENSE.txt file in the root directory
+# of this source tree.
+
+
+# =======
+# Imports
+# =======
+
+import numpy
+import matplotlib
+
+__all__ = ['rgb_hsv']
+
+
+# =======
+# rgb hsv
+# =======
+
+def rgb_hsv(c, shift=0.0, thresh=numpy.inf, n_mod=12.0, n_ph=12.0, vmin=0.35,
+            vmax=1.0, tile_gamma=1.0, tile_mix=1.0):
+    """
+    Convert complex field c to RGB via HSV domain coloring.
+
+    Parameters
+    ----------
+    c : array_like of complex
+        Complex field.
+
+    shift : float, default 0.0
+        Phase offset in turns (1.0 = full 2*pi rotation). Applied to hue.
+
+    thresh : float, default numpy.inf
+        Optional cap on |c| used for magnitude-related terms. Use to prevent
+        very large magnitudes from dominating the encoding.
+
+    n_mod : float, default 12.0
+        Number of modulus steps per 2*pi in log(|c|). Higher -> more concentric
+        rings. Set to 0.0 to disable modulus stepping.
+
+    n_ph : float, default 12.0
+        Number of phase steps per 2*pi in arg(c). Higher -> more angular
+        sectors. Set to 0.0 to disable phase stepping.
+
+    vmin : float, default 0.35
+        Minimum brightness for the tiling shading (darkest parts of tiles).
+
+    vmax : float, default 1.0
+        Maximum brightness for the tiling shading (brightest parts of tiles).
+        Lowering vmax (e.g. 0.8-0.9) can reduce the "neon" look.
+
+    tile_gamma : float, default 1.0
+        Shapes the within-tile ramp. 1.0 = linear sawtooth. >1.0 makes tiles
+        stay darker longer and brighten sharply near boundaries. <1.0 brightens
+        earlier.
+
+    tile_mix : float in [0, 1], default 1.0
+        Mix between original magnitude brightness and tiling shading:
+          0.0 -> value = 1 - exp(-|c|) (your original, no tiling influence)
+          1.0 -> value = tiling shading only (Wegert-style tiling look)
+        Intermediate values overlay tiling onto the original magnitude shading.
+
+    Notes
+    -----
+
+    The coloring technique is inspired from [1]_.
+
+    References
+    ----------
+
+    [1] Wegert, E. (2015) "Visual Complex Functions: An Introduction +with
+        Phase Portraits", Springer.
+        doi: https://doi.org/10.1007/978-3-0348-0180-5
+    """
+
+    hue = (numpy.angle(c) + numpy.pi) / (2.0 * numpy.pi)
+
+    hue = (hue + shift) % 1.0
+
+    r = numpy.abs(c)
+    if numpy.isfinite(thresh):
+        r = numpy.minimum(r, thresh)
+
+    value0 = 1.0 - numpy.exp(-r)
+
+    eps = 1e-300
+    tau = 2.0 * numpy.pi
+
+    g = numpy.ones_like(hue)
+
+    if n_mod and n_mod > 0.0:
+        x_mod = (n_mod / tau) * numpy.log(r + eps)
+        g_mod = numpy.ceil(x_mod) - x_mod
+        g = g * g_mod
+
+    if n_ph and n_ph > 0.0:
+        theta = (numpy.angle(c) + numpy.pi) % tau
+
+        x_ph = (n_ph / tau) * theta
+        g_ph = numpy.ceil(x_ph) - x_ph
+        g = g * g_ph
+
+    g = numpy.clip(g, 0.0, 1.0)
+    if tile_gamma and tile_gamma != 1.0:
+        g = g ** float(tile_gamma)
+
+    vmin = float(numpy.clip(vmin, 0.0, 1.0))
+    vmax = float(numpy.clip(vmax, 0.0, 1.0))
+    if vmax < vmin:
+        vmin, vmax = vmax, vmin
+
+    value_tile = vmin + (vmax - vmin) * g
+
+    tile_mix = float(numpy.clip(tile_mix, 0.0, 1.0))
+    value = (1.0 - tile_mix) * value0 + tile_mix * value_tile
+
+    saturation = numpy.ones_like(hue)
+    hsv = numpy.stack((hue, saturation, numpy.clip(value, 0.0, 1.0)), axis=-1)
+    rgb = matplotlib.colors.hsv_to_rgb(hsv)
+
+    return rgb

{freealg-0.6.3.dist-info → freealg-0.7.1.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: freealg
-Version: 0.6.3
+Version: 0.7.1
 Summary: Free probability for large matrices
 Home-page: https://github.com/ameli/freealg
 Download-URL: https://github.com/ameli/freealg/archive/main.zip

freealg-0.7.1.dist-info/RECORD

@@ -0,0 +1,47 @@
+freealg/__init__.py,sha256=ApA8Dl6Sm29NN_hrT-YRmLFLhrXWEb5N1mc0e0_tDDE,832
+freealg/__version__.py,sha256=2KJZDSMOG7KS82AxYOrZ4ZihYxX0wjfUjDsIZh3L024,22
+freealg/_util.py,sha256=E254DRTCeST7h1QHrfLeyg9jQntRaOyui1oXD7nMaWQ,1890
+freealg/_algebraic_form/__init__.py,sha256=MIB_jVgw2qI-JW_ypqaFSeNAB6c4GvpjNySnap_a6hg,398
+freealg/_algebraic_form/_continuation_algebraic.py,sha256=y2ZKppTM41SBmhD3AsPN02c-MIuJWwrH1EqG385aOdY,13036
+freealg/_algebraic_form/_decompress.py,sha256=gGtixLOVxlMy5S-NsXgoA7lIrB7u7nUZImQk1mIDo3s,21101
+freealg/_algebraic_form/_edge.py,sha256=7l9QyLJDxaEY4WB6MCUFtfEZSf04wyHwH7YPHFJXSbM,10690
+freealg/_algebraic_form/_sheets_util.py,sha256=6OLzWQKu-gN8rxM2rbpbN8TjNZFmD8UJ-8t9kcZdkCo,4174
+freealg/_algebraic_form/algebraic_form.py,sha256=MFWS-zxacFprUWgFWDNdygIk4GFU6Bp-Ed_ns4wVBmk,28298
+freealg/_freeform/__init__.py,sha256=fxHAqlCBK0AqLrK9VhA3FMeid8nutCy5UwsUSgkIhrI,610
+freealg/_freeform/_chebyshev.py,sha256=zkyVA8NLf7uUKlJdLz4ijd_SurdsqUgkA5nHGWSybaE,6916
+freealg/_freeform/_damp.py,sha256=k2vtBtWOxQBf4qXaWu_En81lQBXbEO4QbxxWpvuVhdE,1802
+freealg/_freeform/_decompress.py,sha256=_i37IToZ6oN9DdLXOM8r4y92EbazzcylhcnWwUOpaj0,32108
+freealg/_freeform/_density_util.py,sha256=C_0lKA5QzKAdxPRLJvajNO9zaw2QliZ-n7SI8g2a53M,6745
+freealg/_freeform/_jacobi.py,sha256=z0X6Ws_BEo_h8EQBzDNHGFhLF9F2PUmnGeBVs0bNL7w,10709
+freealg/_freeform/_linalg.py,sha256=U7b_FJwfjU1itO4R4K831Zm2rnUB5eNOdOzvjPSeH98,13142
+freealg/_freeform/_pade.py,sha256=_y89r7rVc2E5lgiN_ZjnxzW2IaVevWwpq5ISor2NVOo,10310
+freealg/_freeform/_plot_util.py,sha256=GKvmc1wjVGeqoomrULPbzBEt6P86FdoR2idBLYh5EDY,20068
+freealg/_freeform/_sample.py,sha256=rhfd_83TCTvvJh8cG8TzEYO4OR8VbtND2YCNtWEhMa8,3205
+freealg/_freeform/_series.py,sha256=33LLCUe4svmV0eWyzhP_XClfDzccQHTW9WBJlYlLfHY,11475
+freealg/_freeform/_support.py,sha256=nxDa2OFlWBgjD0_1qoSMWG7kub6-GIuxIA04n5bdaYw,6614
+freealg/_freeform/freeform.py,sha256=UaRLYSvA0ib5wy1xnLy3_9ndCJ6gUdhNswSm0m9-7go,43571
+freealg/_geometric_form/__init__.py,sha256=mWsXP0nXs3pY8RfUDhPRTgIfhOigKqw_VmoWnJOw2a0,485
+freealg/_geometric_form/_continuation_genus0.py,sha256=4jiXfQaA6w3IhVkJgtKVVjqqtBmavq78FY_YTGUQyY0,4026
+freealg/_geometric_form/_continuation_genus1.py,sha256=X8NZ1_6PxhJJLXZk5ASeGwxej_KwH3-ftuXkBrOFmgU,6021
+freealg/_geometric_form/_elliptic_functions.py,sha256=Rr_pb1A_FjrJlraYQj2G5shdO6f77aVQN2eQzrvIygI,4109
+freealg/_geometric_form/_sphere_maps.py,sha256=NlhTgWXKWXKdyR2dQxMmePsIwHp7IWyYh6uoxgW5Osc,1465
+freealg/_geometric_form/_torus_maps.py,sha256=7m5QsbmnXTWHJE5rWjKG3_TnErHEEQ41vW-3hsOc3lo,3338
+freealg/_geometric_form/geometric_form.py,sha256=whHKYQdakqShtR-jCEugevnje72JEr9M0HSvZ2BYoKs,33379
+freealg/distributions/__init__.py,sha256=POoSAPVazz83y79pfdKeevPtKqpSATNb8S3ocYMdEzI,798
+freealg/distributions/_chiral_block.py,sha256=NxwD-2qx7-ewg5TwMALnPoesUNIWPDktXhyH8RSKi2w,12350
+freealg/distributions/_deformed_marchenko_pastur.py,sha256=cZN3XnAxK4L2zkteE4XL5rlslRz34ArUlInc9TzQIZA,17568
+freealg/distributions/_deformed_wigner.py,sha256=ZC7dj-ViSxaYdHRs45Yx6vf9Yw0WaTm4qZ7FumE0EEo,7239
+freealg/distributions/_kesten_mckay.py,sha256=Abnpz7cRUPvReqbkfXvOh7onrjjvuLDdUUx3wFpu1yA,20088
+freealg/distributions/_marchenko_pastur.py,sha256=G8nzyeZ2kfNIHfv4D2p0uN_9_0Z98JqKvr9PvV9k4o8,20645
+freealg/distributions/_meixner.py,sha256=oxtXpFhEnehPuGHp8d0XEdGCGwmP_KoNJABZRV5GgcQ,17589
+freealg/distributions/_wachter.py,sha256=qn8-_vBoCcpOtsuxXZK6vGd8XABqAgT5wO5FmoSvrc8,17089
+freealg/distributions/_wigner.py,sha256=C8pTM0LMrU79TYIVRvcgPUACG9mwRzMRajS-98YNHPw,16068
+freealg/visualization/__init__.py,sha256=NLq_zwueF7ytZ8sl8zLPqm-AODxxXNvfMozHGmmklcE,435
+freealg/visualization/_glue_util.py,sha256=2oKnEYjUOS4OZfivmciVLauVr53kyHMwi6c2zRKilTQ,693
+freealg/visualization/_rgb_hsv.py,sha256=rEskxXxSlKKxIrHRslVkgxHtD010L3ge9YtcVsOPl8E,3650
+freealg-0.7.1.dist-info/licenses/AUTHORS.txt,sha256=0b67Nz4_JgIzUupHJTAZxu5QdSUM_HRM_X_w4xCb17o,30
+freealg-0.7.1.dist-info/licenses/LICENSE.txt,sha256=J-EEYEtxb3VVf_Bn1TYfWnpY5lMFIM15iLDDcnaDTPA,1443
+freealg-0.7.1.dist-info/METADATA,sha256=xd4C9NmSiyjjLBVQjxFqzMx8VVpQ2rCK-x4tD55afEs,5516
+freealg-0.7.1.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+freealg-0.7.1.dist-info/top_level.txt,sha256=eR2wrgYwDdnnJ9Zf5PruPqe4kQav0GMvRsqct6y00Q8,8
+freealg-0.7.1.dist-info/RECORD,,