def logm(A,disp=1):
"""Matrix logarithm, inverse of expm."""
# Compute using general funm but then use better error estimator and
# make one step in improving estimate using a rotation matrix.
A = mat(asarray(A))
All other matrix functions use A = asarray(A) instead of mat(...) .
Nils