exp-drd/scripts/run-splash2

#!/bin/bash

########################
# Function definitions #
########################

function log2 {
  local i

  for ((i=0;i<64;i++))
  do
    if [ $((2**i)) = $1 ]; then
      echo $i
      return 0
    fi
  done
  echo ""
  return 1
}

function get_cache_size {
  local s
  s=$(</sys/devices/system/cpu/cpu0/cache/index2/size)
  if [ "${s%M}" != "$s" ]; then
    echo $((${s%M}*1024*1024))
  elif [ "${s%K}" != "$s" ]; then
    echo $((${s%K}*1024))
  else
    echo $s
  fi
}

# Read a stream of numbers from stdin (one per line), and print the average
# and standard deviation.
function avgstddev {
  awk '{n++;sum+=$1;sumsq+=$1*$1}END{d=sumsq/n-sum*sum/n/n;print sum/n,(d>0?sqrt(d):0)}'
}

function run_test {
  local tmp avg1=1 stddev1=1 avg2=1 stddev2=1

  tmp="/tmp/test-timing.$$"
  echo "$@"
  for ((i=0;i<3;i++))
  do
    cat "${test_input:-/dev/null}" | \
      /usr/bin/time --format="%e" "$@" 2>&1 | \
      tail -n 1
  done | avgstddev > "$tmp"
  read avg1 stddev1 < "$tmp"
  echo "Average time: ${avg1} +/- ${stddev1} seconds"

  for p in 1 2 4
  do
    echo "$VG --tool=exp-drd $@ -p$p"
    for ((i=0;i<3;i++))
    do
      cat "${test_input:-/dev/null}" | \
        /usr/bin/time --format="%e" $VG --tool=exp-drd "$@" -p$p 2>&1 | \
        tail -n 1
    done | avgstddev > "$tmp"
    read avg2 stddev2 < "$tmp"
    echo "Average time: ${avg2} +/- ${stddev2} seconds"
    awk "END{print "'"'"Ratio ="'"'", ${avg2}/${avg1}, "'"'"+/-"'"'", ${avg2}/${avg1}*(${stddev1}/${avg1}+${stddev2}/${avg2})}" </dev/null
  done

  echo ''

  rm -f "$tmp"
}


# Script body

DRD_SCRIPTS_DIR="$(dirname $0)"
if [ "${DRD_SCRIPTS_DIR}" = "." ]; then
  DRD_SCRIPTS_DIR="$PWD"
fi

SPLASH2="${DRD_SCRIPTS_DIR}/../splash2"
if [ ! -e "${SPLASH2}" ]; then
  echo "Error: splash2 directory not found (${SPLASH2})."
  exit 1
fi

if [ "$VG" = "" ]; then
  VG="${DRD_SCRIPTS_DIR}/../../vg-in-place"
fi

if [ ! -e "$VG" ]; then
  echo "Could not find $VG."
  exit 1
fi

# Results:                      (-p1) (-p2) (-p3) (-p4) ITC (-p4) ITC (-p4)
#                                                        original w/ filter
# .........................................................................
# Cholesky                        46    59    75    92      239      82
# FFT                             15    19   N/A    47       90      41
# LU, contiguous blocks           40    48    53    55      428     128
# LU, non-contiguous blocks       37    47    55    58      428     128
# Ocean, contiguous partitions    20    26   N/A    32       90      28
# Ocean, non-continguous partns   19    24   N/A    34       90      28
# Radiosity                       99    99    99    99      485     163
# Radix                           11    15     ?    17      222      56
# Raytrace                        75    75     ?    75      172      53
# Water-n2                       290   290     ?   290      189      39
# Water-sp                       288   288     ?   287      183      34

cache_size=$(get_cache_size)
log2_cache_size=$(log2 ${cache_size})

# Cholesky
(
  cd ${SPLASH2}/codes/kernels/cholesky/inputs
  for f in *Z
  do
    gzip -cd <$f >${f%.Z}
  done
  run_test ../CHOLESKY -C${cache_size} -n1024 tk29.O
)

# FFT
run_test ${SPLASH2}/codes/kernels/fft/FFT -t -l${log2_cache_size} -m20

# LU, contiguous blocks.
run_test ${SPLASH2}/codes/kernels/lu/contiguous_blocks/LU -n1024

# LU, non-contiguous blocks.
run_test ${SPLASH2}/codes/kernels/lu/non_contiguous_blocks/LU -n1024

# Ocean
run_test ${SPLASH2}/codes/apps/ocean/contiguous_partitions/OCEAN -n2050
run_test ${SPLASH2}/codes/apps/ocean/non_contiguous_partitions/OCEAN -n258

# Radiosity.
run_test ${SPLASH2}/codes/apps/radiosity/RADIOSITY -batch -room

# Radix
run_test ${SPLASH2}/codes/kernels/radix/RADIX -n$((2**24))

# Raytrace
(
  cd ${SPLASH2}/codes/apps/raytrace/inputs
  rm -f *.env *.geo *.rl
  for f in *Z
  do
    gzip -cd <$f >${f%.Z}
  done
  run_test ../RAYTRACE balls4.env
)
run_test ${SPLASH2}/codes/apps/radiosity/RADIOSITY -batch -room

# Water-n2
rm -f water-n2-input
echo "1.5e-16 32768 3 6 -1 3000 3 0 1 6.212752" > water-n2-input
test_input=water-n2-input \
  run_test ${SPLASH2}/codes/apps/water-nsquared/WATER-NSQUARED

# Water-sp
rm -f water-sp-input
echo "1.5e-16 32768 3 6 -1 3000 3 0 1 6.212752" > water-sp-input
test_input=water-sp-input \
  run_test ${SPLASH2}/codes/apps/water-nsquared/WATER-NSQUARED



# Local variables:
# compile-command: "./run-splash2"
# End: