passagemath-standard-no-symbolics 10.6.31rc3__cp314-cp314-macosx_13_0_arm64.whl

sage/interfaces/read_data.py

@@ -0,0 +1,46 @@
+"""
+An interface to read data files
+"""
+###############################################################################
+#   Sage: Open Source Mathematical Software
+#       Copyright (C) 2010 Paul Zimmermann
+#  Distributed under the terms of the GNU General Public License (GPL),
+#  version 2 or later. The full text of the GPL is available at:
+#                  https://www.gnu.org/licenses/
+###############################################################################
+
+
+def read_data(f, t):
+    r"""
+    Read data from file ``f`` and class ``t`` (one element per line),
+    and returns a list of elements.
+
+    INPUT:
+
+    - ``f`` -- a file name
+    - ``t`` -- a class (objects will be coerced to that class)
+
+    OUTPUT: list of elements of class ``t``
+
+    EXAMPLES::
+
+        sage: indata = tmp_filename()
+        sage: f = open(indata, "w")
+        sage: _ = f.write("17\n42\n")
+        sage: f.close()
+        sage: l = read_data(indata, ZZ); l
+        [17, 42]
+        sage: f = open(indata, "w")
+        sage: _ = f.write("1.234\n5.678\n")
+        sage: f.close()
+        sage: l = read_data(indata, RealField(17)); l
+        [1.234, 5.678]
+    """
+    with open(f) as fp:
+        l = []
+        while True:
+            s = fp.readline().strip()
+            if not s:
+                break
+            l.append(t(s))
+    return l

sage/interfaces/scilab.py

