Remove custom code for Col/Colrev/Vec/Vecrev/Vecsmall

Author: jdemeyer

cypari2/gen.pyx

@@ -2294,112 +2294,6 @@ cdef class Gen(Gen_base):
         sig_on()
         return new_gen(vecmin(x.g))
 
-    def Col(x, long n = 0):
-        """
-        Transform the object `x` into a column vector with minimal size `|n|`.
-
-        INPUT:
-
-        - ``x`` -- gen
-
-        - ``n`` -- Make the column vector of minimal length `|n|`. If `n > 0`,
-          append zeros; if `n < 0`, prepend zeros.
-
-        OUTPUT:
-
-        A PARI column vector (type ``t_COL``)
-
-        Examples:
-
-        >>> from cypari2 import Pari
-        >>> pari = Pari()
-
-        >>> pari(1.5).Col()
-        [1.50000000000000]~
-        >>> pari([1,2,3,4]).Col()
-        [1, 2, 3, 4]~
-        >>> pari('[1,2; 3,4]').Col()
-        [[1, 2], [3, 4]]~
-        >>> pari('"CyPari"').Col()
-        ["C", "y", "P", "a", "r", "i"]~
-        >>> pari('x + 3*x^3').Col()
-        [3, 0, 1, 0]~
-        >>> pari('x + 3*x^3 + O(x^5)').Col()
-        [1, 0, 3, 0]~
-
-        We demonstate the `n` argument:
-
-        >>> pari([1,2,3,4]).Col(2)
-        [1, 2, 3, 4]~
-        >>> pari([1,2,3,4]).Col(-2)
-        [1, 2, 3, 4]~
-        >>> pari([1,2,3,4]).Col(6)
-        [1, 2, 3, 4, 0, 0]~
-        >>> pari([1,2,3,4]).Col(-6)
-        [0, 0, 1, 2, 3, 4]~
-
-        See also :meth:`Vec` (create a row vector) for more examples
-        and :meth:`Colrev` (create a column in reversed order).
-        """
-        sig_on()
-        return new_gen(_Vec_append(gtocol(x.g), gen_0, n))
-
-    def Colrev(x, long n = 0):
-        """
-        Transform the object `x` into a column vector with minimal size `|n|`.
-        The order of the resulting vector is reversed compared to :meth:`Col`.
-
-        INPUT:
-
-        - ``x`` -- gen
-
-        - ``n`` -- Make the vector of minimal length `|n|`. If `n > 0`,
-          prepend zeros; if `n < 0`, append zeros.
-
-        OUTPUT:
-
-        A PARI column vector (type ``t_COL``)
-
-        Examples:
-
-        >>> from cypari2 import Pari
-        >>> pari = Pari()
-
-        >>> pari(1.5).Colrev()
-        [1.50000000000000]~
-        >>> pari([1,2,3,4]).Colrev()
-        [4, 3, 2, 1]~
-        >>> pari('[1,2; 3,4]').Colrev()
-        [[3, 4], [1, 2]]~
-        >>> pari('x + 3*x^3').Colrev()
-        [0, 1, 0, 3]~
-
-        We demonstate the `n` argument:
-
-        >>> pari([1,2,3,4]).Colrev(2)
-        [4, 3, 2, 1]~
-        >>> pari([1,2,3,4]).Colrev(-2)
-        [4, 3, 2, 1]~
-        >>> pari([1,2,3,4]).Colrev(6)
-        [0, 0, 4, 3, 2, 1]~
-        >>> pari([1,2,3,4]).Colrev(-6)
-        [4, 3, 2, 1, 0, 0]~
-        """
-        sig_on()
-        # Create a non-reversed column vector
-        cdef GEN v = _Vec_append(gtocol(x.g), gen_0, n)
-        # Reverse it in-place
-        cdef GEN L = v + 1
-        cdef GEN R = v + (lg(v)-1)
-        cdef long t
-        while (L < R):
-            t = L[0]
-            L[0] = R[0]
-            R[0] = t
-            L += 1
-            R -= 1
-        return new_gen(v)
-
     def Ser(f, v=None, long precision=-1):
         """
         Return a power series or Laurent series in the variable `v`
@@ -2590,160 +2484,6 @@ cdef class Gen(Gen_base):
         sig_on()
         return new_gen(Strtex(x.g))
 
-    def Vec(x, long n = 0):
-        """
-        Transform the object `x` into a vector with minimal size `|n|`.
-
-        INPUT:
-
-        - ``x`` -- gen
-
-        - ``n`` -- Make the vector of minimal length `|n|`. If `n > 0`,
-          append zeros; if `n < 0`, prepend zeros.
-
-        OUTPUT:
-
-        A PARI vector (type ``t_VEC``)
-
-        Examples:
-
-        >>> from cypari2 import Pari
-        >>> pari = Pari()
-
-        >>> pari(1).Vec()
-        [1]
-        >>> pari('x^3').Vec()
-        [1, 0, 0, 0]
-        >>> pari('x^3 + 3*x - 2').Vec()
-        [1, 0, 3, -2]
-        >>> pari([1,2,3]).Vec()
-        [1, 2, 3]
-        >>> pari('[1, 2; 3, 4]').Vec()
-        [[1, 3]~, [2, 4]~]
-        >>> pari('"CyPari"').Vec()
-        ["C", "y", "P", "a", "r", "i"]
-        >>> pari('2*x^2 + 3*x^3 + O(x^5)').Vec()
-        [2, 3, 0]
-        >>> pari('2*x^-2 + 3*x^3 + O(x^5)').Vec()
-        [2, 0, 0, 0, 0, 3, 0]
-
-        Note the different term ordering for polynomials and series:
-
-        >>> pari('1 + x + 3*x^3 + O(x^5)').Vec()
-        [1, 1, 0, 3, 0]
-        >>> pari('1 + x + 3*x^3').Vec()
-        [3, 0, 1, 1]
-
-        We demonstate the `n` argument:
-
-        >>> pari([1,2,3,4]).Vec(2)
-        [1, 2, 3, 4]
-        >>> pari([1,2,3,4]).Vec(-2)
-        [1, 2, 3, 4]
-        >>> pari([1,2,3,4]).Vec(6)
-        [1, 2, 3, 4, 0, 0]
-        >>> pari([1,2,3,4]).Vec(-6)
-        [0, 0, 1, 2, 3, 4]
-
-        See also :meth:`Col` (create a column vector) and :meth:`Vecrev`
-        (create a vector in reversed order).
-        """
-        sig_on()
-        return new_gen(_Vec_append(gtovec(x.g), gen_0, n))
-
-    def Vecrev(x, long n = 0):
-        """
-        Transform the object `x` into a vector with minimal size `|n|`.
-        The order of the resulting vector is reversed compared to :meth:`Vec`.
-
-        INPUT:
-
-        - ``x`` -- gen
-
-        - ``n`` -- Make the vector of minimal length `|n|`. If `n > 0`,
-          prepend zeros; if `n < 0`, append zeros.
-
-        OUTPUT:
-
-        A PARI vector (type ``t_VEC``)
-
-        Examples:
-
-        >>> from cypari2 import Pari
-        >>> pari = Pari()
-
-        >>> pari(1).Vecrev()
-        [1]
-        >>> pari('x^3').Vecrev()
-        [0, 0, 0, 1]
-        >>> pari('x^3 + 3*x - 2').Vecrev()
-        [-2, 3, 0, 1]
-        >>> pari([1, 2, 3]).Vecrev()
-        [3, 2, 1]
-        >>> pari('Col([1, 2, 3])').Vecrev()
-        [3, 2, 1]
-        >>> pari('[1, 2; 3, 4]').Vecrev()
-        [[2, 4]~, [1, 3]~]
-        >>> pari('"CyPari"').Vecrev()
-        ["i", "r", "a", "P", "y", "C"]
-
-        We demonstate the `n` argument:
-
-        >>> pari([1,2,3,4]).Vecrev(2)
-        [4, 3, 2, 1]
-        >>> pari([1,2,3,4]).Vecrev(-2)
-        [4, 3, 2, 1]
-        >>> pari([1,2,3,4]).Vecrev(6)
-        [0, 0, 4, 3, 2, 1]
-        >>> pari([1,2,3,4]).Vecrev(-6)
-        [4, 3, 2, 1, 0, 0]
-        """
-        sig_on()
-        return new_gen(_Vec_append(gtovecrev(x.g), gen_0, -n))
-
-    def Vecsmall(x, long n = 0):
-        """
-        Transform the object `x` into a ``t_VECSMALL`` with minimal size `|n|`.
-
-        INPUT:
-
-        - ``x`` -- gen
-
-        - ``n`` -- Make the vector of minimal length `|n|`. If `n > 0`,
-          append zeros; if `n < 0`, prepend zeros.
-
-        OUTPUT:
-
-        A PARI vector of small integers (type ``t_VECSMALL``)
-
-        Examples:
-
-        >>> from cypari2 import Pari
-        >>> pari = Pari()
-
-        >>> pari([1,2,3]).Vecsmall()
-        Vecsmall([1, 2, 3])
-        >>> pari('"CyPari"').Vecsmall()
-        Vecsmall([67, 121, 80, 97, 114, 105])
-        >>> pari(1234).Vecsmall()
-        Vecsmall([1234])
-        >>> pari('x^2 + 2*x + 3').Vecsmall()
-        Vecsmall([1, 2, 3])
-
-        We demonstate the `n` argument:
-
-        >>> pari([1,2,3]).Vecsmall(2)
-        Vecsmall([1, 2, 3])
-        >>> pari([1,2,3]).Vecsmall(-2)
-        Vecsmall([1, 2, 3])
-        >>> pari([1,2,3]).Vecsmall(6)
-        Vecsmall([1, 2, 3, 0, 0, 0])
-        >>> pari([1,2,3]).Vecsmall(-6)
-        Vecsmall([0, 0, 0, 1, 2, 3])
-        """
-        sig_on()
-        return new_gen(_Vec_append(gtovecsmall(x.g), <GEN>0, n))
-
     def bittest(x, long n):
         """
         bittest(x, long n): Returns bit number n (coefficient of
@@ -4994,48 +4734,3 @@ cpdef Gen objtogen(s):
 
     # Simply use the string representation
     return objtogen(str(s))
-
-
-cdef GEN _Vec_append(GEN v, GEN a, long n):
-    """
-    This implements appending zeros (or another constant GEN ``a``) to
-    the result of :meth:`Vec` and similar functions.
-
-    This is a shallow function, copying ``a`` and entries of ``v`` to
-    the result.  The result is simply stored on the PARI stack.
-
-    INPUT:
-
-    - ``v`` -- GEN of type ``t_VEC`` or ``t_COL``
-
-    - ``a`` -- GEN which will be used for the added entries.
-      Normally, this would be ``gen_0``.
-
-    - ``n`` -- Make the vector of minimal length `|n|`. If `n > 0`,
-      append zeros; if `n < 0`, prepend zeros.
-
-    OUTPUT:
-
-    A GEN of the same type as ``v``.
-    """
-    cdef long lenv = lg(v)-1
-    cdef GEN w
-    cdef long i
-    # Do we need to extend the vector with zeros?
-    if n > lenv:
-        w = cgetg(n+1, typ(v))
-        for i from 1 <= i <= lenv:
-            set_gel(w, i, gel(v, i))
-        for i from 1 <= i <= n-lenv:
-            set_gel(w, i+lenv, a)
-        return w
-    elif n < -lenv:
-        n = -n  # Make n positive
-        w = cgetg(n+1, typ(v))
-        for i from 1 <= i <= lenv:
-            set_gel(w, i+(n-lenv), gel(v, i))
-        for i from 1 <= i <= n-lenv:
-            set_gel(w, i, a)
-        return w
-    else:
-        return v