Fix matrix scaling error.

1 files changed, 1 insertions, 1 deletions

diff --git a/services/nnls.rkt b/services/nnls.rkt
index 16703d2..2c5738f 100644
--- a/services/nnls.rkt
+++ b/services/nnls.rkt
@@ -155,7 +155,7 @@
                  (error 'lawson-hanson-1974 "no valid α in inner loop"))
 
                ;; x ← x + α (s − x)
-               (define new-x (matrix+ x (matrix* α (matrix- s x))))
+               (define new-x (matrix+ x (matrix-scale (matrix- s x) α)))
 
                ;; Move to R all indices j in P with x_j <= 0
                (define to-remove '())