#!/bin/bash
TMP_FILE=`mktemp /tmp/bounds.XXXXXX`; trap 'rm -f $TMP_FILE' 0 1 15
outfile=/tmp/bounds2.dat
rm $outfile
for k in -3 -2 -1 0 1; do
    for j in `seq -w 1 9`; do
        beta="0.$j"E"$k"
        for i in `seq 1 1`; do
            echo $beta
            let "seed = k + 5"
            let "seed = seed * 10 + j"
            let "seed = seed * 100"
            seed=100
            ../../utils/createfg --type grid --seed $seed --subtype sg --N 25 --sigma_w $beta --sigma_th $beta --periodic 0  > $TMP_FILE
            echo -n "$beta    " >> $outfile
            ./errorbounds5 $TMP_FILE 1 MF TRW 0 0 >> $outfile
#            cd ~/bound_propagation_1.0
#            domatlab "fast=0; tol=1e-6; gap=1; filename='$TMP_FILE'; boundpropfg" | tail -1 >> $outfile
#            cd ~/svn/libDAI/tests/errorbounds
            echo >> $outfile
            if [ $? != 0 ]; then
                exit
            fi
        done
    done
done
rm -f $TMP_FILE
# beta BPerr(5x) DTgap DTMgap Cavitygap LUgap BoundPropgap
# 5x = min max avg stddev t
