Fix bug #18739.
[pspp] / tests / stats / moments.sh
1 #! /bin/sh
2
3 # Tests calculation of moments.
4
5 TEMPDIR=/tmp/pspp-tst-$$
6
7 # ensure that top_builddir  are absolute
8 if [ -z "$top_builddir" ] ; then top_builddir=. ; fi
9 if [ -z "$top_srcdir" ] ; then top_srcdir=. ; fi
10 top_builddir=`cd $top_builddir; pwd`
11 PSPP=$top_builddir/src/ui/terminal/pspp
12
13 # ensure that top_srcdir is absolute
14 top_srcdir=`cd $top_srcdir; pwd`
15
16 STAT_CONFIG_PATH=$top_srcdir/config
17 export STAT_CONFIG_PATH
18
19
20 cleanup()
21 {
22      cd /
23      rm -rf $TEMPDIR
24      :
25 }
26
27
28 fail()
29 {
30     echo $activity
31     echo FAILED
32     cleanup;
33     exit 1;
34 }
35
36
37 no_result()
38 {
39     echo $activity
40     echo NO RESULT;
41     cleanup;
42     exit 2;
43 }
44
45 pass()
46 {
47     cleanup;
48     exit 0;
49 }
50
51 mkdir -p $TEMPDIR
52
53 cd $TEMPDIR
54 activity="create one-pass moments list"
55 sed -ne 's/#.*//;/^[    ]*$/!p' > $TEMPDIR/moments-list-1p <<'EOF'
56 # Both the one-pass and two-pass algorithms should be 
57 # able to cope properly with these.
58 1 2 3 4 => W=4.000 M1=2.500 M2=1.667 M3=0.000 M4=-1.200
59 1*5 2*5 3*5 4*5 => W=20.000 M1=2.500 M2=1.316 M3=0.000 M4=-1.401
60 1*1 2*2 3*3 4*4 => W=10.000 M1=3.000 M2=1.111 M3=-0.712 M4=-0.450
61 1*0 => W=0.000 M1=sysmis M2=sysmis M3=sysmis M4=sysmis
62 1*1 => W=1.000 M1=1.000 M2=sysmis M3=sysmis M4=sysmis
63 1*2 => W=2.000 M1=1.000 M2=0.000 M3=sysmis M4=sysmis
64 1*3 => W=3.000 M1=1.000 M2=0.000 M3=sysmis M4=sysmis
65 1*2 3 => W=3.000 M1=1.667 M2=1.333 M3=1.732 M4=sysmis
66 1 1.00000001 => W=2.000 M1=1.000 M2=0.000 M3=sysmis M4=sysmis
67 1000001 1000002 1000003 1000004 => W=4.000 M1=1000002.500 M2=1.667 M3=0.000 M4=-1.200
68 EOF
69 if [ $? -ne 0 ] ; then no_result ; fi
70
71 cp $TEMPDIR/moments-list-1p $TEMPDIR/moments-list-2p
72 sed -ne 's/#.*//;/^[    ]*$/!p' >> $TEMPDIR/moments-list-2p <<'EOF'
73 # We used to have an example for which only the two-pass algorithm
74 # produced reasonable results, but the provisional means algorithm
75 # does better, so there aren't any extra tests here.
76 EOF
77
78 activity="create two-pass input file"
79 sed < $TEMPDIR/moments-list-2p >> $TEMPDIR/moments-2p.stat \
80         -e 's#^\(.*\) => \(.*\)$#DEBUG MOMENTS/\1.#'
81 if [ $? -ne 0 ] ; then no_result ; fi
82
83 activity="run two-pass program"
84 $SUPERVISOR $PSPP --testing-mode -o raw-ascii \
85          $TEMPDIR/moments-2p.stat >$TEMPDIR/moments-2p.err 2> $TEMPDIR/moments-2p.out
86
87 activity="compare two-pass output"
88 perl -pi -e 's/^\s*$//g' $TEMPDIR/moments-list-2p $TEMPDIR/moments-2p.out
89 diff -b $TEMPDIR/moments-list-2p $TEMPDIR/moments-2p.out
90 if [ $? -ne 0 ] ; then fail ; fi
91
92 activity="create input file"
93 sed < $TEMPDIR/moments-list-1p >> $TEMPDIR/moments-1p.stat \
94         -e 's#^\(.*\) => \(.*\)$#DEBUG MOMENTS ONEPASS/\1.#'
95 if [ $? -ne 0 ] ; then no_result ; fi
96
97 activity="run one-pass program"
98 $SUPERVISOR $PSPP --testing-mode -o raw-ascii \
99          $TEMPDIR/moments-1p.stat >$TEMPDIR/moments-1p.err 2> $TEMPDIR/moments-1p.out
100
101 activity="compare one-pass output"
102 perl -pi -e 's/^\s*$//g' $TEMPDIR/moments-list-1p $TEMPDIR/moments-1p.out
103 diff -b $TEMPDIR/moments-list-1p $TEMPDIR/moments-1p.out
104 if [ $? -ne 0 ] ; then fail ; fi
105
106 pass