bugGNU Octave - Bugs: bug #30626, QZ for complex matrix

 
 

bug #30626: QZ for complex matrix

Submitted by:  roucaries <broucaries>
Submitted on:  Mon 02 Aug 2010 07:21:16 AM UTC  
 
Category: InterpreterSeverity: 3 - Normal
Priority: 5 - NormalItem Group: Feature Request
Status: FixedAssigned to: None
Originator Name: Open/Closed: Closed
Release: 3.2.4Operating System: GNU/Linux

Add a New Comment(Rich Markup)
   

You are not logged in

Please log in, so followups can be emailed to you.

 

Mon 02 Aug 2010 11:08:00 PM UTC, comment #1:

I think this is already done in the current development sources:

octave:1> a = rand (3) + i * rand (3);
octave:2> b = rand (3) + i * rand (3);
octave:3> [aa, bb, q, z, v, w, lambda] = qz (a, b)
aa =

1.47018 + 0.66439i -1.35595 + 0.53104i 0.77763 + 0.09018i
0.00000 + 0.00000i 0.13679 - 0.18679i -0.35090 - 0.34734i
0.00000 + 0.00000i 0.00000 + 0.00000i -0.27090 + 0.03004i

bb =

0.86253 + 0.00000i -1.44463 + 0.40366i 0.11613 - 0.19049i
0.00000 + 0.00000i 0.99733 + 0.00000i -0.14376 - 0.12654i
0.00000 + 0.00000i 0.00000 + 0.00000i 0.42306 + 0.00000i

q =

0.32665 - 0.53270i 0.42344 - 0.40556i 0.29795 - 0.42069i
-0.28626 - 0.65696i 0.22930 + 0.63063i -0.09847 + 0.16272i
0.09697 + 0.29425i 0.32063 + 0.32115i -0.52100 - 0.65317i

z =

0.11427 - 0.40264i -0.66828 + 0.01684i 0.59865 + 0.13985i
0.29879 + 0.40570i -0.57728 + 0.43868i -0.46132 + 0.08729i
0.71376 + 0.24798i 0.13971 - 0.08894i 0.31717 - 0.54866i

v =

0.11882 - 0.41865i -0.74094 - 0.25906i 0.48385 - 0.34211i
0.31068 + 0.42184i -0.17372 + 0.48285i -0.85294 + 0.14706i
0.74215 + 0.25785i 0.63622 - 0.28382i 0.50845 - 0.04131i

w =

0.050783 + 0.949217i 0.018515 + 0.515701i 0.082583 - 0.250600i
0.583229 - 0.324435i 0.258734 - 0.034069i 0.273073 - 0.273515i
0.158342 + 0.528676i -0.360012 - 0.639988i -0.443721 + 0.556279i

lambda =

1.70450 + 0.77029i
0.13716 - 0.18729i
-0.64033 + 0.07101i

John W. Eaton <jwe>
Project Administrator
Mon 02 Aug 2010 07:21:16 AM UTC, original submission:

Hi,

Qz is not implemented for complex matrix. I could help you to implement it, i will soon post patch

It is really important in linear algebra to have qz working for complex case

roucaries <broucaries>

 

(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

Attach File(s):
   
   
Comment:
   

Attached Files
file #21126:  diff added by broucaries (45KiB - application/octet-stream - Patch to use complex using testing version)

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -unavailable- added by jwe (Posted a comment)
  • -unavailable- added by broucaries (Submitted the item)
  •  

    Do you think this task is very important?
    If so, you can click here to add your encouragement to it.
    This task has 0 encouragements so far.

    Only project members can vote.

     

    Please enter the title of George Orwell's famous dystopian book (it's a date):

     

     

    Follow 3 latest changes.

    Date Changed By Updated Field Previous Value => Replaced By
    Mon 02 Aug 2010 11:08:00 PM UTCjweStatusNone=>Fixed
      Open/ClosedOpen=>Closed
    Mon 02 Aug 2010 12:45:07 PM UTCbroucariesAttached File-=>Added diff, #21126

    Back to the top


    Powered by Savane 3.1-cleanup1