Early release

This is a mostly feature-complete implementation of LDshrink. Sparse matrix support is still lacking, and performance can still be improved.