8697e3710c1701f41841093afd7b5fa56ea51b3a
[ctsim.git] / scripts / test-single-openmp.sh
1 #/bin/sh
2
3 TEST=medium
4 if test $# -ge 1; then
5     TEST=$1
6 fi
7
8 if test "x$TEST" = "xshort"; then
9     PR=165;  PV=144;  IX=115;  IY=$IX
10 elif test "x$TEST" = "xmedium"; then
11     PR=901;  PV=600;  IX=500;  IY=$IX
12 elif test "x$TEST" = "xlong"; then
13     PR=1651;  PV=1440;  IX=1151;  IY=$IX
14 else
15     echo "usage: $0 [test-size]"
16     echo "test-size -- either short, medium, or long"
17     exit
18 fi
19 TIME=
20
21 DIR=/tmp
22 BIN1=$DIR/omp/ctsimtext
23 BIN2=$DIR/sng/ctsimtext
24
25 for PHAN in herman; do
26     fbase=${DIR}/${PHAN}
27     rm -f ${fbase}-*
28     
29     echo -n "$BIN1 Prj: "; $TIME $BIN1 phm2pj $fbase-omp.pj $PR $PV --phantom $PHAN --verbose | tail -1  | sed 's/Run time: (.+) /$1/'
30     echo -n "$BIN2 Prj: "; $TIME $BIN2 phm2pj $fbase-sng.pj $PR $PV --phantom $PHAN --verbose | tail -1 | sed 's/Run time: (.+) /$1/'
31
32     for BP in idiff diff; do
33         for INT in nearest linear cubic; do
34             recbase=${fbase}-${BP}-${INT}
35             echo -n "$BIN1 Rec $BP $INT: "; $TIME $BIN1 pjrec $fbase-omp.pj $recbase-omp.if $IX $IY --verbose --interp $INT --backproj $BP | tail -1
36             echo -n "$BIN2 Rec $BP $INT: "; $TIME $BIN2 pjrec $fbase-sng.pj $recbase-sng.if $IX $IY --verbose --interp $INT --backproj $BP| tail -1
37             echo -n "--> Difference Rec $BP $INT: "; $BIN1 if2 $recbase-omp.if $recbase-sng.if --comp
38             echo
39         done
40     done
41
42     phmbase=${fbase}-phm
43     echo -n "$BIN1 Phm: "; $TIME $BIN1 phm2if ${phmbase}-omp.if $IX $IX --phantom $PHAN --nsample 5 --verbose | tail -1
44     echo -n "$BIN2 Phm: "; $TIME $BIN2 phm2if ${phmbase}-sng.if $IX $IX --phantom $PHAN --nsample 5 --verbose | tail -1
45     echo -n "--> Image diff Phantom: "; $BIN1 if2 ${phmbase}-omp.if ${phmbase}-sng.if --comp
46     echo
47 done