blob: d4335bd4e2d697ab53d444dfa7fb83e49ce0283c (
plain)
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
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
|
/* libguestfs - guestfish shell
* Copyright (C) 2010 Red Hat Inc.
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License along
* with this program; if not, write to the Free Software Foundation, Inc.,
* 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
*/
#ifndef FISH_RMSD_H
#define FISH_RMSD_H
/* Compute the running mean and standard deviation from the
* series of estimated values.
*
* Method:
* http://en.wikipedia.org/wiki/Standard_deviation#Rapid_calculation_methods
* Checked in a test program against answers given by Wolfram Alpha.
*/
struct rmsd {
double a; /* mean */
double i; /* number of samples */
double q;
};
static void
rmsd_init (struct rmsd *r)
{
r->a = 0;
r->i = 1;
r->q = 0;
}
static void
rmsd_add_sample (struct rmsd *r, double x)
{
double a_next, q_next;
a_next = r->a + (x - r->a) / r->i;
q_next = r->q + (x - r->a) * (x - a_next);
r->a = a_next;
r->q = q_next;
r->i += 1.0;
}
static double
rmsd_get_mean (const struct rmsd *r)
{
return r->a;
}
static double
rmsd_get_standard_deviation (const struct rmsd *r)
{
return sqrt (r->q / (r->i - 1.0));
}
#endif /* FISH_RMSD_H */
|
