Fiz o seguinte lcg (gerador congruencial linear) para sortear unsigned int de 0x00000000
a 0xFFFFFFFF
uniformemente. Só para teste, usei a semente igual a zero e sorteando os cinco primeiros números para ver o que sai.
# include <stdio.h>
unsigned long long int state = 0 ;
unsigned int rnd() {
state *= 0x60BDC8431972EFA5ull ;
state += 0x4D268FBC9E53A107ull ;
return (unsigned int)( state >> 32 ) ;
}
int main(void) {
printf( "%u\n" , rnd() ) ;
printf( "%u\n" , rnd() ) ;
printf( "%u\n" , rnd() ) ;
printf( "%u\n" , rnd() ) ;
printf( "%u\n" , rnd() ) ;
return 0 ;
}
O resultado é esse.
1294372796
3162674411
1090754839
3480286764
3255825657
Para sortear um número em formato single floating point dentro do intervalo [a,b]
eu fiz a seguinte função que simplesmente divide o inteiro sorteado por 0xFFFFFFFFu
para funcionar como sorteio de 0.0f
a 1.0f
e usa o resultado para fazer lerp de a
a b
.
float random( float a , float b ){
float f0t1 = rnd()/(float)(0xFFFFFFFFu) ;
return a * ( 1-f0t1 ) + f0t1 * b ;
}
Testando assim,
int main(void) {
printf( "%f\n" , random(0,1) ) ;
printf( "%f\n" , random(0,1) ) ;
printf( "%f\n" , random(0,1) ) ;
printf( "%f\n" , random(0,1) ) ;
printf( "%f\n" , random(0,1) ) ;
return 0 ;
}
o resultado é esse.
0.301370
0.736368
0.253961
0.810317
0.758056
Quer dizer, faz uma divisão que custa muito. Até dá para converter em multiplicação, mas não há um meio mais rápido de calcular isso? Imagino que com números mágicos daria para manusear o código binário do float de maneira eficiente e calcular, correto?