Linalg: function for diagonalize hermitian commuting matrices simultaneously

Dear SciPy Developers, I like to contribute a function for simultaneously diagonalize a set of hermitian commuting matrices. Especially for matrices with degenerate eigenvalues. Is there need for something like that in SciPy ? I think it could be helpful for someone else and I didn't find it in recent SciPy versions or on the mailing list. Hope to hear from you.

Hi Lars-Hendrik This does not answer your question and has nothing to do with scipy, but we have a somewhat specialized implementation of simultaneous diagonalization of hermitian matrices in qutip (a package for quantum dynamics) https://github.com/qutip/qutip/blob/master/qutip/simdiag.py Perhaps it could be useful for you when you work on your implementation. Also, if you have, or find, a better or more efficient way of implementing this function we (qutip project) would be interested in hearing from you as well. Best regards Robert

Thank you Robert for the note. I have seen this implementation but had some trouble in using it. I going to mail you after I figured out what is the problem and maybe we can discuss an alternative implementation. But I am sure this implementation will help me a lot. Regards Lars-Hendrik

Lars-Hendrik Frahm <larshendrikfrahm <at> googlemail.com> writes:
Dear SciPy Developers,
I like to contribute a function for simultaneously diagonalize a set of hermitian commuting matrices. Especially for matrices with degenerate eigenvalues.
Is there need for something like that in SciPy ? I think it could be helpful for someone else and I didn't find it in recent SciPy versions or on the mailing list.
Hope to hear from you.
I do have a need for exactly that capability in python. Have you developed such a module? Best, Murray Daw
participants (3)
-
jrjohansson@gmail.com
-
Lars-Hendrik Frahm
-
Murray Daw