Hello, in my current work, I have a bunch of matrices that I need to simultaneously triangularise. More precisely, given a set of matrices [M_1, ..., M_n] of End(V), I would like to obtain _if possible_ a base v_1, ..., v_k of V such that all M_1, ... M_n are triangular in this basis. Has anyone have some piece of code for it ? I already began to implement it, but I am probably not the only one who needs this feature.
The algorithm is pretty much basic: 1) Find a common eigenvector v_1 (if there is not we stop here: the set is not simultaneously triangularisable) 2) Complete into a basis of V: v_1, v'_2, ..., v'_k 3) Write all M's in this new basis 4) Continue on the (n-1) x (n-1) bottom left submatrix. The main issue is with 1), I do not see how to make this step efficient. Best, Aladin -- You received this message because you are subscribed to the Google Groups "sage-devel" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-devel+unsubscr...@googlegroups.com. To post to this group, send email to sage-devel@googlegroups.com. Visit this group at http://groups.google.com/group/sage-devel. For more options, visit https://groups.google.com/d/optout.