initial release

The code is working now on arbitrary-size matrices. I still have not checked for memory leaks. The code appears to be fast, but I have not profiled it yet.