crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-26 05:05:30 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Address #939 comments

Browse source
This commit is contained in:
Sterling Harper 2017-12-13 18:49:12 -05:00
parent ec302b52d6
commit f6d6b926b1
2 changed files with 7 additions and 6 deletions
Download patch file Download diff file Expand all files Collapse all files

2
CMakeLists.txt
View file

@ -417,7 +417,7 @@ set(LIBOPENMC_FORTRAN_SRC
)
set(LIBOPENMC_CXX_SRC
src/random_lcg.h
src/random_lcg.C)
src/random_lcg.cpp)
add_library(libopenmc SHARED ${LIBOPENMC_FORTRAN_SRC} ${LIBOPENMC_CXX_SRC})
set_target_properties(libopenmc PROPERTIES OUTPUT_NAME openmc)
add_executable(${program} src/main.F90)

11
src/random_lcg.C → src/random_lcg.cpp
View file

@ -21,7 +21,7 @@ constexpr uint64_t prn_mod {0x8000000000000000}; // 2^63
constexpr uint64_t prn_mask {0x7fffffffffffffff}; // 2^63 - 1
constexpr uint64_t prn_stride {152917LL}; // stride between
// particles
constexpr double prn_norm {pow(2, -63)}; // 2^-63
constexpr double prn_norm {1.0 / prn_mod}; // 2^-63
// Current PRNG state
uint64_t prn_seed[N_STREAMS]; // current seed
@ -135,10 +135,11 @@ openmc_set_seed(int64_t new_seed)
{
seed = new_seed;
#pragma omp parallel
for (int i = 0; i < N_STREAMS; i++) {
prn_seed[i] = seed + i;
{
for (int i = 0; i < N_STREAMS; i++) {
prn_seed[i] = seed + i;
}
prn_set_stream(STREAM_TRACKING);
}
prn_set_stream(STREAM_TRACKING);
#pragma end omp parallel
return 0;
}