passagemath-lcalc 10.6.31rc1__cp311-cp311-macosx_13_0_arm64.whl

passagemath_lcalc-10.6.31rc1.dist-info/METADATA

@@ -0,0 +1,92 @@
+Metadata-Version: 2.4
+Name: passagemath-lcalc
+Version: 10.6.31rc1
+Summary: passagemath: L-function calculations with lcalc
+Author-email: The Sage Developers <sage-support@googlegroups.com>
+Maintainer: Matthias Köppe, passagemath contributors
+License-Expression: GPL-2.0-or-later
+Project-URL: release notes, https://github.com/passagemath/passagemath/releases
+Project-URL: repo (upstream), https://github.com/sagemath/sage
+Project-URL: repo, https://github.com/passagemath/passagemath
+Project-URL: documentation, https://passagemath.org/docs/latest
+Project-URL: homepage (upstream), https://www.sagemath.org
+Project-URL: discourse, https://passagemath.discourse.group
+Project-URL: tracker (upstream), https://github.com/sagemath/sage/issues
+Project-URL: tracker, https://github.com/passagemath/passagemath/issues
+Classifier: Development Status :: 6 - Mature
+Classifier: Intended Audience :: Education
+Classifier: Intended Audience :: Science/Research
+Classifier: Operating System :: POSIX
+Classifier: Operating System :: POSIX :: Linux
+Classifier: Operating System :: MacOS :: MacOS X
+Classifier: Programming Language :: Python :: 3 :: Only
+Classifier: Programming Language :: Python :: 3.10
+Classifier: Programming Language :: Python :: 3.11
+Classifier: Programming Language :: Python :: 3.12
+Classifier: Programming Language :: Python :: 3.13
+Classifier: Programming Language :: Python :: Implementation :: CPython
+Classifier: Topic :: Scientific/Engineering :: Mathematics
+Requires-Python: <3.14,>=3.10
+Description-Content-Type: text/x-rst
+Requires-Dist: passagemath-objects==10.6.31rc1
+Requires-Dist: passagemath-categories==10.6.31rc1
+Requires-Dist: passagemath-modules==10.6.31rc1
+Requires-Dist: memory_allocator
+Requires-Dist: cysignals<1.12.4; sys_platform == "win32"
+Requires-Dist: cysignals!=1.12.0,>=1.11.2
+Provides-Extra: test
+
+============================================================================
+ passagemath: L-function calculations with lcalc
+============================================================================
+
+`passagemath <https://github.com/passagemath/passagemath>`__ is open
+source mathematical software in Python, released under the GNU General
+Public Licence GPLv2+.
+
+It is a fork of `SageMath <https://www.sagemath.org/>`__, which has been
+developed 2005-2025 under the motto “Creating a Viable Open Source
+Alternative to Magma, Maple, Mathematica, and MATLAB”.
+
+The passagemath fork uses the motto "Creating a Free Passage Between the
+Scientific Python Ecosystem and Mathematical Software Communities."
+It was created in October 2024 with the following goals:
+
+-  providing modularized installation with pip,
+-  establishing first-class membership in the scientific Python
+   ecosystem,
+-  giving `clear attribution of upstream
+   projects <https://groups.google.com/g/sage-devel/c/6HO1HEtL1Fs/m/G002rPGpAAAJ>`__,
+-  providing independently usable Python interfaces to upstream
+   libraries,
+-  offering `platform portability and integration testing
+   services <https://github.com/passagemath/passagemath/issues/704>`__
+   to upstream projects,
+-  inviting collaborations with upstream projects,
+-  `building a professional, respectful, inclusive
+   community <https://groups.google.com/g/sage-devel/c/xBzaINHWwUQ>`__,
+-  `empowering Sage users to participate in the scientific Python ecosystem
+   <https://github.com/passagemath/passagemath/issues/248>`__ by publishing packages,
+-  developing a port to `Pyodide <https://pyodide.org/en/stable/>`__ for
+   serverless deployment with Javascript,
+-  developing a native Windows port.
+
+`Full documentation <https://passagemath.org/docs/latest/html/en/index.html>`__ is
+available online.
+
+passagemath attempts to support and provides binary wheels suitable for
+all major Linux distributions and recent versions of macOS.
+
+Binary wheels for native Windows (x86_64) are are available for a subset of
+the passagemath distributions. Use of the full functionality of passagemath
+on Windows currently requires the use of Windows Subsystem for Linux (WSL)
+or virtualization.
+
+The supported Python versions in the passagemath 10.6.x series are 3.10.x-3.13.x.
+
+
+About this pip-installable distribution package
+-----------------------------------------------
+
+This pip-installable distribution ``passagemath-lcalc`` provides
+an interface to Michael Rubinstein's L-function calculator.

