add matrix logarithm

[grayj2]

A matrix log is needed to solve A in e^A=B (i.e. not an element wise log). It doesn't look like this method is available under any of the current dependencies, but is available under jeigen.

https://github.com/hughperkins/jeigen#matrix-exponential-matrix-logarithm