statist.cc

TCP westwood code, download

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

// in units of 0.1 milliseconds
#define MAX_DELAY 5000 

int main(int argc, char* argv[])
{
	FILE *finput;
	double time;
	double delay;
	unsigned int station;
	unsigned int traffic_class;
	double avg;
	double sigma;
	long total_packets;

	double ignore_time=0.0;
	if (argc==3) {
		ignore_time=atof(argv[2]);
		if (ignore_time<0.0) return 1;
	}
	total_packets=0; avg=0.0; sigma=0.0;
	printf("Stat for %s (from time %lf): ",argv[1],ignore_time);
	finput=fopen(argv[1],"r");
	while (feof(finput)==0) {
		fscanf(finput,"%lf %lf %d:%d\n",&time,&delay,&station,&traffic_class);
		if (time<ignore_time) continue;
		avg+=delay;
		total_packets+=1;
	}
	fclose(finput);
	avg=avg/total_packets;
	//printf("\t..total number of packets processed for %d:%d is %d\n",station,traffic_class,total_packets);
	printf("avg delay / stnd dev = %2.9lf ",avg);
	finput=fopen(argv[1],"r");
	while (feof(finput)==0) {
		fscanf(finput,"%lf %lf %d:%d\n",&time,&delay,&station,&traffic_class);
		if (time<ignore_time) continue;
		sigma+=pow(delay-avg,2.0);
	}
	fclose(finput);
	sigma=sqrt(sigma/total_packets);
	printf("%2.9lf\n",sigma);
	return 0;
}