dSFMT_init Subroutine

    subroutine dSFMT_init(seed)

        ! Initialise the dSFMT RNG and fill random_store with
        ! a block of random numbers in interval [0,1).
        !
        ! In:
        !    seed: seed for the RNG.

        integer, intent(in) :: seed

        call init_gen_rand_fwrapper(int(seed,c_int32_t))

        call fill_array_close_open_fwrapper(random_store, &
                                            int(random_store_size, c_int))

        current_element = 1

    end subroutine dSFMT_init