diff options
Diffstat (limited to 'libcxx/test/numerics/rand/rand.dis/rand.dist.bern')
8 files changed, 119 insertions, 119 deletions
diff --git a/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.bernoulli/eval.pass.cpp b/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.bernoulli/eval.pass.cpp index 6764efaa56f..d1927b2d751 100644 --- a/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.bernoulli/eval.pass.cpp +++ b/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.bernoulli/eval.pass.cpp @@ -59,10 +59,10 @@ int main() double x_var = d.p()*(1-d.p()); double x_skew = (1 - 2 * d.p())/std::sqrt(x_var); double x_kurtosis = (6 * sqr(d.p()) - 6 * d.p() + 1)/x_var; - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.01); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.02); } { typedef std::bernoulli_distribution D; @@ -95,9 +95,9 @@ int main() double x_var = d.p()*(1-d.p()); double x_skew = (1 - 2 * d.p())/std::sqrt(x_var); double x_kurtosis = (6 * sqr(d.p()) - 6 * d.p() + 1)/x_var; - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.01); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.02); } } diff --git a/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.bernoulli/eval_param.pass.cpp b/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.bernoulli/eval_param.pass.cpp index bb4be855243..5754a8dc6d6 100644 --- a/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.bernoulli/eval_param.pass.cpp +++ b/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.bernoulli/eval_param.pass.cpp @@ -61,10 +61,10 @@ int main() double x_var = p.p()*(1-p.p()); double x_skew = (1 - 2 * p.p())/std::sqrt(x_var); double x_kurtosis = (6 * sqr(p.p()) - 6 * p.p() + 1)/x_var; - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.01); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.02); } { typedef std::bernoulli_distribution D; @@ -99,9 +99,9 @@ int main() double x_var = p.p()*(1-p.p()); double x_skew = (1 - 2 * p.p())/std::sqrt(x_var); double x_kurtosis = (6 * sqr(p.p()) - 6 * p.p() + 1)/x_var; - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.01); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.02); } } diff --git a/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.bin/eval.pass.cpp b/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.bin/eval.pass.cpp index 71bbf198b5e..a919c9dedae 100644 --- a/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.bin/eval.pass.cpp +++ b/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.bin/eval.pass.cpp @@ -31,10 +31,10 @@ int main() { { typedef std::binomial_distribution<> D; - typedef std::minstd_rand G; + typedef std::mt19937_64 G; G g; D d(5, .75); - const int N = 100000; + const int N = 1000000; std::vector<D::result_type> u; for (int i = 0; i < N; ++i) { @@ -64,10 +64,10 @@ int main() double x_var = x_mean*(1-d.p()); double x_skew = (1-2*d.p()) / std::sqrt(x_var); double x_kurtosis = (1-6*d.p()*(1-d.p())) / x_var; - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.01); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.04); } { typedef std::binomial_distribution<> D; @@ -104,10 +104,10 @@ int main() double x_var = x_mean*(1-d.p()); double x_skew = (1-2*d.p()) / std::sqrt(x_var); double x_kurtosis = (1-6*d.p()*(1-d.p())) / x_var; - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.01); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.01); } { typedef std::binomial_distribution<> D; @@ -144,10 +144,10 @@ int main() double x_var = x_mean*(1-d.p()); double x_skew = (1-2*d.p()) / std::sqrt(x_var); double x_kurtosis = (1-6*d.p()*(1-d.p())) / x_var; - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.03); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.01); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.03); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.3); } { typedef std::binomial_distribution<> D; @@ -260,10 +260,10 @@ int main() double x_var = x_mean*(1-d.p()); double x_skew = (1-2*d.p()) / std::sqrt(x_var); double x_kurtosis = (1-6*d.p()*(1-d.p())) / x_var; - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); assert(std::abs(skew - x_skew) < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.01); + assert(std::abs(kurtosis - x_kurtosis) < 0.01); } { typedef std::binomial_distribution<> D; @@ -300,10 +300,10 @@ int main() double x_var = x_mean*(1-d.p()); double x_skew = (1-2*d.p()) / std::sqrt(x_var); double x_kurtosis = (1-6*d.p()*(1-d.p())) / x_var; - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); assert(std::abs(skew - x_skew) < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.01); } { typedef std::binomial_distribution<> D; diff --git a/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.bin/eval_param.pass.cpp b/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.bin/eval_param.pass.cpp index 2632ab5f3ae..264a4a2f11e 100644 --- a/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.bin/eval_param.pass.cpp +++ b/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.bin/eval_param.pass.cpp @@ -32,11 +32,11 @@ int main() { typedef std::binomial_distribution<> D; typedef D::param_type P; - typedef std::minstd_rand G; + typedef std::mt19937_64 G; G g; D d(16, .75); P p(5, .75); - const int N = 100000; + const int N = 1000000; std::vector<D::result_type> u; for (int i = 0; i < N; ++i) { @@ -66,10 +66,10 @@ int main() double x_var = x_mean*(1-p.p()); double x_skew = (1-2*p.p()) / std::sqrt(x_var); double x_kurtosis = (1-6*p.p()*(1-p.p())) / x_var; - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.01); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.04); } { typedef std::binomial_distribution<> D; @@ -108,10 +108,10 @@ int main() double x_var = x_mean*(1-p.p()); double x_skew = (1-2*p.p()) / std::sqrt(x_var); double x_kurtosis = (1-6*p.p()*(1-p.p())) / x_var; - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.01); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.01); } { typedef std::binomial_distribution<> D; @@ -120,7 +120,7 @@ int main() G g; D d(16, .75); P p(40, .25); - const int N = 100000; + const int N = 1000000; std::vector<D::result_type> u; for (int i = 0; i < N; ++i) { @@ -150,9 +150,9 @@ int main() double x_var = x_mean*(1-p.p()); double x_skew = (1-2*p.p()) / std::sqrt(x_var); double x_kurtosis = (1-6*p.p()*(1-p.p())) / x_var; - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.03); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.01); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.04); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.3); } } diff --git a/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.geo/eval.pass.cpp b/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.geo/eval.pass.cpp index 48a6c57aac8..452bb8224c3 100644 --- a/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.geo/eval.pass.cpp +++ b/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.geo/eval.pass.cpp @@ -64,10 +64,10 @@ int main() double x_var = x_mean / d.p(); double x_skew = (2 - d.p()) / std::sqrt((1 - d.p())); double x_kurtosis = 6 + sqr(d.p()) / (1 - d.p()); - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.01); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.01); } { typedef std::geometric_distribution<> D; @@ -104,10 +104,10 @@ int main() double x_var = x_mean / d.p(); double x_skew = (2 - d.p()) / std::sqrt((1 - d.p())); double x_kurtosis = 6 + sqr(d.p()) / (1 - d.p()); - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.02); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.03); } { typedef std::geometric_distribution<> D; @@ -144,10 +144,10 @@ int main() double x_var = x_mean / d.p(); double x_skew = (2 - d.p()) / std::sqrt((1 - d.p())); double x_kurtosis = 6 + sqr(d.p()) / (1 - d.p()); - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.02); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.02); } { typedef std::geometric_distribution<> D; @@ -184,10 +184,10 @@ int main() double x_var = x_mean / d.p(); double x_skew = (2 - d.p()) / std::sqrt((1 - d.p())); double x_kurtosis = 6 + sqr(d.p()) / (1 - d.p()); - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.01); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.02); } { typedef std::geometric_distribution<> D; @@ -224,10 +224,10 @@ int main() double x_var = x_mean / d.p(); double x_skew = (2 - d.p()) / std::sqrt((1 - d.p())); double x_kurtosis = 6 + sqr(d.p()) / (1 - d.p()); - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.02); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.02); } { typedef std::geometric_distribution<> D; @@ -264,9 +264,9 @@ int main() double x_var = x_mean / d.p(); double x_skew = (2 - d.p()) / std::sqrt((1 - d.p())); double x_kurtosis = 6 + sqr(d.p()) / (1 - d.p()); - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.01); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.02); } } diff --git a/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.geo/eval_param.pass.cpp b/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.geo/eval_param.pass.cpp index bea0159faa7..184124df8ce 100644 --- a/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.geo/eval_param.pass.cpp +++ b/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.geo/eval_param.pass.cpp @@ -66,10 +66,10 @@ int main() double x_var = x_mean / p.p(); double x_skew = (2 - p.p()) / std::sqrt((1 - p.p())); double x_kurtosis = 6 + sqr(p.p()) / (1 - p.p()); - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.01); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.01); } { typedef std::geometric_distribution<> D; @@ -108,10 +108,10 @@ int main() double x_var = x_mean / p.p(); double x_skew = (2 - p.p()) / std::sqrt((1 - p.p())); double x_kurtosis = 6 + sqr(p.p()) / (1 - p.p()); - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.03); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.03); } { typedef std::geometric_distribution<> D; @@ -150,9 +150,9 @@ int main() double x_var = x_mean / p.p(); double x_skew = (2 - p.p()) / std::sqrt((1 - p.p())); double x_kurtosis = 6 + sqr(p.p()) / (1 - p.p()); - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.01); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.02); } } diff --git a/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.negbin/eval.pass.cpp b/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.negbin/eval.pass.cpp index 96b4c180261..6bf3d3788ef 100644 --- a/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.negbin/eval.pass.cpp +++ b/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.negbin/eval.pass.cpp @@ -64,10 +64,10 @@ int main() double x_var = x_mean / d.p(); double x_skew = (2 - d.p()) / std::sqrt(d.k() * (1 - d.p())); double x_kurtosis = 6. / d.k() + sqr(d.p()) / (d.k() * (1 - d.p())); - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.01); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.02); } { typedef std::negative_binomial_distribution<> D; @@ -104,10 +104,10 @@ int main() double x_var = x_mean / d.p(); double x_skew = (2 - d.p()) / std::sqrt(d.k() * (1 - d.p())); double x_kurtosis = 6. / d.k() + sqr(d.p()) / (d.k() * (1 - d.p())); - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.01); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.01); } { typedef std::negative_binomial_distribution<> D; @@ -144,10 +144,10 @@ int main() double x_var = x_mean / d.p(); double x_skew = (2 - d.p()) / std::sqrt(d.k() * (1 - d.p())); double x_kurtosis = 6. / d.k() + sqr(d.p()) / (d.k() * (1 - d.p())); - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.03); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.03); } { typedef std::negative_binomial_distribution<> D; @@ -222,10 +222,10 @@ int main() double x_var = x_mean / d.p(); double x_skew = (2 - d.p()) / std::sqrt(d.k() * (1 - d.p())); double x_kurtosis = 6. / d.k() + sqr(d.p()) / (d.k() * (1 - d.p())); - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.04); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.05); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.04); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.05); } { typedef std::negative_binomial_distribution<> D; @@ -262,9 +262,9 @@ int main() double x_var = x_mean / d.p(); double x_skew = (2 - d.p()) / std::sqrt(d.k() * (1 - d.p())); double x_kurtosis = 6. / d.k() + sqr(d.p()) / (d.k() * (1 - d.p())); - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.02); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.03); } } diff --git a/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.negbin/eval_param.pass.cpp b/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.negbin/eval_param.pass.cpp index b8906721e61..eec909ec7b8 100644 --- a/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.negbin/eval_param.pass.cpp +++ b/libcxx/test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.negbin/eval_param.pass.cpp @@ -66,10 +66,10 @@ int main() double x_var = x_mean / p.p(); double x_skew = (2 - p.p()) / std::sqrt(p.k() * (1 - p.p())); double x_kurtosis = 6. / p.k() + sqr(p.p()) / (p.k() * (1 - p.p())); - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.01); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.01); } { typedef std::negative_binomial_distribution<> D; @@ -108,10 +108,10 @@ int main() double x_var = x_mean / p.p(); double x_skew = (2 - p.p()) / std::sqrt(p.k() * (1 - p.p())); double x_kurtosis = 6. / p.k() + sqr(p.p()) / (p.k() * (1 - p.p())); - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.01); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.01); } { typedef std::negative_binomial_distribution<> D; @@ -150,9 +150,9 @@ int main() double x_var = x_mean / p.p(); double x_skew = (2 - p.p()) / std::sqrt(p.k() * (1 - p.p())); double x_kurtosis = 6. / p.k() + sqr(p.p()) / (p.k() * (1 - p.p())); - assert(std::abs(mean - x_mean) / x_mean < 0.01); - assert(std::abs(var - x_var) / x_var < 0.01); - assert(std::abs(skew - x_skew) / x_skew < 0.01); - assert(std::abs(kurtosis - x_kurtosis) / x_kurtosis < 0.03); + assert(std::abs((mean - x_mean) / x_mean) < 0.01); + assert(std::abs((var - x_var) / x_var) < 0.01); + assert(std::abs((skew - x_skew) / x_skew) < 0.01); + assert(std::abs((kurtosis - x_kurtosis) / x_kurtosis) < 0.03); } } |
