Sage Days 96
August 13-16 2018
Sage days 96 to be held at the Fields Institute for the program Teichmüller Theory and its Connections to Geometry, Topology and Dynamics. The schedule is available on the Fields website. This wiki page will be updated with all documents relative to the workshop (worksheets, tutorials, demos, nice pictures, etc).
In order to assist to the workshop you are asked to come with your laptop (though you are also welcome without laptop). It would also be great if you could make sure that Sage is installed on it (see "Sage Installation" below).
Workshop dinner on thursday (Sky Dragon on Dundas and Spadina).
Contents
Worksheets
Introduction with the picture S_2_1.svg to be downloaded in the same folder (the worksheet is also available as intro.en.rst and intro.en.pdf)
- Programming
- Thematic tutorials
Python dictionaries and graph theory (you also need to download the images euler.png and graph0.png)
flipper tutorial (for the last exercises you will need flipper_nf_conversion.py)
Participant projects
Projects for the week
Projects that were done during the week:
- Aaron (Toronto)
- use flatsurf to improve old code
- given a triangulated flat surface and a direction, find a representation as the suspension of an IET
- shear coordinates of hyperbolic surfaces
- Ali (Teheran)
- given a flat surface, compute the algebraic equation of the underlying Riemann surface
- given a triangulated quadratic differential, decide if it admits a Strebel direction.
- equations for covers of Riemann surfaces (and for quotients!)
rewrite code that was previously done in Mathematica (because Sage is OpenSource)
- Anja (Toronto)
- study some infinite iet on [0,1] with rational discontinuities
- order of certain periodic IET's presented in nice way
- Chris (Indiana)
- find the Iso-Delaunay decomposition (in L^2) of Teichmüller discs
- collect evidence the Iso-Delaunay decomposition (in L^2) of Q(S) - are the "cells" contractible?
- Ian (Toronto)
- Delaunay problem: Given a triangulated quadratic differential, compute the L^infinity Delaunay triangulation
- Build the horocycle automaton. (Saul: Once you have the veering automaton, you just need to find the horocyclic edges....)
- Build all of the veering automata for strata of dimension at most four.
- Rex (Toronto)
- work with homogeneous spaces (integration over H^2 / Gamma and possibly working in higher rank)
- Saul (Warwick)
- layout problem: given a triangulated quadratic differential, draw it (ie unfold the triangulation in the plane without crashes)
- automate producing automata - how can we lay out the graphs in a sensible fashion? What are their automorphisms?
- Ser-Wei
- redo some old examples in Sage
- Suppose that D is a (Veech!) Teichmuller disk in M_{0,5}. Fix \tau a track that carries all of the boundary laminations for D. What is the image of \bdy D in P(\tau) (the polytope of projectivised measured laminations)? Draw it in Sage!
- Vincent (Bordeaux)
- given a translation surface, decompose the vertical into components
given a stratum Q(\kappa) find a geometric veering triangulation in it. (aka a seed for the automaton) -> done!!
- given a geometric veering triangulation, decide if it is regular or not.
- generate cylinder diagrams for quadratic differentials
- layout of surface, loaded from flipper
- Adriano (Toronto), Christine (New York), Dan (Toronto), Derrick (Toronto), Paul (Oklahoma), Pratyush (Yale), Yvon (Toronto)
- worksheets
Sage Installation
The installation procedure depends on your operating system:
Debian (versions >= 9.0) and Ubuntu (versions >= 18.04): install the packages sagemath, sagemath-doc-en and sagemath-jupyter
archlinux: install the packages sagemath, sagemath-doc and sagemath-jupyter
gentoo: you can give a try to sage on gentoo
all other systems (GNU/Linux, OSX, Windows): go to the main SageMath website, go to the download section, choose a download server close to your location and download the Sage binary that corresponds to your operating system.
If you run GNU/Linux or OSX you can also try to compile the software from sources. The procedure is described in the documentation.
Installing flipper, snappy, flatsurf and surface_dynamics
We will be interested in these Sage days to four Python-Sage modules that provide many tools to work with surfaces and 3-dimensional manifolds
flipper: mapping classes (homeomorphisms of surfaces)
flatsurf: translation surfaces (affine transformation, linear flow, etc)
surface_dynamics: interval exchange transformations, origamis and more
snappy: 3-d hyperbolic geometry
These packages are installed through pip. Depending how you did install Sage you might run in different kind of complications.
Generic Installation
This installation would work in most situations. Though if you happen to be in one of the specific situations listed below, just jump to that section.
- To install flipper, it should work out of the box with the following command (if you do not know how to open a command line, ask one of the organizer)
$ sage -pip install flipper --user
If the above command does not work, ask for help. If it does work, go to next step.
- To install surface_dynamics and snappy you need a working version of gcc and g++ (that should be the case for most of you). The commands are similar to the one before
$ sage -pip install surface_dynamics --user $ sage -pip install snappy --user
Similarly, report any error or weird behavior.
Lastly, to install flatsurf, look at the instructions at https://github.com/videlec/sage-flatsurf
Mac OS X specifics
If you are using the binaries for Mac OS X, the computer might complain with could not fetch URL .... In that situation, you can try the two following solutions
- You do have Xcode installed with the command line tools. In this case you can recompile the Python shipped with Sage with SSL support with
$ sage -i openssl $ sage -f python2
Then go back to the instructions "Installing flipper, snappy, flatsurf and surface_dynamics". - The alternative to option 1 is to download each package independently: and then install them via
$ sage -pip install flipper-0.12.3.tar.gz $ sage -pip install surface_dynamics-0.4.0.tar.gz $ sage -pip install snappy-2.6.tar.gz $ sage -pip master.zip
Ubuntu 16.04 with ppa
If you run Ubuntu 16.04 and you install Sage with the ppa then do
$ sage -sh (sage-sh)$ pip install flipper --user (sage-sh)$ pip install surface_dynamics --user (sage-sh)$ pip install snappy --user (sage-sh)$ exit
(The sage -sh command brings you to a special environment where you have access to Sage Python)
Ubuntu >= 18.04 and Debian >= 8.0 packages
If you happen to have installed the packages sagemath and sagemath-jupyter from your package manager in Debian/Ubuntu then you would need to do
$ source /usr/share/sagemath/bin/sage-env $ pip install flipper --user $ pip install surface_dynamics --user $ pip install snappy --user
For the installation of flatsurf, follows the instruction at https://github.com/videlec/sage-flatsurf and replace any instance of sage -pip by pip.
Archlinux
Do
$ sage -sh (sage-sh)$ pip2 install flipper --user (sage-sh)$ pip2 install surface_dynamics --user (sage-sh)$ pip2 install snappy --user (sage-sh)$ exit $
For the installation of flatsurf, follows the instruction at https://github.com/videlec/sage-flatsurf and replace any instance of sage -pip by pip.
Windows
Before following the generic instruction do
(sage-sh) $ apt-cyg install libcrypt-devel (sage-sh) $ apt-cyg install liblapack-devel
Organizers and participants
See the Fields Institute webpage.
Questions and answers
Saul: Why does snake_case have that name?
Vincent: See the discussion here - https://en.wikipedia.org/wiki/Snake_case.
Derrick: If I'm running Sage on a Windows installation, how to I make Jupyter launch in the correct directory?
Derrick: [question 41348 on ask](https://ask.sagemath.org/question/41348/sagemath-windows-81-launching-jupyter-from-another-directory-or-drive/)
Aaron: Does Python have a built-in way to iterate over the product of two iterators that don't necessarily terminate?
Vincent: No, but Sage does have a [cantor_product](https://doc.sagemath.org/html/en/reference/misc/sage/misc/mrange.html). Example:
produces (0, 0), (1, 0), (0, 1), (2, 0), (1, 1), (0, 2)... Although the feature isn't documented, iterating up to (0, n) only requires n calls to NN.next().
Rex: Let L be a list. To get a new list K with the elements of L in reverse order, we can write K = L[::-1] What is each the colons and the -1 telling Sage to do exactly?
Vincent: These are [Python slices](https://docs.python.org/2.3/whatsnew/section-slices.html)
Aaron: How do we specify a quadratic differential?
Vincent: triangulated surface, where the lengths live in some arithmetic. (Floating point numbers are not numbers!)
Saul: is there a way to start with floating point surface, and quickly rule out "small" components (eg no short cylinders)?
Rex: I'm running something that's printing a list of integers as it runs. Is there any way to tell Sage to truncate the old output (I'm only interested in seeing how far it has gone along)?
Vincent: If you do
from IPython.display import clear_output
then you can use the function clear_output() to remove everything in the output. Though this slows down a lot the computation.
The following previous answer does not work (in Jupyter)
from sys import stdout from time import sleep stdout.write('hello\r') sleep(1) # wait for 1 sec stdout.write('bouh \r') sleep(1) # wait for 1 sec stdout.write('next \r')
The character '\r' means carriage return and go back the the begining of the line so that you can erase what you just wrote. But with this method, it only works with one line being updated.
Rex: Do different packages ever have any dependencies such that uninstalling one package causes problems for others? For example, if I uninstall flipper at some later point in time (say to reinstall an updated version), could this cause issues with other packages that might rely on flipper?
Saul: Yes, this could happen. Package management is tricky.
Rex: Also, is there a simpler way to update a package than uninstalling and reinstalling it?
Saul:Something like
pip install <PACKAGE> --user --upgrade
Troubles and bug reports
SageMath documentation on archlinux
operating system: archlinux sage version: 8.2 installation: from the package manager pacman problem: the documentation with the question mark is not working properly in the console. This has been reported to https://bbs.archlinux.org/viewtopic.php?id=239525
Archlinux pip installation problem
The local installation (done via "pip install X --user") are different within Sage and within Python Namely
$ pip2 install X --user $ sage -sh -c "pip2 install X --user"
The reason is that PYTHONUSERBASE is set by sage-env before Sage is launched. A patch has been sent to Antonio Rojas (who is the maintainer). And it is now solved!
Misleading comment in cygwin version
Reported at https://groups.google.com/forum/#!topic/sage-devel/wdjj8ur5m9c
Numbering mismatch in surface_dynamics
Rex: the labeling of the square tiles in an origami is mismatched with the labelings of the vertices in its associated graph. The squares start with 1,2,3,... but the graph vertices begin with 0,1,2,...
Documentation requests
The example for Graph.vertex_cover makes it look like the method outputs a plot of the graph with the vertex cover colored. https://doc.sagemath.org/html/en/reference/graphs/sage/graphs/graph.html#sage.graphs.graph.Graph.vertex_cover
There's no unambiguous description of what the EllipticCurve constructor does. https://doc.sagemath.org/html/en/constructions/elliptic_curves.html
The documentation for installing Sage packages under Windows is very misleading. The page
opens with the sentence, "The major known caveat to the current version of the SageMath Windows install is that it is not currently possible to install Sage optional packages," with no indication that it is possible to install Sage pip packages. The Windows Sage shell has a similarly misleading warning: "Do not use this for installing Sage packages using 'sage -i' ", with no indication that sage pip is usable.
The documentation for flatsurf.TranslationSurface.delaunay_single_join is somewhat lacking.
Link to Pad
We used to coordinate using a pad on https://mensuel.framapad.org/ that is now closed.