@@ -0,0 +1,576 @@
+r"""
+Interface to Scilab
+
+Scilab is a scientific software package for numerical computations
+providing a powerful open computing environment for engineering and
+scientific applications.  Scilab includes hundreds of mathematical
+functions with the possibility to add interactively programs from
+various languages (C, C++, Fortran...).  It has sophisticated data
+structures (including lists, polynomials, rational functions, linear
+systems...), an interpreter and a high level programming language.
+
+The commands in this section only work if you have the "scilab"
+interpreter installed and available in your PATH.  It's not necessary
+to install any special Sage packages.
+
+EXAMPLES::
+
+    sage: # optional - scilab
+    sage: scilab.eval('2+2')
+    'ans  =\n \n    4.'
+    sage: scilab('2+2')
+    4.
+    sage: a = scilab(10)
+    sage: a**10
+    1.000D+10
+
+Tutorial based the MATLAB interface tutorial:
+
+EXAMPLES::
+
+    sage: # optional - scilab
+    sage: scilab('4+10')
+    14.
+    sage: scilab('date')
+    15-Feb-2010
+    sage: scilab('5*10 + 6')
+    56.
+    sage: scilab('(6+6)/3')
+    4.
+    sage: scilab('9')^2
+    81.
+    sage: a = scilab(10); b = scilab(20); c = scilab(30)
+    sage: avg = (a+b+c)/3
+    sage: avg
+    20.
+    sage: parent(avg)
+    Scilab
+
+    sage: # optional - scilab
+    sage: my_scalar = scilab('3.1415')
+    sage: my_scalar
+    3.1415
+    sage: my_vector1 = scilab('[1,5,7]')
+    sage: my_vector1
+    1.    5.    7.
+    sage: my_vector2 = scilab('[1;5;7]')
+    sage: my_vector2
+    1.
+    5.
+    7.
+    sage: my_vector1 * my_vector2
+    75.
+
+    sage: # optional - scilab
+    sage: row_vector1 = scilab('[1 2 3]')
+    sage: row_vector2 = scilab('[3 2 1]')
+    sage: matrix_from_row_vec = scilab('[%s; %s]'%(row_vector1.name(), row_vector2.name()))
+    sage: matrix_from_row_vec
+    1.    2.    3.
+    3.    2.    1.
+
+    sage: # optional - scilab
+    sage: column_vector1 = scilab('[1;3]')
+    sage: column_vector2 = scilab('[2;8]')
+    sage: matrix_from_col_vec = scilab('[%s %s]'%(column_vector1.name(), column_vector2.name()))
+    sage: matrix_from_col_vec
+    1.    2.
+    3.    8.
+
+    sage: my_matrix = scilab('[8, 12, 19; 7, 3, 2; 12, 4, 23; 8, 1, 1]')    # optional - scilab
+    sage: my_matrix                                      # optional - scilab
+        8.     12.    19.
+        7.     3.     2.
+        12.    4.     23.
+        8.     1.     1.
+
+    sage: combined_matrix = scilab('[%s, %s]'%(my_matrix.name(), my_matrix.name()))                                        # optional - scilab
+    sage: combined_matrix                               # optional - scilab
+    8.     12.    19.    8.     12.    19.
+    7.     3.     2.     7.     3.     2.
+    12.    4.     23.    12.    4.     23.
+    8.     1.     1.     8.     1.     1.
+
+    sage: tm = scilab('0.5:2:10')                       # optional - scilab
+    sage: tm                                            # optional - scilab
+    0.5    2.5    4.5    6.5    8.5
+
+    sage: # optional - scilab
+    sage: my_vector1 = scilab('[1,5,7]')
+    sage: my_vector1(1)
+    1.
+    sage: my_vector1(2)
+    5.
+    sage: my_vector1(3)
+    7.
+
+Matrix indexing works as follows::
+
+    sage: my_matrix = scilab('[8, 12, 19; 7, 3, 2; 12, 4, 23; 8, 1, 1]')     # optional - scilab
+    sage: my_matrix(3,2)                                # optional - scilab
+    4.
+
+One can also use square brackets::
+
+    sage: my_matrix[3,2]                                # optional - scilab
+    4.
+
+
+Setting using parenthesis cannot work (because of how the Python
+language works).  Use square brackets or the set function::
+
+    sage: # optional - scilab
+    sage: my_matrix = scilab('[8, 12, 19; 7, 3, 2; 12, 4, 23; 8, 1, 1]')
+    sage: my_matrix.set(2,3, 1999)
+    sage: my_matrix
+           8.         12.         19.
+           7.          3.       1999.
+          12.          4.         23.
+           8.          1.          1.
+    sage: my_matrix[2,3] = -126
+    sage: my_matrix
+           8.         12.         19.
+           7.          3.      - 126.
+          12.          4.         23.
+           8.          1.          1.
+
+TESTS::
+
+    sage: # optional - scilab
+    sage: M = scilab(x)
+    Traceback (most recent call last):
+    ...
+    TypeError: ..._interface_init_() takes exactly one argument (0 given)
+    sage: M = scilab(matrix(3,range(9))); M
+        0.    1.    2.
+        3.    4.    5.
+        6.    7.    8.
+    sage: M(10)
+    Traceback (most recent call last):
+    ...
+    TypeError: Error executing code in Scilab
+    ...
+    Invalid index.
+    sage: M[10]
+    Traceback (most recent call last):
+    ...
+    TypeError: Error executing code in Scilab
+    ...
+    Invalid index.
+    sage: M(4,2)
+    Traceback (most recent call last):
+    ...
+    TypeError: Error executing code in Scilab
+    ...
+    Invalid index.
+    sage: M[2,4]
+    Traceback (most recent call last):
+    ...
+    TypeError: Error executing code in Scilab
+    ...
+    Invalid index.
+    sage: M(9) = x
+    Traceback (most recent call last):
+    ...
+    SyntaxError: can...t assign to function call (..., line 1)
+
+AUTHORS:
+
+- Ronan Paixao (2008-11-26), based on the MATLAB tutorial by
+  William Stein (2006-10-11)
+"""
+# ****************************************************************************
+#       Copyright (C) 2006 William Stein <wstein@gmail.com>
+#       Copyright (C) 2008 Ronan Paixao <ronanpaixao@yahoo.com.br>
+#
+#  Distributed under the terms of the GNU General Public License (GPL).
+#
+#    This code is distributed in the hope that it will be useful,
+#    but WITHOUT ANY WARRANTY; without even the implied warranty of
+#    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+#    General Public License for more details.
+#
+#  The full text of the GPL is available at:
+#
+#                  https://www.gnu.org/licenses/
+# ****************************************************************************
+
+import os
+
+from .expect import Expect, ExpectElement
+from sage.misc.instancedoc import instancedoc
+
+
+class Scilab(Expect):
+    """
+    Interface to the Scilab interpreter.
+
+    EXAMPLES::
+
+        sage: # optional - scilab
+        sage: a = scilab('[ 1, 1, 2; 3, 5, 8; 13, 21, 33 ]')
+        sage: b = scilab('[ 1; 3; 13]')
+        sage: c = a * b
+        sage: print(c)
+          30.
+          122.
+          505.
+    """
+    def __init__(self, maxread=None, script_subdirectory=None,
+                 logfile=None, server=None, server_tmpdir=None,
+                 seed=None):
+        """
+        Initialize the Scilab class.
+
+        EXAMPLES::
+
+            sage: from sage.interfaces.scilab import Scilab
+            sage: sci_obj = Scilab()
+            sage: del sci_obj
+        """
+        Expect.__init__(self,
+                        name='scilab',
+                        prompt='-->',
+                        command="scilab -nw",
+                        server=server,
+                        server_tmpdir=server_tmpdir,
+                        script_subdirectory=script_subdirectory,
+                        restart_on_ctrlc=False,
+                        verbose_start=False,
+                        logfile=logfile,
+                        eval_using_file_cutoff=100)
+        self._seed = seed
+
+    def set_seed(self, seed=None):
+        """
+        Set the seed for gp interpreter.
+
+        The seed should be an integer.
+
+        EXAMPLES::
+
+            sage: # optional - scilab
+            sage: from sage.interfaces.scilab import Scilab
+            sage: s = Scilab()
+            sage: s.set_seed(1)
+            1
+            sage: [s.rand() for i in range(5)]
+            [
+            <BLANKLINE>
+                 0.6040239,
+            <BLANKLINE>
+                 0.0079647,
+            <BLANKLINE>
+                 0.6643966,
+            <BLANKLINE>
+                 0.9832111,
+            <BLANKLINE>
+                 0.5321420]
+        """
+        if seed is None:
+            seed = self.rand_seed()
+        self.eval("rand('seed',%d)" % seed)
+        self._seed = seed
+        return seed
+
+    def _quit_string(self):
+        """
+        Return the string used to quit the pexpect interface.
+
+        EXAMPLES::
+
+            sage: scilab._quit_string()                 # optional - scilab
+            'quit;'
+        """
+        return 'quit;'
+
+    def _install_hints(self):
+        """
+        Hints for installing Scilab.
+
+        EXAMPLES::
+
+            sage: print(scilab._install_hints())       # optional - scilab
+            You must ...
+        """
+        return """
+        You must obtain the Scilab program in order to use Scilab
+        from Sage.   You can read all about Scilab at
+                  http://www.scilab.org/
+        The executable must be accessible system-wide.
+        """
+
+    def _start(self):
+        """
+        Starts Scilab and sets some options.
+
+        EXAMPLES::
+
+            sage: scilab._start()                       # optional - scilab
+        """
+        Expect._start(self)
+        self.eval("mode(0)")
+
+        # set random seed
+        self.set_seed(self._seed)
+
+    def eval(self, command, *args, **kwds):
+        """
+        Evaluates commands.
+
+        EXAMPLES::
+
+            sage: scilab.eval("5")                      # optional - scilab
+            'ans  =\n \n    5.'
+            sage: scilab.eval("d=44")                   # optional - scilab
+            'd  =\n \n    44.'
+        """
+        s = Expect.eval(self, command, **kwds)
+        return s.replace("\x1b[?1l\x1b>", "").strip()
+
+    def whos(self, name=None, typ=None):
+        """
+        Return information about current objects.
+        Arguments:
+        nam: first characters of selected names
+        typ: name of selected Scilab variable type
+
+        EXAMPLES::
+
+            sage: scilab.whos("core")                   # optional - scilab
+            'Name                     Type           Size           Bytes...'
+            sage: scilab.whos(typ='function')           # optional - scilab
+            'Name                     Type           Size           Bytes...'
+        """
+        parameters = ""
+        if name:
+            parameters += " -name %s" % (str(name))
+        if typ:
+            parameters += " -type %s" % (str(typ))
+        return self.eval('whos' + parameters)
+
+    def set(self, var, value):
+        """
+        Set the variable var to the given value.
+
+        EXAMPLES::
+
+            sage: scilab.set('a', 123)        # optional - scilab
+            sage: scilab.get('a')               # optional - scilab
+            '\n \n    123.'
+        """
+        cmd = '%s=%s;' % (var, value)
+        out = self.eval(cmd)
+        if out.find("error") != -1:
+            raise TypeError("Error executing code in Scilab\nCODE:\n\t%s\nScilab ERROR:\n\t%s" % (cmd, out))
+
+    def get(self, var):
+        """
+        Get the value of the variable ``var``.
+
+        EXAMPLES::
+
+            sage: scilab.eval('b=124;')                 # optional - scilab
+            ''
+            sage: scilab.get('b')                       # optional - scilab
+            '\n \n    124.'
+        """
+        s = self.eval(f'{var}')
+        i = s.find('=')
+        return s[i+1:]
+
+    def console(self):
+        """
+        Starts Scilab console.
+
+        EXAMPLES::
+
+            sage: scilab.console()          # optional - scilab; not tested
+        """
+        scilab_console()
+
+    def version(self):
+        """
+        Return the version of the Scilab software used.
+
+        EXAMPLES::
+
+            sage: scilab.version()                      # optional - scilab
+            'scilab-...'
+        """
+        return scilab_version()
+
+    def sage2scilab_matrix_string(self, A):
+        """
+        Return a Scilab matrix from a Sage matrix.
+
+        INPUT:
+
+        - ``A`` -- Sage matrix with entries in the rationals or reals
+
+        OUTPUT: string that evaluates to a Scilab matrix
+
+        EXAMPLES::
+
+            sage: M33 = MatrixSpace(QQ,3,3)             # optional - scilab
+            sage: A = M33([1,2,3,4,5,6,7,8,0])          # optional - scilab
+            sage: scilab.sage2scilab_matrix_string(A)   # optional - scilab
+            '[1, 2, 3; 4, 5, 6; 7, 8, 0]'
+        """
+        s = str(A.rows())
+        return s.replace('), (', '; ').replace('(', '').replace(')', '')
+
+    def _object_class(self):
+        """
+        Return the class of the object.
+
+        EXAMPLES::
+
+            sage: scilab._object_class()                # optional - scilab
+            <class 'sage.interfaces.scilab.ScilabElement'>
+            sage: type(scilab(2))                       # optional - scilab
+            <class 'sage.interfaces.scilab.ScilabElement'>
+        """
+        return ScilabElement
+
+
+@instancedoc
+class ScilabElement(ExpectElement):
+    def __getitem__(self, n):
+        """
+        Use parenthesis for Scilab matrices instead.
+
+        EXAMPLES::
+
+            sage: # optional - scilab
+            sage: M = scilab('[1,2,3;4,5,6;7,8,9]')
+            sage: M[1]
+            1.
+            sage: M[7]
+            3.
+            sage: M[3,2]
+            8.
+        """
+        if isinstance(n, tuple):
+            index = str(n)[1:-1]
+        else:
+            index = str(n)
+        return self.parent()('%s(%s)' % (self._name, index))
+
+    def __setitem__(self, n, value):
+        """
+        Set an element of a matrix.
+
+        EXAMPLES::
+
+            sage: # optional - scilab
+            sage: M = scilab('[1,2,3;4,5,6;7,8,9]')
+            sage: M[6] = 0
+            sage: M
+            1.    2.    3.
+            4.    5.    6.
+            7.    0.    9.
+            sage: M[3,2] = 10
+            sage: M
+            1.    2.    3.
+            4.    5.    6.
+            7.    10.    9.
+        """
+        if isinstance(n, tuple):
+            index = str(n)[1:-1]
+        else:
+            index = str(n)
+        self.parent().eval('%s(%s) = %s' % (self._name, index, value))
+
+    def _matrix_(self, R):
+        r"""
+        Return \sage matrix from this scilab element.
+
+        EXAMPLES::
+
+            sage: # optional - scilab
+            sage: A = scilab('[1,2;3,4]')
+            sage: matrix(ZZ, A)
+            [1 2]
+            [3 4]
+            sage: A = scilab('[1,2;3,4.5]')
+            sage: matrix(RR, A)
+            [1.00000000000000 2.00000000000000]
+            [3.00000000000000 4.50000000000000]
+        """
+        from sage.matrix.matrix_space import MatrixSpace
+        s = str(self).strip()
+        v = s.split('\n ')
+        nrows = len(v)
+        if nrows == 0:
+            return MatrixSpace(R, 0, 0)(0)
+        ncols = len(v[0].split())
+        M = MatrixSpace(R, nrows, ncols)
+        v = sum([[x.rstrip('.') for x in w.split()] for w in v], [])
+        return M(v)
+
+    def set(self, i, j, x):
+        """
+        Set the variable var to the given value.
+
+        EXAMPLES::
+
+            sage: scilab.set('c', 125)          # optional - scilab
+            sage: scilab.get('c')               # optional - scilab
+            '\n \n    125.'
+        """
+        P = self._check_valid()
+        z = P(x)
+        P.eval('%s(%s,%s) = %s' % (self.name(), i, j, z.name()))
+
+
+# An instance
+scilab = Scilab()
+
+
+def scilab_console():
+    """
+    This requires that the optional Scilab program be installed and in
+    your PATH, but no optional Sage packages need to be installed.
+
+    EXAMPLES::
+
+        sage: from sage.interfaces.scilab import scilab_console  # optional - scilab
+        sage: scilab_console()                               # optional - scilab; not tested
+                ___________________________________________
+                               scilab-5.0.3
+
+                         Consortium Scilab (DIGITEO)
+                       Copyright (c) 1989-2008 (INRIA)
+                       Copyright (c) 1989-2007 (ENPC)
+                ___________________________________________
+
+
+        Startup execution:
+          loading initial environment
+
+        -->2+3
+        ans  =
+
+           5.
+
+        -->quit
+
+    Typing quit exits the Scilab console and returns you to Sage.
+    Scilab, like Sage, remembers its history from one session to
+    another.
+    """
+    os.system('scilab -nw')
+
+
+def scilab_version():
+    """
+    Return the version of Scilab installed.
+
+    EXAMPLES::
+
+        sage: from sage.interfaces.scilab import scilab_version  # optional - scilab
+        sage: scilab_version()    # optional - scilab
+        'scilab-...'
+    """
+    return str(scilab('getversion()')).strip()

