00001 // Copyright (C) 2009 EPITA Research and Development Laboratory (LRDE) 00002 // 00003 // This file is part of Olena. 00004 // 00005 // Olena is free software: you can redistribute it and/or modify it under 00006 // the terms of the GNU General Public License as published by the Free 00007 // Software Foundation, version 2 of the License. 00008 // 00009 // Olena is distributed in the hope that it will be useful, 00010 // but WITHOUT ANY WARRANTY; without even the implied warranty of 00011 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00012 // General Public License for more details. 00013 // 00014 // You should have received a copy of the GNU General Public License 00015 // along with Olena. If not, see <http://www.gnu.org/licenses/>. 00016 // 00017 // As a special exception, you may use this file as part of a free 00018 // software project without restriction. Specifically, if other files 00019 // instantiate templates or use macros or inline functions from this 00020 // file, or you compile this file and link it with other files to produce 00021 // an executable, this file does not by itself cause the resulting 00022 // executable to be covered by the GNU General Public License. This 00023 // exception does not however invalidate any other reasons why the 00024 // executable file might be covered by the GNU General Public License. 00025 00026 #include <cstdlib> 00027 #include <mln/accu/stat/var.hh> 00028 00029 00030 float my_rand(int c) 00031 { 00032 return (1 + c) * float(std::rand()) / RAND_MAX; 00033 } 00034 00035 00036 int main() 00037 { 00038 using namespace mln; 00039 00040 typedef algebra::vec<3,float> vec3f; 00041 00042 enum { n = 1000 }; 00043 vec3f v[n]; 00044 00045 for (int i = 0; i < n; ++i) 00046 { 00047 v[i][0] = my_rand(0); 00048 v[i][1] = my_rand(1); 00049 v[i][2] = my_rand(2); 00050 } 00051 00052 accu::stat::var<vec3f> a; 00053 for (int i = 0; i < n; ++i) 00054 a.take(v[i]); 00055 00056 mln_assertion(a.n_items() == n); 00057 00058 vec3f m = a.mean(); 00059 mln_assertion(m[0] > 0.4 && m[0] < 0.6); 00060 mln_assertion(m[1] > 0.9 && m[1] < 1.1); 00061 mln_assertion(m[2] > 1.4 && m[2] < 1.6); 00062 00063 algebra::mat<3,3,float> s_1 = a.variance()._1(); 00064 mln_assertion(s_1(0,0) > 11 && s_1(0,0) < 13); 00065 mln_assertion(s_1(1,1) > 2 && s_1(1,1) < 4); 00066 mln_assertion(s_1(2,2) > 1.1 && s_1(2,2) < 1.5); 00067 }