Compute the matrix logarithm of the square matrix A.
The implementation utilizes a Padé approximant and the identity
logm (A) = 2^k * logm (A^(1 / 2^k))
The optional input opt_iters is the maximum number of square roots to compute and defaults to 100.
The optional output iters is the number of square roots actually computed.
See also: expm, sqrtm.
Package: octave