generic-library/cxx
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

[rand.dist.samp.pconst] plus some bug fixes in the tests of the other distributions

Browse Source 
git-svn-id: https://llvm.org/svn/llvm-project/libcxx/trunk@104224 91177308-0d34-0410-b5e6-96231b3b80d8
This commit is contained in:
Howard Hinnant
2010-05-20 15:11:46 +00:00
parent 551d8e4ddb
commit d6d1171f2c
56 changed files with 2593 additions and 393 deletions
Download Patch File Download Diff File Expand all files Collapse all files

30
test/numerics/rand/rand.dis/rand.dist.bern/rand.dist.bern.bin/eval_param.pass.cpp
View File

@@ -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);
}
}