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 version 3 as
  13  *  published by the Free Software Foundation.
  14  *
  15  *  folded-ctf is distributed in the hope that it will be useful, but
  16  *  WITHOUT ANY WARRANTY; without even the implied warranty of
  17  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  18  *  General Public License for more details.
  19  *
  20  *  You should have received a copy of the GNU General Public License
  21  *  along with folded-ctf.  If not, see <http://www.gnu.org/licenses/>.
  22  *
  23  */
  24
  25 #include "gaussian.h"
  26
  27 Gaussian::Gaussian() {
  28   _nb_samples = 0;
  29   _sum = 0.0;
  30   _sum_sq = 0.0;
  31 }
  32
  33 void Gaussian::add_sample(scalar_t x) {
  34   _nb_samples++;
  35   _sum += x;
  36   _sum_sq += x * x;
  37 }
  38
  39 scalar_t Gaussian::expectation() {
  40   return _sum / scalar_t(_nb_samples);
  41 }
  42
  43 scalar_t Gaussian::variance() {
  44   scalar_t e = _sum / scalar_t(_nb_samples);
  45   return (_sum_sq - _sum * e) / scalar_t(_nb_samples - 1);
  46 }
  47
  48 scalar_t Gaussian::standard_deviation() {
  49   return sqrt(variance());
  50 }
  51