#!/bin/bash
Rscript gen_dat.r
rm timefile.txt
echo "N_warmup","N_samples","max_tree","method","time_in_ms" >> timefile.txt
for T in 5 6 7 8 9 10 11 12 13; do
  n_warmup=$((2000/$T))
  n_samples=$((2000/$T))
  start=`date +%s%N`
  STAN_NUM_THREADS=-1 ./minttu_horseshoe_par \
    method=sample algorithm=hmc engine=nuts max_depth=$T \
    adapt delta=.99999	\
	num_warmup=$n_warmup num_samples=$n_samples \
	data file=data.Rdump
  end=`date +%s%N`
  runtime=$(($((end-start))/1000000))
  echo $n_warmup,$n_samples,$T,"par",$runtime >> timefile.txt
  start=`date +%s%N`
  ./minttu_horseshoe \
    method=sample algorithm=hmc engine=nuts max_depth=$T \
    adapt delta=.99999	\
	num_warmup=$n_warmup num_samples=$n_samples \
	data file=data.Rdump
  end=`date +%s%N`
  runtime=$(($((end-start))/1000000))
  echo $n_warmup,$n_samples,$T,"lin",$runtime >> timefile.txt
done;
Rscript post.r