NetBSD/dist/ntp/util/hist.c

/*	$NetBSD: hist.c,v 1.2 2003/12/04 16:23:38 drochner Exp $	*/

/*
 * This program can be used to calibrate the clock reading jitter of a
 * particular CPU and operating system. It first tickles every element
 * of an array, in order to force pages into memory, then repeatedly calls
 * gettimeofday() and, finally, writes out the time values for later
 * analysis. From this you can determine the jitter and if the clock ever
 * runs backwards.
 */

#if 0
#ifdef HAVE_CONFIG_H
# include <config.h>
#endif

#include "ntp_types.h"
#endif

#include <stdio.h>
#include <stdlib.h>
#include <time.h>

#define NBUF 100001		/* size of basic histogram */
#define NSRT 20000		/* size of overflow histogram */
#define NCNT (600 * 1000000)	/* sample interval (us) */

extern int col(hrtime_t *, hrtime_t *);
extern hrtime_t gethrtime(void);

int
main(
	int argc,
	char *argv[]
	)
{
	int i, j, n;
	hrtime_t t, u, v, w, gtod[NBUF], ovfl[NSRT];

	/*
	 * Force pages into memory
	 */
	for (i = 0; i < NBUF; i++)
	    gtod[i] = 0;
	for (i = 0; i < NSRT; i++)
	    ovfl[i] = 0;

	/*
	 * Construct histogram
	 */
	n = 0;
	t = gethrtime();
	v = t;
	while (1) {
		u = gethirestime();
		if (u - v > NCNT)
		    break;
		w = u - t;
		if (w <= 0) {
/*
			printf("error <= 0 %ld %d %d, %d %d\n", w, ts.tv_sec,
			       ts.tv_usec, tr.tv_sec, tr.tv_usec);
*/
		} else if (w > NBUF - 1) {
			ovfl[n] = w;
			if (n < NSRT - 1)
			    n++;
		} else {
			gtod[w]++;
		}
		t = u;
	}

	/*
	 * Write out histogram
	 */
	for (i = 0; i < NBUF - 1; i++) {
		if (gtod[i] > 0)
		    printf("%ld %ld\n", i, gtod[i]);
	}
	if (n == 0)
	    return;
	qsort((char *)ovfl, (size_t)n, sizeof(hrtime_t), col);
	w = 0;
	j = 0;
	for (i = 0; i < n; i++) {
		if (ovfl[i] != w) {
			if (j > 0)
			    printf("%ld %ld\n", w, j);
			w = ovfl[i];
			j = 1;
		} else
		    j++;
	}
	if (j > 0)
	    printf("%ld %ld\n", w, j);
 
	exit(0);
}

int
col(
	hrtime_t *x,
	hrtime_t *y
	)
{
	return (*x - *y);
}