]> rtime.felk.cvut.cz Git - lisovros/iproute2_canprio.git/commitdiff
Rename: tc/normal.c -> netem/normal.c
authorosdl.net!shemminger <osdl.net!shemminger>
Wed, 9 Feb 2005 22:05:41 +0000 (22:05 +0000)
committerosdl.net!shemminger <osdl.net!shemminger>
Wed, 9 Feb 2005 22:05:41 +0000 (22:05 +0000)
(Logical change 1.141)

netem/normal.c

index e69de29bb2d1d6434b8b29ae775ad8c2e48c5391..e6683db8946680910815b05ac8be96b80449e819 100644 (file)
@@ -0,0 +1,56 @@
+/*
+ * Normal distribution table generator
+ * Taken from the uncopyrighted NISTnet code.
+ */
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+#include <string.h>
+#include <limits.h>
+
+#include <linux/types.h>
+#include <linux/pkt_sched.h>
+
+#define TABLESIZE 16384
+#define TABLEFACTOR NETEM_DIST_SCALE
+
+static double
+normal(double x, double mu, double sigma)
+{
+       return .5 + .5*erf((x-mu)/(sqrt(2.0)*sigma));
+}
+
+int
+main(int argc, char **argv)
+{
+       double x, *table;
+       int i, n;
+
+       table = calloc(sizeof(double), TABLESIZE+1);
+       if (!table) {
+               fprintf(stderr, "Not enough memory\n");
+               return 1;
+       }
+
+
+       for (x = -10.0; x < 10.05; x += .00005) {
+               i = (int)rint(TABLESIZE*normal(x, 0.0, 1.0));
+               table[i] = x;
+       }
+
+       
+       printf("# This is the distribution table for the normal distribution.\n");
+       for (i = n = 0; i < TABLESIZE; i += 4) {
+               int value = (int) rint(table[i]*TABLEFACTOR);
+               if (value < SHRT_MIN) value = SHRT_MIN;
+               if (value > SHRT_MAX) value = SHRT_MAX;
+
+               printf(" %d", value);
+               if (++n == 8) {
+                       putchar('\n');
+                       n = 0;
+               }
+       }
+       free(table);
+       return 0;
+}