export ALABB=$HOME/Desktop/bmpru
export SRC=$ALABB/src
export INTLIB=$ALABB/intlib
export PROBS=$ALABB/problems
export BLAS=$ALABB/blas
export LAPACK=$ALABB/lapack
export QP=$ALABB/qp
export GENLIN=$ALABB/genlin
export SIMPLEX=$ALABB/nr-simplx
export MINOS=$ALABB/minos
export CURREN=$ALABB/experiments

gfortran -O4 $INTLIB/*.o $BLAS/*.o $LAPACK/*.o $MINOS/minos.o \
             $QP/*.o $GENLIN/genlin.o $SIMPLEX/*.o $SRC/*.o \
             $PROBS/ballswithinellipse.f -o alabbma

if [ ! -e $CURREN/alabb-table.out ] ; then
  touch $CURREN/alabb-table.out ;
fi

for (( nballs=6; nballs <= 6; nballs++ ))
  do
    let a=4 ;
    let b=2 ;
    echo "$a $b $nballs" > ellipses.dat ;
    $CURREN/alabbma < ellipses.dat ;

    printf "%.1f %.1f %1d" $a $b $nballs >> $CURREN/alabb-table.out ;
    cat $CURREN/alabb-tabline.out        >> $CURREN/alabb-table.out ;

    mv solution.txt solution-$a-$b-$nballs.txt ;

    let a=3 ;
    let b=2 ;
    echo "$a $b $nballs" > ellipses.dat ;
    $CURREN/alabbma < ellipses.dat ;

    printf "%.1f %.1f %1d" $a $b $nballs >> $CURREN/alabb-table.out ;
    cat $CURREN/alabb-tabline.out        >> $CURREN/alabb-table.out ;

    mv solution.txt solution-$a-$b-$nballs.txt ;

    let a=2 ;
    let b=2 ;
    echo "$a $b $nballs" > ellipses.dat ;
    $CURREN/alabbma < ellipses.dat ;

    printf "%.1f %.1f %1d" $a $b $nballs >> $CURREN/alabb-table.out ;
    cat $CURREN/alabb-tabline.out        >> $CURREN/alabb-table.out ;

    mv solution.txt solution-$a-$b-$nballs.txt ;

    let a=2 ;
    let b=1 ;
    echo "$a $b $nballs" > ellipses.dat ;
    $CURREN/alabbma < ellipses.dat ;

    printf "%.1f %.1f %1d" $a $b $nballs >> $CURREN/alabb-table.out ;
    cat $CURREN/alabb-tabline.out        >> $CURREN/alabb-table.out ;

    mv solution.txt solution-$a-$b-$nballs.txt ;
  done