// UniformRandomBitGenerator compient interface:
rnd::minstd_rand::minstd_rand(seed = default_seed);
uint32_t rnd::minstd_rand::max();
uint32_t rnd::minstd_rand::min();
uint32_t rnd::minstd_rand::operator();
// extended interface for generating chuncks of random uniformed distributed float values:
std::unique_ptr<float> rnd::minstd_rand::operator()(size_t num);
void rnd::minstd_rand::operator()(size_t num, float* gen_val);

size_t points_inside_cirl = 0;

std::minstd_rand rgen{rseed};
std::uniform_real_distribution<float> udistr{0.0, 1.0};

constexpr size_t vec_size = 8;
constexpr size_t chuck = 10'000;
constexpr size_t sub_iters = kIters / chuck;

float* val_ptr = new float[2 * chuck];

for (size_t outer_ind = 0; outer_ind < sub_iters; outer_ind++) {
    for (size_t i = 0; i < 2 * chuck; i++) {
        val_ptr[i] = udistr(rgen);
    }

    for (size_t i = 0; i < 2 * chuck; i += 2 * vec_size) {
        Vec8x32f vec_x{&val_ptr[i]};
        Vec8x32f vec_y{&val_ptr[i + vec_size]};

        Vec8x32f sqr = vec_x * vec_x + vec_y * vec_y;
        Vec8x32f one_vec = Vec8x32f{1.0f};

        uint32_t cmp_mask = sqr <= one_vec;
        points_inside_cirl += CountOnes(cmp_mask);
    }
}

delete[] val_ptr;

float pi = points_inside_cirl / total_num_points;

size_t points_inside_cirl = 0;

rnd::minstd_rand rgen{rseed};

constexpr size_t vec_size = 8;
constexpr size_t chuck = 10'000;
constexpr size_t sub_iters = kIters / chuck;

float* val_ptr = new float[2 * chuck];

for (size_t outer_ind = 0; outer_ind < sub_iters; outer_ind++) {
    rgen(2 * chuck, val_ptr);

    for (size_t i = 0; i < 2 * chuck; i += 2 * vec_size) {
        Vec8x32f vec_x{&val_ptr[i]};
        Vec8x32f vec_y{&val_ptr[i + vec_size]};

        Vec8x32f sqr = vec_x * vec_x + vec_y * vec_y;
        Vec8x32f one_vec = Vec8x32f{1.0f};

        uint32_t cmp_mask = sqr <= one_vec;
        points_inside_cirl += CountOnes(cmp_mask);
    }
}

delete[] val_ptr;

float pi = points_inside_cirl / total_num_points;