#include <iostream>
#include "lxm_random/lxm_random.hpp"

int main() {
  // Initialize an RNG with a random seed
  uint64_t seed = 42;
  lxm_random::lxm_random rng(seed);

  // Generate a random number
  std::cout << rng() << std::endl;
}

// Create a root random number generator with a seed
lxm_random::lxm_random rng(seed);

std::vector<std::thread> threads;
for (int i = 0; i < n_threads; i++) {
  // Fork a new thread and pass a new random number generator to it
  // by splitting the root random number generator
  threads.emplace_back([=, rng = rng.split()]() mutable {
    ...
    rng(); // Generate a random number independently from other threads
    ...
  });
}

// Create a root random number generator with a seed
lxm_random::lxm_random rng(seed);

std::vector<std::thread> threads;
for (int i = 0; i < n_threads; i++) {
  // Fork a new thread
  threads.emplace_back([=]() mutable {
    // Split the copied random number generator by giving a unique thread ID
    // so that threads with different IDs will have independent streams
    rng = rng.split(i);
    ...
    rng(); // Generate a random number independently from other threads
    ...
  });
}