#!/bin/bash
TMP_FILE=`mktemp /tmp/trw.XXXXXX`; trap 'rm -f $TMP_FILE' 0 1 15
rm /tmp/trw.dat
echo "# beta     time_MF errlogZ_MF maxdiff_MF     time_TRWUNI errlogZ_TRWUNI maxdiff_TRWUNI     time_TRWOPT errlogZ_TRWOPT maxdiff_TRWOPT     time_TRWOPTbetter errlogZ_TRWOPTbetter maxdiff_TRWOPTbetter     time_JJ errlogZ_JJ     time_JJopt errlogZ_JJopt     time_LK errorlogZ_LK" > /tmp/trw.dat
#echo "# beta     time_MF errlogZ_MF maxdiff_MF     time_JJ errlogZ_JJ" > /tmp/trw.dat
for k in -4 -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
            ../../utils/createfg --type grid --seed "$i"00 --subtype sg --N 25 --sigma_w $beta --sigma_th $beta --periodic 1 > $TMP_FILE
            echo -n "$beta " >> /tmp/trw.dat
            ../../tests/test --aliases ../../tests/aliases.conf --filename $TMP_FILE --methods JTREE_HUGIN MF_SEQRND | tail -1 | awk '{printf("    %d %e %e ", $2, $5, $6);}' >> /tmp/trw.dat
            ../../tests/test --aliases ../../tests/aliases.conf --filename $TMP_FILE --methods JTREE_HUGIN TRWUNI_SEQFIX | tail -1 | awk '{printf("    %d %e %e ", $2, $5, $6);}' >> /tmp/trw.dat
            ../../tests/test --aliases ../../tests/aliases.conf --filename $TMP_FILE --methods JTREE_HUGIN TRWOPT_SEQFIX | tail -1 | awk '{printf("    %d %e %e ", $2, $5, $6);}' >> /tmp/trw.dat
            ../../tests/test --aliases ../../tests/aliases.conf --filename $TMP_FILE --methods JTREE_HUGIN TRWOPT_SEQFIX[outertol=1e-9] | tail -1 | awk '{printf("    %d %e %e", $2, $5, $6);}' >> /tmp/trw.dat
            echo -n "    " >> /tmp/trw.dat
            ./errorbounds3 $TMP_FILE 0 0 5 0 0 >> /tmp/trw.dat
            echo >> /tmp/trw.dat
            if [ $? != 0 ]; then
                exit
            fi
        done
    done
done
rm -f $TMP_FILE
