Strip out new sampling code: simpler to just transpose everything
Re-do indexing throughout the library - untested
Working single loc version
Change sampling scheme, more amenable to future parallelization.
Todo: depricate and remove old sampling code
Clean up matrix multiplication for transposed case. I've really got to document this stuff soon...