Differences between revisions 45 and 46
Revision 45 as of 2021-03-23 20:38:05
Size: 21836
Editor: mkoeppe
Comment: system Python 3.6 deprecated
Revision 46 as of 2021-03-23 20:39:11
Size: 21887
Editor: mkoeppe
Comment: add ticket link
Deletions are marked like this. Additions are marked like this.
Line 211: Line 211:
SageMath 9.3 continues to support system Python versions 3.6.x, 3.7.x, 3.8.x, and 3.9.x; support for Python 3.6.x is deprecated and will be removed in the SageMath 9.4 development cycle. If no suitable system Python 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.8.5 to 3.9.2. [[https://trac.sagemath.org/ticket/30589|#30589]], [[https://trac.sagemath.org/ticket/31419|#31419]] SageMath 9.3 continues to support system Python versions 3.6.x, 3.7.x, 3.8.x, and 3.9.x; support for Python 3.6.x is deprecated and will be removed in the SageMath 9.4 development cycle. If no suitable system Python 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.8.5 to 3.9.2. [[https://trac.sagemath.org/ticket/30589|#30589]], [[https://trac.sagemath.org/ticket/31318|#31318]], [[https://trac.sagemath.org/ticket/31419|#31419]]

Sage 9.3 Release Tour

in progress (2021)

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:

   1 sage: fcube = polytopes.hypercube(4)
   2 sage: tfcube = fcube.face_truncation(fcube.faces(0)[0])
   3 sage: tfcube.facets()[-1]
   4 A 3-dimensional face of a Polyhedron in QQ^4 defined as the convex hull of 8 vertices
   5 sage: sp = tfcube.schlegel_projection(tfcube.facets()[-1])
   6 sage: sp.plot()  # The proper Schlegel diagram is shown

A different values of position changes the projection:

   1 sage: sp = tfcube.schlegel_projection(tfcube.facets()[4],1/2)
   2 sage: sp.plot()
   3 Graphics3d Object
   4 sage: sp = tfcube.schlegel_projection(tfcube.facets()[4],4)
   5 sage: sp.plot()
   6 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

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

Manifolds

Spheres added to the manifold catalog

Spheres of arbitrary dimension have been added to the manifold catalog, so that they can be initialized via manifolds.Sphere(n) (#30804). For instance:

sage: S3 = manifolds.Sphere(3)
sage: S3
3-sphere S^3 of radius 1 smoothly embedded in the 4-dimensional Euclidean space E^4

By default, a single chart is initialized: that of spherical coordinates:

sage: S3.atlas()
[Chart (A, (chi, theta, phi))]
sage: S3.default_chart() is S3.spherical_coordinates()
True

Note that A stands for the domain of spherical coordinates, which is a strict subset of S3:

sage: S3.spherical_coordinates().domain()                                                           
Open subset A of the 3-sphere S^3 of radius 1 smoothly embedded in the 4-dimensional Euclidean space E^4

The sphere S3 is considered as an embedded submanifold of the Euclidean space E4:

sage: S3.embedding().display()
iota: S^3 --> E^4
on A: (chi, theta, phi) |--> (x1, x2, x3, x4) = (cos(phi)*sin(chi)*sin(theta), sin(chi)*sin(phi)*sin(theta), cos(theta)*sin(chi), cos(chi))

S3 is a Riemannian manifold, with the metric induced by the Euclidean metric on E4:

sage: g = S3.induced_metric()
sage: g.display()
gamma = dchi*dchi + sin(chi)^2 dtheta*dtheta + sin(chi)^2*sin(theta)^2 dphi*dphi

Beside spherical coordinates, stereographic coordinates are available too. Those from the North pole (NP) are obtained by

sage: S3.stereographic_coordinates()
Chart (S^3-{NP}, (y1, y2, y3))

The above call has augmented the manifold atlas with various charts and subcharts, corresponding to stereographic coordinates from the North pole and from the South pole:

sage: S3.atlas()
[Chart (A, (chi, theta, phi)),
 Chart (S^3-{NP}, (y1, y2, y3)),
 Chart (S^3-{SP}, (yp1, yp2, yp3)),
 Chart (S^3-{NP,SP}, (y1, y2, y3)),
 Chart (S^3-{NP,SP}, (yp1, yp2, yp3)),
 Chart (A, (y1, y2, y3)),
 Chart (A, (yp1, yp2, yp3))]

All relevant change-of-coordinate formulas have been initialized. For instance, those relating spherical coordinates to stereographic ones are obtained as follows:

sage: spher = S3.spherical_coordinates()
sage: stereo = S3.stereographic_coordinates()
sage: A = spher.domain()
sage: stereoA = stereo.restrict(A)
sage: S3.coord_changes()[(spher, stereoA)].display()
y1 = -cos(phi)*sin(chi)*sin(theta)/(cos(chi) - 1)
y2 = -sin(chi)*sin(phi)*sin(theta)/(cos(chi) - 1)
y3 = -cos(theta)*sin(chi)/(cos(chi) - 1)

The change-of-coordinate formulas are automatically invoked to compute the coordinates of a point from previously known ones in a different chart. For instance:

sage: p = S3((pi/2, pi/3, pi/4), chart=spher)                                                       
sage: p                                                                                             
Point on the 3-sphere S^3 of radius 1 smoothly embedded in the 4-dimensional Euclidean space E^4
sage: spher(p)                                                                                      
(1/2*pi, 1/3*pi, 1/4*pi)
sage: stereo(p)                                                                                     
(1/4*sqrt(3)*sqrt(2), 1/4*sqrt(3)*sqrt(2), 1/2)

Use of ambient metric by default for the Hodge dual of a differential form

When asking for the Hodge dual of a differential form on a pseudo-Riemannian manifold without specifying the metric, the manifold's default metric is assumed (#31322). For example:

sage: E.<x,y,z> = EuclideanSpace()
sage: E.cartesian_coordinates().coframe()
Coordinate coframe (E^3, (dx,dy,dz))
sage: dx = E.cartesian_coordinates().coframe()[1]
sage: dx
1-form dx on the Euclidean space E^3
sage: dx.hodge_dual()
2-form *dx on the Euclidean space E^3
sage: dx.hodge_dual().display()
*dx = dy/\dz

Internal code improvements and bug fixes

Various improvements/refactoring of the code have been performed in this release:

* topological part: #30310, #31243

* differentiable part: #31215, #31255, #31273, #31323.

In addition, various bugs and typos have been fixed: #30174, #30830, #31202.

Numerics

Automatic conversions from from floating-point numbers and symbolic expressions to real and complex intervals are deprecated and will be removed in the future. Code that relies on the existence of coercions

  • from RealField to RealIntervalField,

  • from SymbolicRing to ComplexIntervalField, or

  • from Python floats to ComplexIntervalField,

like

sage: RIF(1/3) + RR(1)

or

sage: 1.1*pi*CIF(i)

needs to be updated. (Other similar conversions, e.g. from RealField to RealBallField or from float to RealField, already needs to be performed explicitly.) Beside leading to more consistent behavior, the removal of these coercions will make it easier to use interval arithmetic reliably. (#15114)

Package upgrades, system packages, user packages

Python

SageMath 9.3 continues to support system Python versions 3.6.x, 3.7.x, 3.8.x, and 3.9.x; support for Python 3.6.x is deprecated and will be removed in the SageMath 9.4 development cycle. If no suitable system Python 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.8.5 to 3.9.2. #30589, #31318, #31419

Other package upgrades

Upgrade tickets, milestone 9.3

SageMath on repology.org

The SageMath distribution is now listed as a repository on repology. The list of packages allows developers to see which packages are in need of updating.

System package information for OpenSUSE

In particular, users of OpenSUSE Tumbleweed will notice the information about system packages that the Sage distribution can use.

Chapter on packages in the Sage reference manual

We have consolidated the information about packages into a single source, build/pkgs/ in the source tree, from which we generate package lists for the reference manual #29655, to which also the Sage website points.

A number of user packages have been added as optional packages: ore_algebra, sage_flatsurf, admcycles, slabbe, snappy, surface_dynamics. Meta-ticket #31164 tracks the task of adding user packages to build/pkgs/.

For developers and packagers: Version constraints for Python packages

For all Python packages in the Sage distribution, build/pkgs/SPKG/install-requires.txt now encodes version constraints (such as lower and upper bounds) #30719. The constraints are in the format of the install_requires key of setup.cfg or setup.py.

The files may include comments (starting with #) that explain why a particular lower bound is warranted or why we wish to include or reject certain versions. For example:

$ cat build/pkgs/sphinx/package-version.txt 
3.1.2.p0
$ cat build/pkgs/sphinx/install-requires.txt
# gentoo uses 3.2.1
sphinx >=3, <3.3

The comments may include links to Trac tickets, as in the following example:

$ cat build/pkgs/packaging/install-requires.txt
packaging >=18.0
# Trac #30975: packaging 20.5 is known to work but we have to silence "DeprecationWarning: Creating a LegacyVersion"

The currently encoded version constraints are merely a starting point. Downstream packagers and developers are invited to refine the version constraints based on their experience and tests. When a package update is made in order to pick up a critical bug fix from a newer version, then the lower bound should be adjusted.

In Sage 9.3, the new files are included only for documentation purposes and to facilitate development and packaging. An effort is underway in #30913 to use the new files to generate dependency metadata for the Sage library. There is no mechanism yet to use system Python packages; see Meta-ticket #29013 for this task.

For developers: Setting up Python packages from PyPI as Sage packages

Setting up Python packages from PyPI as Sage packages has become easier #30974.

$ ./sage -package create scikit_spatial --pypi
Downloading tarball to /Users/mkoeppe/s/sage/sage-rebasing/worktree-algebraic-2018-spring/upstream/scikit-spatial-5.0.0.tar.gz
[......................................................................]
$ ls build/pkgs/scikit_spatial/
SPKG.rst                dependencies            package-version.txt     type
checksums.ini           install-requires.txt    spkg-install.in
$ cat build/pkgs/scikit_spatial/SPKG.rst 
scikit_spatial: Spatial objects and computations based on NumPy arrays
======================================================================

Description
-----------

Spatial objects and computations based on NumPy arrays.

License
-------

BSD license

Upstream Contact
----------------

https://pypi.org/project/scikit-spatial/

Note that this new command does not check the dependencies of the Python package; this still needs to be done manually.

The above command creates a "normal" package (with tarball information in checksums.ini). To create a package as a "pip" package (which will be obtained directly from PyPI on an install), use ./sage -package create scikit_spatial --pypi --source pip instead.

Cleaning of the Sage codebase to conform to best practices

... Status and tools ...

Modularization and packaging of sagelib

A source tarball of the Sage distribution (or a worktree of the Sage git repository after running the ./bootstrap script) contains several self-contained source trees of Python distribution packages in build/pkgs/*/src/, which can be packaged, built and installed by standard Python procedures:

  • setup.py sdist builds a source distribution (it is invoked by the convenience script build/pkgs/*/spkg-src), which can afterwards be installed using pip.

  • setup.py install installs the package directly (this is legacy use of setuptools and is not recommended).

  • setup.py bdist_wheel builds a wheel, which can afterwards be installed using pip.

  • Note, however, that pip install . does not work. In order to keep the monolithic structure of the SAGE_ROOT/src tree unchanged (for the convenience of Sage developers), the source trees of the Python distribution packages make use of symlinks. These symlinks are not compatible with pip install .. However, setup.py sdist follows the symlinks - the resulting source distribution contains ordinary files only and is therefore pip-installable.

The following Python distribution packages exist in Sage 9.3:

In the course of the modularization effort of Meta-ticket #29705, in the Sage 9.4 series we expect to add many more distribution packages in the same format.

  • Splitting out sage_setup (the build system of the Sage library) as a separate distribution package (currently it is shipped as part of sagemath-standard).

  • Providing alternative implementations of sage_conf for different use cases.

  • Splitting sagemath-standard into many namespace packages.

Configuration changes

Editable ("in-place", "develop") installs of the Sage library

Use ./configure --enable-editable to configure the Sage distribution to install the Sage library in "develop" ("editable", "in-place") mode instead of using the Sage library's custom incremental build system. #31377

It has the benefit that to try out changes to Python files, one does not need to run ./sage -b any more; restarting Sage is enough. It may also have benefits in certain develop environments that get confused by sagelib's custom build system.

Note that in an editable install, the source directory will be cluttered with build artifacts (but they are .gitignored). This is normal.

Use of system Jupyter notebook / JupyterLab

The Sage Installation Manual now has instructions on how to run the SageMath Jupyter kernel in a system Jupyter notebook or JupyterLab #30476

If you intend to (a) only use Sage in a terminal, or (b) only use the SageMath kernel with your system Jupyter notebook or JupyterLab, you can now build Sage without the notebook components, by using ./configure --disable-notebook. #30383 #31278

Tickets

Availability of Sage 9.3 and installation help

A release candidate is expected to be available soon.

  • See sage-devel for development discussions and sage-release for announcements of beta versions and release candidates.

ReleaseTours/sage-9.3 (last edited 2022-05-14 16:52:08 by mkoeppe)