From 6b34d9abdb5a44fb2bddb2656ae522c61c2ce1f8 Mon Sep 17 00:00:00 2001
From: NunoSempere <nuno.sempere@protonmail.com>
Date: Fri, 2 Jun 2023 12:50:51 -0600
Subject: [PATCH] feat: add more threads, document rand_r in code.

---
 C-optimized/makefile    |   4 ++++
 C-optimized/out/samples | Bin 22392 -> 22392 bytes
 C-optimized/samples.c   |   4 +++-
 3 files changed, 7 insertions(+), 1 deletion(-)

diff --git a/C-optimized/makefile b/C-optimized/makefile
index dc96b32a..ded7c6bb 100644
--- a/C-optimized/makefile
+++ b/C-optimized/makefile
@@ -46,6 +46,8 @@ multi:
 	OMP_NUM_THREADS=1 ./$(OUTPUT) && echo
 	OMP_NUM_THREADS=2 ./$(OUTPUT) && echo
 	OMP_NUM_THREADS=4 ./$(OUTPUT) && echo
+	OMP_NUM_THREADS=8 ./$(OUTPUT) && echo
+	OMP_NUM_THREADS=16 ./$(OUTPUT) && echo
 	./$(OUTPUT_ONE_THREAD) && echo
 
 time-linux:
@@ -53,6 +55,8 @@ time-linux:
 	OMP_NUM_THREADS=1 /bin/time -f "Time: %es" ./$(OUTPUT) && echo
 	OMP_NUM_THREADS=2 /bin/time -f "Time: %es" ./$(OUTPUT) && echo
 	OMP_NUM_THREADS=4 /bin/time -f "Time: %es" ./$(OUTPUT) && echo
+	OMP_NUM_THREADS=8 /bin/time -f "Time: %es" ./$(OUTPUT) && echo
+	OMP_NUM_THREADS=16 /bin/time -f "Time: %es" ./$(OUTPUT) && echo
 	/bin/time -f "Time: %es" ./$(OUTPUT_ONE_THREAD) && echo
 
 debian-install-dependencies:
diff --git a/C-optimized/out/samples b/C-optimized/out/samples
index 98384ccc85fcb36edf73f64cda3a7c998494fe59..e77b23001169a9830817da5010aeb0b897fdf288 100755
GIT binary patch
delta 74
zcmeydj`7Dj#tj+FB3HKrE#75w%R=1UY^L-nUPeh5L%q#4%&Pn>4bx^YOkOHf3}$SW
X5#GiKVO<lck_NIQyf=p!sHg(~I{g~s

delta 74
zcmeydj`7Dj#tj+FBHZ<s{rkl~KX7_-)|t07`)0aYy29ofW>tO`hKFBSCodH$1~WFx
Y2ybJAu&#+zNwX;I+2^x4#6U$I0F)6OfB*mh

diff --git a/C-optimized/samples.c b/C-optimized/samples.c
index 134c7dc2..d96bf519 100644
--- a/C-optimized/samples.c
+++ b/C-optimized/samples.c
@@ -6,7 +6,7 @@
 
 const float PI = 3.14159265358979323846;
 
-#define N 100000000
+#define N 10000000
 
 //Array helpers
 
@@ -50,6 +50,8 @@ void array_cumsum(float* array_to_sum, float* array_cumsummed, int length)
 float rand_float(float to, unsigned int* seed)
 {
     return ((float)rand_r(seed) / (float)RAND_MAX) * to;
+		// See: <https://stackoverflow.com/questions/43151361/how-to-create-thread-safe-random-number-generator-in-c-using-rand-r>
+		// rand() is not thread-safe, as it relies on (shared) hidden state.
 }
 
 float ur_normal(unsigned int* seed)