come generare il numero casuale intero in fortran 90 nell'intervallo [0,5]?

Question

Sono una specie di nuovo nel fortran proramming. Qualcuno può aiutarmi con la soluzione. Sto avendo un problema di generare numero intero casuale nel range [0,5] in numero casuale FORTRAN utilizzando random_seed e rand

Answer 1

Che dire:

program rand_test 
    use,intrinsic :: ISO_Fortran_env 
    real(REAL32) :: r(5) 
    integer  :: i(5) 

    ! call init_random_seed() would go here 

    call random_number(r) 

    ! Uniform distribution requires floor: Thanks to @francescalus 
    i = floor(r*6._REAL32) 

    print *, i 
end program 

Answer 2

Per sostenere la answer da Alexander Vogt , Generalizzerò.

L'intrinseco random_number(u) restituisce un numero reale u (o un array di tale) dalla distribuzione uniforme nell'intervallo [0,1). [Ovvero, include 0 ma non 1.]

Per avere una distribuzione uniforme discreta sugli interi {n, n + 1, ..., m-1, m} intagliare la distribuzione continua in m + 1-n pezzi di dimensioni uguali, mappando ogni pezzo in un numero intero. Un modo potrebbe essere:

call random_number(u) 
j = n + FLOOR((m+1-n)*u) ! We want to choose one from m-n+1 integers 

Come si può vedere, per la domanda iniziale per {0, 1, 2, 3, 4, 5} Questo riduce al

call random_number(u) 
j = FLOOR(6*u)   ! n=0 and m=5 

e per l'altro caso nel tuo commento {-1, 0, 1}

call random_number(u) 
j = -1 + FLOOR(3*u)  ! n=-1 and m=1 

Naturalmente, altre trasformazioni saranno necessari per i set di numeri interi non contigui, e si dovrebbe prestare attenzione ai problemi numerici.