#include "dataproc.h"
#include <iostream>
#include <math.h>
#include <vector>
#include <cstdlib>
using namespace std;


void Dataproc::add_value(double d)
{
  v.push_back(d);
}

int Dataproc::nentries()
{
  return v.size();
}

double Dataproc::min()
{
  double m;
  int i;

  if (v.size() == 0) {
    cerr << "Dataproc::min() - Empty Dataproc\n";
    exit (1);
  }

  m = v[0];
  for (i = 1; i < v.size(); i++) if(v[i] < m) m = v[i];
  return m;
}

double Dataproc::max()
{
  double m;
  int i;

  if (v.size() == 0) {
    cerr << "Dataproc::max() - Empty Dataproc\n";
    exit (1);
  }

  m = v[0];
  for (i = 1; i < v.size(); i++) if(v[i] > m) m = v[i];
  return m;
}

double Dataproc::mean()
{
  int i;
  double n, total;

  n = v.size();
  if (n == 0) {
    cerr << "Dataproc::mean() - Empty Dataproc\n";
    exit (1);
  }
  total = 0;
  for (i = 0; i < v.size(); i++) total += v[i];
  return total/n;
}

double Dataproc::variance()
{
  double n, vtotal, m, diff;
  int i;

  m = mean();
  n  = v.size();
  vtotal = 0;
  for (i = 0; i < v.size(); i++) {
    diff = (v[i] - m);
    vtotal += (diff*diff);
  }
  return vtotal/n;
}

double Dataproc::stddev()
{
  return sqrt(variance());
}