Alpha release

This is an early alpha release. All functions seem to work so far.
The functions should all use sparse matrices when possible.