diff --git a/sw/blas/gemm/scripts/datagen.py b/sw/blas/gemm/scripts/datagen.py
index 2eb6e2f4db..bead5121d5 100755
--- a/sw/blas/gemm/scripts/datagen.py
+++ b/sw/blas/gemm/scripts/datagen.py
@@ -98,9 +98,9 @@ def emit_header(self, **kwargs):
 
         ctype = du.ctype_from_precision_t(prec)
 
-        a = du.generate_random_array((M, K), prec)
-        b = du.generate_random_array((K, N), prec)
-        c = du.generate_random_array((M, N), prec)
+        a = du.generate_random_array((M, K), prec, seed=42)
+        b = du.generate_random_array((K, N), prec, seed=42)
+        c = du.generate_random_array((M, N), prec, seed=42)
         result = self.exact_golden_model(1, a, b, kwargs['beta'], c)
 
         # Store matrices in transposed form if requested
diff --git a/sw/blas/gemm/scripts/verify.py b/sw/blas/gemm/scripts/verify.py
index 353ea13284..a89b4dbc34 100755
--- a/sw/blas/gemm/scripts/verify.py
+++ b/sw/blas/gemm/scripts/verify.py
@@ -18,7 +18,7 @@ class GemmVerifier(Verifier):
     OUTPUT_UIDS = ['c']
     ERR_THRESHOLD = {
         1: 1e-4,
-        2: 8e-2,
+        2: 5e-1,
         4: 1e-3,
         8: 1e-3
     }
diff --git a/util/sim/data_utils.py b/util/sim/data_utils.py
index 3b732c5ccd..fce39a9b55 100644
--- a/util/sim/data_utils.py
+++ b/util/sim/data_utils.py
@@ -118,7 +118,7 @@ def ctype_from_precision_t(prec):
     return precision_t_to_ctype_map[_integer_precision_t(prec)]
 
 
-def generate_random_array(size, prec='FP64'):
+def generate_random_array(size, prec='FP64', seed=None):
     """Consistent random array generation for Snitch experiments.
 
     Samples values between -1 and 1 from a uniform distribution and
@@ -133,7 +133,7 @@ def generate_random_array(size, prec='FP64'):
             (e.g. "FP64") and integer enumeration values (e.g. 8).
     """
     # Generate in 64b precision and then cast down
-    rand = np.random.default_rng().random(size=size, dtype=np.float64) * 2 - 1
+    rand = np.random.default_rng(seed=seed).random(size=size, dtype=np.float64) * 2 - 1
     # Generate FlexFloat array for 8b floats, casted from 16b Numpy array
     if _integer_precision_t(prec) == 1:
         return ff.array(rand.astype(np.float16), ff_desc_from_precision_t(prec))