bugGNU Octave - Bugs: bug #55831, Incorrect sparse QR factorization

 
 

bug #55831: Incorrect sparse QR factorization

Submitter:  None
Submitted:  Tue 05 Mar 2019 07:37:57 AM UTC
   
 
Category:  Octave Function Severity:  3 - Normal
Priority:  5 - Normal Item Group:  Incorrect Result
Status:  Duplicate 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
   

Wed 06 Mar 2019 08:03:46 PM UTC, comment #2: 

Yes, marking as a duplicate and closing report.

Rik <rik5>
Group administrator
Wed 06 Mar 2019 02:41:02 PM UTC, comment #1: 

I think it is a duplicate (or related to) bug #51950.

Marco Caliari <caliari>
Group Member
Tue 05 Mar 2019 07:37:57 AM UTC, original submission:  

The following code:
[Q,R] = qr(sparse([1,0;0,0]));
outputs
Q = [-1 0
      0 0]

obviously Q is not unitary
it seems to me that each column with all zeros in a sparse input return a column with all zeros in the Q matrix

the size of Q is also incorrect when the input is sparse and an economy size factorization is requested. e.g.:
[Q,R] = qr(sparse(5,2),0);
outputs Q with size 5x5 rather than 5x2. Even though Q is 5x5, R is 2x2.

btw, is there a reason why Octave does not return a sparse Q like Matlab does when the input is sparse? I know Q is usually dense, but in some cases it may be sparse (e.g., economy factorization of a tall and sparse matrix). This could be useful to avoid running out of memory, for instance.

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 caliari (Posted a comment)
  •  

    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 4 latest changes.

    Date Changed by Updated Field Previous Value => Replaced by
    2019-03-06 rik5 StatusNone Duplicate
        Open/ClosedOpen Closed
        Dependencies- Depends on bugs #51950
    2019-03-06 caliari Operating SystemMicrosoft Windows Any

    Back to the top

    Powered by Savane 3.13-758e.
    Corresponding source code