Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
1998
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 1998

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

Search the Archive

Re: JordanDecomposition trouble


  • To: mathgroup@smc.vnet.net
  • Subject: [mg10703] Re: JordanDecomposition trouble
  • From: Paul Abbott <paul@physics.uwa.edu.au>
  • Date: Fri, 30 Jan 1998 04:24:38 -0500
  • Organization: University of Western Australia
  • References: <6ahmmp$pfa@smc.vnet.net>

Josip Loncaric wrote:
 
> Let {s,j} denote Jordan decomposition of the following exact matrix:
> 
> m = {{0,1,0,1},{0,0,0,0},{0,1,0,1},{0,0,0,0}} 
> {s,j} = JordanDecomposition[m]
> 
> We should have m == s.j.Inverse[s], but in Mathematica 3.01 this is not
> the case.
> Instead, we get a rather different result:
> 
> In[3]:= s.j.Inverse[s]
> Out[3]= {{0,1,0,2},{0,0,0,0},{0,1,0,2},{0,0,0,0}}
> 
> The problem is not in matrix inversion since s.Inverse[s] returns
> identity.  The source of trouble appears to be in the
> JordanDecomposition result.

Certainly looks that way.  Since Jordan decomposition works with
symbolic parameters, looking at

In[1]:= m = {{0, 1, 0, 1}, {0, 0, 0, 0}, {0, 1, 0, 1}, {0, 0, 0, a}}; 

In[2]:= {s, j} = JordanDecomposition[m]

In[3]:= s.j.Inverse[s]==m
Out[3]= True

In[4]:= s
Out[4]=
	           1                            1
	{{0, 1, 0, -}, {0, 0, 1, 0}, {-1, 1, 0, -}, {0, 0, 0, 1}}
	           a                            a

may give a hint as to why the Jordan decomposition failed.

Cheers,
	Paul 

____________________________________________________________________ 
Paul Abbott                                   Phone: +61-8-9380-2734
Department of Physics                           Fax: +61-8-9380-1014
The University of Western Australia            Nedlands WA  6907       
mailto:paul@physics.uwa.edu.au  AUSTRALIA                            
http://www.pd.uwa.edu.au/~paul

            God IS a weakly left-handed dice player
____________________________________________________________________



  • Prev by Date: Bug? Can suddenly not edit with LINUX Frontend
  • Next by Date: Re: Lists and Plus
  • Prev by thread: Re: JordanDecomposition trouble
  • Next by thread: Molecular dynamics with Mathematica?