Major changes
* Added local function spaces (`locfunspace`) to `locfun/`
* Class containing a basis of the local Poisson space $V_p(K)$
* Added edge spaces (`edgespace`) to `locfun/`
* Added support for high order (including $p > 3$, unstable for some edges)
* Added barycentric coordinates and Legendre polynomials to `poly/`
* High order edge spaces obtained by traces of integrated Legendre polynomials
* Overhaul of `locfun` traces
* Local function traces are by default considered a list of `polynomial`s, one for each edge
* This list of `polynomial`s is stored in `trace_polys` attribute
* Trace values are computed only as needed
* Trace values can be set manually as before
* The flag `has_poly_trace` can be set to `False` if working functions with functions that do not have a trace that is the trace of a polynomial on each edge
* `locfun` objects are now initialized with the Laplacian polynomial, the list of polynomial traces, and an option
* Overhaul of Nystrom solver
* Replaced `nystrom` module with `nystrom_solver` class
* Overhead of constructing the Nystrom matrix is now consolidated into the constructor
* Harmonic conjugate computations are now handled as method of `nystrom_solver`
* Added `vert` class
* Vertices of a mesh
* Overhaul of `edge` class
* Added topological properties needed for mesh construction
* Parameterization values computed and stored only as needed
* Removed dependence on `copy` package
* Added `integrate_over_edge` method
* Added `global_orientation` attribute
* Overhaul of `cell` class
* Added topological properties needed for mesh construction
* Edge list replaced with `closed_contour` list
* Added `PlanarMesh` class
* Initialized as a collection of `edge` and `vert` objects
* `cell` objects are constructed in situ using topological information from `edge` objects
* Added `solver` class
* Handles all aspects of solving a PDE on a mesh
* Initialized with a `PlanarMesh` object and a `bilinear_form` object
* `solve` method solves the PDE on the mesh
* `plot` method plots the solution
* Added `bilinear_form` class
* Stores info about the PDE to be solved
* Also stores the right-hand side of the PDE as a polynomial
* Added `global_function_space` class to manage global function space
* Added `global_key` class to manage continuity across cells
Minor changes:
* Moved `quad` module to `mesh/` from `mesh/quad/`
* Added `polynomial` functionality
* Added division of a polynomial by a scalar
* Added powers of polynomials
* Added polynomial composition
* Changed directory management in Jupyter notebook examples
* Added `devtools/` directory
* Added `doc/` directory
* Added `.py` versions of examples
* Thanks to `nbconvert` for making this easy!
* Added `ex2-pacman-fem` example
* Updated tests and examples