This is an extremely basic question and I have to be missing something, but when trying to normalize matrix A below so that rows sum to one, some (small) differences remain: A = rand(100,4) A = A./sum(A, dims = 2) sum(sum(A, dims =2) .== 1) # fewer than 100! Why does this happen? Is there a wa…

The result is due to floating point inaccuracies. The following holds: A = rand(100,4) A = A./sum(A, dims = 2) @assert sum(sum(A, dims =2) .≈ 1) == 100 Alternate name for ≈ is isapprox.

Normalizing a vector to sum to one does not work?

General Usage

Sukera October 11, 2021, 12:41pm 8

I think this is relevant here:

5 Likes

Topic		Replies	Views
Custom error to check stochastic matrix General Usage	5	268	October 25, 2022
Rant: I hate that it's possible to iterate over an integer New to Julia	41	3669	March 13, 2022
X[a,1]'b[a]==sum(x[a,1].b[a])? General Usage	2	279	June 7, 2021
Sum of float64 vector gives slightly incorrect answer Performance question	44	4986	April 29, 2019
Discussion about integer overflow Internals & Design numbers , integer-overflow	189	7261	October 15, 2021

Normalizing a vector to sum to one does not work?

Related topics