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