passagemath_lcalc-10.6.31rc1.dist-info/RECORD

@@ -0,0 +1,19 @@
+passagemath_lcalc-10.6.31rc1.dist-info/RECORD,,
+passagemath_lcalc-10.6.31rc1.dist-info/WHEEL,sha256=_2eWE0XUrAnNpxPyXffVnsvXkR9SkDAS3wnEDHyQ9s4,136
+passagemath_lcalc-10.6.31rc1.dist-info/top_level.txt,sha256=Kmzulf9WsphADFQuqgvdy5mvTLDj_V2zkFHU2s3UXos,6
+passagemath_lcalc-10.6.31rc1.dist-info/METADATA,sha256=jJieUmUwbCRyaCF6P0QriJqqcKtENXbgxn-cGNQ7tUQ,4369
+passagemath_lcalc.dylibs/libgf2x.3.dylib,sha256=nb16UZWmpvYh0YMgbx-jLA7SzSwltQ-edWN74AN17wo,137056
+passagemath_lcalc.dylibs/libntl.44.dylib,sha256=8jHzjynRpmMwuhIOq9smjTBkLKbnPVMzakF1E4naKpc,2170096
+passagemath_lcalc.dylibs/libpari-gmp-tls.dylib,sha256=8b_DslZNOGsiLdglj_EmZpYWfLk6TDVaQXQyDQwwWHg,9062448
+passagemath_lcalc.dylibs/libgmp.10.dylib,sha256=DgL3lBxg9OcC_p7PaNkTvNxKG63-zrCIMQ4eEiab3UY,464688
+passagemath_lcalc.dylibs/libmpfr.6.dylib,sha256=UIeJHDY9IKWFhqDYXAn182xRB0aBaHHz_yXjuuX1Ojw,466160
+passagemath_lcalc.dylibs/libLfunction.2.0.0.dylib,sha256=pEVlV8MrcHQhkNNfT2QkVQNZK-SO5KJy0RFMLqIkUy0,979376
+sage/all__sagemath_lcalc.py,sha256=D9iHi3nLdsQCdZPBeiqScQWx038s3f7AAakrOTZp6kE,64
+sage/lfunctions/all__sagemath_lcalc.py,sha256=SoqrpDKr69vZp_s5EhcbUzvMZ1xWT7kGPVFzwriXHdA,155
+sage/lfunctions/lcalc.py,sha256=r8AAsv6jZPEoTB7nr1YJ4l012erjpxUNBLsyZsnEgck,12904
+sage/libs/all__sagemath_lcalc.py,sha256=KQWtVG_LW-OnPa5fWIkPbhXFcSDwkvAILxVmAzMNJyw,44
+sage/libs/lcalc/lcalc_Lfunction.cpython-311-darwin.so,sha256=QEoLWxE2F-2CCsa_wcLS9voAJUfqtFSiZI9B35VDA7w,580768
+sage/libs/lcalc/__init__.py,sha256=KQWtVG_LW-OnPa5fWIkPbhXFcSDwkvAILxVmAzMNJyw,44
+sage/libs/lcalc/lcalc_sage.h,sha256=gcR7C1PK7mt_xFnDKlBjCvmBMfy2O_Om4fTNZ5Xg6Po,1365
+sage/libs/lcalc/lcalc_Lfunction.pyx,sha256=JHuW-AID04L0mUbwIjY7PSqQEgzNklrKd6CyFM0vL9c,37362
+sage/libs/lcalc/lcalc_Lfunction.pxd,sha256=8NThJ5A6o1ASc5Y5pZ6ofA1AjlRM6dX-ypv4imz5MXM,5509

