#include <math.h>
#include <stats.h>
#include <stdlib.h>
#include <dbg.h>
Stats *Stats_recreate(double sum, double sumsq, unsigned long n,
double min, double max)
{
Stats *st = malloc(sizeof(Stats));
check_mem(st);
st->sum = sum;
st->sumsq = sumsq;
st->n = n;
st->min = min;
st->max = max;
return st;
error:
return NULL;
}
Stats *Stats_create()
{
return Stats_recreate(0.0, 0.0, 0L, 0.0, 0.0);
}
void Stats_sample(Stats *st, double s)
{
st->sum += s;
st->sumsq += s * s;
if (st->n == 0) {
st->min = s;
st->max = s;
} else {
if (st->min > s)
st->min = s;
if (st->max < s)
st->max = s;
}
st->n += 1;
}
void Stats_dump(Stats *st)
{
fprintf(stderr,
"sum: %f, sumsq: %f, n: %ld, "
"min: %f, max: %f, mean: %f, stddev: %f",
st->sum, st->sumsq, st->n, st->min, st->max, Stats_mean(st),
Stats_stddev(st));
}