schubmult 2.0.0__py3-none-any.whl

schubmult-2.0.0.dist-info/METADATA

@@ -0,0 +1,455 @@
+Metadata-Version: 2.4
+Name: schubmult
+Version: 2.0.0
+Summary: Package for algebraic computation with single, double, and quantum Schubert polynomials
+Author-email: Matt Samuel <schubmult@gmail.com>
+Maintainer-email: Matt Samuel <schubmult@gmail.com>
+Project-URL: Homepage, http://schubmult.org
+Project-URL: Repository, https://github.com/matthematics/schubmult
+Keywords: Schubert,polynomial,double,algebra
+Classifier: Programming Language :: Python
+Requires-Python: >=3.9
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: numpy
+Requires-Dist: PuLP>=2.7.0
+Requires-Dist: symengine>=0.14.0
+Requires-Dist: sympy>=1.12
+Requires-Dist: psutil
+Requires-Dist: cachetools
+Requires-Dist: sortedcontainers
+Requires-Dist: setuptools>=71.0.0
+Provides-Extra: sage
+Requires-Dist: sagemath-standard; extra == "sage"
+Dynamic: license-file
+
+# schubmult
+
+## Program and package for rapid computation of Littlewood-Richardson coefficients of Schubert polynomials, with optional Sage integration
+
+The main purpose of this python package is for executing scripts to compute coefficients of products of various types of Schubert polynomials. Coproducts can also be computed, as well as substitution of commuting difference operators for quantum double Schubert polynomials. Quantum multiplication also has parabolic subgroup support, computed via the Peterson-Woodward comparison theorem. **Note that except for quantum Schubert polynomial multiplication with the --basic-pieri option, the methodology for quantum/quantum double Schubert polynomials is conjectural at this time.**
+
+[Docs to be hosted on Wiki](https://github.com/matthematics/schubmult/wiki/schubmult-home)
+
+## Basic script command lines, one-line notation
+
+### schubmult_py - ordinary Schubert polynmoials
+
+```
+usage: schubmult_py [-h] [-np] [--code] [--mult MULT [MULT ...]] [--coprod] [--display-mode {basic,pretty,latex}] 
+hyphen-separated list of perms [hyphen-separated list of perms ...]
+
+Compute products of ordinary Schubert polynomials
+
+positional arguments:
+  hyphen-separated list of perms
+                        Space-delimited permutations separated by hyphens, e. g. 3 4 1 2 - 5 1 2 4 3
+
+options:
+  -h, --help            show this help message and exit
+  -np, --no-print       Compute the result but do not print it
+  --code                Permutations represented by the Lehmer code
+  --mult MULT [MULT ...]
+                        Some additional terms in the ring to multiply by
+  --coprod              Compute the coproduct (different syntax: one permutation, then a hyphen, the variable index 
+positions to split on)
+  --display-mode {basic,pretty,latex}
+                        Method of displaying the output. Default basic
+
+Example:
+        schubmult_py 5 1 7 3 2 6 4 - 2 1 6 3 5 4
+            or equivalently
+        schubmult_py --code 4 0 4 1 0 1 - 1 0 3 0 1
+            or alternatively
+        schubmult_py --coprod --code 2 0 3 0 1 - 2 4
+```
+
+
+<!-- ## Quantum commuting difference operators
+
+schubmult_q_double has a feature for displaying the coefficients of the divided difference operators in the evaluation of the quantum double Schubert polynomials on the commuting difference operators of Fomin, Gelfand, and Postnikov. It is necessary to cap the value of n in the group S_n we are working in because as n increases the expression does not stabilize.
+```bash
+schubmult_q_double --nil-hecke 6 --code 2 2 --display-positive
+``` -->
+
+### schubmult_double - Double Schubert polynmoials
+
+```
+usage: schubmult_double [-h] [-np] [--code] [--mult MULT [MULT ...]] [--coprod] [--display-positive] [--optimizer-message] [--down] [-nc] 
+[--mixed-var] [--expand] [--display-mode {basic,pretty,latex}]
+hyphen-separated list of perms [hyphen-separated list of perms ...]
+
+Compute coefficients of product of double Schubert polynomials in the same or different sets of coefficient variables
+
+positional arguments:
+  hyphen-separated list of perms
+                        Space-delimited permutations separated by hyphens, e. g. 3 4 1 2 - 5 1 2 4 3
+
+options:
+  -h, --help            show this help message and exit
+  -np, --no-print       Compute the result but do not print it
+  --code                Permutations represented by the Lehmer code
+  --mult MULT [MULT ...]
+                        Some additional terms in the ring to multiply by
+  --coprod              Compute the coproduct (different syntax: one permutation, then a hyphen, the variable index 
+  positions to split on)
+  --display-positive    Display the result in terms of the positive roots, or if mixed variable attempt to display the 
+  result as a positive algebraic combination of terms of the form y_i - z_j
+  --optimizer-message   Display debug output during integer optimization for --display-positive
+  --down                Reverse multiplication
+  -nc, --no-check       Do not check if positive result matches the original
+  --mixed-var           Used mixed variables y and z
+  --expand              Expand the output rather than leaving it as originally computed (slow)
+  --display-mode {basic,pretty,latex}
+                        Method of displaying the output. Default basic
+
+Example:
+        schubmult_double 5 1 7 3 2 6 4 - 2 1 6 3 5 4 [ --display-positive]
+            or equivalently
+        schubmult_double --code 4 0 4 1 0 1 - 1 0 3 0 1 [ --display-positive]
+            or alternatively
+        schubmult_double --coprod --code 2 0 3 0 1 - 2 4  [ --display-positive]
+```
+
+### schubmult_q - Quantum Schubert polynomials
+
+```
+usage: schubmult_q [-h] [-np] [--code] [--mult MULT [MULT ...]] [--parabolic PARABOLIC [PARABOLIC ...]] [--basic-pieri] 
+[--display-mode {basic,pretty,latex}]
+                   hyphen-separated list of perms [hyphen-separated list of perms ...]
+
+Compute products of quantum Schubert polynomials
+
+positional arguments:
+  hyphen-separated list of perms
+                        Space-delimited permutations separated by hyphens, e. g. 3 4 1 2 - 5 1 2 4 3
+
+options:
+  -h, --help            show this help message and exit
+  -np, --no-print       Compute the result but do not print it
+  --code                Permutations represented by the Lehmer code
+  --mult MULT [MULT ...]
+                        Some additional terms in the ring to multiply by
+  --parabolic PARABOLIC [PARABOLIC ...]
+                        Generators of the parabolic subgroup to compute quantum coeffs for
+  --basic-pieri         Do not apply conjectural computation optimization to quantum
+  --display-mode {basic,pretty,latex}
+                        Method of displaying the output. Default basic
+
+Example:
+        schubmult_q 5 1 7 3 2 6 4 - 2 1 6 3 5 4
+            or equivalently
+        schubmult_q --code 4 0 4 1 0 1 - 1 0 3 0 1
+```
+
+### schubmult_q_double - Quantum double Schubert polynomials
+
+```
+usage: schubmult_q_double [-h] [-np] [--code] [--mult MULT [MULT ...]] [--display-positive] [--optimizer-message] [--down] 
+[-nc] [--mixed-var] [--expand] [--parabolic PARABOLIC [PARABOLIC ...]] [--basic-pieri]
+                          [--nil-hecke N] [--nil-hecke-apply N] [--display-mode {basic,pretty,latex}]
+                          hyphen-separated list of perms [hyphen-separated list of perms ...]
+
+Compute coefficients of products of quantum double Schubert polynomials in the same or different sets of coefficient variables
+
+positional arguments:
+  hyphen-separated list of perms
+                        Space-delimited permutations separated by hyphens, e. g. 3 4 1 2 - 5 1 2 4 3
+
+options:
+  -h, --help            show this help message and exit
+  -np, --no-print       Compute the result but do not print it
+  --code                Permutations represented by the Lehmer code
+  --mult MULT [MULT ...]
+                        Some additional terms in the ring to multiply by
+  --display-positive    Display the result in terms of the positive roots, or if mixed variable attempt to display the 
+result as a positive algebraic combination of terms of the form y_i - z_j
+  --optimizer-message   Display debug output during integer optimization for --display-positive
+  --down                Reverse multiplication
+  -nc, --no-check       Do not check if positive result matches the original
+  --mixed-var           Used mixed variables y and z
+  --expand              Expand the output rather than leaving it as originally computed (slow)
+  --parabolic PARABOLIC [PARABOLIC ...]
+                        Generators of the parabolic subgroup to compute quantum coeffs for
+  --basic-pieri         Do not apply conjectural computation optimization to quantum
+  --nil-hecke N         Substitute up to N of Fomin-Gelfand-Postnikov commuting difference operators
+  --nil-hecke-apply N   Substitute commuting difference operators for perm1, then apply to Schub indexed by perm2
+  --display-mode {basic,pretty,latex}
+                        Method of displaying the output. Default basic
+
+Example:
+        schubmult_q_double 5 1 7 3 2 6 4 - 2 1 6 3 5 4 [ --display-positive]
+            or equivalently
+        schubmult_q_double --code 4 0 4 1 0 1 - 1 0 3 0 1 [ --display-positive]
+```
+
+## Diplaying the result positively
+
+The command line argument `--display-positive `is available in schubmult_double and schubmult_q_double, which displays the result positively (if possible, this is still only always possible conjecturally). It will fail and print out the offending case if it finds a counterexample. This is highly processor intensive.
+
+
+Runtime will vary tremendously by case. The general problem is #P-hard. Though the result is always nonnegative (which at least is known for schubmult_py, schubmult_q, schubmult_double, and schubmult_q_double) and the problem is in GapP, it is not known to be in #P at this time.
+
+schubmult_py is for multiplying ordinary Schubert polynomials. schubmult_double is for multiplying double Schubert polynomials in different sets of coefficient variables (`--mixed-var`) or in the same set of coefficient variables (by default). Similarly, schubmult_q is for multiplying quantum Schubert polynomials, schubmult_q_double is for multiplying quantum double Schubert polynomials (in different sets of coefficient variables with the `--mixed-var` option, or the same set), or in other words it computes the Gromov-Witten invariants, equivariant Gromov-Witten invariants, and (mixed?) equivariant Gromov-Witten invariants of the complete flag variety. All have the same command line syntax as schubmult, except when using the --code option. schubmult_double/schubmult_q_double display the result with nonnegative coefficients (and the q variables) with the `--display-positive` option (either in the negative simple roots r_i, or in y_i - z_j for `--mixed-var`).
+
+schubmult_xx --coprod allows you to split (double) Schubert polynomials along certain indices (not available for quantum). It takes one permutation as an argument, followed by a dash -, then the set of indices you would like to split on.
+
+# Sage integration (as of version 2.0.0)
+
+[SageMath](https://www.sagemath.org/) is a computer algebra system that, while wonderful, is monstrously large and only works on posix-based operating systems (including WSL VMs, so it is still usable on Windows). This is why Sage support is provided optionally in schubmult. The syntax to install the Sage dependencies is
+
+```
+pip install schubmult[sage]
+```
+
+This will install the [sagemath-standard](https://pypi.org/project/sagemath-standard/) python package in addition to the other dependencies. **Again, this only works on Linux, MacOS, or WSL.** To use with a currently installed SageMath distribution, use sage's python interpreter to install the package (the `[sage]` piece is not required in that case).
+
+## Basic sage example
+
+```python
+sage: from schubmult.sage_integration import FastSchubertPolynomial, FastDoubleSchubertPolynomialRing,
+FastQuantumSchubertPolynomial, FastQuantumDoubleSchubertPolynomialRing
+sage: SingleRing = FastSchubertPolynomialRing(ZZ, 100, "x")
+sage: SingleRing([3,4,1,2])
+Sx[3, 4, 1, 2]
+sage: SingleRing([3,4,1,2]) * SingleRing([5,1,4,2,3])
+Sx[7, 3, 4, 1, 2, 5, 6] + Sx[7, 4, 2, 1, 3, 5, 6] + Sx[7, 5, 1, 2, 3, 4,
+6]
+```
+
+## Mixed variable (Molev-Sagan) type products
+
+```python
+sage: DoubleRing = FastDoubleSchubertPolynomialRing(ZZ, 100, "x", ("y", "z"))
+sage: DoubleRing([3,4,1,2])*DoubleRing([1,4,2,3])
+(y1*y2-y1*y4-y2*y4+y4^2)*Sx([3, 4, 1, 2], 'y') + (-y1-y2+y4+y5)*Sx([3, 5, 1, 2, 4], 'y') + Sx([3, 6, 1, 
+2, 4, 5], 'y')
+sage: DoubleRing([3,4,1,2]) * DoubleRing([1,4,2,3],"z")
+(y3^2+y3*y4+y4^2-y3*z1-y4*z1-y3*z2-y4*z2+z1*z2-y3*z3-y4*z3+z1*z3+z2*z3)*Sx([3, 4, 1, 2], 'y') + (y3+y4+
+y5-z1-z2-z3)*Sx([3, 5, 1, 2, 4], 'y') + Sx([3, 6, 1, 2, 4, 5], 
+'y')
+```
+
+## expand()
+
+```python
+sage: SingleRingQ = FastQuantumSchubertPolynomialRing(ZZ, 100, "x")
+sage: SingleRingQ([2,3,1,4]).expand()
+x1*x2 + q_1
+```
+
+## Coercion
+
+Coercion was implemented as widely as possible.
+```python
+sage: DoubleRing([1,4,2,3],"z") * SingleRing([3,4,1,2])
+z1^2*z4^2*Sx([1, 4, 2, 3], 'z') + (z1^2*z4+z1^2*z5)*Sx([1, 5, 2, 3, 4], 'z') + z1^2*Sx([1, 6, 2, 3, 4, 
+5], 'z') + (z1*z4^2+z2*z4^2)*Sx([2, 4, 1, 3], 'z') + (z1*z4+z2*z4+z1*z5+z2*z5)*Sx([2, 5, 1, 3, 4], 'z') +
+ (z1+z2)*Sx([2, 6, 1, 3, 4, 5], 'z') + z4^2*Sx([3, 4, 1, 2], 'z') + (z4+z5)*Sx([3, 5, 1, 2, 4], 'z') + 
+Sx([3, 6, 1, 2, 4, 5], 'z')
+sage: SingleRingQ([2,3,1,4]) * SingleRing([4,1,3,2])
+(-2*q_1^2*q_2+q_1*q_2*q_3)*QSx[1] + q_1*q_2*QSx[1, 3, 4, 2] + (-q_1^2)*QSx[2, 3, 1] + q_1*QSx[2, 4, 3, 
+1] + (-q_1*q_2)*QSx[3, 1, 2] + q_1*QSx[3, 2, 4, 1] + q_1*QSx[3, 4, 1, 2] + (-q_1)*QSx[4, 2, 1, 3] + 
+QSx[5, 2, 3, 1, 4] + QSx[5, 3, 1, 2, 4]
+sage: R.<x1, x2> = PolynomialRing(ZZ, 2)
+sage: SingleRing([1,3,2]) - x1 - x2 == 0
+True
+```
+
+## Coproducts
+
+FastSchubertPolynomialRing and FastDoubleSchubertPolynomialRings are bialgebras and each element implements the `coproduct()` member function. `set_coproduct_indices()` on the base ring will determine the variables to partition on.
+```ada
+sage: DoubleRing.set_coproduct_indices((1,3))
+sage: DoubleRing([4,1,5,2,3], "z").coproduct()
+(y1^2-y1*z2-y1*z3+z2*z3)*Sx([4, 1, 2, 3], 'z') # Sx([1], 'y') + (y1+y2-z2-z3)*Sx([4, 1, 2, 3], 'z') # 
+Sx([2, 1], 'y') + Sx([4, 1, 2, 3], 'z') # Sx([3, 1, 2], 'y') + (y1-z3)*Sx([4, 2, 1, 3], 'z') # Sx([1], 
+'y') + Sx([4, 2, 1, 3], 'z') # Sx([2, 1], 'y') + Sx([4, 3, 1, 2], 'z') # Sx([1], 
+'y')
+```
+
+## Demonstration of quantum double mixed products
+
+```python
+sage: QuantumDoubleRing = FastQuantumDoubleSchubertPolynomialRing(ZZ, 100, "x", ("y","z"))
+sage: QuantumDoubleRing([4,1,3,2])*QuantumDoubleRing([5,1,3,2,4], "z")
+(q_1*q_2*q_3*y1^3+q_1*q_2*q_3*y1^2*y4+q_1*q_2*q_3*y1*y4^2+q_1*q_2*q_3*y4^3-q_1*q_2*q_3*y1^2*z1-
+q_1*q_2*q_3*y1*y4*z1-q_1*q_2*q_3*y4^2*z1-q_1*q_2*q_3*y1^2*z2-q_1*q_2*q_3*y1*y4*z2-q_1*q_2*q_3*y4^2*z2+
+q_1*q_2*q_3*y1*z1*z2+q_1*q_2*q_3*y4*z1*z2-q_1*q_2*q_3*y1^2*z3-q_1*q_2*q_3*y1*y4*z3-q_1*q_2*q_3*y4^2*z3+
+q_1*q_2*q_3*y1*z1*z3+q_1*q_2*q_3*y4*z1*z3+q_1*q_2*q_3*y1*z2*z3+q_1*q_2*q_3*y4*z2*z3-q_1*q_2*q_3*z1*z2*z3-
+q_1*q_2*q_3*y1^2*z4-q_1*q_2*q_3*y1*y4*z4-q_1*q_2*q_3*y4^2*z4+q_1*q_2*q_3*y1*z1*z4+q_1*q_2*q_3*y4*z1*z4+
+q_1*q_2*q_3*y1*z2*z4+q_1*q_2*q_3*y4*z2*z4-q_1*q_2*q_3*z1*z2*z4+q_1*q_2*q_3*y1*z3*z4+q_1*q_2*q_3*y4*z3*z4-
+q_1*q_2*q_3*z1*z3*z4-q_1*q_2*q_3*z2*z3*z4)*QSx([1], 'y') + (q_1*q_2*q_3*y1^2+q_1*q_2*q_3*y1*y4+
+q_1*q_2*q_3*y4^2+q_1*q_2*q_3*y1*y5+q_1*q_2*q_3*y4*y5+q_1*q_2*q_3*y5^2-q_1*q_2*q_3*y1*z1-
+q_1*q_2*q_3*y4*z1-q_1*q_2*q_3*y5*z1-q_1*q_2*q_3*y1*z2-q_1*q_2*q_3*y4*z2-q_1*q_2*q_3*y5*z2+
+q_1*q_2*q_3*z1*z2-q_1*q_2*q_3*y1*z3-q_1*q_2*q_3*y4*z3-q_1*q_2*q_3*y5*z3+q_1*q_2*q_3*z1*z3+
+q_1*q_2*q_3*z2*z3-q_1*q_2*q_3*y1*z4-q_1*q_2*q_3*y4*z4-q_1*q_2*q_3*y5*z4+q_1*q_2*q_3*z1*z4+
+q_1*q_2*q_3*z2*z4+q_1*q_2*q_3*z3*z4)*QSx([1, 2, 3, 5, 4], 'y') + (q_1*q_2*q_3*y1+q_1*q_2*q_3*y4+
+q_1*q_2*q_3*y5+q_1*q_2*q_3*y6-q_1*q_2*q_3*z1-q_1*q_2*q_3*z2-q_1*q_2*q_3*z3-q_1*q_2*q_3*z4)*QSx([1, 2, 3, 
+6, 4, 5], 'y') + q_1*q_2*q_3*QSx([1, 2, 3, 7, 4, 5, 6], 'y') + (q_1*q_3*y1^3+q_1*q_3*y1^2*y4+
+q_1*q_3*y1*y4^2+q_1*q_3*y4^3-q_1*q_3*y1^2*z1-q_1*q_3*y1*y4*z1-q_1*q_3*y4^2*z1-q_1*q_3*y1^2*z2-
+q_1*q_3*y1*y4*z2-q_1*q_3*y4^2*z2+q_1*q_3*y1*z1*z2+q_1*q_3*y4*z1*z2-q_1*q_3*y1^2*z3-q_1*q_3*y1*y4*z3-
+q_1*q_3*y4^2*z3+q_1*q_3*y1*z1*z3+q_1*q_3*y4*z1*z3+q_1*q_3*y1*z2*z3+q_1*q_3*y4*z2*z3-q_1*q_3*z1*z2*z3-
+q_1*q_3*y1^2*z4-q_1*q_3*y1*y4*z4-q_1*q_3*y4^2*z4+q_1*q_3*y1*z1*z4+q_1*q_3*y4*z1*z4+q_1*q_3*y1*z2*z4+
+q_1*q_3*y4*z2*z4-q_1*q_3*z1*z2*z4+q_1*q_3*y1*z3*z4+q_1*q_3*y4*z3*z4-q_1*q_3*z1*z3*z4-
+q_1*q_3*z2*z3*z4)*QSx([1, 4, 2, 3], 'y') + (q_1*y1^3*y3+q_1*y1^3*y4+q_1*y1^2*y3*y4+q_1*y1^2*y4^2+
+q_1*y1*y3*y4^2+q_1*y1*y4^3+q_1*y3*y4^3-q_1*y1^3*z1-q_1*y1^2*y3*z1-2*q_1*y1^2*y4*z1-q_1*y1*y3*y4*z1-
+2*q_1*y1*y4^2*z1-q_1*y3*y4^2*z1-q_1*y4^3*z1+q_1*y1^2*z1^2+q_1*y1*y4*z1^2+q_1*y4^2*z1^2-q_1*y1^3*z2-
+q_1*y1^2*y3*z2-2*q_1*y1^2*y4*z2-q_1*y1*y3*y4*z2-2*q_1*y1*y4^2*z2-q_1*y3*y4^2*z2-q_1*y4^3*z2+
+2*q_1*y1^2*z1*z2+q_1*y1*y3*z1*z2+3*q_1*y1*y4*z1*z2+q_1*y3*y4*z1*z2+2*q_1*y4^2*z1*z2-q_1*y1*z1^2*z2-
+q_1*y4*z1^2*z2+q_1*y1^2*z2^2+q_1*y1*y4*z2^2+q_1*y4^2*z2^2-q_1*y1*z1*z2^2-q_1*y4*z1*z2^2-q_1*y1^2*y3*z3-
+q_1*y1^2*y4*z3-q_1*y1*y3*y4*z3-q_1*y1*y4^2*z3-q_1*y3*y4^2*z3+q_1*y1^2*z1*z3+q_1*y1*y3*z1*z3+
+2*q_1*y1*y4*z1*z3+q_1*y3*y4*z1*z3+q_1*y4^2*z1*z3-q_1*y1*z1^2*z3-q_1*y4*z1^2*z3+q_1*y1^2*z2*z3+
+q_1*y1*y3*z2*z3+2*q_1*y1*y4*z2*z3+q_1*y3*y4*z2*z3+q_1*y4^2*z2*z3-2*q_1*y1*z1*z2*z3-q_1*y3*z1*z2*z3-
+2*q_1*y4*z1*z2*z3+q_1*z1^2*z2*z3-q_1*y1*z2^2*z3-q_1*y4*z2^2*z3+q_1*z1*z2^2*z3-q_1*y1^2*y3*z4-
+q_1*y1^2*y4*z4-q_1*y1*y3*y4*z4-q_1*y1*y4^2*z4-q_1*y3*y4^2*z4+q_1*y1^2*z1*z4+q_1*y1*y3*z1*z4+
+2*q_1*y1*y4*z1*z4+q_1*y3*y4*z1*z4+q_1*y4^2*z1*z4-q_1*y1*z1^2*z4-q_1*y4*z1^2*z4+q_1*y1^2*z2*z4+
+q_1*y1*y3*z2*z4+2*q_1*y1*y4*z2*z4+q_1*y3*y4*z2*z4+q_1*y4^2*z2*z4-2*q_1*y1*z1*z2*z4-q_1*y3*z1*z2*z4-
+2*q_1*y4*z1*z2*z4+q_1*z1^2*z2*z4-q_1*y1*z2^2*z4-q_1*y4*z2^2*z4+q_1*z1*z2^2*z4+q_1*y1*y3*z3*z4+
+q_1*y1*y4*z3*z4+q_1*y3*y4*z3*z4-q_1*y1*z1*z3*z4-q_1*y3*z1*z3*z4-q_1*y4*z1*z3*z4+q_1*z1^2*z3*z4-
+q_1*y1*z2*z3*z4-q_1*y3*z2*z3*z4-q_1*y4*z2*z3*z4+q_1*z1*z2*z3*z4+q_1*z2^2*z3*z4)*QSx([1, 4, 3, 2], 'y') + 
+(q_1*y1^3+q_1*y1^2*y4+q_1*y1*y4^2+q_1*y4^3-q_1*y1^2*z1-q_1*y1*y4*z1-q_1*y4^2*z1-q_1*y1^2*z2-q_1*y1*y4*z2-
+q_1*y4^2*z2+q_1*y1*z1*z2+q_1*y4*z1*z2-q_1*y1^2*z3-q_1*y1*y4*z3-q_1*y4^2*z3+q_1*y1*z1*z3+q_1*y4*z1*z3+
+q_1*y1*z2*z3+q_1*y4*z2*z3-q_1*z1*z2*z3-q_1*y1^2*z4-q_1*y1*y4*z4-q_1*y4^2*z4+q_1*y1*z1*z4+q_1*y4*z1*z4+
+q_1*y1*z2*z4+q_1*y4*z2*z4-q_1*z1*z2*z4+q_1*y1*z3*z4+q_1*y4*z3*z4-q_1*z1*z3*z4-q_1*z2*z3*z4)*QSx([1, 4, 
+5, 2, 3], 'y') + (q_1*q_3*y1^2+q_1*q_3*y1*y4+q_1*q_3*y4^2+q_1*q_3*y1*y5+q_1*q_3*y4*y5+q_1*q_3*y5^2-
+q_1*q_3*y1*z1-q_1*q_3*y4*z1-q_1*q_3*y5*z1-q_1*q_3*y1*z2-q_1*q_3*y4*z2-q_1*q_3*y5*z2+q_1*q_3*z1*z2-
+q_1*q_3*y1*z3-q_1*q_3*y4*z3-q_1*q_3*y5*z3+q_1*q_3*z1*z3+q_1*q_3*z2*z3-q_1*q_3*y1*z4-q_1*q_3*y4*z4-
+q_1*q_3*y5*z4+q_1*q_3*z1*z4+q_1*q_3*z2*z4+q_1*q_3*z3*z4)*QSx([1, 5, 2, 3, 4], 'y') + (q_1*y1^3+
+q_1*y1^2*y3+q_1*y1^2*y4+q_1*y1*y3*y4+q_1*y1*y4^2+q_1*y3*y4^2+q_1*y1^2*y5+q_1*y1*y3*y5+q_1*y1*y4*y5+
+q_1*y3*y4*y5+q_1*y1*y5^2+q_1*y3*y5^2-2*q_1*y1^2*z1-q_1*y1*y3*z1-2*q_1*y1*y4*z1-q_1*y3*y4*z1-q_1*y4^2*z1-
+2*q_1*y1*y5*z1-q_1*y3*y5*z1-q_1*y4*y5*z1-q_1*y5^2*z1+q_1*y1*z1^2+q_1*y4*z1^2+q_1*y5*z1^2-2*q_1*y1^2*z2-
+q_1*y1*y3*z2-2*q_1*y1*y4*z2-q_1*y3*y4*z2-q_1*y4^2*z2-2*q_1*y1*y5*z2-q_1*y3*y5*z2-q_1*y4*y5*z2-
+q_1*y5^2*z2+3*q_1*y1*z1*z2+q_1*y3*z1*z2+2*q_1*y4*z1*z2+2*q_1*y5*z1*z2-q_1*z1^2*z2+q_1*y1*z2^2+
+q_1*y4*z2^2+q_1*y5*z2^2-q_1*z1*z2^2-q_1*y1^2*z3-q_1*y1*y3*z3-q_1*y1*y4*z3-q_1*y3*y4*z3-q_1*y1*y5*z3-
+q_1*y3*y5*z3+2*q_1*y1*z1*z3+q_1*y3*z1*z3+q_1*y4*z1*z3+q_1*y5*z1*z3-q_1*z1^2*z3+2*q_1*y1*z2*z3+
+q_1*y3*z2*z3+q_1*y4*z2*z3+q_1*y5*z2*z3-2*q_1*z1*z2*z3-q_1*z2^2*z3-q_1*y1^2*z4-q_1*y1*y3*z4-q_1*y1*y4*z4-
+q_1*y3*y4*z4-q_1*y1*y5*z4-q_1*y3*y5*z4+2*q_1*y1*z1*z4+q_1*y3*z1*z4+q_1*y4*z1*z4+q_1*y5*z1*z4-q_1*z1^2*z4+
+2*q_1*y1*z2*z4+q_1*y3*z2*z4+q_1*y4*z2*z4+q_1*y5*z2*z4-2*q_1*z1*z2*z4-q_1*z2^2*z4+q_1*y1*z3*z4+
+q_1*y3*z3*z4-q_1*z1*z3*z4-q_1*z2*z3*z4)*QSx([1, 5, 3, 2, 4], 'y') + (q_1*y1^2+q_1*y1*y4+q_1*y4^2+
+q_1*y1*y5+q_1*y4*y5+q_1*y5^2-q_1*y1*z1-q_1*y4*z1-q_1*y5*z1-q_1*y1*z2-q_1*y4*z2-q_1*y5*z2+q_1*z1*z2-
+q_1*y1*z3-q_1*y4*z3-q_1*y5*z3+q_1*z1*z3+q_1*z2*z3-q_1*y1*z4-q_1*y4*z4-q_1*y5*z4+q_1*z1*z4+q_1*z2*z4+
+q_1*z3*z4)*QSx([1, 5, 4, 2, 3], 'y') + (q_1*q_3*y1+q_1*q_3*y4+q_1*q_3*y5+q_1*q_3*y6-q_1*q_3*z1-
+q_1*q_3*z2-q_1*q_3*z3-q_1*q_3*z4)*QSx([1, 6, 2, 3, 4, 5], 'y') + (q_1*y1^2+q_1*y1*y3+q_1*y1*y4+q_1*y3*y4+
+q_1*y1*y5+q_1*y3*y5+q_1*y1*y6+q_1*y3*y6-2*q_1*y1*z1-q_1*y3*z1-q_1*y4*z1-q_1*y5*z1-q_1*y6*z1+q_1*z1^2-
+2*q_1*y1*z2-q_1*y3*z2-q_1*y4*z2-q_1*y5*z2-q_1*y6*z2+2*q_1*z1*z2+q_1*z2^2-q_1*y1*z3-q_1*y3*z3+q_1*z1*z3+
+q_1*z2*z3-q_1*y1*z4-q_1*y3*z4+q_1*z1*z4+q_1*z2*z4)*QSx([1, 6, 3, 2, 4, 5], 'y') + (q_1*y1+q_1*y4+q_1*y5+
+q_1*y6-q_1*z1-q_1*z2-q_1*z3-q_1*z4)*QSx([1, 6, 4, 2, 3, 5], 'y') + q_1*q_3*QSx([1, 7, 2, 3, 4, 5, 6], 
+'y') + (q_1*y1+q_1*y3-q_1*z1-q_1*z2)*QSx([1, 7, 3, 2, 4, 5, 6], 'y') + q_1*QSx([1, 7, 4, 2, 3, 5, 6], 
+'y') + (q_1*q_2*q_3*y1^2+q_1*q_2*q_3*y1*y2+q_1*q_2*q_3*y2^2+q_1*q_2*q_3*y1*y4+q_1*q_2*q_3*y2*y4+
+q_1*q_2*q_3*y4^2-q_1*q_2*q_3*y1*z1-q_1*q_2*q_3*y2*z1-q_1*q_2*q_3*y4*z1-q_1*q_2*q_3*y1*z2-
+q_1*q_2*q_3*y2*z2-q_1*q_2*q_3*y4*z2+q_1*q_2*q_3*z1*z2-q_1*q_2*q_3*y1*z3-q_1*q_2*q_3*y2*z3-
+q_1*q_2*q_3*y4*z3+q_1*q_2*q_3*z1*z3+q_1*q_2*q_3*z2*z3-q_1*q_2*q_3*y1*z4-q_1*q_2*q_3*y2*z4-
+q_1*q_2*q_3*y4*z4+q_1*q_2*q_3*z1*z4+q_1*q_2*q_3*z2*z4+q_1*q_2*q_3*z3*z4)*QSx([2, 1], 'y') + 
+(q_1*q_2*q_3*y1+q_1*q_2*q_3*y2+q_1*q_2*q_3*y4+q_1*q_2*q_3*y5-q_1*q_2*q_3*z1-q_1*q_2*q_3*z2-
+q_1*q_2*q_3*z3-q_1*q_2*q_3*z4)*QSx([2, 1, 3, 5, 4], 'y') + q_1*q_2*q_3*QSx([2, 1, 3, 6, 4, 5], 'y') + 
+(q_1*q_3*y1^2+q_1*q_3*y1*y2+q_1*q_3*y2^2+q_1*q_3*y1*y4+q_1*q_3*y2*y4+q_1*q_3*y4^2-q_1*q_3*y1*z1-
+q_1*q_3*y2*z1-q_1*q_3*y4*z1-q_1*q_3*y1*z2-q_1*q_3*y2*z2-q_1*q_3*y4*z2+q_1*q_3*z1*z2-q_1*q_3*y1*z3-
+q_1*q_3*y2*z3-q_1*q_3*y4*z3+q_1*q_3*z1*z3+q_1*q_3*z2*z3-q_1*q_3*y1*z4-q_1*q_3*y2*z4-q_1*q_3*y4*z4+
+q_1*q_3*z1*z4+q_1*q_3*z2*z4+q_1*q_3*z3*z4)*QSx([2, 4, 1, 3], 'y') + (q_1*y1^2*y3+q_1*y1*y2*y3+
+q_1*y2^2*y3+q_1*y1^2*y4+q_1*y1*y2*y4+q_1*y2^2*y4+q_1*y1*y3*y4+q_1*y2*y3*y4+q_1*y1*y4^2+q_1*y2*y4^2+
+q_1*y3*y4^2+q_1*y4^3-q_1*y1^2*z1-q_1*y1*y2*z1-q_1*y2^2*z1-q_1*y1*y3*z1-q_1*y2*y3*z1-2*q_1*y1*y4*z1-
+2*q_1*y2*y4*z1-q_1*y3*y4*z1-2*q_1*y4^2*z1+q_1*y1*z1^2+q_1*y2*z1^2+q_1*y4*z1^2-q_1*y1^2*z2-q_1*y1*y2*z2-
+q_1*y2^2*z2-q_1*y1*y3*z2-q_1*y2*y3*z2-2*q_1*y1*y4*z2-2*q_1*y2*y4*z2-q_1*y3*y4*z2-2*q_1*y4^2*z2+
+2*q_1*y1*z1*z2+2*q_1*y2*z1*z2+q_1*y3*z1*z2+3*q_1*y4*z1*z2-q_1*z1^2*z2+q_1*y1*z2^2+q_1*y2*z2^2+
+q_1*y4*z2^2-q_1*z1*z2^2-q_1*y1*y3*z3-q_1*y2*y3*z3-q_1*y1*y4*z3-q_1*y2*y4*z3-q_1*y3*y4*z3-q_1*y4^2*z3+
+q_1*y1*z1*z3+q_1*y2*z1*z3+q_1*y3*z1*z3+2*q_1*y4*z1*z3-q_1*z1^2*z3+q_1*y1*z2*z3+q_1*y2*z2*z3+q_1*y3*z2*z3+
+2*q_1*y4*z2*z3-2*q_1*z1*z2*z3-q_1*z2^2*z3-q_1*y1*y3*z4-q_1*y2*y3*z4-q_1*y1*y4*z4-q_1*y2*y4*z4-
+q_1*y3*y4*z4-q_1*y4^2*z4+q_1*y1*z1*z4+q_1*y2*z1*z4+q_1*y3*z1*z4+2*q_1*y4*z1*z4-q_1*z1^2*z4+q_1*y1*z2*z4+
+q_1*y2*z2*z4+q_1*y3*z2*z4+2*q_1*y4*z2*z4-2*q_1*z1*z2*z4-q_1*z2^2*z4+q_1*y3*z3*z4+q_1*y4*z3*z4-
+q_1*z1*z3*z4-q_1*z2*z3*z4)*QSx([2, 4, 3, 1], 'y') + (q_1*y1^2+q_1*y1*y2+q_1*y2^2+q_1*y1*y4+q_1*y2*y4+
+q_1*y4^2-q_1*y1*z1-q_1*y2*z1-q_1*y4*z1-q_1*y1*z2-q_1*y2*z2-q_1*y4*z2+q_1*z1*z2-q_1*y1*z3-q_1*y2*z3-
+q_1*y4*z3+q_1*z1*z3+q_1*z2*z3-q_1*y1*z4-q_1*y2*z4-q_1*y4*z4+q_1*z1*z4+q_1*z2*z4+q_1*z3*z4)*QSx([2, 4, 5, 
+1, 3], 'y') + (q_1*q_3*y1+q_1*q_3*y2+q_1*q_3*y4+q_1*q_3*y5-q_1*q_3*z1-q_1*q_3*z2-q_1*q_3*z3-
+q_1*q_3*z4)*QSx([2, 5, 1, 3, 4], 'y') + (q_1*y1^2+q_1*y1*y2+q_1*y2^2+q_1*y1*y3+q_1*y2*y3+q_1*y1*y4+
+q_1*y2*y4+q_1*y3*y4+q_1*y4^2+q_1*y1*y5+q_1*y2*y5+q_1*y3*y5+q_1*y4*y5+q_1*y5^2-2*q_1*y1*z1-2*q_1*y2*z1-
+q_1*y3*z1-2*q_1*y4*z1-2*q_1*y5*z1+q_1*z1^2-2*q_1*y1*z2-2*q_1*y2*z2-q_1*y3*z2-2*q_1*y4*z2-2*q_1*y5*z2+
+3*q_1*z1*z2+q_1*z2^2-q_1*y1*z3-q_1*y2*z3-q_1*y3*z3-q_1*y4*z3-q_1*y5*z3+2*q_1*z1*z3+2*q_1*z2*z3-q_1*y1*z4-
+q_1*y2*z4-q_1*y3*z4-q_1*y4*z4-q_1*y5*z4+2*q_1*z1*z4+2*q_1*z2*z4+q_1*z3*z4)*QSx([2, 5, 3, 1, 4], 'y') + 
+(q_1*y1+q_1*y2+q_1*y4+q_1*y5-q_1*z1-q_1*z2-q_1*z3-q_1*z4)*QSx([2, 5, 4, 1, 3], 'y') + q_1*q_3*QSx([2, 6, 
+1, 3, 4, 5], 'y') + (q_1*y1+q_1*y2+q_1*y3+q_1*y4+q_1*y5+q_1*y6-2*q_1*z1-2*q_1*z2-q_1*z3-q_1*z4)*QSx([2, 
+6, 3, 1, 4, 5], 'y') + q_1*QSx([2, 6, 4, 1, 3, 5], 'y') + q_1*QSx([2, 7, 3, 1, 4, 5, 6], 'y') + 
+(q_1*q_2*q_3*y1+q_1*q_2*q_3*y2+q_1*q_2*q_3*y3+q_1*q_2*q_3*y4-q_1*q_2*q_3*z1-q_1*q_2*q_3*z2-
+q_1*q_2*q_3*z3-q_1*q_2*q_3*z4)*QSx([3, 1, 2], 'y') + q_1*q_2*q_3*QSx([3, 1, 2, 5, 4], 'y') + 
+(q_1*y1^2*y3+q_1*y1*y3^2+q_1*y1^2*y4+2*q_1*y1*y3*y4+q_1*y3^2*y4+q_1*y1*y4^2+q_1*y3*y4^2-q_1*y1^2*z1-
+2*q_1*y1*y3*z1-q_1*y3^2*z1-2*q_1*y1*y4*z1-2*q_1*y3*y4*z1-q_1*y4^2*z1+q_1*y1*z1^2+q_1*y3*z1^2+q_1*y4*z1^2-
+q_1*y1^2*z2-2*q_1*y1*y3*z2-q_1*y3^2*z2-2*q_1*y1*y4*z2-2*q_1*y3*y4*z2-q_1*y4^2*z2+2*q_1*y1*z1*z2+
+2*q_1*y3*z1*z2+2*q_1*y4*z1*z2-q_1*z1^2*z2+q_1*y1*z2^2+q_1*y3*z2^2+q_1*y4*z2^2-q_1*z1*z2^2-q_1*y1*y3*z3-
+q_1*y1*y4*z3-q_1*y3*y4*z3+q_1*y1*z1*z3+q_1*y3*z1*z3+q_1*y4*z1*z3-q_1*z1^2*z3+q_1*y1*z2*z3+q_1*y3*z2*z3+
+q_1*y4*z2*z3-q_1*z1*z2*z3-q_1*z2^2*z3-q_1*y1*y3*z4-q_1*y1*y4*z4-q_1*y3*y4*z4+q_1*y1*z1*z4+q_1*y3*z1*z4+
+q_1*y4*z1*z4-q_1*z1^2*z4+q_1*y1*z2*z4+q_1*y3*z2*z4+q_1*y4*z2*z4-q_1*z1*z2*z4-q_1*z2^2*z4+q_1*q_3*y1+
+q_1*q_3*y2+q_1*q_3*y3+q_1*q_3*y4-q_1*q_3*z1-q_1*q_3*z2-q_1*q_3*z3-q_1*q_3*z4)*QSx([3, 4, 1, 2], 'y') + 
+(q_1*y1*y3+q_1*y2*y3+q_1*y3^2+q_1*y1*y4+q_1*y2*y4+2*q_1*y3*y4+q_1*y4^2-q_1*y1*z1-q_1*y2*z1-2*q_1*y3*z1-
+2*q_1*y4*z1+q_1*z1^2-q_1*y1*z2-q_1*y2*z2-2*q_1*y3*z2-2*q_1*y4*z2+2*q_1*z1*z2+q_1*z2^2-q_1*y3*z3-
+q_1*y4*z3+q_1*z1*z3+q_1*z2*z3-q_1*y3*z4-q_1*y4*z4+q_1*z1*z4+q_1*z2*z4)*QSx([3, 4, 2, 1], 'y') + (q_1*y1+
+q_1*y2+q_1*y3+q_1*y4-q_1*z1-q_1*z2-q_1*z3-q_1*z4)*QSx([3, 4, 5, 1, 2], 'y') + (q_1*y1^2+2*q_1*y1*y3+
+q_1*y3^2+q_1*y1*y4+q_1*y3*y4+q_1*y1*y5+q_1*y3*y5-2*q_1*y1*z1-2*q_1*y3*z1-q_1*y4*z1-q_1*y5*z1+q_1*z1^2-
+2*q_1*y1*z2-2*q_1*y3*z2-q_1*y4*z2-q_1*y5*z2+2*q_1*z1*z2+q_1*z2^2-q_1*y1*z3-q_1*y3*z3+q_1*z1*z3+q_1*z2*z3-
+q_1*y1*z4-q_1*y3*z4+q_1*z1*z4+q_1*z2*z4+q_1*q_3)*QSx([3, 5, 1, 2, 4], 'y') + (q_1*y1+q_1*y2+2*q_1*y3+
+q_1*y4+q_1*y5-2*q_1*z1-2*q_1*z2-q_1*z3-q_1*z4)*QSx([3, 5, 2, 1, 4], 'y') + q_1*QSx([3, 5, 4, 1, 2], 'y') 
++ (q_1*y1+q_1*y3-q_1*z1-q_1*z2)*QSx([3, 6, 1, 2, 4, 5], 'y') + q_1*QSx([3, 6, 2, 1, 4, 5], 'y') + 
+(q_3*y4^4-q_3*y4^3*z1-q_3*y4^3*z2+q_3*y4^2*z1*z2-q_3*y4^3*z3+q_3*y4^2*z1*z3+q_3*y4^2*z2*z3-
+q_3*y4*z1*z2*z3-q_3*y4^3*z4+q_3*y4^2*z1*z4+q_3*y4^2*z2*z4-q_3*y4*z1*z2*z4+q_3*y4^2*z3*z4-q_3*y4*z1*z3*z4-
+q_3*y4*z2*z3*z4+q_3*z1*z2*z3*z4)*QSx([4, 1, 2, 3], 'y') + (y1*y4^4+y3*y4^4-y1*y4^3*z1-y3*y4^3*z1-y4^4*z1+
+y4^3*z1^2-y1*y4^3*z2-y3*y4^3*z2-y4^4*z2+y1*y4^2*z1*z2+y3*y4^2*z1*z2+2*y4^3*z1*z2-y4^2*z1^2*z2+y4^3*z2^2-
+y4^2*z1*z2^2-y1*y4^3*z3-y3*y4^3*z3+y1*y4^2*z1*z3+y3*y4^2*z1*z3+y4^3*z1*z3-y4^2*z1^2*z3+y1*y4^2*z2*z3+
+y3*y4^2*z2*z3+y4^3*z2*z3-y1*y4*z1*z2*z3-y3*y4*z1*z2*z3-2*y4^2*z1*z2*z3+y4*z1^2*z2*z3-y4^2*z2^2*z3+
+y4*z1*z2^2*z3-y1*y4^3*z4-y3*y4^3*z4+y1*y4^2*z1*z4+y3*y4^2*z1*z4+y4^3*z1*z4-y4^2*z1^2*z4+y1*y4^2*z2*z4+
+y3*y4^2*z2*z4+y4^3*z2*z4-y1*y4*z1*z2*z4-y3*y4*z1*z2*z4-2*y4^2*z1*z2*z4+y4*z1^2*z2*z4-y4^2*z2^2*z4+
+y4*z1*z2^2*z4+y1*y4^2*z3*z4+y3*y4^2*z3*z4-y1*y4*z1*z3*z4-y3*y4*z1*z3*z4-y4^2*z1*z3*z4+y4*z1^2*z3*z4-
+y1*y4*z2*z3*z4-y3*y4*z2*z3*z4-y4^2*z2*z3*z4+y1*z1*z2*z3*z4+y3*z1*z2*z3*z4+2*y4*z1*z2*z3*z4-z1^2*z2*z3*z4+
+y4*z2^2*z3*z4-z1*z2^2*z3*z4)*QSx([4, 1, 3, 2], 'y') + (y4^4-y4^3*z1-y4^3*z2+y4^2*z1*z2-y4^3*z3+
+y4^2*z1*z3+y4^2*z2*z3-y4*z1*z2*z3-y4^3*z4+y4^2*z1*z4+y4^2*z2*z4-y4*z1*z2*z4+y4^2*z3*z4-y4*z1*z3*z4-
+y4*z2*z3*z4+z1*z2*z3*z4)*QSx([4, 1, 5, 2, 3], 'y') + (y4^4-y4^3*z1-y4^3*z2+y4^2*z1*z2-y4^3*z3+y4^2*z1*z3+
+y4^2*z2*z3-y4*z1*z2*z3-y4^3*z4+y4^2*z1*z4+y4^2*z2*z4-y4*z1*z2*z4+y4^2*z3*z4-y4*z1*z3*z4-y4*z2*z3*z4+
+z1*z2*z3*z4)*QSx([4, 2, 3, 1], 'y') + (q_1*y1+q_1*y3+q_1*y4+q_1*y5-q_1*z1-q_1*z2-q_1*z3-q_1*z4)*QSx([4, 
+5, 1, 2, 3], 'y') + q_1*QSx([4, 5, 2, 1, 3], 'y') + q_1*QSx([4, 6, 1, 2, 3, 5], 'y') + (q_3*y4^3+
+q_3*y4^2*y5+q_3*y4*y5^2+q_3*y5^3-q_3*y4^2*z1-q_3*y4*y5*z1-q_3*y5^2*z1-q_3*y4^2*z2-q_3*y4*y5*z2-
+q_3*y5^2*z2+q_3*y4*z1*z2+q_3*y5*z1*z2-q_3*y4^2*z3-q_3*y4*y5*z3-q_3*y5^2*z3+q_3*y4*z1*z3+q_3*y5*z1*z3+
+q_3*y4*z2*z3+q_3*y5*z2*z3-q_3*z1*z2*z3-q_3*y4^2*z4-q_3*y4*y5*z4-q_3*y5^2*z4+q_3*y4*z1*z4+q_3*y5*z1*z4+
+q_3*y4*z2*z4+q_3*y5*z2*z4-q_3*z1*z2*z4+q_3*y4*z3*z4+q_3*y5*z3*z4-q_3*z1*z3*z4-q_3*z2*z3*z4)*QSx([5, 1, 
+2, 3, 4], 'y') + (y1*y4^3+y3*y4^3+y1*y4^2*y5+y3*y4^2*y5+y1*y4*y5^2+y3*y4*y5^2+y1*y5^3+y3*y5^3-y1*y4^2*z1-
+y3*y4^2*z1-y4^3*z1-y1*y4*y5*z1-y3*y4*y5*z1-y4^2*y5*z1-y1*y5^2*z1-y3*y5^2*z1-y4*y5^2*z1-y5^3*z1+y4^2*z1^2+
+y4*y5*z1^2+y5^2*z1^2-y1*y4^2*z2-y3*y4^2*z2-y4^3*z2-y1*y4*y5*z2-y3*y4*y5*z2-y4^2*y5*z2-y1*y5^2*z2-
+y3*y5^2*z2-y4*y5^2*z2-y5^3*z2+y1*y4*z1*z2+y3*y4*z1*z2+2*y4^2*z1*z2+y1*y5*z1*z2+y3*y5*z1*z2+2*y4*y5*z1*z2+
+2*y5^2*z1*z2-y4*z1^2*z2-y5*z1^2*z2+y4^2*z2^2+y4*y5*z2^2+y5^2*z2^2-y4*z1*z2^2-y5*z1*z2^2-y1*y4^2*z3-
+y3*y4^2*z3-y1*y4*y5*z3-y3*y4*y5*z3-y1*y5^2*z3-y3*y5^2*z3+y1*y4*z1*z3+y3*y4*z1*z3+y4^2*z1*z3+y1*y5*z1*z3+
+y3*y5*z1*z3+y4*y5*z1*z3+y5^2*z1*z3-y4*z1^2*z3-y5*z1^2*z3+y1*y4*z2*z3+y3*y4*z2*z3+y4^2*z2*z3+y1*y5*z2*z3+
+y3*y5*z2*z3+y4*y5*z2*z3+y5^2*z2*z3-y1*z1*z2*z3-y3*z1*z2*z3-2*y4*z1*z2*z3-2*y5*z1*z2*z3+z1^2*z2*z3-
+y4*z2^2*z3-y5*z2^2*z3+z1*z2^2*z3-y1*y4^2*z4-y3*y4^2*z4-y1*y4*y5*z4-y3*y4*y5*z4-y1*y5^2*z4-y3*y5^2*z4+
+y1*y4*z1*z4+y3*y4*z1*z4+y4^2*z1*z4+y1*y5*z1*z4+y3*y5*z1*z4+y4*y5*z1*z4+y5^2*z1*z4-y4*z1^2*z4-y5*z1^2*z4+
+y1*y4*z2*z4+y3*y4*z2*z4+y4^2*z2*z4+y1*y5*z2*z4+y3*y5*z2*z4+y4*y5*z2*z4+y5^2*z2*z4-y1*z1*z2*z4-
+y3*z1*z2*z4-2*y4*z1*z2*z4-2*y5*z1*z2*z4+z1^2*z2*z4-y4*z2^2*z4-y5*z2^2*z4+z1*z2^2*z4+y1*y4*z3*z4+
+y3*y4*z3*z4+y1*y5*z3*z4+y3*y5*z3*z4-y1*z1*z3*z4-y3*z1*z3*z4-y4*z1*z3*z4-y5*z1*z3*z4+z1^2*z3*z4-
+y1*z2*z3*z4-y3*z2*z3*z4-y4*z2*z3*z4-y5*z2*z3*z4+2*z1*z2*z3*z4+z2^2*z3*z4)*QSx([5, 1, 3, 2, 4], 'y') + 
+(y4^3+y4^2*y5+y4*y5^2+y5^3-y4^2*z1-y4*y5*z1-y5^2*z1-y4^2*z2-y4*y5*z2-y5^2*z2+y4*z1*z2+y5*z1*z2-y4^2*z3-
+y4*y5*z3-y5^2*z3+y4*z1*z3+y5*z1*z3+y4*z2*z3+y5*z2*z3-z1*z2*z3-y4^2*z4-y4*y5*z4-y5^2*z4+y4*z1*z4+y5*z1*z4+
+y4*z2*z4+y5*z2*z4-z1*z2*z4+y4*z3*z4+y5*z3*z4-z1*z3*z4-z2*z3*z4)*QSx([5, 1, 4, 2, 3], 'y') + (y4^3+
+y4^2*y5+y4*y5^2+y5^3-y4^2*z1-y4*y5*z1-y5^2*z1-y4^2*z2-y4*y5*z2-y5^2*z2+y4*z1*z2+y5*z1*z2-y4^2*z3-
+y4*y5*z3-y5^2*z3+y4*z1*z3+y5*z1*z3+y4*z2*z3+y5*z2*z3-z1*z2*z3-y4^2*z4-y4*y5*z4-y5^2*z4+y4*z1*z4+y5*z1*z4+
+y4*z2*z4+y5*z2*z4-z1*z2*z4+y4*z3*z4+y5*z3*z4-z1*z3*z4-z2*z3*z4)*QSx([5, 2, 3, 1, 4], 'y') + (q_3*y4^2+
+q_3*y4*y5+q_3*y5^2+q_3*y4*y6+q_3*y5*y6+q_3*y6^2-q_3*y4*z1-q_3*y5*z1-q_3*y6*z1-q_3*y4*z2-q_3*y5*z2-
+q_3*y6*z2+q_3*z1*z2-q_3*y4*z3-q_3*y5*z3-q_3*y6*z3+q_3*z1*z3+q_3*z2*z3-q_3*y4*z4-q_3*y5*z4-q_3*y6*z4+
+q_3*z1*z4+q_3*z2*z4+q_3*z3*z4)*QSx([6, 1, 2, 3, 4, 5], 'y') + (y1*y4^2+y3*y4^2+y1*y4*y5+y3*y4*y5+y1*y5^2+
+y3*y5^2+y1*y4*y6+y3*y4*y6+y1*y5*y6+y3*y5*y6+y1*y6^2+y3*y6^2-y1*y4*z1-y3*y4*z1-y4^2*z1-y1*y5*z1-y3*y5*z1-
+y4*y5*z1-y5^2*z1-y1*y6*z1-y3*y6*z1-y4*y6*z1-y5*y6*z1-y6^2*z1+y4*z1^2+y5*z1^2+y6*z1^2-y1*y4*z2-y3*y4*z2-
+y4^2*z2-y1*y5*z2-y3*y5*z2-y4*y5*z2-y5^2*z2-y1*y6*z2-y3*y6*z2-y4*y6*z2-y5*y6*z2-y6^2*z2+y1*z1*z2+y3*z1*z2+
+2*y4*z1*z2+2*y5*z1*z2+2*y6*z1*z2-z1^2*z2+y4*z2^2+y5*z2^2+y6*z2^2-z1*z2^2-y1*y4*z3-y3*y4*z3-y1*y5*z3-
+y3*y5*z3-y1*y6*z3-y3*y6*z3+y1*z1*z3+y3*z1*z3+y4*z1*z3+y5*z1*z3+y6*z1*z3-z1^2*z3+y1*z2*z3+y3*z2*z3+
+y4*z2*z3+y5*z2*z3+y6*z2*z3-2*z1*z2*z3-z2^2*z3-y1*y4*z4-y3*y4*z4-y1*y5*z4-y3*y5*z4-y1*y6*z4-y3*y6*z4+
+y1*z1*z4+y3*z1*z4+y4*z1*z4+y5*z1*z4+y6*z1*z4-z1^2*z4+y1*z2*z4+y3*z2*z4+y4*z2*z4+y5*z2*z4+y6*z2*z4-
+2*z1*z2*z4-z2^2*z4+y1*z3*z4+y3*z3*z4-z1*z3*z4-z2*z3*z4)*QSx([6, 1, 3, 2, 4, 5], 'y') + (y4^2+y4*y5+y5^2+
+y4*y6+y5*y6+y6^2-y4*z1-y5*z1-y6*z1-y4*z2-y5*z2-y6*z2+z1*z2-y4*z3-y5*z3-y6*z3+z1*z3+z2*z3-y4*z4-y5*z4-
+y6*z4+z1*z4+z2*z4+z3*z4)*QSx([6, 1, 4, 2, 3, 5], 'y') + (y4^2+y4*y5+y5^2+y4*y6+y5*y6+y6^2-y4*z1-y5*z1-
+y6*z1-y4*z2-y5*z2-y6*z2+z1*z2-y4*z3-y5*z3-y6*z3+z1*z3+z2*z3-y4*z4-y5*z4-y6*z4+z1*z4+z2*z4+z3*z4)*QSx([6, 
+2, 3, 1, 4, 5], 'y') + (q_3*y4+q_3*y5+q_3*y6+q_3*y7-q_3*z1-q_3*z2-q_3*z3-q_3*z4)*QSx([7, 1, 2, 3, 4, 5, 
+6], 'y') + (y1*y4+y3*y4+y1*y5+y3*y5+y1*y6+y3*y6+y1*y7+y3*y7-y1*z1-y3*z1-y4*z1-y5*z1-y6*z1-y7*z1+z1^2-
+y1*z2-y3*z2-y4*z2-y5*z2-y6*z2-y7*z2+2*z1*z2+z2^2-y1*z3-y3*z3+z1*z3+z2*z3-y1*z4-y3*z4+z1*z4+
+z2*z4)*QSx([7, 1, 3, 2, 4, 5, 6], 'y') + (y4+y5+y6+y7-z1-z2-z3-z4)*QSx([7, 1, 4, 2, 3, 5, 6], 'y') + (y4+
+y5+y6+y7-z1-z2-z3-z4)*QSx([7, 2, 3, 1, 4, 5, 6], 'y') + q_3*QSx([8, 1, 2, 3, 4, 5, 6, 7], 'y') + (y1+y3-
+z1-z2)*QSx([8, 1, 3, 2, 4, 5, 6, 7], 'y') + QSx([8, 1, 4, 2, 3, 5, 6, 7], 'y') + QSx([8, 2, 3, 1, 4, 5, 
+6, 7], 'y')
+```
+This output of roughly 17,000 characters took about 2 seconds to compute on my laptop. **Again note that quantum double computations are technically conjectural, but a proof is likely forthcoming.**
+
+[Homepage of schubmult](http://schubmult.org/)

schubmult-2.0.0.dist-info/RECORD

@@ -0,0 +1,36 @@
+schubmult/__init__.py,sha256=jfMYUkKiruYd9cM9tY-_mLygW0fRj75GOhW2tEhKhOI,21
+schubmult/_base_argparse.py,sha256=Gpr2rKcjmXH4Ufl84e0s7mX-ezDrMUZSUiZ5PPg2ndY,5523
+schubmult/perm_lib.py,sha256=Q2Z_8fIDzBE2Y-f9BMEqmxyxRWw-Auu2vuyjO02v5-k,31815
+schubmult/sage_integration/__init__.py,sha256=_EeNYfhM92qVbabAPJ9jdwedhD750D--EpIhWzkAP1o,797
+schubmult/sage_integration/_fast_double_schubert_polynomial_ring.py,sha256=mtlcjrqrOuJm0jXPlkDL49Kt3y4zalBZ7YcSgWBtsUw,20067
+schubmult/sage_integration/_fast_schubert_polynomial_ring.py,sha256=m2CY9mp8WLdcvP8byoeud_V1piXt61vL_zsWRqiuni8,12836
+schubmult/sage_integration/_indexing.py,sha256=29zYqbylKJZz9XqYGfJXFfeaYa3ifXtxpEQIvAQEfI4,1576
+schubmult/schubmult_double/__init__.py,sha256=2-mZz_moyqcq1ncx7_YB6FNtztNycHZvGiJjtRlLPiE,256
+schubmult/schubmult_double/__main__.py,sha256=qhN6z5Evi4PakpzTM4wXHsYCSwfA0T-M2RYDzVcmNdE,91
+schubmult/schubmult_double/_funcs.py,sha256=XMf3b1RZEus10zM0zE1gQU8iHJIreYeAhEmNCVTVtds,51551
+schubmult/schubmult_double/_script.py,sha256=tfkKNREFUkaIkMiSHTf4Op_PIl4wsv6ZzMmgKfrpPlM,14808
+schubmult/schubmult_double/_vars.py,sha256=vGd7uBcGGBgFGQ7AeDqzEm8kXzFXiM84Mx8jraHqKX0,263
+schubmult/schubmult_py/__init__.py,sha256=WCH_afeTqjjozAFkqJDyO1t62U0PRXsZF8r8JKNjs7s,102
+schubmult/schubmult_py/__main__.py,sha256=qhN6z5Evi4PakpzTM4wXHsYCSwfA0T-M2RYDzVcmNdE,91
+schubmult/schubmult_py/_funcs.py,sha256=G9nw3GNLIiY09Xq6L0nZVOwWx9C8q1FN2-7jhn-3_0M,3440
+schubmult/schubmult_py/_script.py,sha256=w41UoMb5UXX65b45l45QlC3q5-wmyjzbWVjTiUIpirk,3919
+schubmult/schubmult_py/_vars.py,sha256=10b7uZhrKg65eMdInI_74CMnhGOXMxInxhMqKbp_hco,68
+schubmult/schubmult_q/__init__.py,sha256=9RA1-qRIe1A5EumnXTygFrasdeJZHh9o-WESZJ9X3tw,140
+schubmult/schubmult_q/__main__.py,sha256=qhN6z5Evi4PakpzTM4wXHsYCSwfA0T-M2RYDzVcmNdE,91
+schubmult/schubmult_q/_funcs.py,sha256=RQZG-zKQtpjWmlxcSmvUzOBA6Vh0zWPq4muDqc2uTQE,10453
+schubmult/schubmult_q/_script.py,sha256=lRdScEcdSFPa8oJvxJkkfpnXPLTL3X6acbuGaRhvEEg,4870
+schubmult/schubmult_q/_vars.py,sha256=dU1f9CgMersgRpP9dld7ViRkyLA76sqx7ejq0JM1ECs,295
+schubmult/schubmult_q_double/__init__.py,sha256=lkRcgbX8kFCaOB2hBmzH5n4DUNPRN0ZClJITikoVgCU,206
+schubmult/schubmult_q_double/__main__.py,sha256=qhN6z5Evi4PakpzTM4wXHsYCSwfA0T-M2RYDzVcmNdE,91
+schubmult/schubmult_q_double/_funcs.py,sha256=cmxKSCkmBAan0W9FNXJA1NIZiK4iC7mrQJjDMXv8gbM,18388
+schubmult/schubmult_q_double/_script.py,sha256=Vntga-N_0x57gd8IW-Gh5sAs_y3WMnsKi1l5Mphr5bI,14419
+schubmult/schubmult_q_double/_vars.py,sha256=JPL4kUPRtxkBJU9FkjA4mThNMbNVev5w_SmTcEsGvoM,329
+schubmult-2.0.0.dist-info/licenses/LICENSE,sha256=OXLcl0T2SZ8Pmy2_dmlvKuetivmyPd5m1q-Gyd-zaYY,35149
+tests/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+tests/test_fast_double_schubert.py,sha256=TGzjZ4kbsGpRncwHI-KNQlwOVo76ESQdCWuQxM3z0v4,4980
+tests/test_fast_schubert.py,sha256=cF1Wckg-FBycCxv0KYASIpAPQ7AOYqmtGJl9Q_puKXA,1315
+schubmult-2.0.0.dist-info/METADATA,sha256=PVM_5KxFdoEJUw8sxIp7aTtxKKh2pLrWOH4vRziu5nI,32373
+schubmult-2.0.0.dist-info/WHEEL,sha256=1tXe9gY0PYatrMPMDd6jXqjfpz_B-Wqm32CPfRC58XU,91
+schubmult-2.0.0.dist-info/entry_points.txt,sha256=FN_a2Ct5ZWBWWjST2GKx1h-VjU_c0UIkpnNeyevHPuo,244
+schubmult-2.0.0.dist-info/top_level.txt,sha256=zYdC9Gl5JX_FHOKDElw9g21WKiJTL2-cGBAXbC92Lo8,16
+schubmult-2.0.0.dist-info/RECORD,,

schubmult-2.0.0.dist-info/WHEEL

@@ -0,0 +1,5 @@
+Wheel-Version: 1.0
+Generator: setuptools (77.0.3)
+Root-Is-Purelib: true
+Tag: py3-none-any
+

schubmult-2.0.0.dist-info/entry_points.txt

@@ -0,0 +1,5 @@
+[console_scripts]
+schubmult_double = schubmult.schubmult_double.__main__:main
+schubmult_py = schubmult.schubmult_py.__main__:main
+schubmult_q = schubmult.schubmult_q.__main__:main
+schubmult_q_double = schubmult.schubmult_q_double.__main__:main