gaussian.cc

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