Issues about the module "m_random"

[leonfoks]

Issues about the module "m_random"

[leonfoks]

I just confirmed this, thanks for the report! I'll see what I can come up with!

[leonfoks]

Okay, it seems I have an issue with my "global prng" and how that is being set. If I use a local instantiation of the Prng class, it works fine for me (see below).

Using a local instance is in my opinion preferred as it gives you access to two different generators, and the seed for reproducibility purposes. The Prng class is also thread safe if you wish to use it in parallel and you can jump the state to ensure that each core (or rank in MPI) is getting a different set of random numbers.

I will work on fixing the other approach using a global prng, so you can either switch to this class approach or wait!

Thanks again!

Program main
    ! Import modules
    Use variableKind, only: r64, i32
    use Prng_Class, only: Prng

    Implicit None

    ! Data dictionary: declare parameters

    real(r64) :: r
    real(r64), dimension(10) :: rA

    integer(i32) :: i
    integer(i32), dimension(10) :: iA

    type(Prng) :: rng

    rng = Prng(big = .true., display = .true.)

    call rng%rngUniform(r, 0.d0, 1.d0)
    write(*,*) r
    call rng%rngUniform(rA, 0.d0, 1.d0)
    write(*,*) rA

    call rng%rngNormal(r)
    write(*,*) r
    call rng%rngNormal(rA)
    write(*,*) rA

    call rng%rngInteger(i, 0, 10)
    write(*,*) i

    call rng%rngInteger(iA, 0, 10)
    write(*,*) iA
End Program

[WangYun1995]

You are welcome. BTW, 'coretran' is so cool!

[leonfoks]

Thank you!!