# Last edited on 2022-04-14 17:42:13 by stolfi

# Processing the results for the HotFill tests for the paper.

# Creating the Delta x Rfab plots:

  # Directory with the data files:

    datadir=2021-12-11-143000

  # Got a bunch of files from Elis "${datadir}/plot_data_filter_maxband${mu}.csv"

  # Generate the plots of Delta x Rfab for various maxband mu:

    # Datasets eliminated to reduce text size
    # adfoot_0 bkwren_0 hlatch_90 nkhang_90 

    datasets=( \
      dstop2_90 runleg_0 \
      tkwren_90 dstop1_0 unwren_90 \
      grille_0 cthang_90 flange_0 \
      runleg_90 \
    )
    fabmus=( 5 15 20 25 80 )
    plotdir=${datadir}/Rcon-graphs
    mkdir -p ${plotdir}
    for mu in ${fabmus[@]}; do 
      plot_hotfill_performance.sh \
        ${datadir} \
        plot_data_filter_maxband${mu}.csv \
        ${plotdir}/D008000-mb${mu}.pdf \
        ${datasets[@]}
    done

  Files for upload are "${plotdir}/D008000-mb${mu}.png"

# Viewing the solutions:

    datadir=2021-12-03-210000
    mu="20"
    for ds in ${datasets[@]}; do
      show_hotfill_paths.sh ${datadir} ${ds} ${mu} 
    done

# Plotting the Tcpu dependency on max bnd width:

# Obtaining the data:

  # Got from Elis/Rodrigo a zip file with the "out.txt" files

    datadir=2021-12-11-143000
    ( cd ${datadir} && unzip outfiltered.zip )
    
  # Had to rename the files to remove the slice index.
  # Created by hand a file "${datadir}/rename_datasets.sh"
  
    (cd ${datadir} && rename_datasets.sh )

  # Extracting the Tcpu raw data from the "out.txt" files:

    datadir=2021-12-11-143000
    plotdir=${datadir}/Tcpu-graphs
    mkdir -p ${plotdir}
    cpumus=( 10 20 30 40 50 60 70 80 )
    Delta="64.0"
    extract_tcpus.sh "${datadir}" "${plotdir}" "${Delta}" "${cpumus[@]}" ${datasets[@]}
    
  # Plotting the Tcpus as function of mu:
  
    mumax=${cpumus[-1]}
    plot_tcpus.sh "${plotdir}" "${Delta}" ${mumax} ${datasets[@]}

# ======================================================================
# JUNK

  # Get (Tcoolmax, Rfab) data points for SCN, SCA, RP3, SLIC3R.
  # Did not look good because those points did not
  # easily connect to the curves.

    datadir=
    for alg in SLIC3R RP3 SCANLINE ALTSCANLINE ; do
      extract_ref_data.sh ${datadir} ${alg} "${datasets[@]}"
    done

# Estimating the Tcpu dependency on mu by polynominal regression.

  # Creating the input files for regression:
  
    regmus=( 40 50 60 70 80 )
    regsets=( \
      bkwren_0 cthang_90 dstop1_0 flange_0 hlatch_90 \
      nkhang_90 runleg_0 runleg_90 tkwren_90 \
    )
    for ds in ${regsets[@]}; do 
      dsregdir=${regdir}/${ds}_D${Delta}
      rawfile=${dsregdir}/raw.txt
      kukfile=${dsregdir}/cooked.txt
      gawk \
          ' /[0-9]/{
              x=($1+0)/100; y=$2; k++;
              printf "%03d %8.0f %10.8f %10.8f %10.8f %10.8f\n", k, y, 1, x, x*x, x*x*x
            }
           ' \
        < ${rawfile} \
        > ${kukfile}
    done

  # Performing the regression:
  
    for ds in ${regsets[@]}; do 
      dsregdir=${regdir}/${ds}_D${Delta}
      kukfile=${dsregdir}/cooked.txt
      fitfile=${dsregdir}/fitted.txt
      formfile=${dsregdir}/form.txt
      linear_fit \
          -terms 4 \
          -writeFormula ${formfile} \
        < ${kukfile} \
        > ${fitfile}
    done

    nev 'x=10; 0.783*x; 0.0692*x^2; 0.000977*x^3'
    nev 'x=20; 0.783*x; 0.0692*x^2; 0.000977*x^3'
    nev 'x=70; 0.783*x; 0.0692*x^2; 0.000977*x^3'
    nev 'x=80; 0.783*x; 0.0692*x^2; 0.000977*x^3'

    nev 'x=80; 0.850*x; 0.104*x^2; 0.000829*x^3'