#!/usr/bin/env bash

BASE="ftp://math.nist.gov/pub/MatrixMarket2/Harwell-Boeing/"
                
MATRIX_SOLVERS_BENCHMARK="tnl-matrix-solvers-benchmark-dbg"
STOP_TIME="1"

SOLVERS="sor cg bicgstab gmres"
#SOLVERS="cg"
SOR_OMEGAS="0.4 0.8 1.2 1.6 2.0"
GMRES_RESTARTINGS="8 16 32 64 128 256"

LOG_FILE="matrix-solvers-benchmark.html"

export CUDA_PROFILE=0
export CUDA_PROFILE_CONFIG="$TNL_SOURCE_DIR/tests/cuda-profiler.conf"
PROCESS_CUDA_PROFILE="$TNL_SOURCE_DIR/tests/process-cuda-profile.pl"
source ./matrix-market
source ./florida-matrix-market

export TNL_SPARSE_MATRIX_CHECK_CFG_DESC_FILE="$TNL_SOURCE_DIR/tests/tnl-sparse-matrix-check.cfg.desc"

#MM_MATRICES=""
#FLORIDA_MM_MATRICES=""


benchmark_matrix()
# $1 input file
# $2 input mtx file
# $3 complete log file
# $4 verbose
{
   echo "          <tr>" >> $3
   echo "             <td>$2</td>" >> $3
   $MATRIX_SOLVERS_BENCHMARK --input-file $1 \
                             --input-mtx-file $2 \
                             --solver-name none\
                             --matrix-stats-file matrix-stats.html \
                             --verbose $4
   cat matrix-stats.html >> $3
   for solver in $SOLVERS;
   do     
      echo "Testing with $solver solver." 
      if [ x$solver == xcg ] || [ x$solver == xbicgstab ];
      then
         logfile="$2.results.$solver.log"
         if test ! -e $logfile;
         then
            $MATRIX_SOLVERS_BENCHMARK --input-file $1 \
                                      --input-mtx-file $2 \
                                      --solver-name $solver \
                                      --log-file $logfile \
                                      --verbose $4
            echo ""
         else
            echo "Test has already passed."
         fi
         cat $logfile >> $3
      fi
      if test x$solver == xsor;
      then
         for omega in $SOR_OMEGAS;
         do
            echo "Setting sor-omega to $omega"
            logfile="$2.results.sor.omega-$omega.log"
            if test ! -e $logfile;
            then
               $MATRIX_SOLVERS_BENCHMARK --input-file $1 \
                                         --input-mtx-file $2 \
                                         --solver-name $solver \
                                         --sor-omega $omega \
                                         --log-file $logfile \
                                         --verbose $4
               echo ""
            else
               echo "Test has already passed."
            fi
            cat $logfile >> $3
         done
      fi
      if test x$solver == xgmres;
      then
         for gmres_restarting in $GMRES_RESTARTINGS;
         do
            echo "Setting GMRES restarting to $gmres_restarting"
            logfile="$2.results.gmres.restarting-$gmres_restarting.log"
            if test ! -e $logfile;
            then
               $MATRIX_SOLVERS_BENCHMARK --input-file $1 \
                                         --input-mtx-file $2 \
                                         --solver-name $solver \
                                         --gmres-restarting $gmres_restarting \
                                         --log-file $logfile \
                                         --verbose $4
               echo ""
            else
               echo "Test has already passed."
            fi
            cat $logfile >> $3
         done
      fi
   done
   echo "          </tr>" >> $3
}

