Function matrix_log 

pub fn matrix_log<T: Float + Sum + FloatConst>(r: &M33<T>) -> M33<T>

Brief.

Compute the matrix logarithm of a rotation matrix

Function arguments: r: &M33 rotation matrix

Output: M33: representing the matrix logarithm of the input