Commit fdea8abe authored by Rusty Russell's avatar Rusty Russell

tally: simplify, add Y axis to graph.

parent 1c7178a0
...@@ -15,12 +15,11 @@ struct tally { ...@@ -15,12 +15,11 @@ struct tally {
ssize_t min, max; ssize_t min, max;
size_t total[2]; size_t total[2];
/* This allows limited frequency analysis. */ /* This allows limited frequency analysis. */
size_t buckets; unsigned buckets, step_bits;
size_t step_bits;
size_t counts[1 /* [buckets] */ ]; size_t counts[1 /* [buckets] */ ];
}; };
struct tally *tally_new(size_t buckets) struct tally *tally_new(unsigned buckets)
{ {
struct tally *tally; struct tally *tally;
...@@ -29,20 +28,8 @@ struct tally *tally_new(size_t buckets) ...@@ -29,20 +28,8 @@ struct tally *tally_new(size_t buckets)
return NULL; return NULL;
tally = malloc(sizeof(*tally) + sizeof(tally->counts[0])*buckets); tally = malloc(sizeof(*tally) + sizeof(tally->counts[0])*buckets);
if (tally) { if (tally) {
/* SSIZE_MAX isn't portable, so make it one of these types. */ tally->max = ((size_t)1 << (SIZET_BITS - 1));
BUILD_ASSERT(sizeof(tally->min) == sizeof(int) tally->min = ~tally->max;
|| sizeof(tally->min) == sizeof(long)
|| sizeof(tally->min) == sizeof(long long));
if (sizeof(tally->min) == sizeof(int)) {
tally->min = INT_MAX;
tally->max = INT_MIN;
} else if (sizeof(tally->min) == sizeof(long)) {
tally->min = LONG_MAX;
tally->max = LONG_MIN;
} else if (sizeof(tally->min) == sizeof(long long)) {
tally->min = (ssize_t)LLONG_MAX;
tally->max = (ssize_t)LLONG_MIN;
}
tally->total[0] = tally->total[1] = 0; tally->total[0] = tally->total[1] = 0;
/* There is always 1 bucket. */ /* There is always 1 bucket. */
tally->buckets = buckets+1; tally->buckets = buckets+1;
...@@ -444,7 +431,7 @@ char *tally_histogram(const struct tally *tally, ...@@ -444,7 +431,7 @@ char *tally_histogram(const struct tally *tally,
memcpy(tmp->counts, tally->counts, memcpy(tmp->counts, tally->counts,
sizeof(tally->counts[0]) * tmp->buckets); sizeof(tally->counts[0]) * tmp->buckets);
while ((max_bucket = get_max_bucket(tmp)) >= height) while ((max_bucket = get_max_bucket(tmp)) >= height)
renormalize(tmp, tmp->min, tmp->max *= 2); renormalize(tmp, tmp->min, tmp->max * 2);
/* Restore max */ /* Restore max */
tmp->max = tally->max; tmp->max = tally->max;
tally = tmp; tally = tmp;
...@@ -463,23 +450,29 @@ char *tally_histogram(const struct tally *tally, ...@@ -463,23 +450,29 @@ char *tally_histogram(const struct tally *tally,
free(tmp); free(tmp);
return NULL; return NULL;
} }
for (i = 0; i < height; i++) { for (i = 0; i < height; i++) {
unsigned covered = 0; unsigned covered = 1;
count = (double)tally->counts[i] / largest_bucket * width; count = (double)tally->counts[i] / largest_bucket * (width-1)+1;
if (i == 0) if (i == 0)
covered = snprintf(p, width, "%zi", tally->min); covered = snprintf(p, width, "%zi", tally->min);
else if (i == height - 1) else if (i == height - 1)
covered = snprintf(p, width, "%zi", tally->max); covered = snprintf(p, width, "%zi", tally->max);
if (covered) { else if (i == bucket_of(tally->min, tally->step_bits, 0))
if (covered > width) *p = '+';
covered = width; else
p += covered; *p = '|';
if (count > covered)
count -= covered; if (covered > width)
else covered = width;
count = 0; p += covered;
}
if (count > covered)
count -= covered;
else
count = 0;
memset(p, '*', count); memset(p, '*', count);
p += count; p += count;
*p = '\n'; *p = '\n';
......
...@@ -10,9 +10,11 @@ struct tally; ...@@ -10,9 +10,11 @@ struct tally;
* *
* This allocates a tally structure using malloc(). The greater the value * This allocates a tally structure using malloc(). The greater the value
* of @buckets, the more accurate tally_approx_median() and tally_approx_mode() * of @buckets, the more accurate tally_approx_median() and tally_approx_mode()
* and tally_graph() will be, but more memory is consumed. * and tally_histogram() will be, but more memory is consumed. If you want
* to use tally_histogram(), the optimal bucket value is the same as that
* @height argument.
*/ */
struct tally *tally_new(size_t buckets); struct tally *tally_new(unsigned int buckets);
/** /**
* tally_add - add a value. * tally_add - add a value.
......
...@@ -8,7 +8,7 @@ int main(void) ...@@ -8,7 +8,7 @@ int main(void)
char *graph, *p; char *graph, *p;
bool trunc; bool trunc;
plan_tests(100 + 1 + 10 + 1 + 100 + 1 + 10 + 1 + 10 + 2 + 1); plan_tests(100 + 1 + 10 + 1 + 100 + 1 + 10 + 1 + 10 * 2 + 1);
/* Uniform distribution, easy. */ /* Uniform distribution, easy. */
tally = tally_new(100); tally = tally_new(100);
...@@ -100,11 +100,16 @@ int main(void) ...@@ -100,11 +100,16 @@ int main(void)
/* Check min/max labels. */ /* Check min/max labels. */
if (i == 0) if (i == 0)
ok1(strncmp(p, "-5*", 3) == 0); ok1(strncmp(p, "-5*", 3) == 0);
if (i == 9) else if (i == 9)
ok1(strncmp(p, "4*", 2) == 0); ok1(strncmp(p, "4*", 2) == 0);
else if (i == 5)
ok1(p[0] == '+'); /* 0 marker */
else
ok1(p[0] == '|');
p = eol + 1; p = eol + 1;
} }
ok1(!*p); ok1(!*p);
diag("Here's the pretty: %s", graph);
free(graph); free(graph);
free(tally); free(tally);
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment