[Racket] Ferti hydroponic nutrient solver, redux.
Fix matrix scaling error.
services/nnls.rkt
@@ -155,7 +155,7 @@
155
155
(error 'lawson-hanson-1974 "no valid α in inner loop"))
156
156
157
157
;; x ← x + α (s − x)
158
Removed:
(define new-x (matrix+ x (matrix* α (matrix- s x))))
158
Added:
(define new-x (matrix+ x (matrix-scale (matrix- s x) α)))
159
159
160
160
;; Move to R all indices j in P with x_j <= 0
161
161
(define to-remove '())