bugGNU Octave - Bugs: bug #60650, [octave forge] (linear-algebra)...

 
 

bug #60650: [octave forge] (linear-algebra) nmf_bpas fails for rank k=1

Submitter:  Tasos Papastylianou <tpapastylianou>
Submitted:  Thu 20 May 2021 05:41:15 PM UTC
   
 
Category:  Octave Package Severity:  3 - Normal
Priority:  5 - Normal Item Group:  Unexpected Error or Warning
Status:  None Assigned to:  None
Originator Name:  Tasos Papastylianou Open/Closed:  * Open
Release:  * 6.2.0 Operating System:  * Any
Fixed Release:  None Planned Release:  None
* Mandatory Fields

Add a New Comment Rich Markup
   

Thu 20 May 2021 05:41:15 PM UTC, original submission:  

To reproduce:


pkg load linear-algebra
V=rand(4, 20);
[w, h, i, s ]=nmf_bpas(V,1);


see https://stackoverflow.com/q/67618375/4183191 for context

I have narrowed down the culprit to line 373:


pGrad = [gradW(gradW<0|W>0); gradH(gradH<0|H>0)];


which should change to something like:


pGrad = [gradW(gradW<0|W>0); gradH(gradH<0|H>0)(:)];


This is because when the matrices involved are one-dimensional, logical indexing may result in a horizontal vector instead of a vertical one, which breaks the concatenation there.

Tasos Papastylianou <tpapastylianou>

 

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

Attach Files:
   
   
Comment:
   

No files currently attached

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -email is unavailable- added by tpapastylianou (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 2 latest changes.

    Date Changed by Updated Field Previous Value => Replaced by
    2021-05-21 mmuetzel Carbon-CopyRemoved 102357 -
    2021-05-21 mmuetzel Summarynmf_bpas from linear-algebra package fails for rank k=1 [octave forge] (linear-algebra) nmf_bpas fails for rank k=1

    Back to the top

    Powered by Savane 3.13-f8d8.
    Corresponding source code