Successfully rewrote func, per issue #3. However ...

It looks like it won't be as efficient for small matrices. It *might* be okay
for larger ones? Maybe even better? (Thinking memory, not speed... need to be
more thorough about TCR, though.)