passagemath-categories 10.6.32__cp314-cp314t-musllinux_1_2_aarch64.whl

sage/ext/fast_eval.pxd

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

sage/ext/fast_eval.pyx

@@ -0,0 +1,102 @@
+# sage_setup: distribution = sagemath-categories
+r"""
+Fast Numerical Evaluation
+
+For many applications such as numerical integration, differential
+equation approximation, plotting a 3d surface, optimization problems,
+monte-carlo simulations, etc., one wishes to pass around and evaluate
+a single algebraic expression many, many times at various floating
+point values. Doing this via recursive calls over a python
+representation of the object (even if Maxima or other outside packages
+are not involved) is extremely inefficient.
+
+The solution implemented in this module, by Robert Bradshaw (2008-10),
+has been superseded by :func:`~sage.ext.fast_callable.fast_callable`.
+All that remains here is a compatible interface function :func:`fast_float`.
+
+AUTHORS:
+
+- Robert Bradshaw (2008-10): Initial version
+"""
+
+#*****************************************************************************
+#       Copyright (C) 2008 Robert Bradshaw <robertwb@math.washington.edu>
+#
+# This program 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, either version 2 of the License, or
+# (at your option) any later version.
+#                  http://www.gnu.org/licenses/
+#*****************************************************************************
+
+from sage.ext.fast_callable import fast_callable, Wrapper
+
+
+def fast_float(f, *vars, expect_one_var=False):
+    """
+    Try to create a function that evaluates f quickly using
+    floating-point numbers, if possible.
+
+    On failure, returns the input unchanged.
+
+    This is an alternative interface to :func:`sage.ext.fast_callable.fast_callable`.
+    :issue:`32268` proposes to deprecate this function.
+
+    INPUT:
+
+    - ``f`` -- an expression
+    - ``vars`` -- the names of the arguments
+    - ``expect_one_var`` -- don't give deprecation warning if ``vars`` is
+      omitted, as long as expression has only one var
+
+    EXAMPLES::
+
+        sage: from sage.ext.fast_eval import fast_float
+        sage: x,y = var('x,y')                                                          # needs sage.symbolic
+        sage: f = fast_float(sqrt(x^2+y^2), 'x', 'y')                                   # needs sage.symbolic
+        sage: f(3,4)                                                                    # needs sage.symbolic
+        5.0
+
+    Specifying the argument names is essential, as fast_float objects
+    only distinguish between arguments by order. ::
+
+        sage: # needs sage.symbolic
+        sage: f = fast_float(x-y, 'x','y')
+        sage: f(1,2)
+        -1.0
+        sage: f = fast_float(x-y, 'y','x')
+        sage: f(1,2)
+        1.0
+    """
+    if isinstance(f, (tuple, list)):
+        return tuple([fast_float(x, *vars, expect_one_var=expect_one_var) for x in f])
+
+    cdef int i
+    for i from 0 <= i < len(vars):
+        if not isinstance(vars[i], str):
+            v = str(vars[i])
+            # inexact generators display as 1.00..0*x
+            if '*' in v:
+                v = v[v.index('*')+1:]
+            vars = vars[:i] + (v,) + vars[i+1:]
+
+    try:
+        return fast_callable(f, vars=vars, domain=float,
+                             expect_one_var=expect_one_var)
+    except AttributeError:
+        pass
+
+    try:
+        from sage.symbolic.ring import SR
+        return fast_float(SR(f), *vars)
+    except TypeError:
+        pass
+
+    if f is None:
+        raise TypeError("no way to make fast_float from None")
+
+    return f
+
+
+def is_fast_float(x):
+    return isinstance(x, Wrapper)

sage/ext/interpreters/__init__.py

@@ -0,0 +1 @@
+# Automatically generated by /tmp/build-env-lslpcmg7/lib/python3.14t/site-packages/sage_setup/autogen/interpreters/internal/generator.py.  Do not edit!

sage/ext/interpreters/all__sagemath_categories.py

@@ -0,0 +1,2 @@
+# sage_setup: distribution = sagemath-categories
+# Automatically generated by /tmp/build-env-lslpcmg7/lib/python3.14t/site-packages/sage_setup/autogen/interpreters/internal/generator.py.  Do not edit!

