a -dom_eigenvects...rangercolsrrrrrsz"dom_eigenvects..domaincs g|]}fdd|DqS)cs g|]}tj|gdjqS)r )r from_listrepritemr!lr rr+rrr rrowr&r rr*scsg|]}fdd|DqS)csg|] }|qSr r r$rr rr-rrr r(r*r rr-rcs&g|]fddtDqS)csg|]}|krnjqSr rrrr rr0rrrrrrrr/s)charpolyshaper!rlenrr Z nullspaceappendrr"rr#to_ddm)Ar'r+rows_Zfactorsrational_eigenvectsalgebraic_eigenvectsbaseexpZEE_itemsZEEZbasisminpolyZAA_itemsZAAr )rr!rrr'rdom_eigenvectss>     r8c  sg}|D]D\}}|j}fdd|D}|||fq|D]\}}|j}jdfdd|D}}tfi|} t| |krʇfddt |D} | D]*fdd|D}|||fqqR|S)Ncs"g|]}fdd|DqS)csg|]}|qSr to_sympyrxr*r rrDr6dom_eigenvects_to_sympy...r rZvect)Matrixrr rrCsz+dom_eigenvects_to_sympy..rcsg|]}fdd|DqS)csg|]}|qSr r9r;r*r rrLrr=r r>r*r rrLrcsg|]}t|qSr r)ridx)r'r7r rrRrcs$g|]}fdd|DqS)csg|]}|qSr )subsr;) eigenvaluer'r rrVrr=r r>)r?rBr'r rrUs) r#r/r:r.ZgensdegreeZas_exprrr-r) r3r4r?kwargsresultZ multiplicityZ eigenvectsZnew_eigenvectsrCZ eigenvalsr )r?rBrr'r7rdom_eigenvects_to_sympy:s.      rFN)__doc__Zsympy.core.symbolrZagca.extensionsrZ factortoolsrZ polyrootsrZ polytoolsrZ rootoftoolsr Z domainmatrixr r8rFr r r rs       )