Include replicate number in output and add merge result.

Author: highperformancecoder

models/SAR.py

@@ -12,6 +12,8 @@
 # we want initialisation to be identical across all processes
 #randomSeed(1)
 
+replicate=sys.argv[1] if len(sys.argv)>=2 else 1
+
 # we want the array operations to have a different seed across processes
 array_urand.seed(int(1024*random()))
 
@@ -33,7 +35,7 @@ def av(x):
 #mut_max=float(sys.argv[3])   # mutation rate
 #migration=float(sys.argv[4]) # initial migration rate
 
-print("Area, Mutation rate, Migration rate, Number of species, Interaction strength^2",flush=True)
+print("Replicate, Area, Mutation rate, Migration rate, Number of species, Interaction strength^2",flush=True)
 for A in [1, 2, 4, 6, 9, 12, 16]:
     for mut_max in [1e-5, 1e-4, 1e-3]:
         for migration in [1e-6, 1e-5]:
@@ -73,7 +75,7 @@ def av(x):
                     step()
 
                 conn=sum([x*x for x in ecolab.interaction.val()])/len(ecolab.species)**2 if len(ecolab.species)>0 else 0
-                print(numX*numY, av(ecolab.mutation()), av(ecolab.migration()), len(ecolab.species), conn, flush=True,sep=',')
+                print(replicate,numX*numY, av(ecolab.mutation()), av(ecolab.migration()), len(ecolab.species), conn, flush=True,sep=',')
 
 
 

models/SAR.sh

@@ -1,6 +1,12 @@
 i=0
 while [ $[i++] -lt 8 ]; do
-    python3 SAR.py > SARResults_$i.csv &
+    python3 SAR.py $i > SARResults_$i.csv &
     sleep 1
 done
 wait
+
+head -1 SARResults_1.csv >SARResults.csv
+i=0
+while [ $[i++] -lt 8 ]; do
+    tail -n +2  SARResults_$i.csv >>SARResults.csv
+done