Is there any ready code which gives an algorithm of expressing a matrix pencil in its Kronecker's Canonical Form?

There is an old result which gives an algorithm for an arbitrary pencil but it is lengthy and looks quite tedious.



GUPTRI by Jim Demmel and Bo Kagstrom computes a triangular decomposition that reveals the Kronecker structure of a pencil. It is Fortran code that can be called from Matlab using a Mex-file interface.

The code is quite old, though. Probably the Fortran part can still be compiled and run without changes, but I can speak for experience that getting old MEX files to run on a modern 64-bit Matlab can require lots of tweaking.

  • $\begingroup$ Thanks! I will try and let you know. Meanwhile I am almost done with mathematica code and I will post a comparison if I manage to run GUPTRI. $\endgroup$
    – Mari Gachi
    Sep 13 '16 at 11:13
  • $\begingroup$ Could you share how you solved the problem, please? I am having the same problem, too. $\endgroup$
    – Misaki
    Oct 17 '16 at 16:24
  • 1
    $\begingroup$ @Misaki I converted your non-answer to a comment, but it's not really clear what your question is. It might be better to ask a new question, linking back to this thread, but you should consider also if this is the right site for your exact concern (which may have more to do with programming languages than mathematics). $\endgroup$
    – Todd Trimble
    Oct 17 '16 at 17:49

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.