#include <random>
#include <cmath>
+#define EPSILON (100*std::numeric_limits<double>::epsilon())
+
TEST_CASE("xbt::random: Random Number Generation")
{
SECTION("Using XBT_RNG_xbt")
{
simgrid::xbt::random::set_mersenne_seed(12345);
-
- REQUIRE(fabs(simgrid::xbt::random::exponential(25) - 0.00291934351538427348) <= 1e-14);
+ REQUIRE(simgrid::xbt::random::exponential(25) == Approx(0.00291934351538427348).epsilon(EPSILON));
REQUIRE(simgrid::xbt::random::uniform_int(1, 6) == 4);
- REQUIRE(fabs(simgrid::xbt::random::uniform_real(0, 1) - 0.31637556043369124970) <= 1e-14);
- REQUIRE(fabs(simgrid::xbt::random::normal(0, 2) - 1.62746784745133976635) <= 1e-14);
+ REQUIRE(simgrid::xbt::random::uniform_real(0, 1) == Approx(0.31637556043369124970).epsilon(EPSILON));
+ REQUIRE(simgrid::xbt::random::normal(0, 2) == Approx(1.62746784745133976635).epsilon(EPSILON));
}
SECTION("Using XBT_RNG_std")
std::uniform_real_distribution<> distC(0, 1);
std::normal_distribution<> distD(0, 2);
- REQUIRE(fabs(simgrid::xbt::random::exponential(25) - distA(gen)) <= 1e-14);
+ REQUIRE(simgrid::xbt::random::exponential(25) == Approx(distA(gen)).epsilon(EPSILON));
REQUIRE(simgrid::xbt::random::uniform_int(1, 6) == distB(gen));
- REQUIRE(fabs(simgrid::xbt::random::uniform_real(0, 1) - distC(gen)) <= 1e-14);
- REQUIRE(fabs(simgrid::xbt::random::normal(0, 2) - distD(gen)) <= 1e-14);
+ REQUIRE(simgrid::xbt::random::uniform_real(0, 1) == Approx(distC(gen)).epsilon(EPSILON));
+ REQUIRE(simgrid::xbt::random::normal(0, 2) == Approx(distD(gen)).epsilon(EPSILON));
}
}