tests/variable.c


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42

/* -*- mode: C; c-basic-offset: 2; indent-tabs-mode: nil; -*- */

#include <math.h>
#include "variable.h"

void
variable_init (Variable *variable)
{
  variable->weight = 0.0;
  variable->sum = 0.0;
  variable->sum2 = 0.0;
}

void
variable_add_weighted (Variable *variable,
                       double    value,
                       double    weight)
{
  variable->weight += weight;
  variable->sum += weight * value;
  variable->sum2 += weight * value * value;
}

void
variable_add (Variable *variable,
              double    value)
{
  variable_add_weighted (variable, value, 1.);
}

double
variable_mean (Variable *variable)
{
  return variable->sum / variable->weight;
}

double
variable_standard_deviation (Variable *variable)
{
  double mean = variable_mean (variable);
  return sqrt (variable->sum2 / variable->weight - mean * mean);
}