bugGNU Octave - Bugs: bug #47095, Wrong generalized eigenvalue...

 
 

bug #47095: Wrong generalized eigenvalue solutions for hermitian matrices

Submitter:  None
Submitted:  Tue 09 Feb 2016 12:13:37 AM UTC
   
 
Category:  Octave Function Severity:  3 - Normal
Priority:  5 - Normal Item Group:  Incorrect Result
Status:  Duplicate Assigned to:  None
Originator Name:  Mohammad Elahi Originator Email:  -email is unavailable-
Open/Closed:  * Closed Release:  * 4.0.0
Operating System:  * Any Fixed Release:  None
Planned Release:  None
* Mandatory Fields

Add a New Comment Rich Markup
   

Tue 09 Feb 2016 12:47:18 AM UTC, comment #1: 

Duplicate of bug #45511.

This has already been fixed on the development branch.  Closing report.

Rik <rik5>
Group administrator
Tue 09 Feb 2016 12:13:37 AM UTC, original submission:  

The results of generalized eigenvalue problem solved by 'eig' function are incorrect when the input matrices are hermitian.
A=hermitian matrix
B=hermitian matrix
eig(A,B).

The bug exists in 'libinterp/corefcn/EIG.cc':

octave_idx_type
SEIG::init (const ComplexMatrix& a, const ComplexMatrix& b, bool calc_ev)
{
...
  if (a.is_hermitian () && b.is_hermitian () && info == 0)
    return hermitian_init (a, calc_ev);
...
}

So the hermitian_init is called without 'b' which is incorrect. The trivial fix is the following:

octave_idx_type
SEIG::init (const ComplexMatrix& a, const ComplexMatrix& b, bool calc_ev)
{
...
  if (a.is_hermitian () && b.is_hermitian () && info == 0)
    return hermitian_init (a, b, calc_ev);
...
}


Anonymous

 

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

Attach Files:
   
   
Comment:
   

No files currently attached

 

Digest:
   bug dependencies.

Items that depend on this one: None found

 

Carbon-Copy List
  • -email is unavailable- added by rik5 (Posted a comment)
  • -email is unavailable- added by None (Submitted 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 3 latest changes.

    Date Changed by Updated Field Previous Value => Replaced by
    2016-02-09 rik5 StatusNone Duplicate
        Open/ClosedOpen Closed
    2016-02-09 rik5 Dependencies- Depends on bugs #45511

    Back to the top

    Powered by Savane 3.13-02a9.
    Corresponding source code