[julia-users] Julia non-negative least squares nnls is 100x slower than R nnls

[David C Cohen]

Hi,

I understand that Julia and R nnls do not use the same algorithm. A 400x300 
problem on my machine takes about 0.02 s on R nnls, and about 2 s on Julia 
nnls.

The documentations on nnls is rather poor:

?nnls

nnls (generic function with 2 methods)


Am I missing something here, or is this implementation is meant to be slow?