sage/interfaces/tests.py

@@ -0,0 +1,81 @@
+"""
+TESTS:
+
+We test coercions::
+
+    sage: 2 * gp('2')
+    4
+    sage: a = 2 * gp('2'); a
+    4
+    sage: parent(a)
+    PARI/GP interpreter
+    sage: a = 2 * gap('2'); a
+    4
+    sage: parent(a)
+    Gap
+    sage: a = 2 * maxima('2'); a                                                        # needs sage.symbolic
+    4
+    sage: parent(a)                                                                     # needs sage.symbolic
+    Maxima
+    sage: a = 2 * singular('2'); a
+    4
+    sage: parent(a)
+    Singular
+
+Test that write errors to stderr are handled gracefully by GAP
+(see :issue:`13211`) and other interfaces::
+
+    sage: import subprocess
+    sage: try:
+    ....:     f = open('/dev/full', 'w')
+    ....: except IOError:
+    ....:     f = open('/dev/null', 'w')
+    sage: kwds = dict(shell=True, stdout=f, stderr=f)
+    sage: subprocess.call("echo syntax error | gap", **kwds) in (0, 1)
+    True
+    sage: subprocess.call("echo syntax error | gp", **kwds)
+    0
+    sage: subprocess.call("echo syntax error | ipython", **kwds) in (0, 1, 120)
+    True
+    sage: subprocess.call("echo syntax error | Singular", **kwds)
+    0
+    sage: f.close()
+"""
+import sys
+
+from .all import *
+from sage.misc.timing import cputime, walltime
+
+
+def manyvars(s, num=70000, inlen=1, step=2000):
+    """
+    Test that > 65,000 variable names works in each system.
+    """
+    print(f"Testing -- {s}")
+    t = '"%s"' % ('9' * int(inlen))
+    try:
+        t = cputime()
+        w = walltime()
+        v = []
+        for i in range(num):
+            if not i % step:
+                sys.stdout.write(f'{i} ')
+                sys.stdout.flush()
+            v.append(s(t))
+        print('\nsuccess -- time = cpu: %s, wall: %s' % (cputime(t),
+                                                         walltime(w)))
+    except Exception:
+        print("%s -- failed!" % s)
+
+
+def manyvars_all(num=70000):
+    # for s in [gap, gp, singular, kash, magma, octave, maxima, mathematica]:
+    for s in [kash, magma, octave, maxima, mathematica]:
+        manyvars(s, num)
+
+# bad: maple -- infinite loop -- exception pexpect.EOF: <pexpect.EOF instance at 0xb091250c> in  ignored
+
+
+def manyvars_all2(num=70000):
+    for s in [singular, maxima, mathematica, octave]:
+        manyvars(s, num)

sage/libs/all.py

@@ -0,0 +1,11 @@
+from sage.libs.all__sagemath_combinat import *
+from sage.libs.all__sagemath_gap import *
+from sage.libs.all__sagemath_flint import *
+from sage.libs.all__sagemath_ntl import *
+from sage.libs.all__sagemath_pari import *
+from sage.libs.all__sagemath_eclib import *
+
+try:
+    from sage.libs.all__sagemath_symbolics import *
+except ImportError:
+    pass