Why do you convert a single-precision random number to double precision
when you can generate a double-precision random number directly?
!---------------------------------------------!
function drand()
implicit none
real(8) :: r
call random_number(r)
drand = r
end function drand
!---------------------------------------------!
Cheers, Wes
On Wed, 16 Nov 2011, Adriano Amaricci wrote:
> !---------------------------------------------!
> function drand()
> implicit none
> real(8) :: drand
> real(4) :: r
> call random_number(r)
> drand=dble(r)
> end function drand
> !---------------------------------------------!
|