export ALGENCAN=$HOME/algencan-4.0.0

gfortran -c -O3 covering-america.f90
gfortran -c -O3 modamerica.f90
gfortran -c -O3 vorintpols.f90
gfortran -c -O3 drand.f90
gfortran -c -O3 geometry.f90
gfortran -c -O3 geompack2.f90

gfortran covering-america.o modamerica.o vorintpols.o drand.o geometry.o geompack2.o \
         -L$ALGENCAN/sources/algencan/lib -L$ALGENCAN/sources/hsl/lib \
	 -L$ALGENCAN/sources/blas/lib -lalgencan -lhsl -lblas -o covering-america

for nballs in `seq 1 1 100`
do
  export bestrad=1000000.0d0
  for itrial in `seq 1 1 10000`
  do
    echo $nballs $itrial $bestrad | timeout --preserve-status 30 ./covering-america
    export statdtris=$?
    if [ -e tabline.txt ] ; then
      sed -e 's/$/  S/' -i tabline.txt
      cat tabline.txt >> table-covering-america.txt ;
      bestrad=`cat bestrad.txt`
    else
      if [ $statdtris -eq 0 ] && [ -e alsolver-interrupted-tabline.txt ]; then
        sed -e 's/$/  S/' -i alsolver-interrupted-tabline.txt
        cat alsolver-interrupted-tabline.txt >> table-covering-america.txt ;
      elif [ $statdtris != 0 ] && [ -e alsolver-interrupted-tabline.txt ]; then
        sed -e 's/$/  F/' -i alsolver-interrupted-tabline.txt
        cat alsolver-interrupted-tabline.txt >> table-covering-america.txt ;
      else
        printf " no information\n" >> table-covering-america.txt ;
      fi
    fi
    rm -f tabline.txt alsolver-interrupted-tabline.txt \
          solver-interrupted-tabline.txt newtkkt-interrupted-tabline.txt
  done
done
