From: sgunderson@bigfoot.com <>
Date: Wed, 31 Jan 2007 20:13:08 +0000 (+0100)
Subject: Add a Pareto random number generator.
X-Git-Url: https://git.sesse.net/?a=commitdiff_plain;h=c7b34c15c91b9552b41463be7c2e598922874d4e;p=jam

Add a Pareto random number generator.
---

diff --git a/jam.c b/jam.c
index 3b38d02..e81e166 100644
--- a/jam.c
+++ b/jam.c
@@ -1,4 +1,5 @@
 #include <stdio.h>
+#include <math.h>
 #include <string.h>
 #include <stdlib.h>
 #include <getopt.h>
@@ -28,6 +29,12 @@ const static struct option longopts[] = {
 	{ NULL, 0, NULL, 0 }
 };
 
+double gen_pareto_random(double min, double k)
+{
+	double u = rand() / (RAND_MAX+1.0);
+	return min * pow(u, -1.0 / k);
+}
+
 void read_ip_list(char *filename, struct in_addr **addr_list, unsigned *num, unsigned *room)
 {
 	char buf[256];
@@ -113,7 +120,15 @@ void parse_options(int argc, char **argv)
 
 void *sender_worker(void *arg)
 {
+	unsigned i;
+
 	printf("Dummy sender worker\n");
+
+	for (i = 0; i < 1000; ++i) {
+		unsigned bytes = (unsigned)gen_pareto_random(1048576.0, 1.0);
+		printf("Sending %u bytes\n", bytes);
+	}
+
 	pthread_exit(0);
 }