passagemath-pari 10.6.32__cp314-cp314-musllinux_1_2_x86_64.whl

sage/schemes/elliptic_curves/all__sagemath_pari.py

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

sage/schemes/elliptic_curves/descent_two_isogeny_pari.pyx

@@ -0,0 +1,46 @@
+# sage_setup: distribution = sagemath-pari
+
+from sage.libs.gmp.mpz cimport mpz_t
+
+from cypari2.paridecl cimport (GEN, cgetg, t_POL, set_gel, gel, stoi, lg,
+                               evalvarn, evalsigne, Z_issquare,
+                               hyperellratpoints)
+from cypari2.stack cimport clear_stack
+from cysignals.signals cimport sig_on
+
+from sage.libs.pari.convert_gmp cimport _new_GEN_from_mpz_t
+
+
+cdef bint ratpoints_mpz_exists_only(mpz_t *coeffs, long degree, long H) except -1:
+    """
+    Search for projective points on the hyperelliptic curve
+    ``y^2 = P(x)``.
+
+    INPUT:
+
+    - ``coeffs`` -- an array of length ``degree + 1`` giving the
+      coefficients of ``P``, starting with the constant coefficient
+
+    - ``degree`` -- degree of ``P``
+
+    - ``H`` -- bound on the naive height for search
+
+    OUTPUT: boolean, whether or not a projective point was found
+    """
+    sig_on()
+    cdef GEN pol = cgetg(degree + 3, t_POL)
+    pol[1] = evalvarn(0) + evalsigne(1)
+    cdef long i
+    for i in range(degree + 1):
+        set_gel(pol, i+2, _new_GEN_from_mpz_t(coeffs[i]))
+
+    # PARI checks only for affine points, so we manually check for
+    # points at infinity (of the smooth model)
+    cdef int r
+    if degree % 2 or Z_issquare(gel(pol, degree+2)):
+        r = 1
+    else:
+        R = hyperellratpoints(pol, stoi(H), 1)
+        r = (lg(R) > 1)
+    clear_stack()
+    return r

sage_wheels/bin/gp2c-run

@@ -0,0 +1,57 @@
+#!/bin/sh -e
+#Copyright (C) 2001-2013  The PARI group.
+#
+#This file is part of the GP2C package.
+#
+#PARI/GP is free software; you can redistribute it and/or modify it under the
+#terms of the GNU General Public License as published by the Free Software
+#Foundation. It is distributed in the hope that it will be useful, but WITHOUT
+#ANY WARRANTY WHATSOEVER.
+#
+#Check the License for details. You should have received a copy of it, along
+#with the package; see the file 'COPYING'. If not, write to the Free Software
+#Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
+
+name="stdin";
+for p in "$@"; do
+case "$p" in
+  -h|--help) cat <<EOF
+gp2c-run <options> <file.gp>
+Compile file.gp with gp2c and launch a gp session with functions of file.gp
+added.
+<options> are passed to gp2c
+
+gp2c-run <file.c>
+Compile file.c and launch a gp session with functions of file.c added.
+
+ENVIRONMENT:
+GP2C: path to the gp2c compiler, default: /host/sage-musllinux_1_2_x86_64/bin/gp2c
+GP: path to the gp calculator, default: /host/sage-musllinux_1_2_x86_64/bin/gp
+EOF
+exit
+;;
+  -l|-t|-v)  exec $GP2C "$@" ;;
+  -o) echo "gp2c-run does not support option -o" 1>&2; exit 1;;
+  -*) ;; #We discard options
+  *) name=$p;;
+esac
+done
+
+if test "x$GP2C" = "x"; then
+  GP2C="/host/sage-musllinux_1_2_x86_64/bin/gp2c"
+fi
+
+
+if test "x$GP" = "x"; then
+  GP="/host/sage-musllinux_1_2_x86_64/bin/gp"
+fi
+
+case $name in
+*.c) name=${name%.c};;
+*) $GP2C -o $name.c "$@" ;;
+esac
+
+command=`echo "gcc -c -o %s.o -O3 -Wall -fno-strict-aliasing -g -O2 -fPIC -I\"/host/sage-musllinux_1_2_x86_64/include\" %s.c && gcc -o %s.so -shared -O3 -Wall -fno-strict-aliasing -g -O2 -fPIC -Wl,-shared -Wl,-rpath-link,/host/sage-musllinux_1_2_x86_64/lib -L/host/sage-musllinux_1_2_x86_64/lib -Wl,-rpath,/host/sage-musllinux_1_2_x86_64/lib -L/host/sage-musllinux_1_2_x86_64/lib -Wl,-rpath-link,/host/sage-musllinux_1_2_x86_64/lib -L/host/sage-musllinux_1_2_x86_64/lib -Wl,-rpath,/host/sage-musllinux_1_2_x86_64/lib %s.o -Wl,-rpath,/host/sage-musllinux_1_2_x86_64/lib -lc -lm -L/host/sage-musllinux_1_2_x86_64/lib -lpari" | sed -e "s %s $name g"`
+eval $command || exit $?
+grep "^GP;" $name.c | sed 's/^GP;//' >$name.run
+$GP $name.run