Function FillRandom

Src/Base/AMReX_Random.cpp:263–296

}

void FillRandom (Real* p, Long N)
{
#ifdef AMREX_USE_CUDA

#  ifdef BL_USE_FLOAT
    AMREX_CURAND_SAFE_CALL(curandGenerateUniform(gpu_rand_generator, p, N));
#  else
    AMREX_CURAND_SAFE_CALL(curandGenerateUniformDouble(gpu_rand_generator, p, N));
#  endif
    Gpu::synchronize();

#elif defined(AMREX_USE_HIP)

#  ifdef BL_USE_FLOAT
    AMREX_HIPRAND_SAFE_CALL(hiprandGenerateUniform(gpu_rand_generator, p, N));
#  else
    AMREX_HIPRAND_SAFE_CALL(hiprandGenerateUniformDouble(gpu_rand_generator, p, N));
#  endif
    Gpu::synchronize();

#elif defined(AMREX_USE_SYCL)

    oneapi::mkl::rng::uniform<Real> distr;
    auto event = oneapi::mkl::rng::generate(distr, *gpu_rand_generator, N, p);
    event.wait();

#else
    std::uniform_real_distribution<Real> distribution(Real(0.0), Real(1.0));
    auto& gen = generators[OpenMP::get_thread_num()];
    for (Long i = 0; i < N; ++i) {
        p[i] = distribution(gen);
    }
#endif
}

void FillRandomNormal (Real* p, Long N, Real mean, Real stddev)
{