> your approach gives a non uniform distribution Unless I made some significant error, it is a uniform distribution over multiples of 2^-52. It's also a discre

> your approach gives a non uniform distribution Unless I made some significant error, it is a uniform distribution over multiples of 2^-52. It's also a discrete approximation of the real uniform distribution in that for any subdivision of [0,1) into N equal (half-open) intervals, the most-populated interval has at most one more inhabitant the least populated one. It is however strictly inferior to rand53()*2^-53 (NextDoubleInner at your github link) unless your int64->double or float-multiply operations are quite bad.