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

 
 

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

Submitted by:  Tasos Papastylianou <tpapastylianou>
Submitted on:  Thu 20 May 2021 05:41:15 PM UTC  
 
Category:  Octave Forge 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

Add a New Comment (Rich Markup)
   

You are not logged in

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

 

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)
  •  

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

    Only project 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.6