fltx/examples/example_random.cpp at main · bitloop-dev/fltx · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
#include <iostream>
#include <iomanip>
#include <limits>
#include <array>

#include <fltx.h>

using namespace bl;
using namespace bl::literals;

template<std::size_t N, typename T>
consteval T avg(const std::array<T, N>& values)
{
    T total = 0;
    for (const auto& v : values)
        total += v;
    return total / N;
}

int main()
{
    // fqd: random uniform arrays (min=0, max=100)
    constexpr auto a = bl::uniform_real_array<100>(0_qd, 100_qd, bl::mt19937_64{ 0x123ull }); // explicit engine / custom seed
    constexpr auto b = bl::uniform_real_array<100>(0_qd, 100_qd, 0x123ull);                   // default engine  / custom seed
    constexpr auto c = bl::uniform_real_array<100>(0_qd, 100_qd);                             // default engine  / default seed

    // fdd: random normal array (mean=50, stddev=50)
    constexpr auto d = bl::normal_array<100>(50_dd, 50_dd, bl::mt19937_64{ 0x123ull }); // explicit engine / custom seed
    constexpr auto e = bl::normal_array<100>(50_dd, 50_dd, 0x123ull);                   // default engine  / custom seed
    constexpr auto f = bl::normal_array<100>(50_dd, 50_dd);                             // default engine  / default seed

    constexpr fqd a_avg = avg(a);
    constexpr fqd b_avg = avg(b);
    constexpr fqd c_avg = avg(c);

    constexpr fdd d_avg = avg(d);
    constexpr fdd e_avg = avg(e);
    constexpr fdd f_avg = avg(f);

    std::cout
        << std::fixed
        << std::setprecision(std::numeric_limits<fqd>::digits10)
        << "a_avg: " << a_avg << "\n"
        << "b_avg: " << b_avg << "\n"
        << "c_avg: " << c_avg << "\n\n"
        << "d_avg: " << d_avg << "\n"
        << "e_avg: " << e_avg << "\n"
        << "f_avg: " << f_avg << "\n";
}