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, (C) IDIAP                                //
  16 // Contact <francois.fleuret@idiap.ch> for comments & bug reports        //
  17 ///////////////////////////////////////////////////////////////////////////
  18
  19 #include "gaussian.h"
  20
  21 Gaussian::Gaussian() {
  22   _nb_samples = 0;
  23   _sum = 0.0;
  24   _sum_sq = 0.0;
  25 }
  26
  27 void Gaussian::add_sample(scalar_t x) {
  28   _nb_samples++;
  29   _sum += x;
  30   _sum_sq += x * x;
  31 }
  32
  33 scalar_t Gaussian::expectation() {
  34   return _sum / scalar_t(_nb_samples);
  35 }
  36
  37 scalar_t Gaussian::variance() {
  38   scalar_t e = _sum / scalar_t(_nb_samples);
  39   return (_sum_sq - _sum * e) / scalar_t(_nb_samples - 1);
  40 }
  41
  42 scalar_t Gaussian::standard_deviation() {
  43   return sqrt(variance());
  44 }
  45