Releases: sagemath/sage
10.2
SageMath version 10.2 (released Dec 3, 2023)
Release Notes
- major package upgrades, Cython 3
- removes support for Python 3.8 and GCC older than 8.4
What's Changed
- upgrade msolve to 0.5.0, update interface by @mezzarobba in #36103
- using mute variables in combinat by @fchapoton in #36107
- some details in quadratic forms by @fchapoton in #36108
- using
# needs
in posets folder by @fchapoton in #36114 - fix cardinality, #36119, #36118, #36116 by @mantepse in #36124
- python-style loops in
partn_ref.data_structures
(partial) + some pep8 by @fchapoton in #36140 - Drop support for GCC < 8.4, drop testing of
debian-buster
andfedora-29
by @mkoeppe in #35810 sage.categories
: Update# needs
by @mkoeppe in #35945- Initial implementation of dynamical semigroups by @MatheMagicianPi in #35947
- Implemented five new operations for dynamical semigroups by @MatheMagicianPi in #35988
- Working MCDS implementation using connectivity constraints by @MrBanananator in #35990
- Implement Goss polynomials of Drinfeld modules by @DavidAyotte in #35991
- New functions
is_chebyshev
andis_Lattes
for one dimensional projective dynamical systems by @guojing0 in #36000 - Fix tests with singular 4.3.2p4 by @antonio-rojas in #36018
- Implement algebra_containment from Singular (issue #34502) by @jhpalmieri in #36030
- Implemented two new methods for dynamical semigroups by @MatheMagicianPi in #36041
- Fix sqrt(sqrt(2)) memory leak in ginac numeric.cpp by @vbraun in #36046
- New GitHub workflow to continue our images on Docker Hub by @soehms in #36047
- normaliz: Update to 3.10.1 by @mkoeppe in #36050
- cython-lint : add note about unused imports by @fchapoton in #36053
sage.arith
: Update# needs
by @mkoeppe in #36057- Speed up the creation of submatrices of
Matrix_modn_dense_template
matrices by @marizee in #36059 sage.{cpython,data_structures,databases,ext,parallel,structure}
: Update# needs
by @mkoeppe in #36067sage.{matrix,modules,probability,stats,quadratic_forms,tensor}
: Update# needs
, modularization fixes by @mkoeppe in #36070- Correct typo in class name of
Polynomial_ZZ_pEX
in.pxd
file by @GiacomoPope in #36071 sage.{interacts,misc,repl,typeset}
: Update# needs
, modularization fixes by @mkoeppe in #36078- fix and activate pycodestyle E301 in pyx files by @fchapoton in #36079
- build/pkgs/jupyter_core: Update to 4.12.0 by @mkoeppe in #36083
- some details in coxeter groups by @fchapoton in #36085
- using python-style for-loops in modular folder by @fchapoton in #36086
- using python-style for loops in calculus by @fchapoton in #36087
sage.geometry.polyhedron
: fix doctest error due to modularization by @dcoudert in #36089- Speed-up the creation of a zero matrix of type
Matrix_modn_dense_template
by @marizee in #36093 - Height function for projective subvarieties by @MatheMagicianPi in #36094
- add ordered set partitions as supported findstat collection by @mantepse in #36096
- Add additional bindings from NTL to
Polynomial_ZZ_pEX
by @GiacomoPope in #36097 - Detect and use system MathJax-3.x by @orlitzky in #36098
sage.crypto
: Update# needs
, modularization fixes by @mkoeppe in #36106- build/pkgs/cython: Update to 0.29.36 by @mkoeppe in #35084
- libatomic_ops: Update to 7.8.0 by @ivmai in #35399
- Drop support for Python 3.8 by @tobiasdiez in #35404
sage.functions
: Decouple fromsage.symbolic.expression
andmpmath
, update# needs
by @mkoeppe in #35716- gc: update to 8.2.4 by @ivmai in #35743
- Implement cocharacter and primitive Eulerian polynomials for hyperplane arrangements by @tscrim in #35914
- Delete .zenodo.json by @dimpase in #36132
- add some links to python error classes in the doc by @fchapoton in #36088
- run GitHub actions locally using gh act by @dimpase in #36117
- Bandaid for polynomial evaluation by @fchapoton in #36127
sage.combinat.root_system
: Modularization fixes,# needs
by @mkoeppe in #36136- fix Conrey number of trivial Dirichlet character by @fchapoton in #36139
- python-style loops in
sage.graphs
by @dcoudert in #36142 - Fix two "Computational Math.." book doctests by @orlitzky in #36143
- build/pkgs/mathjax/distros/gentoo.txt: Fix by @mkoeppe in #36145
- Fixing many e305 by @fchapoton in #36147
- fix the 2 broken linters by @fchapoton in #36148
- Fix doctest in
src/sage/algebras/quantum_groups/quantum_group_gap.py
by @dcoudert in #36149 - random elt for colored permutations by @fchapoton in #36155
- detect the zero series when returning an exact stream by @mantepse in #36156
- some pep for E30 and more in quadratic forms by @fchapoton in #36157
- Fix documentation previews using mathjax cdns by @kwankyu in #36160
- build/pkgs/singular: Upgrade to 4.3.2p7, reject system Singular without FLINT by @mkoeppe in #35934
- Implementing a generic one method for unital algebras by @tscrim in #36095
sage.rings.finite_rings.element_*
: Remove unused imports by @mkoeppe in #36105- Prepare for updating Cython to 3.0.0 by @infmagic2047 in #36109
- small improvements for functorial composition by @mantepse in #36115
- improve performance of stream cauchy invert by @mantepse in #35338
- Update Zariski-van Kampen functions by @enriqueartal in #35376
- using PyLong API instead of legacy PyInt API by @fchapoton in #36111
- Upgrades: scipy 1.11 and dependencies by @mkoeppe in #36112
- Coxeter arrangement by @fchapoton in #36151
sage.rings
: Update# needs
by @mkoeppe in #36152- minor details in graph_plot by @fchapoton in #36162
- Fix the ascii art for the Heisenberg Lie algebras by @tscrim in #36170
- Fix quadratic memory of berlekamp_massey by @grhkm21 in #36173
- fix and activate pycodestyle E305 in py files by @fchapoton in #36177
- Improve methods in
sage.misc.rest_index_of_methods
to fix the documentation of graphs. by @dcoudert in #36179 - more fixes in quadratic forms by @fchapoton in #36180
- fix E228 and E225 in algebras/ by @fchapoton in #36182
sage.misc.sageinspect
: fix doctest warnings due to modularisation by @dcoudert in #36183- src/sage/tests: fix another "Computational Math..." doctest by @orlitzky in #36185
- use semi-primitive root when checking kernel polynomials of isogenies by @yyyyx4 in #36187
- fix E228 and E225 in schemes/ by @fchapoton in #36191
- fix E228 then E225 in rings/ by @fchapoton in #36192
- schemes/elliptic_curves: disable a test that takes too long by @tornaria in #36195
- Fixing the iterator of SemistandardMultiSkewTableaux by @tscrim in #36200
- some micro-details in Dyck words and Parking functions by @fchapoton in #36202
- Remove self-
cimport
s by @mkoeppe in htt...
10.1
SageMath version 10.1 (released Aug 20, 2023)
Release Notes
- major package upgrades; new development tools
What's Changed
- Some fixes for Specht modules and diagrams by @tscrim in #35432
- Implement down-up algebras and their Verma modules by @tscrim in #35484
sage.rings.factorint
: Modularization fixes by @mkoeppe in #35502- Deprecate
is_SymbolicVariable
,is_SymbolicEquation
by @mkoeppe in #35505 - Some cython-linting in
modular/
folder by @fchapoton in #35516 - Make
rank
return anInteger
object by @kryzar in #35519 src/sage_docbuild/conf.py
: Do not fail ifpplpy
doc is not installed by @mkoeppe in #35536- For groups: generic
centralizer
,subgroup
methods; improvingcenter
by @tscrim in #35540 - Review of #35540 by @dwbump in #35547
- Fvector for fans by @fchapoton in #35548
- Fixing some pycodestyle warnings E251 by @fchapoton in #35549
- Removing some unused imports in cython files by @fchapoton in #35550
- Enable access to some methods in
sage/graphs/graph_decompositions
by @dcoudert in #35551 - cython-lint cleaning in
combinat/
by @fchapoton in #35553 - Make doc preview message from github bot one liner by @kwankyu in #35563
- Replace obsolescent
egrep
by @fchapoton in #35565 - Partial cython-lint in
plot
folder by @fchapoton in #35566 4ti2
: Update to 1.6.10, fixesegrep
issue, adds GCC 13 support by @mkoeppe in #35570- Some details in
posets
by @fchapoton in #35573 - Support
tox
4.x for testing modularized distributions by @mkoeppe in #35208 - Implement the logarithm and the exponential of a Drinfeld module by @DavidAyotte in #35260
- Add a few operations from linear symplectic geometry by @tobiasdiez in #35354
openblas
: Update to 0.3.23 by @mkoeppe in #35371- The Fusion Ring of the Drinfeld Double of a Finite Group by @dwbump in #35387
- Remove tag fetching in docker workflow by @tobiasdiez in #35417
sage.categories
: Add# optional
for modularization; reformat doctests by @mkoeppe in #35422- Adding bigraded Betti numbers functionality by @OP5642 in #35430
- Add ability to generate graphs based on correlations of sequences #25933 by @Bruno-TT in #35009
src/doc/en/developer/portability_testing.rst
: Update after migration by @mkoeppe in #35108- Add parallel algorithm to Graph
chromatic_number
by @dcoudert in #35145 - Some typing and code tweaks in MZV by @fchapoton in #35567
- Details on complex reflection groups by @fchapoton in #35574
- Some details in commutative dga by @fchapoton in #35575
- Simplify some trivial iterations by @fchapoton in #35576
- Refactor: outsource methods in combinatorial polyhedron by @kliem in #35577
- Refactor: improve initialization of combinatorial polyhedron by @kliem in #35578
sage.{topology,homology}
: Modularization fixes by @mkoeppe in #35581- Fix discrepancy in partitions between provided number and starting partition by @trevorkarn in #35582
- Support
networkx
3.1 by @tornaria in #35584 sage.geometry
: Add some# optional
, reformat doctests by @mkoeppe in #35586bootstrap-conda
: Skip the dummy package_sagemath
by @mkoeppe in #35589- Implement twisted affine Lie algebras by @tscrim in #35591
- Fix Max=Min warning in coin interacts by @tobiasdiez in #35596
- Remove global
sage.all
import in pytest by @tobiasdiez in #35598 src/doc/common/python3.inv
: Update by @mkoeppe in #35608- Check for empty string before trying to evaluate in
_element_constructor
by @nbruin in #35610 - Support
linbox
1.7.0 and 1.6.3 at the same time by @tornaria in #35612 - Upgrade Maxima to 5.46.0 by @dimpase in #35619
Animation.ffmpeg
: Shell-quote filenames by @mkoeppe in #35622- cython-linting of unused imports in
modules
anddynamics
by @fchapoton in #35623 - Remove jupyter-packaging for conda by @tobiasdiez in #35627
- Implement the octonions by @tscrim in #35630
- Simplify some
all(all(...))
incombinat
andgeometry
by @fchapoton in #35631 - Use fixed random-seed in the "Build & Test" CI by @tobiasdiez in #35632
- Activating more pycodestyle checks by @fchapoton in #35634
- Support
sympy
1.12 by @tornaria in #35635 - Fix typo in documentation by @antonmosich in #35641
- Avoid redundant computations in finite field
.conjugate()
method by @remyoudompheng in #35643 sage.graphs
: More# optional
and other modularization fixes by @mkoeppe in #35647- pkg-config should be
_bootstrap
, not_prereq
, on all distros by @mkoeppe in #35407 sage.rings
: Reformat doctests, add# optional
annotations by @mkoeppe in #35457- Clear up left right confusion in
ore_polynomial_element.pyx
by @kwankyu in #35562 - Make
EllipticCurve.lift_x()
deterministic by @yyyyx4 in #35644 - Remove abuse of predefined
x
in doctests by @mkoeppe in #35648 - Fix bug caused by cached representation in cohomology classes by @miguelmarco in #35650
- Installation guide: Document workaround for Jupyter server port clash on WSL by @mkoeppe in #35651
- Cleanup and details in
finite_word.py
by @fchapoton in #35657 - Fix build with
sphinx
7 by @antonio-rojas in #35658 - CI Linux: Update platforms by @mkoeppe in #35660
- SymPy 1.12 upgrade by @mkoeppe in #35662
- Changing the knot theory PD-code convention by @soehms in #35665
- cython-lint: removing unused imports in
polynomials
by @fchapoton in #35672 - cylint: remove unused imports in
pyx
files in foldersa*-geo*
by @fchapoton in #35673 - cylint: remove unused imports in
graphs
,groups
,interacts
by @fchapoton in #35674 - cylint: remove unused imports in
matrix
folder by @fchapoton in #35675 - Some pep8 in elliptic curves by @fchapoton in #35677
- Fix building the Singular docstring dict when Singular info is built with recent texinfo by @antonio-rojas in #35678
- Register magic modes for lazy imported interfaces by @miguelmarco in #35680
- Deprecation warning for #34880 by @yyyyx4 in #35683
- Unify input of matrices and vectors by @dkrenn in #35685
- Fixing a lot of pycostyle warnings in
calculus
folder by @fchapoton in #35686 - Some fixes for pycodestyle E275 by @fchapoton in #35687
- cython-lint: removed many unused imports in
rings/
by @fchapoton in #35688 - cylint: remove unused imports in folders
[s-z]*
by @fchapoton in #35689 - Correct spelling in
permgroup_named
: Diyclic => Dicyclic by @dwbump in #35694 - Fix incorrect linear branch number for non-invertible SBox by @rusydi in #35695
- cylint: remove imports in
misc/
by @fchapoton in #35697 - Fix corner case of ordered set partitions iteration by @videlec in #35700
- GAP: fix definition...
10.0
SageMath version 10.0 (released May 20, 2023)
Release Notes
- adds support for GCC 13; major package upgrades
- new optimization solvers via CVXpy, major new features in combinatorics and algebra
What's Changed
- Fix deployment of docs to netlify by @tobiasdiez in #34984
- Fix W391 linter issues by @tobiasdiez in #35109
tox.ini
(fedora-33): Do not setIGNORE_MISSING_SYSTEM_PACKAGES=no
by @mkoeppe in #35126- Mark unstable tests in
klyachko.py
due to #32773 by @tornaria in #35058 - ci-linux: Remove unmaintained local-conda runs by @mkoeppe in #35125
- Fix Cygwin CI after #32841 by @mkoeppe in #35052
- Add default devcontainer using conda by @tobiasdiez in #35072
- Improve camera positioning for threejs by @kwankyu in #34960
- Show explicit formulas in documentation of
WeierstrassIsomorphism
by @yyyyx4 in #34967 - Compute matrix kernels modulo composites by @yyyyx4 in #34968
- Cross-link matrix methods
image
andcolumn_space
by @Sandstorm831 in #34972 SignedPermutation
should allow iterables as input by @Sandstorm831 in #34974- Add
.torsion_basis()
method toEllipticCurve_finite_field
by @yyyyx4 in #34982 - Add construction of strongly regular digraph by @MatteoCati in #34986
- Trac role to GitHub by @roed314 in #35015
- Add doctest for #20847 by @mezzarobba in #34963
- Add
.twists()
toEllipticCurve_finite_field
by @yyyyx4 in #34981 - Add new skew Hadamard matrices by @MatteoCati in #34985
- Fix doctests to support numpy 1.24 by @tornaria in #34994
- Fix bug due to UB in conversion from python
int
toZZ
(python 3.11, 32 bit, gcc12) by @tornaria in #34997 - Enable dependabot for github action updates by @tobiasdiez in #34999
- Add doctest for #20846 by @mezzarobba in #35001
- Add doctest for #16031 by @mezzarobba in #35003
- Speed up some methods for quasimodular forms ring elements by @DavidAyotte in #35004
- Fix some "cannot" by @fchapoton in #35007
- Add doctest for #13569 by @mezzarobba in #35010
sage.{categories,matrix,structure}
: Replace imports fromsage.*.all
for namespace packages by @mkoeppe in #35012- Add
limit=
argument toInteger.prime_divisors()
by @yyyyx4 in #35019 - Use NTL's
MinPolyMod()
for more rings by @yyyyx4 in #34966 qqbar.clear_denominators
: crude but fast alternative algorithm by @mezzarobba in #34961- Support for tachyon >= 0.99.2 by @tornaria in #34995
- Implement check for Lorentzian polynomials #28252 by @fchapoton in #35021
- Remove upstream urls from
checksums.ini
that point to trac by @roed314 in #35023 - Fix usage of
verbose
with positional argument by @alexjbest in #35024 - Implement
__getitem__
and alias weight methods for quasimodular forms by @DavidAyotte in #35025 - gitignore another (temporary) autoconf file by @alexjbest in #35027
sage.schemes
: Replace imports fromsage.*.all
for namespace packages by @mkoeppe in #35033- Make
_multiple_x_*()
methods work for all n≠0 by @yyyyx4 in #35035 - Document argument
is_open
inManifoldSubset.complement
anddifference
by @egourgoulhon in #35040 - Add Hadamard matrices up to order 1000 by @MatteoCati in #35059
- A bijectionist's toolkit by @mantepse in #35060
- Add notes about recent changes on
KnotInfo
by @soehms in #35063 - Removing some unused imports by @fchapoton in #35064
- Upgrade eclib to 20221012 by @alexjbest in #35067
- Conform to doc requirements so that the docs look more beautiful by @alexjbest in #35069
- Document sage installation method with pip by @zovelsanj in #35070
- Combinatorial polyhedron: move list of pairs to dedicated class by @kliem in #35073
- Add ABCs
CommutativePolynomial
,MPolynomial_libsingular
,InfinitePolynomial
; deprecateis_Polynomial
,is_MPolynomial
by @mkoeppe in #35076 - Remove 'docker' from the names of the Docker images published on ghcr.io by @mkoeppe in #35079
- Add
is_supergreedy()
to linear extensions by @Sandstorm831 in #34970 - Deprecate constructing number-field fractional ideals via orders'
.ideal()
method by @yyyyx4 in #34979 is_prime
for ideals uses factorization, can be VERY slow by @yyyyx4 in #34980- Integer-valued polynomial ring by @fchapoton in #34988
- Drinfeld modules by @kryzar in #35026
- Fix the method
monomials_of_degree
by @DavidAyotte in #35043 - Convert result of multivariate polynomial evaluation into correct parent by @mezzarobba in #35045
sage.{coding,groups}
: Replace imports from sage.*.all for namespace packages by @mkoeppe in #35099- Remove direct use of
setup.py sdist
, add targetsmake SPKG-sdist
by @mkoeppe in #35104 sage.{functions,interfaces,symbolic}
: Replace imports fromsage.*.all
for namespace packages by @mkoeppe in #35105sage.{arith,crypto,databases,dynamics,lfunctions,quadratic_forms}
: Replace imports fromsage.*.all
for namespace packages by @mkoeppe in #35106sage.{finance,interacts,libs,numerical,stats,tests}
: Replace imports fromsage.*.all
for namespace packages by @mkoeppe in #35107- Deprecate
is_FiniteField
etc., makesage.rings.finite_rings
a namespace package by @mkoeppe in #35119 - Fix use of
sig_on()
/sig_off()
inCBF.integral()
by @mezzarobba in #35044 - Implementing q-commuting Laurent polynomials. by @tscrim in #35054
- Upgrade scipy to 1.10.1 by @mkoeppe in #35082
build/pkgs/pcre
: Remove by @mkoeppe in #35083- Upgrade polymake to 4.9 by @mkoeppe in #35085
- Move
single_valued
method of MZV to auxiliary F ring by @fchapoton in #35086 - Use enum for face iterator status by @kliem in #35087
sage.{algebras,combinat,matroids}
: Replace imports from sage.*.all for namespace packages by @mkoeppe in #35090- Fix and activate W391 in
pyx
files by @fchapoton in #35091 - Update GAP to 4.12.2 by @dimpase in #35093
- To determine
GAP_SO
,sage.env
looks forlibgap.so
but it should look forlibgap.so*
by @dimpase in #35094 - Simplicial set group by @miguelmarco in #35097
sage.{topology,homology}
: Replace imports fromsage.*.all
for namespace packages by @mkoeppe in #35098- New ABC
sage.structure.element.NumberFieldElement
,deprecate is_NumberFieldElement
by @mkoeppe in #35100 scip_backend
: Remove use of deprecatedsage.ext.memory_allocator
by @mkoeppe in #35101- Fix all pycodestyle E303 warnings in all folders
c*
by @fchapoton in #35151 - Update msolve to 0.4.9 by @dimpase in #35159
- Fix a slow doctest in
matrix_integer_dense_hnf.py
by @tornaria in #35162 - Fix the documentation for the so matrix Lie algebra by @tscrim in #35163
- Compute the matrix of an isogeny on a given n-torsion subgroup by @yyyyx4 in #35164
- Activate W293 and E714 in
pyx
files by @fchapoton in #35165 - Fix bug in
is_eulerian
by @dcoudert in #35170 - Speedup of Poset cha...
9.8
SageMath version 9.8 (released Feb 11, 2023)
Release Notes
- first release after our move from Trac to GitHub
- adds support for Python 3.11, removes support for GCC older than 8
- major new solvers, major package upgrades
- 402 tickets/PRs merged
What's Changed
Merged in sage-9.8.beta0:
#12696: Lorenz Panny: rename .rational_reconstruct() to .rational_reconstruction() for consistency [Reviewed by Kwankyu Lee]
#17965: Frédéric Chapoton: Uniformize the API to compute the inverse of an element [Reviewed by Travis Scrimshaw]
#20467: Frédéric Chapoton: stackoverflow in Maxima integral [Reviewed by Nils Bruin]
#29619: Matthias Koeppe: Matrix and Components should have a sparse iterator [Reviewed by Eric Gourgoulhon]
#30229: Matthias Koeppe: Submodules of TensorFreeModule defined by the symmetries of a Components object [Reviewed by Eric Gourgoulhon]
#30241: Matthias Koeppe: New implementation class FiniteRankDualFreeModule [Reviewed by Eric Gourgoulhon]
#32297: Linden Disney-Hogg: Abel-Jacobi map on Riemann surfaces [Reviewed by Nils Bruin, Travis Scrimshaw]
#32669: Dennis Jahn: Adding upper and lower Bruhat cones of M. Dyer to sage/combinat/root_system/reflection_group_real.py [Reviewed by Frédéric Chapoton, Travis Scrimshaw]
#32921: Gabriel F. Lipnik: k-recursive sequences with inhomogeneities [Reviewed by Clemens Heuberger]
#33002: Sébastien Labbé: Method tikz of polyhedron class can now return an object of type TikzPicture [Reviewed by Laith Rastanawi]
#33062: Matthias Koeppe: GH Actions (docker): Run a job for "make build-local" first, cache image for job "make build" [Reviewed by Dima Pasechnik]
#33093: John Palmieri: 5 doctests failing related to external program octave [Reviewed by Frédéric Chapoton]
#33431: Matthias Koeppe: Category of submodules of sage.modules.free_module.FreeModule_ambient_pid [Reviewed by Travis Scrimshaw]
#33613: Matthias Koeppe: Add pip packages to conda environments, add missing conda.txt files [Reviewed by Dima Pasechnik]
#33640: Frédéric Chapoton: sage fails to factor some easy expressions [Reviewed by David Lowry-Duda]
#33981: Matthias Koeppe: latte_int: ExponentialSubst.h:62:3: error: ISO C++17 does not allow dynamic exception specifications [Reviewed by John Palmieri]
#34118: Oscar Benjamin, Matthias Koeppe: Update sympy to 1.11.1 [Reviewed by Matthias Koeppe, Antonio Rojas, Eric Gourgoulhon]
#34195: Matthias Koeppe: sage.geometry.polyhedron: More # optional - sage.rings.number_field [Reviewed by Jonathan Kliem]
#34203: Matthias Koeppe: Find and uninstall broken installed SPKGs and wheels [Reviewed by John Palmieri]
#34228: Matthias Koeppe: tox -e docker-...-incremental [Reviewed by Dima Pasechnik]
#34235: Matthias Koeppe: Error in comparing two Symbolic Ring elements [Reviewed by Travis Scrimshaw]
#34242: Julian Rüth: Update docker build [Reviewed by MartinPepin]
#34266: Matthias Koeppe: Drop support for GCC < 8 in Sage 9.8 [Reviewed by Dima Pasechnik]
#34282: Sebastian Oehms, Kwankyu Lee: Make feature TeXFile check latex first [Reviewed by Kwankyu Lee, Sebastian Oehms]
#34287: Frédéric Chapoton: modernize super in crystals and asymptotics [Reviewed by Travis Scrimshaw]
#34290: Frédéric Chapoton: fix some W391 [Reviewed by Travis Scrimshaw]
#34313: David Coudert: pycodestyle cleanup in strongly_regular_db.pyx (part 3) [Reviewed by Kwankyu Lee]
#34320: Frédéric Chapoton: fix E251 in schemes [Reviewed by David Coudert]
#34323: Frédéric Chapoton: fix E251 in groups [Reviewed by David Coudert]
#34342: Frédéric Chapoton: fix various suggestions by lgtm [Reviewed by David Coudert]
#34345: Trevor K. Karn: Method to get cells in horizontal/vertical border strip [Reviewed by Travis Scrimshaw]
#34350: Travis Scrimshaw: speed up powers of lazy Taylor series [Reviewed by Martin Pepin]
#34379: Travis Scrimshaw: Implement a hook to access free (graded) resolutions [Reviewed by Kwankyu Lee]
#34393: Frédéric Chapoton: add method "tensor_factors" to tensor products [Reviewed by Matthias Koeppe, Travis Scrimshaw]
#34412: Travis Scrimshaw: q-commuting polynomials [Reviewed by Frédéric Chapoton]
#34427: Matthias Koeppe: TensorFreeModule.isomorphism_with_fixed_basis [Reviewed by Eric Gourgoulhon]
#34449: Matthias Koeppe: Add package pyproject_metadata [Reviewed by François Bissey]
#34451: Matthias Koeppe: sage.tensor: Canonicalize sym, antisym [Reviewed by Eric Gourgoulhon]
#34453: Travis Scrimshaw: Tensor products of commutative algebras do not know they are commutative rings [Reviewed by Frédéric Chapoton]
#34466: Lorenz Panny: fix various linter errors [Reviewed by Frédéric Chapoton]
#34467: Lorenz Panny: fix random doctest failure in EllipticCurveHom_velusqrt [Reviewed by Kwankyu Lee]
#34474: Eric Gourgoulhon: Make FiniteRankFreeModule.tensor_module(0, 1) return the dual [Reviewed by Matthias Koeppe]
#34478: Matthias Koeppe: pynormaliz: Add distros information for conda, arch [Reviewed by Isuru Fernando]
#34484: Matthias Koeppe: Methods quotient, quotient_module are not the same for some modules from sage.modules [Reviewed by John Palmieri, Travis Scrimshaw]
#34490: Thierry Monteil: sage --package clean
command to remove outdated source tarballs [Reviewed by David Coudert]
#34491: Matthias Koeppe: Upgrade igraph to 0.9.10 [Reviewed by David Coudert]
#34494: Frédéric Chapoton: better monomials in shuffle algebras [Reviewed by Travis Scrimshaw]
#34495: Matthias Koeppe: Construction for invariant/equivariant submodules [Reviewed by Travis Scrimshaw]
#34498: Matthias Koeppe: Update igraph, python_igraph to 0.10.x [Reviewed by David Coudert]
#34505: Frédéric Chapoton: using items in indexed_elements.pyx [Reviewed by Travis Scrimshaw]
#34509: Matthias Koeppe: Make IndexedFreeModuleElement compatible with collections.abc, change method support to return a SupportView [Reviewed by Travis Scrimshaw]
#34510: Trevor K. Karn: Add to check for fixed total/length in IntegerVector [Reviewed by Travis Scrimshaw]
#34514: David Roe: Incorrect decomposition type for relative number fields [Reviewed by Travis Scrimshaw]
#34519: Marc Mezzarobba: Fix + extend msolve interface [Reviewed by Matthias Koeppe]
#34527: Trevor K. Karn: Make Composition
s into a collections.abc.Sequence
[Reviewed by Travis Scrimshaw, Matthias Köppe]
#34535: Trevor K. Karn: Segfault when floor dividing in infinite polynomial ring [Reviewed by Travis Scrimshaw]
#34541: Frédéric Chapoton: enhance free dendriform algebras [Reviewed by Travis Scrimshaw]
#34561: Jan Groenewald: Documentation: There is no package in debian/ubuntu called openssh [Reviewed by Matthias Koeppe]
#34568: Frédéric Chapoton: rename example to QuaternionicProjectivePlane [Reviewed by John Palmieri]
Merged in sage-9.8.beta1:
#21129: Paul Fili, Holly Krieger, Jing Guo: Arakelov-Zhang pairing of rational maps [Reviewed by Alexander Galarraga]
#25046: Lorenz Panny: Add warning for is_prime(n) that produce correct but unexpected output [Reviewed by Kwankyu Lee]
#31664: Matthias Koeppe: Add package msolve 0.4.4 (multivariate polynomial system solver) [Reviewed by Marc Mezzarobba]
#32483: Matthias Koeppe: Remove workarounds for packages that use flit_core or poetry_core as build system [Reviewed by John Palmieri]
#33805: Dima Pasechnik: use Homebrew's primecount [Reviewed by Matthias Koeppe]
#33812: Matthias Koeppe: Refactor distributions sagemath-{objects,categories} through sagemath-{environment,repl} [Reviewed by Kwankyu Lee]
#33850: Sophia Elia: Unstable doctests in equivariant Ehrhart methods [Reviewed by Matthias Koeppe]
#33999: Xavier Caruso: Reduced charpoly for skew polynomials [Reviewed by Travis Scrimshaw, Frédéric Chapoton]
#34110: Matthias Koeppe: numpy 1.23.x [Reviewed by Dima Pasechnik, John Palmieri]
#34212: Jing Guo: log_embedding
for number field elements [Reviewed by Alexander Galarraga]
#34272: Matthias Koeppe: Update pip to 22.2.2 [Reviewed by Kwankyu Lee]
#34314: David Coudert: pycodestyle cleanup in strongly_regular_db.pyx (part 4) [Reviewed by Matthias Koeppe]
#34344: Matthias Koeppe: Update setuptools 63.4.3, setuptools_wheel 65.x, tomlkit 0.11.4, importlib_resources 5.9.0, importlib_metadata 4.12.0, distlib 0.3.6, hatchling 1.10.0, pathspec 0.10.1 [Reviewed by John Palmieri]
#34358: David Coudert: pycodestyle cleanup in src/sage/graphs/generic_graph.py (part 4) [Reviewed by Kwankyu Lee]
#34375: Travis Scrimshaw: Replace sage.algebras.yangian.GeneratorIndexingSet with cartesian_product [Reviewed by Matthias Koeppe]
#34383: Martin Rubey: revert for LazyTaylorSeries and LazySymmetricFunction is missing [Reviewed by Travis Scrimshaw]
#34392: David Coudert: pycodestyle cleanup in src/sage/graphs/generic_graph.py (part 5) [Reviewed by Matthias Koeppe]
#34409: Lorenz Panny: remove experimental warning for composite elliptic-curve isogenies [Reviewed by Kwankyu Lee]
#34413: Martin Rubey: implement derivatives of lazy series [Reviewed by Travis Scrimshaw]
#34422: Martin Rubey: implement functorial composition of lazy symmetric functions [Reviewed by Travis Scrimshaw]
#34423: Martin Rubey: implement arithmetic product of lazy symmetric functions [Reviewed by Travis Scrimshaw]
#34432: Samuel Lelièvre: Upgrade: jupyter-packaging 0.12.3 [Reviewed by Matthias Koeppe]
#34435: Trevor K. Karn: Add method to trim trailing zeros from IntegerVector [Reviewed by Travis Scrimshaw]
#34448: Matthias Koeppe: Put tensor modules of FiniteRankFreeModule in Modules().TensorProducts() [Reviewed by Eric Gourgoulhon]
#34454: Matthias Koeppe: sage --package create: Add option --source wheel [Reviewed by Kwankyu Lee]
#34462: Kwankyu Lee: Uniformize the headline: finite rings [Reviewed by Matthias Koeppe]
#34465: John Palmieri: Invalid escape sequence in special.py [Reviewed by Dima Pasechnik]
#34468: Fr...
9.7
SageMath version 9.7 (released Sep 19, 2022)
Release Notes
- removes support for GCC older than 6.3, removes support for Python 3.7
- new code in algebra, number theory, elliptic curves; major package upgrades; new developer tools
- 549 tickets (PRs) merged, 92 contributors
9.6
Release Tour
Sage 9.6 was released on May 15, 2022 (changelog)
- adds support for GCC 12; major package upgrades
- major improvements to user interface and graphics; new developer tools
- 497 tickets (PRs) merged
A total of 83 people were involved as authors or reviewers of code contributions to Sage 9.6.
Here is an overview of some of the main changes in this version.
User interface, plotting and graphics
JupyterLab 3.3
JupyterLab, the latest web-based interactive development environment for notebooks, code, and data, is slated to replace the now-classic Jupyter notebook interface. The version of JupyterLab in the Sage distribution has been upgraded to the major new version 3.3. #32069, #33607
After ./sage -i jupyterlab_widgets
, you can run it using
./sage -n jupyterlab
Also two new interface variants are provided:
./sage -n nbclassic
and
./sage -n retrolab
LaTeX displays in JupyterLab
Users of Sage in JupyterLab got used to expressions displayed at center in the LaTeX display mode. For compatibility with displays in classic Jupyter, we decided to change the behavior so that now expressions are displayed aligned left by default.
If you belong to the minority preferring centered displays, you can set your preference by
dm = get_display_manager()
dm.preferences.align_latex = 'center' # or 'left'
in the ~/.sage/init.sage
script.
Interactive graph editing with phitigra
With the new optional package phitigra (use ./sage -i phitigra
to install), graphs can be edited by interactively placing vertices, edges, etc. This works both in the classic Jupyter notebook and in JupyterLab. It can also be used to make animations (see the demo notebook at https://github.com/jfraymond/phitigra for examples). Done in #30540 and #33639.
Hyperbolic plots
- Added the ability to choose the hyperbolic model for hyperbolic plots. #22081
Graphics with TikZ
The TikzPicture
module which was developed in the slabbe package for more than 5 years is now in Sage. This was done in ticket #20343. The module is within the new file sage/misc/latex_standalone.py
and its documentation in the reference manual is available here: https://doc.sagemath.org/html/en/reference/misc/sage/misc/latex_standalone.html. Below are some usage examples.
First example shows that it takes any tikz picture string as input:
sage: from sage.misc.latex_standalone import TikzPicture
sage: s = '\\begin{tikzpicture}\n\\draw[->,green,very thick](0,0) -- (1,1);\\end{tikzpicture}'
sage: t = TikzPicture(s)
sage: t # in Jupyter, rich representation will show the image instead
\documentclass[tikz]{standalone}
\begin{document}
\begin{tikzpicture}
\draw[->,green,very thick](0,0) -- (1,1);\end{tikzpicture}
\end{document}
sage: path_to_file = t.pdf() # and opens the image in a viewer
Of course, conversion to pdf format necessitates pdflatex
or lualatex
. If lualatex
is available it uses it in preference to pdflatex
because it handles better the very big pictures in terms of memory limits.
One can provide a local filename to save to, or convert the image to other formats (using pdftocairo or imagemagick external packages):
sage: path_to_file = t.pdf('file.pdf') # when providing a filename, it just saves
# the file locally, does not open in a viewer
sage: path_to_file = t.png() # conversion to png
sage: path_to_file = t.svg() # to svg
sage: path_to_file = t.tex() # print the tex source to a file
Another example with graphs where additional usepackage are necessary to compile the image correctly:
sage: from sage.misc.latex_standalone import TikzPicture
sage: g = graphs.PetersenGraph()
sage: t = TikzPicture(latex(g), standalone_config=["border=4mm"], usepackage=['tkz-graph'])
sage: t # in Jupyter, rich representation will show the image instead
\documentclass[tikz]{standalone}
\standaloneconfig{border=4mm}
\usepackage{tkz-graph}
\begin{document}
\begin{tikzpicture}
\definecolor{cv0}{rgb}{0.0,0.0,0.0}
\definecolor{cfv0}{rgb}{1.0,1.0,1.0}
\definecolor{clv0}{rgb}{0.0,0.0,0.0}
\definecolor{cv1}{rgb}{0.0,0.0,0.0}
---
65 lines not printed (3695 characters in total).
Use print to see the full content.
---
\Edge[lw=0.1cm,style={color=cv6v8,},](v6)(v8)
\Edge[lw=0.1cm,style={color=cv6v9,},](v6)(v9)
\Edge[lw=0.1cm,style={color=cv7v9,},](v7)(v9)
%
\end{tikzpicture}
\end{document}
sage: _ = t.pdf() # or t.png() or t.svg()
sage: from sage.misc.latex_standalone import TikzPicture
sage: V = [[1,0,1],[1,0,0],[1,1,0],[0,0,-1],[0,1,0],[-1,0,0],[0,1,1],[0,0,1],[0,-1,0]]
sage: P = Polyhedron(vertices=V).polar()
sage: s = P.projection().tikz([674,108,-731],112)
sage: t = TikzPicture(s)
sage: _ = t.pdf() # or t.png() or t.svg()
The module also contains a class Standalone
, from which the class TikzPicture
inherits:
sage: from sage.misc.latex_standalone import Standalone
sage: s = Standalone('Hello World', usepackage=['amsmath'], standalone_config=['beamer=true','border=1mm'])
sage: s # in Jupyter, rich representation will show the image instead
\documentclass{standalone}
\standaloneconfig{beamer=true}
\standaloneconfig{border=1mm}
\usepackage{amsmath}
\begin{document}
Hello World
\end{document}
sage: _ = s.pdf() # or s.png() or s.svg()
Another example using Standalone
with a tableau:
sage: P = Permutations(10)
sage: p = P.random_element()
sage: p
[3, 10, 1, 9, 5, 6, 7, 2, 8, 4]
sage: t = p.to_tableau_by_shape([3,3,3,1])
sage: t
[[2, 8, 4], [5, 6, 7], [10, 1, 9], [3]]
sage: s = Standalone(latex(t), standalone_config=["border=1mm"])
sage: s
\documentclass{standalone}
\standaloneconfig{border=1mm}
\begin{document}
{\def\lr#1{\multicolumn{1}{|@{\hspace{.6ex}}c@{\hspace{.6ex}}|}{\raisebox{-.3ex}{$#1$}} }
\raisebox{-.6ex}{$\begin{array}[b]{*{3}c}\cline{1-3}
\lr{2}&\lr{8}&\lr{4}\\\cline{1-3}
\lr{5}&\lr{6}&\lr{7}\\\cline{1-3}
\lr{10}&\lr{1}&\lr{9}\\\cline{1-3}
\lr{3}\\\cline{1-1}
\end{array}$}
}
\end{document}
sage: _ = s.pdf() # or s.png() or s.svg()
In a next step, a method tikz()
will be added to graphs, polytopes, posets, etc. to return an object of type TikzPicture
see #33002.
Complex plots
The complex plotting package phase_mag_plot has been incorporated into Sage. Now complex_plot
allows contouring, tiling, and matplotlib
-compatible colormaps. This was added in ticket #33416.
To use a colormap, one can pass in a string as in
sage: complex_plot((x - 5)*sqrt(x), (-10, 10), (-10, 10), cmap='twilight')
Contouring or tiling are enabled through keyword options. To look smooth, it's typically necessary to plot the function on additional points through the use of plot_points
. This looks like
sage: complex_plot((x - 5)*sqrt(x), (-10, 10), (-10, 10), cmap='twilight', plot_points=500, contoured=True)
sage: complex_plot((x - 5)*sqrt(x), (-10, 10), (-10, 10), cmap='twilight', plot_points=500, tiled=True)
Linear algebra
NumPy integration
The new classes Matrix_numpy_integer_dense
and Vector_numpy_integer_dense
implement matrices and vectors with 64-bit integer entries backed by numpy
arrays. #32465.
As a first application, several methods of GenericGraph
that return matrices, such as adjacency_matrix
, now accept keyword arguments that can select the matrix implementation. #33377, #33387, #33388, #33389
sage: graphs.PathGraph(5).adjacency_matrix(sparse=False, implementation='numpy')
[0 1 0 0 0]
[1 0 1 0 0]
[0 1 0 1 0]
[0 0 1 0 1]
[0 0 0 1 0]
sage: type(_)
<class 'sage.matrix.matrix_numpy_integer_dense.Matrix_numpy_integer_dense'>
CombinatorialFreeModule improvements
Performing sums and similar constructions for CombinatorialFreeModule
have been made faster. #33267
#...
9.5
Release Tour
SageMath version 9.5 was released on Jan 30, 2022 (changelog)
- adds support for Apple M1, Linux distributions using glibc ≥ 2.34, and system Python 3.10
- new code from 2021 Google Summer of Code projects in algebra and number theory
- 663 tickets (PRs) merged, 94 contributors
Symbolics
Changes to symbolic expressions
symbolic_expression
is now able to create vectors and matrices of symbolic expressions for more general inputs. #16761
For example, if the input is a list or tuple of lists/tuples/vectors:
sage: M = symbolic_expression([[1, x, x^2], (x, x^2, x^3), vector([x^2, x^3, x^4])]); M
[ 1 x x^2]
[ x x^2 x^3]
[x^2 x^3 x^4]
sage: M.parent()
Full MatrixSpace of 3 by 3 dense matrices over Symbolic Ring
- Symbolic expressions can no longer be called with positional arguments. #14270
This was deprecated since Sage 4.0, although a bug prevented the deprecation warning from being issued in Sage versions 8.4 to 9.3. #32319
Instead of (x+1)(2)
, write either (x+1)(x=2)
, or (x+1).subs(x=2)
, or ((x+1).function(x))(2)
.
Interface to Mathics, a free implementation of the Wolfram language
Sage now has an optional package providing Mathics, a free (open-source) general-purpose computer algebra system featuring Mathematica-compatible syntax and functions, and an interface to it.#31778
Linear Algebra
- The Cholesky decomposition for sparse RDF/CDF matrices now uses a specialized fast algorithm when cvxopt is available. #13674
- The
is_hermitian()
method for sparse RDF/CDF matrices now uses a small tolerance by default to mitigate numerical issues. This brings it to parity with its dense counterpart. #33031
Manifolds
The full list of changes is available in this changelog.
De Rham cohomology and characteristic classes
The de Rham cohomology has been made an algebra (#32270).
The method characteristic_class
for vector bundles is now outdated and replaced by the method characteristic_cohomology_class
. This change reflects the difference between characteristic classes, seen as natural transformations, and characteristic cohomology classes in a more rigorous way (see #29581). The previous usability and syntax remains intact. Among other things, the following has been changed:
- The performance of computing characteristic forms has been improved significantly by using a Faddeev-LeVerrier-like algorithm.
- The characteristic forms of Pontryagin/Chern/Euler classes w.r.t. to a given connection are cached in order to speed up computations of all characteristic forms w.r.t. the same connection.
Furthermore, new features have been added. For example, characteristic cohomology classes now form an algebra:
sage: M = Manifold(4, 'M')
sage: E = M.vector_bundle(2, 'E', field='complex')
sage: R = E.characteristic_cohomology_class_ring(); R
Algebra of characteristic cohomology classes of the Differentiable complex vector bundle E -> M of rank 2 over the base space 4-dimensional differentiable manifold M
sage: R.gens()
[Characteristic cohomology class (c_1)(E) of the Differentiable complex vector bundle E -> M of rank 2 over the base space 4-dimensional differentiable manifold M,
Characteristic cohomology class (c_2)(E) of the Differentiable complex vector bundle E -> M of rank 2 over the base space 4-dimensional differentiable manifold M]
sage: c_1, c_2 = R.gens()
Therefore, elements can be added and multiplied:
sage: c_1 + c_2
Characteristic cohomology class (c_1 + c_2)(E) of the Differentiable complex vector bundle E -> M of rank 2 over the base space 4-dimensional differentiable manifold M
sage: c_1 * c_1
Characteristic cohomology class (c_1^2)(E) of the Differentiable complex vector bundle E -> M of rank 2 over the base space 4-dimensional differentiable manifold M
Additive, multiplicative and Pfaffian cohomology classes are now related to the generators of the characteristic cohomology class ring via additive/multiplicative sequences:
sage: ch = E.characteristic_cohomology_class('ChernChar'); ch
Characteristic cohomology class ch(E) of the Differentiable complex vector bundle E -> M of rank 2 over the base space 4-dimensional differentiable manifold M
sage: ch == 2 + c_1 + c_1^2 / 2 - c_2 # additive sequence of exp(x)
True
As for the tangent bundle of a manifold, as long as an orientation and a metric is provided, the characteristic form of the Euler class (and therefore all Pfaffian classes) w.r.t. the Levi-Civita connection is now computed automatically (previously, a compatible curvature form matrix had to be provided by the user):
sage: M.<x,y> = manifolds.Sphere(2, coordinates='stereographic')
sage: g = M.metric()
sage: nab = g.connection()
sage: nab.set_immutable()
sage: TM = M.tangent_bundle()
sage: e = TM.characteristic_cohomology_class('Euler'); e
Characteristic cohomology class e(TS^2) of the Tangent bundle TS^2 over the 2-sphere S^2 of radius 1 smoothly embedded in the Euclidean space E^3
sage: e_form = e.get_form(nab)
sage: e_form
Mixed differential form e(TS^2, nabla_g) on the 2-sphere S^2 of radius 1 smoothly embedded in the Euclidean space E^3
sage: e_form.display_expansion()
e(TS^2, nabla_g) = 2/(pi + pi*x^4 + pi*y^4 + 2*pi*x^2 + 2*(pi + pi*x^2)*y^2) dx∧dy
Automatic coordinate change in curve plots
The method plot
of curves is now allowed to perform a coordinate change to make the plot in terms of the coordinates specified via the argument chart
, in case the curve is not known in these coordinates (#32578). For example, a cardioid is defined in terms of polar coordinates:
sage: E.<r, phi> = [[EuclideanSpace]](coordinates='polar')
sage: c = E.curve((1 + cos(phi), phi), (phi, 0, 2*pi))
sage: c.display()
(0, 2*pi) → E^2
phi ↦ (r, phi) = (cos(phi) + 1, phi)
and its plot in terms of Cartesian coordinates can now be obtained simply by
sage: c.plot(chart=E.cartesian_coordinates(), aspect_ratio=1)
The above command has triggered the computation of the curve's expression in terms of Cartesian coordinates:
sage: c.display()
(0, 2*pi) → E^2
phi ↦ (r, phi) = (cos(phi) + 1, phi)
phi ↦ (x, y) = (cos(phi)^2 + cos(phi), (cos(phi) + 1)*sin(phi))
Internal code improvements and bug fixes
Various improvements have been performed in the internal code, some of them in view of SageMath modularization:
- faster generation of non-redundant indices (#32318)
- unnecessary uses of symbolic functions removed from
sage.tensor.modules
(#32415) - doctests involving
SR
marked optional insage.tensor.modules
(#32712) sage.tensor.modules
made independent fromsage.manifolds
(#32708).
Some bugs have been fixed: #31781, #32457, #32355, #32929.
Number theory
Logarithms
- Logarithms modulo composite integers are now dramatically faster in some important cases (such as prime-power moduli). #32375
- Logarithms in binary finite fields now use index calculus instead of generic-group algorithms when appropriate (via PARI's fflog()). #32842
Binary quadratic forms
- Binary quadratic forms' .solve_integer() method now uses PARI's qfbsolve() instead of a brute-force search, which is often exponentially faster. #32782
Prime counting
- Prime counting and related functions, in particular prime_pi, are now implemented using external libraries, primecount and primesieve. This improved performance and fixed a long-standing bug #24960
Elliptic curves and isogenies
- Composite isogenies of elliptic curves have been added as an experimental feature. This includes support for compact smooth-degree isogenies. #32744
- Isomorphisms of elliptic curves now expose the same interface as other isogenies. [#32502](https://trac.sagemath...
9.4
Release Tour
SageMath 9.4 was released on Aug 22, 2021 (changelog), 442 tickets (PRs) merged, 73 contributors.
- adds support for GCC 11, removes support for Python 3.6
- major advances in symbolics, convex and differential geometry, knot theory, coding theory
Symbolics
Extended interface with SymPy
The SymPy package has been updated to version 1.8.
SageMath has a bidirectional interface with SymPy. Symbolic expressions in Sage provide a _sympy_
method, which converts to SymPy; also, Sage attaches _sage_
methods to various SymPy classes, which provide the opposite conversion.
In Sage 9.4, several conversions have been added. Now there is a bidirectional interface as well for matrices and vectors. #31942
sage: M = matrix([[sin(x), cos(x)], [-cos(x), sin(x)]]); M
[ sin(x) cos(x)]
[-cos(x) sin(x)]
sage: sM = M._sympy_(); sM
Matrix([
[ sin(x), cos(x)],
[-cos(x), sin(x)]])
sage: sM.subs(x, pi/4) # computation in SymPy
Matrix([
[ sqrt(2)/2, sqrt(2)/2],
[-sqrt(2)/2, sqrt(2)/2]])
Work is underway to make SymPy's symbolic linear algebra methods available in Sage via this route.
Callable symbolic expressions, such as those created using the Sage preparser's f(...) = ...
syntax, now convert to a SymPy Lambda
. #32130
sage: f(x, y) = x^2 + y^2; f
(x, y) |--> x^2 + y^2
sage: f._sympy_()
Lambda((x, y), x**2 + y**2)
Sage has added a formal set membership function element_of for use in symbolic expressions; it converts to a SymPy's Contains
expression. #24171
Moreover, all sets and algebraic structures (Parent
s) of SageMath are now accessible to SymPy by way of a wrapper class SageSet
, which implements the SymPy Set API. #31938
sage: F = Family([2, 3, 5, 7]); F
Family (2, 3, 5, 7)
sage: sF = F._sympy_(); sF
SageSet(Family (2, 3, 5, 7)) # this is how the wrapper prints
sage: sF._sage_() is F
True # bidirectional
sage: bool(sF)
True
sage: len(sF)
4
sage: sF.is_finite_set # SymPy property
True
Finite or infinite, we can wrap it:
sage: W = WeylGroup(["A",1,1])
sage: sW = W._sympy_(); sW
SageSet(Weyl Group of type ['A', 1, 1] (as a matrix group acting on the root space))
sage: sW.is_finite_set
False
sage: sW.is_iterable
True
sage: sB3 = WeylGroup(["B", 3])._sympy_(); sB3
SageSet(Weyl Group of type ['B', 3] (as a matrix group acting on the ambient space))
sage: len(sB3)
48
Some parents or constructions have a more specific conversion to SymPy. #31931, #32015
sage: ZZ3 = cartesian_product([ZZ, ZZ, ZZ])
sage: sZZ3 = ZZ3._sympy_(); sZZ3
ProductSet(Integers, Integers, Integers)
sage: (1, 2, 3) in sZZ3
sage: NN = NonNegativeIntegers()
sage: NN._sympy_()
Naturals0
sage: (RealSet(1, 2).union(RealSet.closed(3, 4)))._sympy_()
Union(Interval.open(1, 2), Interval(3, 4))
sage: X = Set(QQ).difference(Set(ZZ)); X
Set-theoretic difference of
Set of elements of Rational Field and
Set of elements of Integer Ring
sage: X._sympy_()
Complement(Rationals, Integers)
sage: X = Set(ZZ).difference(Set(QQ)); X
Set-theoretic difference of
Set of elements of Integer Ring and
Set of elements of Rational Field
sage: X._sympy_()
EmptySet
See Meta-ticket #31926: Connect Sage sets to SymPy sets
ConditionSet
Sage 9.4 introduces a class ConditionSet for subsets of a parent (or another set) consisting of elements that satisfy the logical "and" of finitely many predicates. #32089
sage: in_small_oblong(x, y) = x^2 + 3 * y^2 <= 42
sage: SmallOblongUniverse = ConditionSet(QQ^2, in_small_oblong)
sage: SmallOblongUniverse
{ (x, y) ∈ Vector space of dimension 2 over Rational Field : x^2 + 3*y^2 <= 42 }
sage: parity_check(x, y) = abs(sin(pi/2*(x + y))) < 1/1000
sage: EvenUniverse = ConditionSet(ZZ^2, parity_check); EvenUniverse
{ (x, y) ∈ Ambient free module of rank 2 over the principal ideal
domain Integer Ring : abs(sin(1/2*pi*x + 1/2*pi*y)) < (1/1000) }
sage: SmallOblongUniverse & EvenUniverse
{ (x, y) ∈ Free module of degree 2 and rank 2 over Integer Ring
Echelon basis matrix:
[1 0]
[0 1] : x^2 + 3*y^2 <= 42, abs(sin(1/2*pi*x + 1/2*pi*y)) < (1/1000) }
The name ConditionSet
is borrowed from SymPy. In fact, if the given predicates (condition) are symbolic, a ConditionSet
can be converted to a SymPy ConditionSet
; the _sympy_
method falls back to creating a SageSet
wrapper otherwise.
symbolic_expression(lambda x, y: ...)
Sage 9.4 has added a new way to create callable symbolic expressions. #32103
The global function symbolic_expression
now accepts a callable such as those created by lambda
expressions. The result is a callable symbolic expression, in which the formal arguments of the callable are the symbolic arguments.
Example:
symbolic_expression(lambda x,y: x^2+y^2) == (SR.var("x")^2 + SR.var("y")^2).function(SR.var("x"), SR.var("y"))
This provides a convenient syntax in particular in connection to ConditionSet
.
Instead of
sage: predicate(x, y, z) = sqrt(x^2 + y^2 + z^2) < 12 # preparser syntax, creates globals
sage: ConditionSet(ZZ^3, predicate)
one is now able to write
sage: ConditionSet(ZZ^3, symbolic_expression(lambda x, y, z:
....: sqrt(x^2 + y^2 + z^2) < 12))
Convex geometry
ABC for convex sets
Sage 9.4 has added an abstract base class ConvexSet_base (as well as abstract subclasses ConvexSet_closed
, ConvexSet_compact
, ConvexSet_relatively_open
, ConvexSet_open
) for convex subsets of finite-dimensional real vector spaces. The abstract methods and default implementations of methods provide a unifying API to the existing classes Polyhedron_base
, ConvexRationalPolyhedralCone
, LatticePolytope
, and PolyhedronFace
. #31919, #31959, #31990, #31993
Several methods previously only available for Polyhedron_base
instances, are now available for all convex sets. The method an_affine_basis
returns a sequence of points that span by affine linear combinations the affine hull, i.e., the smallest affine subspace in which the convex set lies. The method affine_hull
returns the latter as a polyhedron. The method affine_hull_projection (renamed from affine_hull
in Sage 9.1) computes an affine linear transformation of the convex set to a new ambient vector space, in which the image is full-dimensional. The generalized method also provides additional data: the right inverse (section map) of the projection. #27366, #31963, #31993
As part of the ConvexSet_base
API, there are new methods for point-set topology such as is_open
, relative_interior
, and closure
. For example, taking the relative_interior
of a polyhedron constructs an instance of RelativeInterior, a simple object that provides a __contains__
method and all other methods of the ConvexSet_base
API. #31916
sage: P = Polyhedron(vertices=[(1,0), (-1,0)])
sage: ri_P = P.relative_interior(); ri_P
Relative interior of
a 1-dimensional polyhedron in ZZ^2 defined as the convex hull of 2 vertices
sage: (0, 0) in ri_P
True
sage: (1, 0) in ri_P
False
ConvexSet_base
is a subclass of the new abstract base class Set_base. #32013
This makes various methods that were previously only defined for sets constructed using the Set
constructor available for polyhedra and other convex sets. As an example, we can now do:
sage: polytopes.cube().union(polytopes.tetrahedron()) ...
9.3
Release Tour
SageMath 9.3 was released on May 9, 2021 (changelog).
- adds support for macOS 11 "Big Sur"
- major package upgrades
- 679 tickets (PRs) merged, 112 contributors
Linear and multilinear algebra
Bär--Faddeev--LeVerrier algorithm for the Pfaffian of skew-symmetric matrices
According to https://arxiv.org/abs/2008.04247, the Pfaffian of skew-symmetric matrices over commutative torsion-free rings can be computed with a Faddeev--LeVerrier-like algorithm. This algorithm is now implemented under the weaker assumption of the base ring's fraction field being a QQ
-algebra (#30681). It leads to a significant increase of computational speed in comparison to the definition involving perfect matchings, which has been the only algorithm available in Sage so far.
Using the definition of the Pfaffian:
sage: A = matrix([(0, 0, 1, 0, -1, -2, -1, 0, 2, 1),
(0, 0, 1, -3/2, 0, -1, 1/2, 3, 3/2, -1/2),
(-1, -1, 0, 2, 0, 5/2, 1, 0, -2, 1),
(0, 3/2, -2, 0, 5/2, -1, 2, 0, -1, -3/2),
(1, 0, 0, -5/2, 0, 0, -1, 1/2, 1, -1),
(2, 1, -5/2, 1, 0, 0, 2, 1, 2, 1),
(1, -1/2, -1, -2, 1, -2, 0, 0, -3, -1),
(0, -3, 0, 0, -1/2, -1, 0, 0, 1/2, 1/2),
(-2, -3/2, 2, 1, -1, -2, 3, -1/2, 0, 1),
(-1, 1/2, -1, 3/2, 1, -1, 1, -1/2, -1, 0)])
sage: %%time
....: A.pfaffian(algorithm='definition')
CPU times: user 18.7 ms, sys: 0 ns, total: 18.7 ms
Wall time: 18.6 ms
817/16
With Bär--Faddeev--LeVerrier:
sage: A = matrix([(0, 0, 1, 0, -1, -2, -1, 0, 2, 1),
(0, 0, 1, -3/2, 0, -1, 1/2, 3, 3/2, -1/2),
(-1, -1, 0, 2, 0, 5/2, 1, 0, -2, 1),
(0, 3/2, -2, 0, 5/2, -1, 2, 0, -1, -3/2),
(1, 0, 0, -5/2, 0, 0, -1, 1/2, 1, -1),
(2, 1, -5/2, 1, 0, 0, 2, 1, 2, 1),
(1, -1/2, -1, -2, 1, -2, 0, 0, -3, -1),
(0, -3, 0, 0, -1/2, -1, 0, 0, 1/2, 1/2),
(-2, -3/2, 2, 1, -1, -2, 3, -1/2, 0, 1),
(-1, 1/2, -1, 3/2, 1, -1, 1, -1/2, -1, 0)])
sage: %%time
....: A.pfaffian(algorithm='bfl')
CPU times: user 554 µs, sys: 41 µs, total: 595 µs
Wall time: 599 µs
817/16
Other improvements
- Speedup access items in gf2e dense matrices: 29853
- Speedup conjugation of double dense matrices: 31283
- New Bunch-Kaufman block_ldlt() factorization for possibly indefinite matrices: 10332
- New numerically stable is_positive_semidefinite() method for matrices: 10332
Polyhedral geometry
New features
The Schlegel diagrams are now repaired (they previously broke convexity). Now, one specifies which facet to use to do the projection 30015:
#!python
sage: fcube = polytopes.hypercube(4)
sage: tfcube = fcube.face_truncation(fcube.faces(0)[0])
sage: tfcube.facets()[-1]
A 3-dimensional face of a Polyhedron in QQ^4 defined as the convex hull of 8 vertices
sage: sp = tfcube.schlegel_projection(tfcube.facets()[-1])
sage: sp.plot() # The proper Schlegel diagram is shown
A different values of position
changes the projection:
#!python
sage: sp = tfcube.schlegel_projection(tfcube.facets()[4],1/2)
sage: sp.plot()
Graphics3d Object
sage: sp = tfcube.schlegel_projection(tfcube.facets()[4],4)
sage: sp.plot()
Graphics3d Object
New features:
- 30704: Upgrade to Normaliz 3.8.9 and PyNormaliz 2.13
- 30946: Add "minimal=True" option to affine_hull_projection
- 30954: Implement a proper equality check for polyhedron representation objects
Implementation improvements
The zonotope construction got improved:
Before:
sage: from itertools import combinations
sage: cu = polytopes.cube()
sage: sgmt = [p.vector()-q.vector() for p,q in combinations(cu.vertices(),2)]
sage: sgmt2 = set(tuple(x) for x in sgmt)
sage: # %time polytopes.zonotope(sgmt) # killed due to memory overflow
sage: %time polytopes.zonotope(sgmt2)
CPU times: user 2.06 s, sys: 23.9 ms, total: 2.09 s
Wall time: 2.09 s
A 3-dimensional polyhedron in ZZ^3 defined as the convex hull of 96 vertices
With 31038:
sage: from itertools import combinations
sage: cu = polytopes.cube()
sage: sgmt = [p.vector()-q.vector() for p,q in combinations(cu.vertices(),2)]
sage: sgmt2 = set(tuple(x) for x in sgmt)
sage: %time polytopes.zonotope(sgmt)
CPU times: user 138 ms, sys: 0 ns, total: 138 ms
Wall time: 138 ms
A 3-dimensional polyhedron in ZZ^3 defined as the convex hull of 96 vertices
sage: %time polytopes.zonotope(sgmt2)
CPU times: user 58 ms, sys: 0 ns, total: 58 ms
Wall time: 57.6 ms
A 3-dimensional polyhedron in ZZ^3 defined as the convex hull of 96 vertices
Improvements:
- 30040: Improve face iterator for simple/simplicial polytopes
There are also some bug fixes and other improvements. For more details see the release notes for optimization and polyhedral geometry software interactions in Sage.
Graph theory {#graph_theory}
Major improvements in the backends:
- 30777: Deleting edges
- 30665: Edge iterator and copy
- 30776: Subgraph and equality check
- 30753: Obtaining subgraphs.
- 31117, 31154: Breadth First Search
- 31129: Depth first search
- 31197: Use binary matrix data structure for bitsets.
Algebra
Power Series Ring
-
The method
set_default_prec
is now deprecated since it led to unwanted behavior (see #18416 for details). If another default precision is needed, a new power series ring must be created:sage: R. = [[PowerSeriesRing]](QQ, default_prec=10)
sage: sin(x)
x - 1/6x^3 + 1/120x^5 - 1/5040x^7 + 1/362880x^9 + O(x^10)
sage: R. = [[PowerSeriesRing]](QQ, default_prec=15)
sage: sin(x)
x - 1/6x^3 + 1/120x^5 - 1/5040x^7 + 1/362880x^9 - 1/39916800x^11 + 1/6227020800x^13 + O(x^15)This change does not affect the behavior of its ring elements. Code that relies on this method needs to be updated.
-
Inversion of power series ring elements now provides the correct parent: #8972
Other additions, improvements, and key bug-fixes
-
Implement the symplectic derivation Lie algebra following https://arxiv.org/abs/2006.06064:
sage: lie_algebras.SymplecticDerivation(QQ, 4)
Symplectic derivation Lie algebra of rank 4 over Rational Field -
Implement the *-insertion algorithm from https://arxiv.org/abs/1911.08732:
sage: from sage.combinat.rsk import [[RuleStar]]
sage: p,q = [RuleStar].forward_rule([1,1,2,2,4,4], [1,3,2,4,2,4])
sage: ascii_art(p, q)2 4 1 1 2
4 2 4
3 4
sage: line1,line2...
9.2
Release Tour
SageMath 9.2 was released October 24, 2020 (changelog), 727 tickets (PRs) merged, 134 contributors.
- Python 3 support expanded to 3.6, 3.7, 3.8, 3.9; Python 2 support dropped
- major package upgrades
Python 3 transition completed
SageMath 9.0 was the first version of Sage running on Python 3 by default. SageMath 9.1 continued to support Python 2.
Support for Python 2 removed
Sage 9.2 has removed support for Python 2. The Sage library now makes use of Python language and library features that are only available in Python 3.6 or newer; and large amounts of compatibility code have been removed.
However, note that this is unrelated to the minimal requirements for a source installation of the Sage distribution: Sage 9.2 is still able to build on a system that only provides Python 2.x or Python 3.5 or older. In this case, the SageMath distribution builds its own copy of Python 3.
Support for Python 3.6, 3.8, and 3.9 added
Sage 9.2 has added support for Python 3.8 in #27754 and Python 3.9 in #30184.
Sage 9.2 has also added support for Python 3.6. This allows Sage to use the system Python on some older Linux distributions that are still in widespread use in scientific computing, including centos-8
and fedora-{26,27,28}
(although Python 3.7.x packages are also available for these). See #29033 for more details.
Hence, Sage 9.2 conforms to (and exceeds) NumPy Enhancement Proposal 29 regarding Python version support policies.
If no suitable system Python, versions 3.6.x, 3.7.x, 3.8.x, or 3.9.x is found, Sage installs its own copy of Python 3 from source. The version of Python shipped with the Sage distribution has been upgraded from 3.7.3 to 3.8.5.
For developers: Using Python 3.6+ features in sagelib
Meta-ticket #29756 provides a starting point for a discussion of new features of the Python language and standard library to bring them to systematic use in sagelib. All features provided by Python 3.6 can be used immediately; features introduced in Python 3.7 or later will require backporting or a decision to drop the goal of supporting Python 3.6.
More details
- Trac tickets with keyword/component python3 in milestone 9.2
- See Python3-Switch for more details.
Package upgrades
The removal of support for Python 2 has enabled major package upgrades.
Major user-visible package upgrades below...
matplotlib
Dropping Python 2 support allowed us to make a major jump from matplotlib 2.2.5 to 3.3.1. See matplotlib's release notes for 3.0, 3.1, 3.2,3.3. In addition to improved output, this update will likely enable Sage developers to implement new features for plotting and graphics.
rpy2 and R
The rpy2 Python package is the foundation for SageMath's interface to R. Dropping Python 2 support allowed us to make the major upgrade from 2.8.2 to 3.3.5 in #29441; see the release notes for details.
We only did a minor upgrade of R itself in the Sage distribution, to 3.6.3, the latest in the 3.6.x series. Of course, if R 4.0.x is installed in the system, Sage will use it instead of building its own copy.
The SageMath developers are eager to learn from users how they use the SageMath-R interface, and what needs to be added to it to become more powerful. Let us know at sage-devel.
sphinx
Sage uses Sphinx to build its documentation. Sage 9.2 has updated Sphinx from 1.8.5 to 3.1.2; see Sphinx release notes for more information.
SymPy
SymPy has been updated from 1.5 to 1.6.2 in #29730, #30425. See the Release notes.
IPython, Jupyter notebook, JupyterLab
Dropping support for Python 2 allowed us to upgrade IPython from 5.8.0 to 7.13.0 in #28197. See the release notes for the 6.x and 7.x series.
We have also upgraded the Jupyter notebook from 5.7.6 to 6.1.1 in #26919; see the notebook changelog for more information. Besides, the pdf export of Jupyter notebooks has been fixed, so that LaTeX-typeset outputs are now rendered in the pdf file (#23330).
JupyterLab is now fully supported as an optional, alternative interface #30246, including interacts. To use it, install it first, using the command sage -i jupyterlab_widgets
. Then you can start it using ./sage -n jupyterlab
.
Normaliz
The optional package Normaliz, a tool for computations in affine monoids, vector configurations, lattice polytopes, rational cones, and algebraic polyhedra has been upgraded from 3.7.2 to 3.8.8, and PyNormaliz to version 2.12.
The upgrade adds support for incremental ("dynamic") computations, the computation of automorphism groups and refined triangulations of cones and polyhedra, and limited support for semiopen cones and polyhedra.
To install Normaliz and PyNormaliz, use sage -i pynormaliz
.
SageTeX
Updated to version 3.5, improving Python 3 compatibility, also updated to version 3.5 on CTAN.
sws2rst + usage example
In ticket #28838, the command sage -sws2rst
was resurrected via a new pip-installable package sage-sws2rst. It can be installed in Sage 9.2 using sage -i sage_sws2rst
.
Below is an example of usage. First we download a sage worksheet (.sws
) prepared for Sage Days 20 at CIRM (Marseille, 2010):
$ wget http://slabbe.org/Sage/2010-perpignan/CIRM_Tutorial_1.sws
$ ls
CIRM_Tutorial_1.sws
We translate the sws worksheet into a ReStructuredText syntax file (.rst
) using sage -sws2rst
. This creates also a directory of images:
$ sage -sws2rst CIRM_Tutorial_1.sws
Processing CIRM_Tutorial_1.sws
File at CIRM_Tutorial_1.rst
Image directory at CIRM_Tutorial_1_media
$ ls
CIRM_Tutorial_1_media CIRM_Tutorial_1.rst CIRM_Tutorial_1.sws
Then, we can check that it works properly by looking at the generated rst file. Alternatively, we can translate it to a basic html file using rst2html
:
$ rst2html.py CIRM_Tutorial_1.rst CIRM_Tutorial_1.html
CIRM_Tutorial_1.rst:176: (WARNING/2) Explicit markup ends without a blank line; unexpected unindent.
CIRM_Tutorial_1.rst:334: (WARNING/2) Inline strong start-string without end-string.
As seen above, there are few warnings sometimes because the translation made by sws2rst
is not 100% perfect, but most of it is okay:
$ firefox CIRM_Tutorial_1.html
Moreover, one can use the sage -rst2ipynb
script to translate the rst file obtained above to a Jupyter notebook:
$ sage -rst2ipynb ...