passagemath_lcalc-10.6.31rc1.dist-info/WHEEL

@@ -0,0 +1,6 @@
+Wheel-Version: 1.0
+Generator: setuptools (80.9.0)
+Root-Is-Purelib: false
+Tag: cp311-cp311-macosx_13_0_arm64
+Generator: delocate 0.13.0
+

passagemath_lcalc-10.6.31rc1.dist-info/top_level.txt

@@ -0,0 +1,2 @@
+
+sage

sage/all__sagemath_lcalc.py

@@ -0,0 +1,2 @@
+# sage_setup: distribution = sagemath-lcalc
+# delvewheel: patch

sage/lfunctions/all__sagemath_lcalc.py

@@ -0,0 +1,7 @@
+# sage_setup: distribution = sagemath-lcalc
+
+from sage.misc.lazy_import import lazy_import
+
+lazy_import("sage.lfunctions.lcalc", "lcalc")
+
+del lazy_import

sage/lfunctions/lcalc.py

@@ -0,0 +1,413 @@
+# sage_setup: distribution = sagemath-lcalc
+r"""
+Rubinstein's `L`-function calculator
+
+This interface provides complete
+access to Rubinstein's lcalc calculator with extra PARI
+functionality compiled in
+and is a standard part of Sage.
+
+.. NOTE::
+
+   Each call to ``lcalc`` runs a complete
+   ``lcalc`` process. On a typical Linux system, this
+   entails about 0.3 seconds overhead.
+
+AUTHORS:
+
+- Michael Rubinstein (2005): released under GPL the C++ program lcalc
+
+- William Stein (2006-03-05): wrote Sage interface to lcalc
+"""
+
+# ######################################################################
+#       Copyright (C) 2006 William Stein <wstein@gmail.com>
+#
+#  Distributed under the terms of the GNU General Public License (GPL)
+#
+#                  https://www.gnu.org/licenses/
+# ######################################################################
+
+import os
+
+from sage.structure.sage_object import SageObject
+from sage.misc.lazy_import import lazy_import
+from sage.misc.pager import pager
+from sage.rings.integer_ring import ZZ
+from sage.rings.rational_field import QQ
+
+lazy_import('sage.rings.complex_mpfr', 'ComplexField')
+lazy_import('sage.rings.real_mpfr', 'RealField')
+lazy_import('sage.schemes.elliptic_curves.ell_generic', 'EllipticCurve_generic', as_='EllipticCurve')
+
+prec = 32
+
+
+class LCalc(SageObject):
+    r"""
+    Rubinstein's `L`-functions Calculator.
+
+    Type ``lcalc.[tab]`` for a list of useful commands that
+    are implemented using the command line interface, but return
+    objects that make sense in Sage. For each command the possible
+    inputs for the `L`-function are:
+
+
+    - ``"`` -- (default) the Riemann zeta function
+
+    - ``'tau'`` -- the L function of the Ramanujan delta function
+
+    - ``E`` -- an elliptic curve over `\QQ`; defines `L(E,s)`
+
+    You can also use the complete command-line interface of
+    Rubinstein's `L`-functions calculations program via this
+    class. Type ``lcalc.help()`` for a list of commands and
+    how to call them.
+    """
+    def _repr_(self):
+        return "Rubinsteins L-function Calculator"
+
+    def __call__(self, args):
+        cmd = 'lcalc %s' % args
+        with os.popen(cmd) as f:
+            res = f.read().strip()
+        return res
+
+    def _compute_L(self, L):
+        if isinstance(L, str):
+            if L == 'tau':
+                return '--tau'
+            return L
+        if isinstance(L, EllipticCurve):
+            if L.base_ring() == QQ:
+                L = L.minimal_model()
+                return '-e --a1 %s --a2 %s --a3 %s --a4 %s --a6 %s' % tuple(L.a_invariants())
+        raise TypeError("$L$-function of %s not known" % L)
+
+    def help(self):
+        try:
+            h = self.__help
+        except AttributeError:
+            h = "-" * 70 + '\n'
+            h += "   Call lcalc with one argument, e.g., \n"
+            h += "      sage: lcalc('--tau -z 1000')\n"
+            h += "   is translated into the command line\n"
+            h += "      $ lcalc --tau -z 1000\n"
+            h += "-" * 70 + '\n'
+            h += '\n' + self('--help')
+            self.__help = h
+        pager()(h)
+
+    def zeros(self, n, L=''):
+        """
+        Return the imaginary parts of the first `n` nontrivial
+        zeros of the `L`-function in the upper half plane, as
+        32-bit reals.
+
+        INPUT:
+
+        - ``n`` -- integer
+
+        - ``L`` -- defines `L`-function (default: Riemann zeta function)
+
+        This function also checks the Riemann Hypothesis and makes sure no
+        zeros are missed. This means it looks for several dozen zeros to
+        make sure none have been missed before outputting any zeros at all,
+        so takes longer than
+        ``self.zeros_of_zeta_in_interval(...)``.
+
+        EXAMPLES::
+
+            sage: lcalc.zeros(4)                           # long time
+            [14.1347251, 21.0220396, 25.0108576, 30.4248761]
+            sage: lcalc.zeros(5, L='--tau')                # long time
+            [9.22237940, 13.9075499, 17.4427770, 19.6565131, 22.3361036]
+            sage: lcalc.zeros(3, EllipticCurve('37a'))     # long time
+            [0.000000000, 5.00317001, 6.87039122]
+        """
+        L = self._compute_L(L)
+        RR = RealField(prec)
+        X = self('-z %s %s' % (int(n), L))
+        return [RR(z) for z in X.split()]
+
+    def zeros_in_interval(self, x, y, stepsize, L=''):
+        r"""
+        Return the imaginary parts of (most of) the nontrivial zeros of the
+        `L`-function on the line `\Re(s)=1/2` with positive
+        imaginary part between `x` and `y`, along with a
+        technical quantity for each.
+
+        INPUT:
+
+        - ``x, y, stepsize`` -- positive floating point numbers
+
+        - ``L`` -- defines `L`-function (default: Riemann zeta function)
+
+        OUTPUT: list of pairs (zero, S(T)).
+
+        Rubinstein writes: The first column outputs the imaginary part of
+        the zero, the second column a quantity related to `S(T)`
+        (it increases roughly by 2 whenever a sign change, i.e. pair of
+        zeros, is missed). Higher up the critical strip you should use a
+        smaller stepsize so as not to miss zeros.
+
+        EXAMPLES::
+
+            sage: lcalc.zeros_in_interval(10, 30, 0.1)
+            [(14.1347251, 0.184672916), (21.0220396, -0.0677893290), (25.0108576, -0.0555872781)]
+        """
+        L = self._compute_L(L)
+        RR = RealField(prec)
+        X = self('--zeros-interval -x %s -y %s --stepsize=%s %s' % (
+            float(x), float(y), float(stepsize), L))
+        return [tuple([RR(z) for z in t.split()]) for t in X.split('\n')]
+
+    def value(self, s, L=''):
+        r"""
+        Return `L(s)` for `s` a complex number.
+
+        INPUT:
+
+        - ``s`` -- complex number
+
+        - ``L`` -- defines `L`-function (default: Riemann zeta function)
+
+        EXAMPLES::
+
+            sage: I = CC.0
+            sage: lcalc.value(0.5 + 100*I)
+            2.69261989 - 0.0203860296*I
+
+        Note, Sage can also compute zeta at complex numbers (using the PARI
+        C library)::
+
+            sage: (0.5 + 100*I).zeta()
+            2.69261988568132 - 0.0203860296025982*I
+        """
+        L = self._compute_L(L)
+        CC = ComplexField(prec)
+        s = CC(s)
+        x, y = self('-v -x %s -y %s %s' % (s.real(), s.imag(), L)).split()
+        return CC((float(x), float(y)))
+
+    def values_along_line(self, s0, s1, number_samples, L=''):
+        r"""
+        Return values of `L(s)` at ``number_samples``
+        equally-spaced sample points along the line from `s_0` to
+        `s_1` in the complex plane.
+
+        INPUT:
+
+        - ``s0, s1`` -- complex numbers
+
+        - ``number_samples`` -- integer
+
+        - ``L`` -- defines `L`-function (default: Riemann zeta function)
+
+        OUTPUT: list of pairs `(s, L(s))`, where the `s` are equally spaced
+        sampled points on the line from `s_0` to `s_1`
+
+        EXAMPLES::
+
+            sage: I = CC.0
+            sage: values = lcalc.values_along_line(0.5, 0.5+20*I, 5)
+            sage: values[0][0] # abs tol 1e-8
+            0.5
+            sage: values[0][1] # abs tol 1e-8
+            -1.46035451 + 0.0*I
+            sage: values[1][0] # abs tol 1e-8
+            0.5 + 4.0*I
+            sage: values[1][1] # abs tol 1e-8
+            0.606783764 + 0.0911121400*I
+            sage: values[2][0] # abs tol 1e-8
+            0.5 + 8.0*I
+            sage: values[2][1] # abs tol 1e-8
+            1.24161511 + 0.360047588*I
+            sage: values[3][0] # abs tol 1e-8
+            0.5 + 12.0*I
+            sage: values[3][1] # abs tol 1e-8
+            1.01593665 - 0.745112472*I
+            sage: values[4][0] # abs tol 1e-8
+            0.5 + 16.0*I
+            sage: values[4][1] # abs tol 1e-8
+            0.938545408 + 1.21658782*I
+
+        Sometimes warnings are printed (by lcalc) when this command is
+        run::
+
+            sage: E = EllipticCurve('389a')
+            sage: values = E.lseries().values_along_line(0.5, 3, 5)
+            sage: values[0][0] # abs tol 1e-8
+            0.0
+            sage: values[0][1] # abs tol 1e-8
+            0.209951303  + 0.0*I
+            sage: values[1][0] # abs tol 1e-8
+            0.5
+            sage: values[1][1] # abs tol 1e-8
+            0.0  + 0.0*I
+            sage: values[2][0] # abs tol 1e-8
+            1.0
+            sage: values[2][1] # abs tol 1e-8
+            0.133768433 - 0.0*I
+            sage: values[3][0] # abs tol 1e-8
+            1.5
+            sage: values[3][1] # abs tol 1e-8
+            0.360092864 - 0.0*I
+            sage: values[4][0] # abs tol 1e-8
+            2.0
+            sage: values[4][1] # abs tol 1e-8
+            0.552975867 + 0.0*I
+        """
+        L = self._compute_L(L)
+        CC = ComplexField(prec)
+        s0 = CC(s0)
+        s1 = CC(s1)
+        v = self('--value-line-segment -x %s -y %s -X %s -Y %s --number-samples %s %s' % (
+            (s0.real(), s0.imag(), s1.real(), s1.imag(), int(number_samples), L)))
+        w = []
+        for a in v.split('\n'):
+            try:
+                x0, y0, x1, y1 = a.split()
+                w.append((CC(x0, y0), CC(x1, y1)))
+            except ValueError:
+                print('lcalc: {}'.format(a))
+        return w
+
+    def twist_values(self, s, dmin, dmax, L=''):
+        r"""
+        Return values of `L(s, \chi_k)` for each quadratic
+        character `\chi_k` whose discriminant `d` satisfies
+        `d_{\min} \leq d \leq d_{\max}`.
+
+        INPUT:
+
+        - ``s`` -- complex numbers
+
+        - ``dmin`` -- integer
+
+        - ``dmax`` -- integer
+
+        - ``L`` -- defines `L`-function (default: Riemann zeta function)
+
+        OUTPUT: list of pairs `(d, L(s,\chi_d))`
+
+        EXAMPLES::
+
+            sage: values = lcalc.twist_values(0.5, -10, 10)
+            sage: values[0][0]
+            -8
+            sage: values[0][1] # abs tol 1e-8
+            1.10042141 + 0.0*I
+            sage: values[1][0]
+            -7
+            sage: values[1][1] # abs tol 1e-8
+            1.14658567 + 0.0*I
+            sage: values[2][0]
+            -4
+            sage: values[2][1] # abs tol 1e-8
+            0.667691457 + 0.0*I
+            sage: values[3][0]
+            -3
+            sage: values[3][1] # abs tol 1e-8
+            0.480867558 + 0.0*I
+            sage: values[4][0]
+            5
+            sage: values[4][1] # abs tol 1e-8
+            0.231750947 + 0.0*I
+            sage: values[5][0]
+            8
+            sage: values[5][1] # abs tol 1e-8
+            0.373691713 + 0.0*I
+        """
+        L = self._compute_L(L)
+        CC = ComplexField(prec)
+        s = CC(s)
+        typ = '--twist-quadratic'
+        dmin = int(dmin)
+        dmax = int(dmax)
+        v = self('-v -x %s -y %s %s --start %s --finish %s %s' % (
+            (s.real(), s.imag(), typ, dmin, dmax, L)))
+        w = []
+        if len(v) == 0:
+            return w
+        if len(v) == 0:
+            return w
+        for a in v.split('\n'):
+            d, x, y = a.split()
+            w.append((ZZ(d), CC(x, y)))
+        return w
+
+    def twist_zeros(self, n, dmin, dmax, L=''):
+        r"""
+        Return first `n` real parts of nontrivial zeros for each
+        quadratic character `\chi_k` whose discriminant `d` satisfies
+        `d_{\min} \leq d \leq d_{\max}`.
+
+        INPUT:
+
+        - ``n`` -- integer
+
+        - ``dmin`` -- integer
+
+        - ``dmax`` -- integer
+
+        - ``L`` -- defines `L`-function (default: Riemann zeta function)
+
+        OUTPUT: dictionary; keys are the discriminants `d`, and values are list
+        of corresponding zeros
+
+        EXAMPLES::
+
+            sage: lcalc.twist_zeros(3, -3, 6)
+            {-3: [8.03973716, 11.2492062, 15.7046192], 5: [6.64845335, 9.83144443, 11.9588456]}
+        """
+        L = self._compute_L(L)
+        RR = RealField(prec)
+        typ = '--twist-quadratic'
+        n = int(n)
+        v = self('-z %s %s --start %s --finish %s %s' % (
+            (n, typ, dmin, dmax, L)))
+        w = {}
+        if len(v) == 0:
+            return w
+        for a in v.split('\n'):
+            d, x = a.split()
+            x = RR(x)
+            d = ZZ(d)
+            if d in w:
+                w[d].append(x)
+            else:
+                w[d] = [x]
+        return w
+
+    def analytic_rank(self, L=''):
+        r"""
+        Return the analytic rank of the `L`-function at the central
+        critical point.
+
+        INPUT:
+
+        - ``L`` -- defines `L`-function (default: Riemann zeta function)
+
+        OUTPUT: integer
+
+        .. NOTE::
+
+           Of course this is not provably correct in general, since it
+           is an open problem to compute analytic ranks provably
+           correctly in general.
+
+        EXAMPLES::
+
+            sage: E = EllipticCurve('37a')
+            sage: lcalc.analytic_rank(E)
+            1
+        """
+        L = self._compute_L(L)
+        s = self('--rank-compute %s' % L)
+        i = s.find('equals')
+        return ZZ(s[i + 6:])
+
+
+# An instance
+lcalc = LCalc()

sage/libs/all__sagemath_lcalc.py

@@ -0,0 +1 @@
+# sage_setup: distribution = sagemath-lcalc

sage/libs/lcalc/__init__.py

@@ -0,0 +1 @@
+# sage_setup: distribution = sagemath-lcalc