write_log_header()
# $1 log file name
{
   echo "<html>" > $1
   echo "   <body>" >> $1
   echo "      <table border=1>" >> $1
   echo "          <tr>" >> $1
   echo "             <td rowspan=1 colspan=3 align=center>Matrix</td>" >> $1
   echo "             <td rowspan=1 colspan=15 align=center>SOR</td>" >> $1
   echo "             <td rowspan=2 colspan=3 align=center>CG</td>" >> $1
   echo "             <td rowspan=2 colspan=3 align=center>BICGStab</td>" >> $1
   echo "             <td rowspan=1 colspan=18 align=center>GMRES</td>" >> $1
   echo "          </tr>" >> $1
         
   echo "          <tr>" >> $1
   echo "             <td rowspan=2>Name</td>" >> $1                      # Matrix description
   echo "             <td rowspan=2>Size</td>" >> $1
   echo "             <td rowspan=2>NonZeros No.</td>" >> $1
   echo "             <td colspan=3>0.4</td>" >> $1                      # SOR Omegas
   echo "             <td colspan=3>0.8</td>" >> $1                      # SOR Omegas
   echo "             <td colspan=3>1.2</td>" >> $1                      # SOR Omegas
   echo "             <td colspan=3>1.6</td>" >> $1                      # SOR Omegas
   echo "             <td colspan=3>2.0</td>" >> $1                      # SOR Omegas
#   echo "             <td colspan=3></td>" >> $1                         # CG
#   echo "             <td colspan=3></td>" >> $1                         # BICGStab
   echo "             <td colspan=3>8</td>" >> $1                        # GMRES restartings
   echo "             <td colspan=3>16</td>" >> $1                       # GMRES restartings
   echo "             <td colspan=3>32</td>" >> $1                       # GMRES restartings
   echo "             <td colspan=3>64</td>" >> $1                       # GMRES restartings
   echo "             <td colspan=3>128</td>" >> $1                      # GMRES restartings
   echo "             <td colspan=3>256</td>" >> $1                      # GMRES restartings
   echo "          </tr>" >> $1
   
   echo "          <tr>" >> $1
   echo "             <td>Residue</td>" >> $1                           # SOR Results omega = 0.4
   echo "             <td>Iterations</td>" >> $1
   echo "             <td>CPU</td>" >> $1
   echo "             <td>Residue</td>" >> $1                           # SOR Results omega = 0.8
   echo "             <td>Iterations</td>" >> $1
   echo "             <td>CPU</td>" >> $1
   echo "             <td>Residue</td>" >> $1                           # SOR Results omega = 1.2
   echo "             <td>Iterations</td>" >> $1
   echo "             <td>CPU</td>" >> $1
   echo "             <td>Residue</td>" >> $1                           # SOR Results omega = 1.6
   echo "             <td>Iterations</td>" >> $1
   echo "             <td>CPU</td>" >> $1
   echo "             <td>Residue</td>" >> $1                           # SOR Results omega = 2.0
   echo "             <td>Iterations</td>" >> $1
   echo "             <td>CPU</td>" >> $1
   echo "             <td>Residue</td>" >> $1                           # CG
   echo "             <td>Iterations</td>" >> $1
   echo "             <td>CPU</td>" >> $1
   echo "             <td>Residue</td>" >> $1                           # BICGStab
   echo "             <td>Iterations</td>" >> $1
   echo "             <td>CPU</td>" >> $1
   echo "             <td>Residue</td>" >> $1                           # GMRES Results restarting 8
   echo "             <td>Iterations</td>" >> $1
   echo "             <td>CPU</td>" >> $1
   echo "             <td>Residue</td>" >> $1                           # GMRES Results restarting 16
   echo "             <td>Iterations</td>" >> $1
   echo "             <td>CPU</td>" >> $1
   echo "             <td>Residue</td>" >> $1                           # GMRES Results restarting 32
   echo "             <td>Iterations</td>" >> $1
   echo "             <td>CPU</td>" >> $1
   echo "             <td>Residue</td>" >> $1                           # GMRES Results restarting 64
   echo "             <td>Iterations</td>" >> $1
   echo "             <td>CPU</td>" >> $1
   echo "             <td>Residue</td>" >> $1                           # GMRES Results restarting 128
   echo "             <td>Iterations</td>" >> $1
   echo "             <td>CPU</td>" >> $1
   echo "             <td>Residue</td>" >> $1                           # GMRES Results restarting 256
   echo "             <td>Iterations</td>" >> $1
   echo "             <td>CPU</td>" >> $1
   echo "          </tr>" >> $1
}

close_log()
# $1 log file name
{
      echo "      </table>" >> $1
   echo "   </body>" >> $1
   echo "</html>" >> $1
}

write_log_header $LOG_FILE

for link in $MM_MATRICES;
do
   echo "###############################################################################################"
   matrix=matrices`echo $link | sed 's/ftp:\/\/math.nist.gov\/pub//'`
   unzipped_matrix=`echo $matrix | sed 's/.gz//'`
   if test ! -e $matrix;
   then
      echo "Matrix $matrix is missing !!! Run the script 'get-matrices' first."
      #echo "Matrix $matrix is missing !!! Run the script 'get-matrices' first." >> sparse-matrix-benchmark.log            
   else
      if test ! -e $unzipped_matrix.double.bin.bz2;
      then
         echo "Missing $unzipped_matrix.double.bin.bz2 !!! Run the script 'convert-matrices'."
         #echo "Missing $unzipped_matrix.double.bin.bz2 !!! Run the script 'convert-matrices'."  >> sparse-matrix-benchmark.log
      else
         gunzip -f $matrix
         echo "Checking with the matrix $unzipped_matrix in double precison ..."
         export CUDA_PROFILE_LOG=$unzipped_matrix.double.log 
         benchmark_matrix $unzipped_matrix.double.bin.bz2 $unzipped_matrix $LOG_FILE 1
         gzip $unzipped_matrix
      fi          
   fi
done

for link in $FLORIDA_MM_MATRICES;
do
   matrix=matrices`echo $link | sed 's/http:\/\/www.cise.ufl.edu\/research\/sparse//'`
   if test ! -e $matrix;
   then      
      echo "Matrix $matrix is missing !!! Run the script 'get-matrices' first."
      #echo "Matrix $matrix is missing !!! Run the script 'get-matrices' first." >> sparse-matrix-benchmark.log
   else
     DIRNAME=`dirname $matrix`
     FILENAME=`basename $matrix`
     cd $DIRNAME
     tar zxvf $FILENAME
     cd $IWD
     SUBDIRNAME=`echo $FILENAME | sed 's/.tar.gz//'`
     rm -f $DIRNAME/$SUBDIRNAME/*_b.mtx # these are usualy in array format
     for file in $DIRNAME/$SUBDIRNAME/*.mtx;
     do
         echo "###############################################################################################"
         if test ! -e $file.double.bin.bz2;
         then
            echo "Missing $file.double.bin.bz2 !!! Run the script 'convert-matrices'."
         else   
            echo "Checking with the matrix $file ..."
            benchmark_matrix $file.double.bin.bz2 $file $LOG_FILE 1                        
         fi                           
         rm $file         
      done
   fi
done

close_log matrix-solvers-benchmark.log