sage/ext/interpreters/wrapper_el.pxd

@@ -0,0 +1,18 @@
+# sage_setup: distribution = sagemath-categories
+# Automatically generated by /tmp/build-env-lslpcmg7/lib/python3.14t/site-packages/sage_setup/autogen/interpreters/internal/generator.py.  Do not edit!
+
+from cpython.ref cimport PyObject
+
+from sage.ext.fast_callable cimport Wrapper
+
+cdef class Wrapper_el(Wrapper):
+    cdef int _n_args
+    cdef object _list_constants
+    cdef int _n_constants
+    cdef PyObject** _constants
+    cdef object _list_stack
+    cdef int _n_stack
+    cdef PyObject** _stack
+    cdef object _domain
+    cdef int _n_code
+    cdef int* _code

sage/ext/interpreters/wrapper_el.pyx

@@ -0,0 +1,148 @@
+# sage_setup: distribution = sagemath-categories
+# Automatically generated by /tmp/build-env-lslpcmg7/lib/python3.14t/site-packages/sage_setup/autogen/interpreters/internal/generator.py.  Do not edit!
+from cpython.number cimport PyNumber_TrueDivide
+from sage.structure.element cimport Element
+
+cdef public object el_check_element(object v, parent):
+    cdef Element v_el
+
+    if isinstance(v, Element):
+        v_el = <Element>v
+        if v_el._parent is parent:
+            return v_el
+
+    return parent(v)
+
+
+from cpython.ref cimport PyObject
+cdef extern from "Python.h":
+    void Py_DECREF(PyObject *o)
+    void Py_INCREF(PyObject *o)
+    void Py_CLEAR(PyObject *o)
+
+    object PyList_New(Py_ssize_t len)
+    ctypedef struct PyListObject:
+        PyObject **ob_item
+
+    ctypedef struct PyTupleObject:
+        PyObject **ob_item
+
+from cysignals.memory cimport check_allocarray, sig_free
+
+from sage.ext.fast_callable cimport Wrapper
+
+cdef extern from "interp_el.c":
+    object interp_el(PyObject** args,
+        PyObject** constants,
+        PyObject** stack,
+        PyObject* domain,
+        int* code)
+cdef class Wrapper_el(Wrapper):
+    # attributes are declared in corresponding .pxd file
+
+    def __init__(self, args):
+        Wrapper.__init__(self, args, metadata)
+        cdef int i
+        cdef int count
+        count = args['args']
+        self._n_args = count
+        val = args['constants']
+        self._n_constants = len(val)
+        self._list_constants = PyList_New(self._n_constants)
+        self._constants = (<PyListObject *>self._list_constants).ob_item
+        for i in range(len(val)):
+            self._constants[i] = <PyObject *>val[i]; Py_INCREF(self._constants[i])
+        count = args['stack']
+        self._n_stack = count
+        self._list_stack = PyList_New(self._n_stack)
+        self._stack = (<PyListObject *>self._list_stack).ob_item
+        self._domain = args['domain']
+        val = args['code']
+        self._n_code = len(val)
+        self._code = <int*>check_allocarray(self._n_code, sizeof(int))
+        for i in range(len(val)):
+            self._code[i] = val[i]
+
+    def __dealloc__(self):
+        cdef int i
+        if self._code:
+            sig_free(self._code)
+
+    def __call__(self, *args):
+        if self._n_args != len(args): raise ValueError
+        mapped_args = [self._domain(a) for a in args]
+        try:
+            return interp_el((<PyListObject*>mapped_args).ob_item
+                , self._constants
+                , self._stack
+                , <PyObject*>self._domain
+                , self._code
+                )
+        except BaseException:
+            for i in range(self._n_stack):
+                Py_CLEAR(self._stack[i])
+            raise
+
+
+from sage.ext.fast_callable import CompilerInstrSpec, InterpreterMetadata
+metadata = InterpreterMetadata(by_opname={
+  'load_arg':
+  (CompilerInstrSpec(0, 1, ['args']), 0),
+  'load_const':
+  (CompilerInstrSpec(0, 1, ['constants']), 1),
+  'return':
+  (CompilerInstrSpec(1, 0, []), 2),
+  'py_call':
+  (CompilerInstrSpec(0, 1, ['constants', 'n_inputs']), 3),
+  'add':
+  (CompilerInstrSpec(2, 1, []), 4),
+  'sub':
+  (CompilerInstrSpec(2, 1, []), 5),
+  'mul':
+  (CompilerInstrSpec(2, 1, []), 6),
+  'div':
+  (CompilerInstrSpec(2, 1, []), 7),
+  'floordiv':
+  (CompilerInstrSpec(2, 1, []), 8),
+  'pow':
+  (CompilerInstrSpec(2, 1, []), 9),
+  'ipow':
+  (CompilerInstrSpec(1, 1, ['constants']), 10),
+  'neg':
+  (CompilerInstrSpec(1, 1, []), 11),
+  'invert':
+  (CompilerInstrSpec(1, 1, []), 12),
+  'abs':
+  (CompilerInstrSpec(1, 1, []), 13),
+ },
+ by_opcode=[
+  ('load_arg',
+   CompilerInstrSpec(0, 1, ['args'])),
+  ('load_const',
+   CompilerInstrSpec(0, 1, ['constants'])),
+  ('return',
+   CompilerInstrSpec(1, 0, [])),
+  ('py_call',
+   CompilerInstrSpec(0, 1, ['constants', 'n_inputs'])),
+  ('add',
+   CompilerInstrSpec(2, 1, [])),
+  ('sub',
+   CompilerInstrSpec(2, 1, [])),
+  ('mul',
+   CompilerInstrSpec(2, 1, [])),
+  ('div',
+   CompilerInstrSpec(2, 1, [])),
+  ('floordiv',
+   CompilerInstrSpec(2, 1, [])),
+  ('pow',
+   CompilerInstrSpec(2, 1, [])),
+  ('ipow',
+   CompilerInstrSpec(1, 1, ['constants'])),
+  ('neg',
+   CompilerInstrSpec(1, 1, [])),
+  ('invert',
+   CompilerInstrSpec(1, 1, [])),
+  ('abs',
+   CompilerInstrSpec(1, 1, [])),
+ ],
+ ipow_range=True)

