experiment.sh

#!/bin/sh

echo "trial,block,thread,our_time,their_time" > analysis.csv
for trials in $(seq 500 500 2500)
#trials=$(($trials*1000))
do 
  for i in $(seq 1 10)
    do 
      for jj in $(seq 0 5)
        do
          j=$((2**jj))
          c=$({  time -p  /scratch/hpc_proj_yuz12012_rjm11010/mcmc/seq_gibbs data.txt $trials >/dev/null; } 2>&1 | grep real |awk '{print $2}')
          a=$({  time -p ./mygibbs data.txt $trials $i $j >/dev/null; } 2>&1 | grep real |awk '{print $2}')
          b=$({  time -p ./gibbs data.txt $trials >/dev/null; } 2>&1 | grep real |awk '{print $2}')
	  echo "$trials,$i,$j,$c,$a,$b" >> analysis.csv

        done
    done             
done