bugGNU Octave - Bugs: bug #56026, Incorrect output of generalized...

 
 

bug #56026: Incorrect output of generalized eigenvalue problem

Submitter:  None
Submitted:  Fri 29 Mar 2019 05:06:42 AM UTC
   
 
Category:  Octave Function Severity:  3 - Normal
Priority:  5 - Normal Item Group:  Incorrect Result
Status:  Fixed Assigned to:  None
Originator Name:  Originator Email:  -email is unavailable-
Open/Closed:  * Closed Release:  * 5.1.0
Operating System:  * Any Fixed Release:  None
Planned Release:  None
* Mandatory Fields

Add a New Comment Rich Markup
   

Mon 01 Apr 2019 09:49:22 PM UTC, comment #4: 

I fixed this on the stable branch here (https://hg.savannah.gnu.org/hgweb/octave/rev/fc73dafece57).

Closing report.

 Rik <rik5>
Group administrator
Mon 01 Apr 2019 07:53:22 PM UTC, comment #3: 

That makes sense.  It should follow the same syntax as for calc_rev.  So


  F77_XFCN (zggev, ZGGEV, (F77_CONST_CHAR_ARG2 ("N", 1),
                           F77_CONST_CHAR_ARG2 (calc_rev ? "V" : "N", 1),


should be


  F77_XFCN (zggev, ZGGEV, (F77_CONST_CHAR_ARG2 (calc_lev ? "V" : "N", 1),
                           F77_CONST_CHAR_ARG2 (calc_rev ? "V" : "N", 1),



 Rik <rik5>
Group administrator
Mon 01 Apr 2019 07:38:51 PM UTC, comment #2: 

The problem is in the first argument passed to zggev (EIG.cc). It should be


F77_CONST_CHAR_ARG2 (calc_lev ? "V" : "N", 1)


(two occurrences).

 Marco Caliari <caliari>
Group Member
Fri 29 Mar 2019 03:20:47 PM UTC, comment #1: 

Confirmed.  I've attached a simple test script to illustrate the problem.


A = magic(4);
B = complex (magic(4), magic(4));
[v,l,w] = eig (A,B)


I tested back to Octave version 4.2.1 and this problem has always been present.  Apparently there was not a test for this in the test suite.


(file #46661)

 Rik <rik5>
Group administrator
Fri 29 Mar 2019 05:06:42 AM UTC, original submission:  

Calling eig function with the following form

[V, LAMBDA, W] = eig (A, B)

returns a zero matrix on W when the inputs are complex.

This problem doesn't occur with qz function

 Anonymous

 

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

Attach Files:
   
   
Comment:
   

Attached Files
file #46661:  tst_eig.m added by rik5 (69B - text/x-matlab)

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -email is unavailable- added by caliari (Posted a comment)
  • -email is unavailable- added by rik5 (Updated the item)

    •  

    There are 0 votes so far. Votes easily highlight which items people would like to see resolved in priority, independently of the priority of the item set by tracker managers.

    Only group members can vote.

     

    Follow 5 latest changes.

    Date Changed by Updated Field Previous Value => Replaced by
    2019-04-01 rik5 StatusConfirmed Fixed
        Open/ClosedOpen Closed
    2019-03-29 rik5 Attached File- Added tst_eig.m, #46661
        StatusNone Confirmed
        Operating SystemMicrosoft Windows Any

    Back to the top

    Copyright © 2024 Free Software Foundation, Inc.
    Verbatim copying and distribution of this entire article is permitted in any medium, provided this notice is preserved.
    The Levitating, Meditating, Flute-playing Gnu logo is a GNU GPL'ed image provided by the Nevrax Design Team.
    Page source code

    Powered by Savane 3.14-f13d.
    Corresponding source code