b45ch1/algopy — pull request #10 by argriffing (2012-10-16 17:10:17+00:00)
Title: more sophisticated expm from scipy, but somewhat incomplete and under-tested

I've added expm_pade(A, q) which you can feed an explicit order q in {3,5,7,9,13}, and I've also added an incomplete and under-tested implementation of the full Higham 2005 variable-order Pade approximation with scaling and squaring. More test matrices will be needed for code coverage, and functions like min, max, ceil, linalg.norm(A, 1) may need to be added if they do not already exist and if they are even philosophically compatible with automatic differentiation.