Ticket #7150: TestVar.cpp

#include <iostream>
#include <stdlib.h>
#include <boost/accumulators/accumulators.hpp>
#include <boost/accumulators/statistics/stats.hpp>
#include <boost/accumulators/statistics/mean.hpp>
#include <boost/accumulators/statistics/variance.hpp>

using namespace boost::accumulators;

typedef accumulator_set<double, stats<tag::variance> > Var;
typedef accumulator_set<double, stats<tag::variance, tag::mean> > VarMean;

Var accumulate_Var(double *v, int n) {
  Var acc;
  for(unsigned int i=0; i<n; i++) {
       acc(v[i]);
  }
  return acc;
}

VarMean accumulate_VarMean(double *v, int n) {
  VarMean acc;
  for(unsigned int i=0; i<n; i++) {
       acc(v[i]);
  }
  return acc;
}

int main(int argc, char *argv[]) {
   int n=argc-1;
   double x[n];
   for (int i=1; i<=n; i++) {
     x[i-1] = atoi(argv[i]);
   }

   std::cout << "    var: " << variance(accumulate_Var(x, n)) << std::endl;
   std::cout << "VarMean: " << variance(accumulate_VarMean(x, n)) << std::endl;

   return 0;
}