automatic commit
[folded-ctf.git] / gaussian.cc
1
2 ///////////////////////////////////////////////////////////////////////////
3 // This program is free software: you can redistribute it and/or modify  //
4 // it under the terms of the version 3 of the GNU General Public License //
5 // as published by the Free Software Foundation.                         //
6 //                                                                       //
7 // This program is distributed in the hope that it will be useful, but   //
8 // WITHOUT ANY WARRANTY; without even the implied warranty of            //
9 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU      //
10 // General Public License for more details.                              //
11 //                                                                       //
12 // You should have received a copy of the GNU General Public License     //
13 // along with this program. If not, see <http://www.gnu.org/licenses/>.  //
14 //                                                                       //
15 // Written by Francois Fleuret                                           //
16 // (C) Idiap Research Institute                                          //
17 //                                                                       //
18 // Contact <francois.fleuret@idiap.ch> for comments & bug reports        //
19 ///////////////////////////////////////////////////////////////////////////
20
21 #include "gaussian.h"
22
23 Gaussian::Gaussian() {
24   _nb_samples = 0;
25   _sum = 0.0;
26   _sum_sq = 0.0;
27 }
28
29 void Gaussian::add_sample(scalar_t x) {
30   _nb_samples++;
31   _sum += x;
32   _sum_sq += x * x;
33 }
34
35 scalar_t Gaussian::expectation() {
36   return _sum / scalar_t(_nb_samples);
37 }
38
39 scalar_t Gaussian::variance() {
40   scalar_t e = _sum / scalar_t(_nb_samples);
41   return (_sum_sq - _sum * e) / scalar_t(_nb_samples - 1);
42 }
43
44 scalar_t Gaussian::standard_deviation() {
45   return sqrt(variance());
46 }
47