Implementation of the Cheng and Church algorithm

My in-progress implementation of the Cheng and Church algorithm. Still to do:

- [ ] get base implementation working.
- [ ] get it working on sparse matrices
- [ ] incremental update of MSR
- [ ] tests
- [ ] data sample generator
- [ ] documentation
- [ ] examples