27 #include <mln/accu/stat/var.hh>
28 #include <mln/fun/stat/mahalanobis.hh>
33 return (1 + c) * float(std::rand()) / RAND_MAX;
41 typedef algebra::vec<3,float> vec3f;
46 for (
int i = 0; i < n; ++i)
54 for (
int i = 0; i < n; ++i)
57 fun::stat::mahalanobis<vec3f> f(a.variance(), a.mean());
58 mln_assertion(f(a.mean()) == 0.f);
61 for (
int i = 0; i < n; ++i)
66 mln_assertion(std::abs(sum / n - 3.f) < 0.00002f);