suzanne.soy/racket
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Fix math/matrix for new type of expt.

Browse Source
This commit is contained in:
Vincent St-Amour 2014-11-05 10:21:29 -05:00
parent 4afd6c86d0
commit 8a0c859d14
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
pkgs/math-pkgs/math-lib/math/private/matrix/matrix-basic.rkt
View File

@ -238,8 +238,8 @@
;; Compute this divided by the maximum to avoid underflow and overflow
(define mx (array-all-max M))
(cond [(and (rational? mx) (positive? mx))
(* mx (expt (array-all-sum (inline-array-map (λ (x) (expt (abs (/ x mx)) p)) M))
(/ p)))]
(fl (* mx (expt (array-all-sum (inline-array-map (λ (x) (expt (abs (/ x mx)) p)) M))
(/ p))))]
[else mx]))))
(: matrix-norm (case-> ((Matrix Flonum) -> Nonnegative-Flonum)