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

StefanKarpinski Split this topic October 12, 2021, 2:18pm 10

A post was merged into an existing topic: Discussion about integer overflow

Topic		Replies	Views
Custom error to check stochastic matrix General Usage	5	267	October 25, 2022
Rant: I hate that it's possible to iterate over an integer New to Julia	41	3661	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	4978	April 29, 2019
Discussion about integer overflow Internals & Design numbers , integer-overflow	189	7243	October 15, 2021

Normalizing a vector to sum to one does not work?

Related topics