forked from mpusz/mp-units
refactored random distribution number wrapper and test
This commit is contained in:
Ramzi Sabra
committed by
Mateusz Pusz
Mateusz Pusz
parent
1ce7949cf5
commit
2f67372e17
@@ -1,3 +1,25 @@
|
||||
// The MIT License (MIT)
|
||||
//
|
||||
// Copyright (c) 2018 Mateusz Pusz
|
||||
//
|
||||
// Permission is hereby granted, free of charge, to any person obtaining a copy
|
||||
// of this software and associated documentation files (the "Software"), to deal
|
||||
// in the Software without restriction, including without limitation the rights
|
||||
// to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
|
||||
// copies of the Software, and to permit persons to whom the Software is
|
||||
// furnished to do so, subject to the following conditions:
|
||||
//
|
||||
// The above copyright notice and this permission notice shall be included in all
|
||||
// copies or substantial portions of the Software.
|
||||
//
|
||||
// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
|
||||
// IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
|
||||
// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
|
||||
// AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
|
||||
// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
|
||||
// OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
|
||||
// SOFTWARE.
|
||||
|
||||
#include <numeric>
|
||||
#include <units/random.h>
|
||||
#include <units/physical/si/length.h>
|
||||
@@ -8,120 +30,120 @@ using namespace units::physical::si;
|
||||
|
||||
TEST_CASE("uniform_int_distribution")
|
||||
{
|
||||
using Rep = std::int64_t;
|
||||
using Qty = length<metre, Rep>;
|
||||
using rep = std::int64_t;
|
||||
using q = length<metre, rep>;
|
||||
|
||||
SECTION("default")
|
||||
{
|
||||
auto dist = units::uniform_int_distribution<Qty>();
|
||||
auto dist = units::uniform_int_distribution<q>();
|
||||
|
||||
CHECK(dist.a() == Qty::zero());
|
||||
CHECK(dist.b() == Qty::max());
|
||||
CHECK(dist.a() == q::zero());
|
||||
CHECK(dist.b() == q::max());
|
||||
}
|
||||
|
||||
SECTION ("parametrized") {
|
||||
constexpr Rep a = 5;
|
||||
constexpr Rep b = 2;
|
||||
constexpr rep a = 5;
|
||||
constexpr rep b = 2;
|
||||
|
||||
auto stl_dist = std::uniform_int_distribution(a, b);
|
||||
auto units_dist = units::uniform_int_distribution(Qty(a), Qty(b));
|
||||
auto units_dist = units::uniform_int_distribution(q(a), q(b));
|
||||
|
||||
CHECK(units_dist.a() == Qty(stl_dist.a()));
|
||||
CHECK(units_dist.b() == Qty(stl_dist.b()));
|
||||
CHECK(units_dist.min() == Qty(stl_dist.min()));
|
||||
CHECK(units_dist.max() == Qty(stl_dist.max()));
|
||||
CHECK(units_dist.a() == q(stl_dist.a()));
|
||||
CHECK(units_dist.b() == q(stl_dist.b()));
|
||||
CHECK(units_dist.min() == q(stl_dist.min()));
|
||||
CHECK(units_dist.max() == q(stl_dist.max()));
|
||||
}
|
||||
}
|
||||
|
||||
TEST_CASE("uniform_real_distribution")
|
||||
{
|
||||
using Rep = long double;
|
||||
using Qty = length<metre, Rep>;
|
||||
using rep = long double;
|
||||
using q = length<metre, rep>;
|
||||
|
||||
SECTION("default")
|
||||
{
|
||||
auto dist = units::uniform_real_distribution<Qty>();
|
||||
auto dist = units::uniform_real_distribution<q>();
|
||||
|
||||
CHECK(dist.a() == Qty::zero());
|
||||
CHECK(dist.b() == Qty::one());
|
||||
CHECK(dist.a() == q::zero());
|
||||
CHECK(dist.b() == q::one());
|
||||
}
|
||||
|
||||
SECTION ("parametrized") {
|
||||
constexpr Rep a = 5.0;
|
||||
constexpr Rep b = 2.0;
|
||||
constexpr rep a = 5.0;
|
||||
constexpr rep b = 2.0;
|
||||
|
||||
auto stl_dist = std::uniform_real_distribution(a, b);
|
||||
auto units_dist = units::uniform_real_distribution(Qty(a), Qty(b));
|
||||
auto units_dist = units::uniform_real_distribution(q(a), q(b));
|
||||
|
||||
CHECK(units_dist.a() == Qty(stl_dist.a()));
|
||||
CHECK(units_dist.b() == Qty(stl_dist.b()));
|
||||
CHECK(units_dist.min() == Qty(stl_dist.min()));
|
||||
CHECK(units_dist.max() == Qty(stl_dist.max()));
|
||||
CHECK(units_dist.a() == q(stl_dist.a()));
|
||||
CHECK(units_dist.b() == q(stl_dist.b()));
|
||||
CHECK(units_dist.min() == q(stl_dist.min()));
|
||||
CHECK(units_dist.max() == q(stl_dist.max()));
|
||||
}
|
||||
}
|
||||
|
||||
TEST_CASE("binomial_distribution")
|
||||
{
|
||||
using Rep = std::int64_t;
|
||||
using Qty = length<metre, Rep>;
|
||||
using rep = std::int64_t;
|
||||
using q = length<metre, rep>;
|
||||
|
||||
SECTION("default")
|
||||
{
|
||||
auto dist = units::binomial_distribution<Qty>();
|
||||
auto dist = units::binomial_distribution<q>();
|
||||
|
||||
CHECK(dist.p() == 0.5);
|
||||
CHECK(dist.t() == Qty::one());
|
||||
CHECK(dist.t() == q::one());
|
||||
}
|
||||
|
||||
SECTION ("parametrized") {
|
||||
constexpr Rep t = 5;
|
||||
constexpr rep t = 5;
|
||||
constexpr double p = 0.25;
|
||||
|
||||
auto stl_dist = std::binomial_distribution(t, p);
|
||||
auto units_dist = units::binomial_distribution(Qty(t), p);
|
||||
auto units_dist = units::binomial_distribution(q(t), p);
|
||||
|
||||
CHECK(units_dist.p() == stl_dist.p());
|
||||
CHECK(units_dist.t() == Qty(stl_dist.t()));
|
||||
CHECK(units_dist.min() == Qty(stl_dist.min()));
|
||||
CHECK(units_dist.max() == Qty(stl_dist.max()));
|
||||
CHECK(units_dist.t() == q(stl_dist.t()));
|
||||
CHECK(units_dist.min() == q(stl_dist.min()));
|
||||
CHECK(units_dist.max() == q(stl_dist.max()));
|
||||
}
|
||||
}
|
||||
|
||||
TEST_CASE("negative_binomial_distribution")
|
||||
{
|
||||
using Rep = std::int64_t;
|
||||
using Qty = length<metre, Rep>;
|
||||
using rep = std::int64_t;
|
||||
using q = length<metre, rep>;
|
||||
|
||||
SECTION("default")
|
||||
{
|
||||
auto dist = units::negative_binomial_distribution<Qty>();
|
||||
auto dist = units::negative_binomial_distribution<q>();
|
||||
|
||||
CHECK(dist.p() == 0.5);
|
||||
CHECK(dist.k() == Qty::one());
|
||||
CHECK(dist.k() == q::one());
|
||||
}
|
||||
|
||||
SECTION ("parametrized") {
|
||||
constexpr Rep k = 5;
|
||||
constexpr rep k = 5;
|
||||
constexpr double p = 0.25;
|
||||
|
||||
auto stl_dist = std::negative_binomial_distribution(k, p);
|
||||
auto units_dist = units::negative_binomial_distribution(Qty(k), p);
|
||||
auto units_dist = units::negative_binomial_distribution(q(k), p);
|
||||
|
||||
CHECK(units_dist.p() == stl_dist.p());
|
||||
CHECK(units_dist.k() == Qty(stl_dist.k()));
|
||||
CHECK(units_dist.min() == Qty(stl_dist.min()));
|
||||
CHECK(units_dist.max() == Qty(stl_dist.max()));
|
||||
CHECK(units_dist.k() == q(stl_dist.k()));
|
||||
CHECK(units_dist.min() == q(stl_dist.min()));
|
||||
CHECK(units_dist.max() == q(stl_dist.max()));
|
||||
}
|
||||
}
|
||||
|
||||
TEST_CASE("geometric_distribution")
|
||||
{
|
||||
using Rep = std::int64_t;
|
||||
using Qty = length<metre, Rep>;
|
||||
using rep = std::int64_t;
|
||||
using q = length<metre, rep>;
|
||||
|
||||
SECTION("default")
|
||||
{
|
||||
auto dist = units::geometric_distribution<Qty>();
|
||||
auto dist = units::geometric_distribution<q>();
|
||||
|
||||
CHECK(dist.p() == 0.5);
|
||||
}
|
||||
@@ -129,23 +151,23 @@ TEST_CASE("geometric_distribution")
|
||||
SECTION ("parametrized") {
|
||||
constexpr double p = 0.25;
|
||||
|
||||
auto stl_dist = std::geometric_distribution<Rep>(p);
|
||||
auto units_dist = units::geometric_distribution<Qty>(p);
|
||||
auto stl_dist = std::geometric_distribution<rep>(p);
|
||||
auto units_dist = units::geometric_distribution<q>(p);
|
||||
|
||||
CHECK(units_dist.p() == stl_dist.p());
|
||||
CHECK(units_dist.min() == Qty(stl_dist.min()));
|
||||
CHECK(units_dist.max() == Qty(stl_dist.max()));
|
||||
CHECK(units_dist.min() == q(stl_dist.min()));
|
||||
CHECK(units_dist.max() == q(stl_dist.max()));
|
||||
}
|
||||
}
|
||||
|
||||
TEST_CASE("poisson_distribution")
|
||||
{
|
||||
using Rep = std::int64_t;
|
||||
using Qty = length<metre, Rep>;
|
||||
using rep = std::int64_t;
|
||||
using q = length<metre, rep>;
|
||||
|
||||
SECTION("default")
|
||||
{
|
||||
auto dist = units::poisson_distribution<Qty>();
|
||||
auto dist = units::poisson_distribution<q>();
|
||||
|
||||
CHECK(dist.mean() == 1.0);
|
||||
}
|
||||
@@ -153,23 +175,23 @@ TEST_CASE("poisson_distribution")
|
||||
SECTION ("parametrized") {
|
||||
constexpr double mean = 5.0;
|
||||
|
||||
auto stl_dist = std::poisson_distribution<Rep>(mean);
|
||||
auto units_dist = units::poisson_distribution<Qty>(mean);
|
||||
auto stl_dist = std::poisson_distribution<rep>(mean);
|
||||
auto units_dist = units::poisson_distribution<q>(mean);
|
||||
|
||||
CHECK(units_dist.mean() == stl_dist.mean());
|
||||
CHECK(units_dist.min() == Qty(stl_dist.min()));
|
||||
CHECK(units_dist.max() == Qty(stl_dist.max()));
|
||||
CHECK(units_dist.min() == q(stl_dist.min()));
|
||||
CHECK(units_dist.max() == q(stl_dist.max()));
|
||||
}
|
||||
}
|
||||
|
||||
TEST_CASE("exponential_distribution")
|
||||
{
|
||||
using Rep = long double;
|
||||
using Qty = length<metre, Rep>;
|
||||
using rep = long double;
|
||||
using q = length<metre, rep>;
|
||||
|
||||
SECTION("default")
|
||||
{
|
||||
auto dist = units::exponential_distribution<Qty>();
|
||||
auto dist = units::exponential_distribution<q>();
|
||||
|
||||
CHECK(dist.lambda() == 1.0);
|
||||
}
|
||||
@@ -177,23 +199,23 @@ TEST_CASE("exponential_distribution")
|
||||
SECTION ("parametrized") {
|
||||
constexpr double lambda = 2.0;
|
||||
|
||||
auto stl_dist = std::exponential_distribution<Rep>(lambda);
|
||||
auto units_dist = units::exponential_distribution<Qty>(lambda);
|
||||
auto stl_dist = std::exponential_distribution<rep>(lambda);
|
||||
auto units_dist = units::exponential_distribution<q>(lambda);
|
||||
|
||||
CHECK(units_dist.lambda() == stl_dist.lambda());
|
||||
CHECK(units_dist.min() == Qty(stl_dist.min()));
|
||||
CHECK(units_dist.max() == Qty(stl_dist.max()));
|
||||
CHECK(units_dist.min() == q(stl_dist.min()));
|
||||
CHECK(units_dist.max() == q(stl_dist.max()));
|
||||
}
|
||||
}
|
||||
|
||||
TEST_CASE("gamma_distribution")
|
||||
{
|
||||
using Rep = long double;
|
||||
using Qty = length<metre, Rep>;
|
||||
using rep = long double;
|
||||
using q = length<metre, rep>;
|
||||
|
||||
SECTION("default")
|
||||
{
|
||||
auto dist = units::gamma_distribution<Qty>();
|
||||
auto dist = units::gamma_distribution<q>();
|
||||
|
||||
CHECK(dist.alpha() == 1.0);
|
||||
CHECK(dist.beta() == 1.0);
|
||||
@@ -203,249 +225,249 @@ TEST_CASE("gamma_distribution")
|
||||
constexpr double alpha = 5.0;
|
||||
constexpr double beta = 2.0;
|
||||
|
||||
auto stl_dist = std::gamma_distribution<Rep>(alpha, beta);
|
||||
auto units_dist = units::gamma_distribution<Qty>(alpha, beta);
|
||||
auto stl_dist = std::gamma_distribution<rep>(alpha, beta);
|
||||
auto units_dist = units::gamma_distribution<q>(alpha, beta);
|
||||
|
||||
CHECK(units_dist.alpha() == stl_dist.alpha());
|
||||
CHECK(units_dist.beta() == stl_dist.beta());
|
||||
CHECK(units_dist.min() == Qty(stl_dist.min()));
|
||||
CHECK(units_dist.max() == Qty(stl_dist.max()));
|
||||
CHECK(units_dist.min() == q(stl_dist.min()));
|
||||
CHECK(units_dist.max() == q(stl_dist.max()));
|
||||
}
|
||||
}
|
||||
|
||||
TEST_CASE("weibull_distribution")
|
||||
{
|
||||
using Rep = long double;
|
||||
using Qty = length<metre, Rep>;
|
||||
using rep = long double;
|
||||
using q = length<metre, rep>;
|
||||
|
||||
SECTION("default")
|
||||
{
|
||||
auto dist = units::weibull_distribution<Qty>();
|
||||
auto dist = units::weibull_distribution<q>();
|
||||
|
||||
CHECK(dist.a() == 1.0);
|
||||
CHECK(dist.b() == 1.0);
|
||||
}
|
||||
|
||||
SECTION ("parametrized") {
|
||||
constexpr Rep a = 5.0;
|
||||
constexpr Rep b = 2.0;
|
||||
constexpr rep a = 5.0;
|
||||
constexpr rep b = 2.0;
|
||||
|
||||
auto stl_dist = std::weibull_distribution(a, b);
|
||||
auto units_dist = units::weibull_distribution<Qty>(a, b);
|
||||
auto units_dist = units::weibull_distribution<q>(a, b);
|
||||
|
||||
CHECK(units_dist.a() == stl_dist.a());
|
||||
CHECK(units_dist.b() == stl_dist.b());
|
||||
CHECK(units_dist.min() == Qty(stl_dist.min()));
|
||||
CHECK(units_dist.max() == Qty(stl_dist.max()));
|
||||
CHECK(units_dist.min() == q(stl_dist.min()));
|
||||
CHECK(units_dist.max() == q(stl_dist.max()));
|
||||
}
|
||||
}
|
||||
|
||||
TEST_CASE("extreme_value_distribution")
|
||||
{
|
||||
using Rep = long double;
|
||||
using Qty = length<metre, Rep>;
|
||||
using rep = long double;
|
||||
using q = length<metre, rep>;
|
||||
|
||||
SECTION("default")
|
||||
{
|
||||
auto dist = units::extreme_value_distribution<Qty>();
|
||||
auto dist = units::extreme_value_distribution<q>();
|
||||
|
||||
CHECK(dist.a() == Qty::zero());
|
||||
CHECK(dist.a() == q::zero());
|
||||
CHECK(dist.b() == 1.0);
|
||||
}
|
||||
|
||||
SECTION ("parametrized") {
|
||||
constexpr Rep a = 5.0;
|
||||
constexpr Rep b = 2.0;
|
||||
constexpr rep a = 5.0;
|
||||
constexpr rep b = 2.0;
|
||||
|
||||
auto stl_dist = std::extreme_value_distribution(a, b);
|
||||
auto units_dist = units::extreme_value_distribution<Qty>(Qty(a), b);
|
||||
auto units_dist = units::extreme_value_distribution<q>(q(a), b);
|
||||
|
||||
CHECK(units_dist.a() == Qty(stl_dist.a()));
|
||||
CHECK(units_dist.a() == q(stl_dist.a()));
|
||||
CHECK(units_dist.b() == stl_dist.b());
|
||||
CHECK(units_dist.min() == Qty(stl_dist.min()));
|
||||
CHECK(units_dist.max() == Qty(stl_dist.max()));
|
||||
CHECK(units_dist.min() == q(stl_dist.min()));
|
||||
CHECK(units_dist.max() == q(stl_dist.max()));
|
||||
}
|
||||
}
|
||||
|
||||
TEST_CASE("normal_distribution")
|
||||
{
|
||||
using Rep = long double;
|
||||
using Qty = length<metre, Rep>;
|
||||
using rep = long double;
|
||||
using q = length<metre, rep>;
|
||||
|
||||
SECTION("default")
|
||||
{
|
||||
auto dist = units::normal_distribution<Qty>();
|
||||
auto dist = units::normal_distribution<q>();
|
||||
|
||||
CHECK(dist.mean() == Qty::zero());
|
||||
CHECK(dist.stddev() == Qty::one());
|
||||
CHECK(dist.mean() == q::zero());
|
||||
CHECK(dist.stddev() == q::one());
|
||||
}
|
||||
|
||||
SECTION("parametrized")
|
||||
{
|
||||
constexpr Rep mean = 5.0;
|
||||
constexpr Rep stddev = 2.0;
|
||||
constexpr rep mean = 5.0;
|
||||
constexpr rep stddev = 2.0;
|
||||
|
||||
auto stl_dist = std::normal_distribution(mean, stddev);
|
||||
auto units_dist = units::normal_distribution(Qty(mean), Qty(stddev));
|
||||
auto units_dist = units::normal_distribution(q(mean), q(stddev));
|
||||
|
||||
CHECK(units_dist.mean() == Qty(stl_dist.mean()));
|
||||
CHECK(units_dist.stddev() == Qty(stl_dist.stddev()));
|
||||
CHECK(units_dist.min() == Qty(stl_dist.min()));
|
||||
CHECK(units_dist.max() == Qty(stl_dist.max()));
|
||||
CHECK(units_dist.mean() == q(stl_dist.mean()));
|
||||
CHECK(units_dist.stddev() == q(stl_dist.stddev()));
|
||||
CHECK(units_dist.min() == q(stl_dist.min()));
|
||||
CHECK(units_dist.max() == q(stl_dist.max()));
|
||||
}
|
||||
}
|
||||
|
||||
TEST_CASE("lognormal_distribution")
|
||||
{
|
||||
using Rep = long double;
|
||||
using Qty = length<metre, Rep>;
|
||||
using rep = long double;
|
||||
using q = length<metre, rep>;
|
||||
|
||||
SECTION("default")
|
||||
{
|
||||
auto dist = units::lognormal_distribution<Qty>();
|
||||
auto dist = units::lognormal_distribution<q>();
|
||||
|
||||
CHECK(dist.m() == Qty::zero());
|
||||
CHECK(dist.s() == Qty::one());
|
||||
CHECK(dist.m() == q::zero());
|
||||
CHECK(dist.s() == q::one());
|
||||
}
|
||||
|
||||
SECTION("parametrized")
|
||||
{
|
||||
constexpr Rep m = 5.0;
|
||||
constexpr Rep s = 2.0;
|
||||
constexpr rep m = 5.0;
|
||||
constexpr rep s = 2.0;
|
||||
|
||||
auto stl_dist = std::lognormal_distribution(m, s);
|
||||
auto units_dist = units::lognormal_distribution(Qty(m), Qty(s));
|
||||
auto units_dist = units::lognormal_distribution(q(m), q(s));
|
||||
|
||||
CHECK(units_dist.m() == Qty(stl_dist.m()));
|
||||
CHECK(units_dist.s() == Qty(stl_dist.s()));
|
||||
CHECK(units_dist.min() == Qty(stl_dist.min()));
|
||||
CHECK(units_dist.max() == Qty(stl_dist.max()));
|
||||
CHECK(units_dist.m() == q(stl_dist.m()));
|
||||
CHECK(units_dist.s() == q(stl_dist.s()));
|
||||
CHECK(units_dist.min() == q(stl_dist.min()));
|
||||
CHECK(units_dist.max() == q(stl_dist.max()));
|
||||
}
|
||||
}
|
||||
|
||||
TEST_CASE("chi_squared_distribution")
|
||||
{
|
||||
using Rep = long double;
|
||||
using Qty = length<metre, Rep>;
|
||||
using rep = long double;
|
||||
using q = length<metre, rep>;
|
||||
|
||||
SECTION("default")
|
||||
{
|
||||
auto dist = units::chi_squared_distribution<Qty>();
|
||||
auto dist = units::chi_squared_distribution<q>();
|
||||
|
||||
CHECK(dist.n() == 1.0);
|
||||
}
|
||||
|
||||
SECTION("parametrized")
|
||||
{
|
||||
constexpr Rep n = 5.0;
|
||||
constexpr rep n = 5.0;
|
||||
|
||||
auto stl_dist = std::chi_squared_distribution(n);
|
||||
auto units_dist = units::chi_squared_distribution<Qty>(n);
|
||||
auto units_dist = units::chi_squared_distribution<q>(n);
|
||||
|
||||
CHECK(units_dist.n() == stl_dist.n());
|
||||
CHECK(units_dist.min() == Qty(stl_dist.min()));
|
||||
CHECK(units_dist.max() == Qty(stl_dist.max()));
|
||||
CHECK(units_dist.min() == q(stl_dist.min()));
|
||||
CHECK(units_dist.max() == q(stl_dist.max()));
|
||||
}
|
||||
}
|
||||
|
||||
TEST_CASE("cauchy_distribution")
|
||||
{
|
||||
using Rep = long double;
|
||||
using Qty = length<metre, Rep>;
|
||||
using rep = long double;
|
||||
using q = length<metre, rep>;
|
||||
|
||||
SECTION("default")
|
||||
{
|
||||
auto dist = units::cauchy_distribution<Qty>();
|
||||
auto dist = units::cauchy_distribution<q>();
|
||||
|
||||
CHECK(dist.a() == Qty::zero());
|
||||
CHECK(dist.b() == Qty::one());
|
||||
CHECK(dist.a() == q::zero());
|
||||
CHECK(dist.b() == q::one());
|
||||
}
|
||||
|
||||
SECTION ("parametrized") {
|
||||
constexpr Rep a = 5.0;
|
||||
constexpr Rep b = 2.0;
|
||||
constexpr rep a = 5.0;
|
||||
constexpr rep b = 2.0;
|
||||
|
||||
auto stl_dist = std::cauchy_distribution(a, b);
|
||||
auto units_dist = units::cauchy_distribution(Qty(a), Qty(b));
|
||||
auto units_dist = units::cauchy_distribution(q(a), q(b));
|
||||
|
||||
CHECK(units_dist.a() == Qty(stl_dist.a()));
|
||||
CHECK(units_dist.b() == Qty(stl_dist.b()));
|
||||
CHECK(units_dist.min() == Qty(stl_dist.min()));
|
||||
CHECK(units_dist.max() == Qty(stl_dist.max()));
|
||||
CHECK(units_dist.a() == q(stl_dist.a()));
|
||||
CHECK(units_dist.b() == q(stl_dist.b()));
|
||||
CHECK(units_dist.min() == q(stl_dist.min()));
|
||||
CHECK(units_dist.max() == q(stl_dist.max()));
|
||||
}
|
||||
}
|
||||
|
||||
TEST_CASE("fisher_f_distribution")
|
||||
{
|
||||
using Rep = long double;
|
||||
using Qty = length<metre, Rep>;
|
||||
using rep = long double;
|
||||
using q = length<metre, rep>;
|
||||
|
||||
SECTION("default")
|
||||
{
|
||||
auto dist = units::fisher_f_distribution<Qty>();
|
||||
auto dist = units::fisher_f_distribution<q>();
|
||||
|
||||
CHECK(dist.m() == 1.0);
|
||||
CHECK(dist.n() == 1.0);
|
||||
}
|
||||
|
||||
SECTION ("parametrized") {
|
||||
constexpr Rep m = 5.0;
|
||||
constexpr Rep n = 2.0;
|
||||
constexpr rep m = 5.0;
|
||||
constexpr rep n = 2.0;
|
||||
|
||||
auto stl_dist = std::fisher_f_distribution<Rep>(m, n);
|
||||
auto units_dist = units::fisher_f_distribution<Qty>(m, n);
|
||||
auto stl_dist = std::fisher_f_distribution<rep>(m, n);
|
||||
auto units_dist = units::fisher_f_distribution<q>(m, n);
|
||||
|
||||
CHECK(units_dist.m() == stl_dist.m());
|
||||
CHECK(units_dist.n() == stl_dist.n());
|
||||
CHECK(units_dist.min() == Qty(stl_dist.min()));
|
||||
CHECK(units_dist.max() == Qty(stl_dist.max()));
|
||||
CHECK(units_dist.min() == q(stl_dist.min()));
|
||||
CHECK(units_dist.max() == q(stl_dist.max()));
|
||||
}
|
||||
}
|
||||
|
||||
TEST_CASE("student_t_distribution")
|
||||
{
|
||||
using Rep = long double;
|
||||
using Qty = length<metre, Rep>;
|
||||
using rep = long double;
|
||||
using q = length<metre, rep>;
|
||||
|
||||
SECTION("default")
|
||||
{
|
||||
auto dist = units::student_t_distribution<Qty>();
|
||||
auto dist = units::student_t_distribution<q>();
|
||||
|
||||
CHECK(dist.n() == 1.0);
|
||||
}
|
||||
|
||||
SECTION ("parametrized") {
|
||||
constexpr Rep n = 2.0;
|
||||
constexpr rep n = 2.0;
|
||||
|
||||
auto stl_dist = std::student_t_distribution<Rep>(n);
|
||||
auto units_dist = units::student_t_distribution<Qty>(n);
|
||||
auto stl_dist = std::student_t_distribution<rep>(n);
|
||||
auto units_dist = units::student_t_distribution<q>(n);
|
||||
|
||||
CHECK(units_dist.n() == stl_dist.n());
|
||||
CHECK(units_dist.min() == Qty(stl_dist.min()));
|
||||
CHECK(units_dist.max() == Qty(stl_dist.max()));
|
||||
CHECK(units_dist.min() == q(stl_dist.min()));
|
||||
CHECK(units_dist.max() == q(stl_dist.max()));
|
||||
}
|
||||
}
|
||||
|
||||
TEST_CASE("discrete_distribution")
|
||||
{
|
||||
using Rep = std::int64_t;
|
||||
using Qty = length<metre, Rep>;
|
||||
using rep = std::int64_t;
|
||||
using q = length<metre, rep>;
|
||||
|
||||
SECTION("default")
|
||||
{
|
||||
auto stl_dist = std::discrete_distribution<Rep>();
|
||||
auto units_dist = units::discrete_distribution<Qty>();
|
||||
auto stl_dist = std::discrete_distribution<rep>();
|
||||
auto units_dist = units::discrete_distribution<q>();
|
||||
|
||||
CHECK(units_dist.min() == Qty(stl_dist.min()));
|
||||
CHECK(units_dist.max() == Qty(stl_dist.max()));
|
||||
CHECK(units_dist.min() == q(stl_dist.min()));
|
||||
CHECK(units_dist.max() == q(stl_dist.max()));
|
||||
CHECK(units_dist.probabilities() == stl_dist.probabilities());
|
||||
}
|
||||
|
||||
SECTION ("parametrized_input_it") {
|
||||
constexpr std::array<double, 3> weights = {1.0, 2.0, 3.0};
|
||||
|
||||
auto stl_dist = std::discrete_distribution<Rep>(weights.cbegin(), weights.cend());
|
||||
auto units_dist = units::discrete_distribution<Qty>(weights.cbegin(), weights.cend());
|
||||
auto stl_dist = std::discrete_distribution<rep>(weights.cbegin(), weights.cend());
|
||||
auto units_dist = units::discrete_distribution<q>(weights.cbegin(), weights.cend());
|
||||
|
||||
CHECK(units_dist.probabilities() == stl_dist.probabilities());
|
||||
}
|
||||
@@ -453,8 +475,8 @@ TEST_CASE("discrete_distribution")
|
||||
SECTION ("parametrized_initializer_list") {
|
||||
std::initializer_list<double> weights = {1.0, 2.0, 3.0};
|
||||
|
||||
auto stl_dist = std::discrete_distribution<Rep>(weights);
|
||||
auto units_dist = units::discrete_distribution<Qty>(weights);
|
||||
auto stl_dist = std::discrete_distribution<rep>(weights);
|
||||
auto units_dist = units::discrete_distribution<q>(weights);
|
||||
|
||||
CHECK(units_dist.probabilities() == stl_dist.probabilities());
|
||||
}
|
||||
@@ -463,8 +485,8 @@ TEST_CASE("discrete_distribution")
|
||||
constexpr std::size_t count = 3;
|
||||
constexpr double xmin = 1, xmax = 3;
|
||||
|
||||
auto stl_dist = std::discrete_distribution<Rep>(count, xmin, xmax, [](double val) { return val; });
|
||||
auto units_dist = units::discrete_distribution<Qty>(count, xmin, xmax, [](double val) { return val; });
|
||||
auto stl_dist = std::discrete_distribution<rep>(count, xmin, xmax, [](double val) { return val; });
|
||||
auto units_dist = units::discrete_distribution<q>(count, xmin, xmax, [](double val) { return val; });
|
||||
|
||||
CHECK(units_dist.probabilities() == stl_dist.probabilities());
|
||||
}
|
||||
@@ -472,19 +494,19 @@ TEST_CASE("discrete_distribution")
|
||||
|
||||
TEST_CASE("piecewise_constant_distribution")
|
||||
{
|
||||
using Rep = long double;
|
||||
using Qty = length<metre, Rep>;
|
||||
using rep = long double;
|
||||
using q = length<metre, rep>;
|
||||
|
||||
std::vector<Rep> intervals_rep_vec = {1.0, 2.0, 3.0};
|
||||
std::vector<Qty> intervals_qty_vec = {1.0q_m, 2.0q_m, 3.0q_m};
|
||||
std::vector<rep> intervals_rep_vec = {1.0, 2.0, 3.0};
|
||||
std::vector<q> intervals_qty_vec = {1.0q_m, 2.0q_m, 3.0q_m};
|
||||
|
||||
SECTION("default")
|
||||
{
|
||||
auto stl_dist = std::piecewise_constant_distribution<Rep>();
|
||||
auto units_dist = units::piecewise_constant_distribution<Qty>();
|
||||
auto stl_dist = std::piecewise_constant_distribution<rep>();
|
||||
auto units_dist = units::piecewise_constant_distribution<q>();
|
||||
|
||||
CHECK(units_dist.min() == Qty(stl_dist.min()));
|
||||
CHECK(units_dist.max() == Qty(stl_dist.max()));
|
||||
CHECK(units_dist.min() == q(stl_dist.min()));
|
||||
CHECK(units_dist.max() == q(stl_dist.max()));
|
||||
CHECK(stl_dist.intervals().size() == 2);
|
||||
CHECK(units_dist.intervals().size() == 2);
|
||||
CHECK(stl_dist.densities().size() == 1);
|
||||
@@ -492,12 +514,12 @@ TEST_CASE("piecewise_constant_distribution")
|
||||
}
|
||||
|
||||
SECTION ("parametrized_input_it") {
|
||||
constexpr std::array<Rep, 3> intervals_rep = {1.0, 2.0, 3.0};
|
||||
constexpr std::array<Qty, 3> intervals_qty = {1.0q_m, 2.0q_m, 3.0q_m};
|
||||
constexpr std::array<Rep, 3> weights = {1.0, 2.0, 3.0};
|
||||
constexpr std::array<rep, 3> intervals_rep = {1.0, 2.0, 3.0};
|
||||
constexpr std::array<q, 3> intervals_qty = {1.0q_m, 2.0q_m, 3.0q_m};
|
||||
constexpr std::array<rep, 3> weights = {1.0, 2.0, 3.0};
|
||||
|
||||
auto stl_dist = std::piecewise_constant_distribution<Rep>(intervals_rep.cbegin(), intervals_rep.cend(), weights.cbegin());
|
||||
auto units_dist = units::piecewise_constant_distribution<Qty>(intervals_qty.cbegin(), intervals_qty.cend(), weights.cbegin());
|
||||
auto stl_dist = std::piecewise_constant_distribution<rep>(intervals_rep.cbegin(), intervals_rep.cend(), weights.cbegin());
|
||||
auto units_dist = units::piecewise_constant_distribution<q>(intervals_qty.cbegin(), intervals_qty.cend(), weights.cbegin());
|
||||
|
||||
CHECK(stl_dist.intervals() == intervals_rep_vec);
|
||||
CHECK(units_dist.intervals() == intervals_qty_vec);
|
||||
@@ -505,11 +527,11 @@ TEST_CASE("piecewise_constant_distribution")
|
||||
}
|
||||
|
||||
SECTION ("parametrized_initializer_list") {
|
||||
std::initializer_list<Rep> intervals_rep = {1.0, 2.0, 3.0};
|
||||
std::initializer_list<Qty> intervals_qty = {1.0q_m, 2.0q_m, 3.0q_m};
|
||||
std::initializer_list<rep> intervals_rep = {1.0, 2.0, 3.0};
|
||||
std::initializer_list<q> intervals_qty = {1.0q_m, 2.0q_m, 3.0q_m};
|
||||
|
||||
auto stl_dist = std::piecewise_constant_distribution<Rep>(intervals_rep, [](Rep val) { return val; });
|
||||
auto units_dist = units::piecewise_constant_distribution<Qty>(intervals_qty, [](Qty qty) { return qty.count(); });
|
||||
auto stl_dist = std::piecewise_constant_distribution<rep>(intervals_rep, [](rep val) { return val; });
|
||||
auto units_dist = units::piecewise_constant_distribution<q>(intervals_qty, [](q qty) { return qty.count(); });
|
||||
|
||||
CHECK(units_dist.intervals() == intervals_qty_vec);
|
||||
CHECK(units_dist.densities() == stl_dist.densities());
|
||||
@@ -517,11 +539,11 @@ TEST_CASE("piecewise_constant_distribution")
|
||||
|
||||
SECTION ("parametrized_range") {
|
||||
constexpr std::size_t nw = 2;
|
||||
constexpr Rep xmin_rep = 1.0, xmax_rep = 3.0;
|
||||
constexpr Qty xmin_qty = 1.0q_m, xmax_qty = 3.0q_m;
|
||||
constexpr rep xmin_rep = 1.0, xmax_rep = 3.0;
|
||||
constexpr q xmin_qty = 1.0q_m, xmax_qty = 3.0q_m;
|
||||
|
||||
auto stl_dist = std::piecewise_constant_distribution<Rep>(nw, xmin_rep, xmax_rep, [](Rep val) { return val; });
|
||||
auto units_dist = units::piecewise_constant_distribution<Qty>(nw, xmin_qty, xmax_qty, [](Qty qty) { return qty.count(); });
|
||||
auto stl_dist = std::piecewise_constant_distribution<rep>(nw, xmin_rep, xmax_rep, [](rep val) { return val; });
|
||||
auto units_dist = units::piecewise_constant_distribution<q>(nw, xmin_qty, xmax_qty, [](q qty) { return qty.count(); });
|
||||
|
||||
CHECK(units_dist.intervals() == intervals_qty_vec);
|
||||
CHECK(units_dist.densities() == stl_dist.densities());
|
||||
@@ -530,19 +552,19 @@ TEST_CASE("piecewise_constant_distribution")
|
||||
|
||||
TEST_CASE("piecewise_linear_distribution")
|
||||
{
|
||||
using Rep = long double;
|
||||
using Qty = length<metre, Rep>;
|
||||
using rep = long double;
|
||||
using q = length<metre, rep>;
|
||||
|
||||
std::vector<Rep> intervals_rep_vec = {1.0, 2.0, 3.0};
|
||||
std::vector<Qty> intervals_qty_vec = {1.0q_m, 2.0q_m, 3.0q_m};
|
||||
std::vector<rep> intervals_rep_vec = {1.0, 2.0, 3.0};
|
||||
std::vector<q> intervals_qty_vec = {1.0q_m, 2.0q_m, 3.0q_m};
|
||||
|
||||
SECTION("default")
|
||||
{
|
||||
auto stl_dist = std::piecewise_linear_distribution<Rep>();
|
||||
auto units_dist = units::piecewise_linear_distribution<Qty>();
|
||||
auto stl_dist = std::piecewise_linear_distribution<rep>();
|
||||
auto units_dist = units::piecewise_linear_distribution<q>();
|
||||
|
||||
CHECK(units_dist.min() == Qty(stl_dist.min()));
|
||||
CHECK(units_dist.max() == Qty(stl_dist.max()));
|
||||
CHECK(units_dist.min() == q(stl_dist.min()));
|
||||
CHECK(units_dist.max() == q(stl_dist.max()));
|
||||
CHECK(stl_dist.intervals().size() == 2);
|
||||
CHECK(units_dist.intervals().size() == 2);
|
||||
CHECK(stl_dist.densities().size() == 2);
|
||||
@@ -550,12 +572,12 @@ TEST_CASE("piecewise_linear_distribution")
|
||||
}
|
||||
|
||||
SECTION ("parametrized_input_it") {
|
||||
constexpr std::array<Rep, 3> intervals_rep = {1.0, 2.0, 3.0};
|
||||
constexpr std::array<Qty, 3> intervals_qty = {1.0q_m, 2.0q_m, 3.0q_m};
|
||||
constexpr std::array<Rep, 3> weights = {1.0, 2.0, 3.0};
|
||||
constexpr std::array<rep, 3> intervals_rep = {1.0, 2.0, 3.0};
|
||||
constexpr std::array<q, 3> intervals_qty = {1.0q_m, 2.0q_m, 3.0q_m};
|
||||
constexpr std::array<rep, 3> weights = {1.0, 2.0, 3.0};
|
||||
|
||||
auto stl_dist = std::piecewise_linear_distribution<Rep>(intervals_rep.cbegin(), intervals_rep.cend(), weights.cbegin());
|
||||
auto units_dist = units::piecewise_linear_distribution<Qty>(intervals_qty.cbegin(), intervals_qty.cend(), weights.cbegin());
|
||||
auto stl_dist = std::piecewise_linear_distribution<rep>(intervals_rep.cbegin(), intervals_rep.cend(), weights.cbegin());
|
||||
auto units_dist = units::piecewise_linear_distribution<q>(intervals_qty.cbegin(), intervals_qty.cend(), weights.cbegin());
|
||||
|
||||
CHECK(stl_dist.intervals() == intervals_rep_vec);
|
||||
CHECK(units_dist.intervals() == intervals_qty_vec);
|
||||
@@ -563,11 +585,11 @@ TEST_CASE("piecewise_linear_distribution")
|
||||
}
|
||||
|
||||
SECTION ("parametrized_initializer_list") {
|
||||
std::initializer_list<Rep> intervals_rep = {1.0, 2.0, 3.0};
|
||||
std::initializer_list<Qty> intervals_qty = {1.0q_m, 2.0q_m, 3.0q_m};
|
||||
std::initializer_list<rep> intervals_rep = {1.0, 2.0, 3.0};
|
||||
std::initializer_list<q> intervals_qty = {1.0q_m, 2.0q_m, 3.0q_m};
|
||||
|
||||
auto stl_dist = std::piecewise_linear_distribution<Rep>(intervals_rep, [](Rep val) { return val; });
|
||||
auto units_dist = units::piecewise_linear_distribution<Qty>(intervals_qty, [](Qty qty) { return qty.count(); });
|
||||
auto stl_dist = std::piecewise_linear_distribution<rep>(intervals_rep, [](rep val) { return val; });
|
||||
auto units_dist = units::piecewise_linear_distribution<q>(intervals_qty, [](q qty) { return qty.count(); });
|
||||
|
||||
CHECK(units_dist.intervals() == intervals_qty_vec);
|
||||
CHECK(units_dist.densities() == stl_dist.densities());
|
||||
@@ -575,11 +597,11 @@ TEST_CASE("piecewise_linear_distribution")
|
||||
|
||||
SECTION ("parametrized_range") {
|
||||
constexpr std::size_t nw = 2;
|
||||
constexpr Rep xmin_rep = 1.0, xmax_rep = 3.0;
|
||||
constexpr Qty xmin_qty = 1.0q_m, xmax_qty = 3.0q_m;
|
||||
constexpr rep xmin_rep = 1.0, xmax_rep = 3.0;
|
||||
constexpr q xmin_qty = 1.0q_m, xmax_qty = 3.0q_m;
|
||||
|
||||
auto stl_dist = std::piecewise_linear_distribution<Rep>(nw, xmin_rep, xmax_rep, [](Rep val) { return val; });
|
||||
auto units_dist = units::piecewise_linear_distribution<Qty>(nw, xmin_qty, xmax_qty, [](Qty qty) { return qty.count(); });
|
||||
auto stl_dist = std::piecewise_linear_distribution<rep>(nw, xmin_rep, xmax_rep, [](rep val) { return val; });
|
||||
auto units_dist = units::piecewise_linear_distribution<q>(nw, xmin_qty, xmax_qty, [](q qty) { return qty.count(); });
|
||||
|
||||
CHECK(units_dist.intervals() == intervals_qty_vec);
|
||||
CHECK(units_dist.densities() == stl_dist.densities());
|
||||
|
||||
Reference in New Issue
Block a user