sage/ext/interpreters/wrapper_py.pxd

@@ -0,0 +1,17 @@
+# sage_setup: distribution = sagemath-categories
+# Automatically generated by /tmp/build-env-lslpcmg7/lib/python3.14t/site-packages/sage_setup/autogen/interpreters/internal/generator.py.  Do not edit!
+
+from cpython.ref cimport PyObject
+
+from sage.ext.fast_callable cimport Wrapper
+
+cdef class Wrapper_py(Wrapper):
+    cdef int _n_args
+    cdef object _list_constants
+    cdef int _n_constants
+    cdef PyObject** _constants
+    cdef object _list_stack
+    cdef int _n_stack
+    cdef PyObject** _stack
+    cdef int _n_code
+    cdef int* _code

sage/ext/interpreters/wrapper_py.pyx

@@ -0,0 +1,133 @@
+# sage_setup: distribution = sagemath-categories
+# Automatically generated by /tmp/build-env-lslpcmg7/lib/python3.14t/site-packages/sage_setup/autogen/interpreters/internal/generator.py.  Do not edit!
+from cpython.number cimport PyNumber_TrueDivide
+
+
+from cpython.ref cimport PyObject
+cdef extern from "Python.h":
+    void Py_DECREF(PyObject *o)
+    void Py_INCREF(PyObject *o)
+    void Py_CLEAR(PyObject *o)
+
+    object PyList_New(Py_ssize_t len)
+    ctypedef struct PyListObject:
+        PyObject **ob_item
+
+    ctypedef struct PyTupleObject:
+        PyObject **ob_item
+
+from cysignals.memory cimport check_allocarray, sig_free
+
+from sage.ext.fast_callable cimport Wrapper
+
+cdef extern from "interp_py.c":
+    object interp_py(PyObject** args,
+        PyObject** constants,
+        PyObject** stack,
+        int* code)
+cdef class Wrapper_py(Wrapper):
+    # attributes are declared in corresponding .pxd file
+
+    def __init__(self, args):
+        Wrapper.__init__(self, args, metadata)
+        cdef int i
+        cdef int count
+        count = args['args']
+        self._n_args = count
+        val = args['constants']
+        self._n_constants = len(val)
+        self._list_constants = PyList_New(self._n_constants)
+        self._constants = (<PyListObject *>self._list_constants).ob_item
+        for i in range(len(val)):
+            self._constants[i] = <PyObject *>val[i]; Py_INCREF(self._constants[i])
+        count = args['stack']
+        self._n_stack = count
+        self._list_stack = PyList_New(self._n_stack)
+        self._stack = (<PyListObject *>self._list_stack).ob_item
+        val = args['code']
+        self._n_code = len(val)
+        self._code = <int*>check_allocarray(self._n_code, sizeof(int))
+        for i in range(len(val)):
+            self._code[i] = val[i]
+
+    def __dealloc__(self):
+        cdef int i
+        if self._code:
+            sig_free(self._code)
+
+    def __call__(self, *args):
+        if self._n_args != len(args): raise ValueError
+        try:
+            return interp_py((<PyTupleObject*>args).ob_item
+                , self._constants
+                , self._stack
+                , self._code
+                )
+        except BaseException:
+            for i in range(self._n_stack):
+                Py_CLEAR(self._stack[i])
+            raise
+
+
+from sage.ext.fast_callable import CompilerInstrSpec, InterpreterMetadata
+metadata = InterpreterMetadata(by_opname={
+  'load_arg':
+  (CompilerInstrSpec(0, 1, ['args']), 0),
+  'load_const':
+  (CompilerInstrSpec(0, 1, ['constants']), 1),
+  'return':
+  (CompilerInstrSpec(1, 0, []), 2),
+  'py_call':
+  (CompilerInstrSpec(0, 1, ['constants', 'n_inputs']), 3),
+  'add':
+  (CompilerInstrSpec(2, 1, []), 4),
+  'sub':
+  (CompilerInstrSpec(2, 1, []), 5),
+  'mul':
+  (CompilerInstrSpec(2, 1, []), 6),
+  'div':
+  (CompilerInstrSpec(2, 1, []), 7),
+  'floordiv':
+  (CompilerInstrSpec(2, 1, []), 8),
+  'pow':
+  (CompilerInstrSpec(2, 1, []), 9),
+  'ipow':
+  (CompilerInstrSpec(1, 1, ['constants']), 10),
+  'neg':
+  (CompilerInstrSpec(1, 1, []), 11),
+  'invert':
+  (CompilerInstrSpec(1, 1, []), 12),
+  'abs':
+  (CompilerInstrSpec(1, 1, []), 13),
+ },
+ by_opcode=[
+  ('load_arg',
+   CompilerInstrSpec(0, 1, ['args'])),
+  ('load_const',
+   CompilerInstrSpec(0, 1, ['constants'])),
+  ('return',
+   CompilerInstrSpec(1, 0, [])),
+  ('py_call',
+   CompilerInstrSpec(0, 1, ['constants', 'n_inputs'])),
+  ('add',
+   CompilerInstrSpec(2, 1, [])),
+  ('sub',
+   CompilerInstrSpec(2, 1, [])),
+  ('mul',
+   CompilerInstrSpec(2, 1, [])),
+  ('div',
+   CompilerInstrSpec(2, 1, [])),
+  ('floordiv',
+   CompilerInstrSpec(2, 1, [])),
+  ('pow',
+   CompilerInstrSpec(2, 1, [])),
+  ('ipow',
+   CompilerInstrSpec(1, 1, ['constants'])),
+  ('neg',
+   CompilerInstrSpec(1, 1, [])),
+  ('invert',
+   CompilerInstrSpec(1, 1, [])),
+  ('abs',
+   CompilerInstrSpec(1, 1, [])),
+ ],
+ ipow_range=True)