MathGroup Archive 2012

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: JordanDecompose[]

  • To: mathgroup at
  • Subject: [mg124024] Re: JordanDecompose[]
  • From: danl at
  • Date: Fri, 6 Jan 2012 04:15:19 -0500 (EST)
  • Delivered-to:
  • References: <je3vti$eeh$>
  • Reply-to: comp.soft-sys.math.mathematica at

Here is a mildly edited comment from the code.

"This routine uses the method based on matrix presentations.
   It is a more general method than the one used in [...]
   (the latter is faster but does not always work, for derogatory matrices).
   Initially the matrix is decomposed into subblocks, where each
   subblock has eigenvalues of the same value.  Eigenvalues in
   different subblocks are distinct.  Then elementary row and column
   operations are done on the characteristic presentations of the

This was from 1993, when two of us worked on the implementation. It was based on course notes I took in 1985. (Coincidently, the class was taught by one of the original Mathematica developers, Dan Grayson.) It's been a while since bugs were spotted in that code, so I have to confess I do not remember the details.

Daniel Lichtblau
Wolfram Research

  • Prev by Date: SortBy won't order irrationals
  • Next by Date: DirectoryStack[] not showing top directory stack item
  • Previous by thread: Re: SortBy won't order irrationals
  • Next by thread: Re: JordanDecompose[]