--- a/COBRAxy/flux_simulation_beta.py	Mon Sep 22 16:21:14 2025 +0000
+++ b/COBRAxy/flux_simulation_beta.py	Mon Sep 22 16:29:01 2025 +0000
@@ -212,7 +212,7 @@
     
     for i in range(n_batches):
         batch_filename = f"{ARGS.output_path}/{model_name}_{i}_OPTGP.npy"
-        batch_data = np.load(batch_filename)
+        batch_data = np.load(batch_filename, allow_pickle=True)
         all_samples.append(batch_data)
     
     # Concatenate all batches
@@ -289,13 +289,13 @@
     
     for i in range(n_batches):
         batch_filename = f"{ARGS.output_path}/{model_name}_{i}_CBS.npy"
-        batch_data = np.load(batch_filename)
+        batch_data = np.load(batch_filename, allow_pickle=True)
         all_samples.append(batch_data)
     
     # Concatenate all batches
     samplesTotal_array = np.vstack(all_samples)
     
-    # Convert back to DataFrame with proper column names
+    # Convert back to DataFrame with proper column namesq
     samplesTotal = pd.DataFrame(samplesTotal_array, columns=reaction_ids)
     
     # Save the final merged result as CSV