[Racket] Ferti hydroponic nutrient solver, redux.
Greek letters ftw.
services/nnls.rkt
@@ -214,12 +214,12 @@
214
214
(test-case "Two variables, known solution"
215
215
(define A (matrix [[1 0] [0 1]]))
216
216
(define y (col-matrix [3 4]))
217
Removed:
(define epsilon 1e-6)
217
Added:
(define ε 1e-6)
218
218
219
Removed:
(define result (lawson-hanson-1974 A y epsilon))
219
Added:
(define result (lawson-hanson-1974 A y ε))
220
220
221
Removed:
(check-= (matrix-ref result 0 0) 3.0 epsilon "x1 should be 3")
222
Removed:
(check-= (matrix-ref result 1 0) 4.0 epsilon "x2 should be 4"))
221
Added:
(check-= (matrix-ref result 0 0) 3.0 ε "x1 should be 3")
222
Added:
(check-= (matrix-ref result 1 0) 4.0 ε "x2 should be 4"))
223
223
224
224
(test-case "Overdetermined system"
225
225
(define A (matrix [[1 1] [2 1] [1 2]]))
@@ -235,9 +235,9 @@
235
235
(test-case "Non-negativity enforcement"
236
236
(define A (matrix [[1 -1] [1 1]]))
237
237
(define y (col-matrix [1 3]))
238
Removed:
(define epsilon 1e-6)
238
Added:
(define ε 1e-6)
239
239
240
Removed:
(define result (lawson-hanson-1974 A y epsilon))
240
Added:
(define result (lawson-hanson-1974 A y ε))
241
241
242
242
;; All results should be non-negative
243
243
(check-true (>= (matrix-ref result 0 0) 0) "x1 should be non-negative")