Output-sensitive modular algorithms for polynomial matrix normal forms