Fix parallel rng

[vandenman]

tested with

library(bgms)
x = Wenchuan[1:50, 1:5]
p = ncol(x)

samples1 = bgm(x, iter = 1000, nuts_max_depth = 10, cores = 4, chains = 4, seed = 1234)
samples2 = bgm(x, iter = 1000, nuts_max_depth = 10, cores = 4, chains = 4, seed = 1234)

print(all.equal(samples1$raw_samples$main,     samples2$raw_samples$main))
print(all.equal(samples1$raw_samples$pairwise, samples2$raw_samples$pairwise))

prints TRUE for both (on my machine).

[MaartenMarsman]

Comment