1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
|
#!/bin/bash
set -e
pkg="scrm"
export LC_ALL=C.UTF-8
if [ "${AUTOPKGTEST_TMP}" = "" ] ; then
AUTOPKGTEST_TMP=$(mktemp -d /tmp/${pkg}-test.XXXXXX)
trap "rm -rf ${AUTOPKGTEST_TMP}" 0 INT QUIT ABRT PIPE TERM
fi
test_scrm () {
echo -n " scrm $@ "
for i in `seq 1 10`; do
echo -n "."
# Test using scrm self-checks
scrm $@ -seed $i > /dev/null
if [ $? -ne 0 ]; then
echo ""
echo "Executing \"./scrm $@ -seed $i\" failed."
echo "Debug Call: make -mj2 scrm_dbg && ./scrm_dbg $@ -seed $i 2>&1 | less"
exit 1
fi
done
echo " done."
}
echo "Testing Initial Tree"
test_scrm 5 100 -t 5 || exit 1
test_scrm 3 100 -t 5 -L -oSFS -transpose-segsites || exit 1
test_scrm 50 20 -O || exit 1
test_scrm 25 1 -T || exit 1
test_scrm 10 5 -I 1 5 -eI 1.0 5 -L || exit 1
test_scrm 10 5 -I 2 2 2 0.5 -eI 1.0 3 3 -L || exit 1
echo ""
echo "Testing Recombinations"
test_scrm 4 10 -r 5 100 -l -1 || exit 1
test_scrm 6 10 -r 1 100 -t 5 -L -T -transpose-segsites -l -1 || exit 1
test_scrm 8 10 -r 1 100 -t 5 -oSFS -O -l -1 || exit 1
test_scrm 10 5 -r 5 100 -I 2 2 2 0.5 -eI 1.0 3 3 -L -l -1 || exit 1
echo ""
echo "Testing Pruning"
test_scrm 10 20 -r 10 500 -l 10 -t 5 -oSFS || exit 1
test_scrm 3 50 -r 10 500 -l 0 -t 1 -oSFS || exit 1
test_scrm 10 20 -r 10 500 -l 2r -L || exit 1
test_scrm 10 20 -r 1 500 -l -1 -T || exit 1
test_scrm 10 5 -r 5 100 -I 2 2 2 0.5 -eI 1.0 3 3 -L -l 2r || exit 1
echo ""
echo "Testing Migration"
test_scrm 5 5 -r 5 100 -I 2 3 2 1.2 || exit 1
test_scrm 10 2 -r 20 200 -I 5 2 2 2 2 2 0.75 -l 5 || exit 1
test_scrm 10 2 -r 10 100 -I 2 7 3 0.5 -eM 0.3 1.1 -O || exit 1
test_scrm 10 2 -r 10 100 -I 2 7 3 -m 1 2 0.3 -em 0.5 2 1 0.6 -eM 2.0 1 || exit 1
test_scrm 20 2 -I 3 2 2 2 1.0 -eI 1.0 2 2 2 -eI 2.0 2 3 3 || exit 1
echo ""
echo "Testing Size Change"
test_scrm 10 2 -r 1 100 -I 3 3 3 4 0.5 -eN 0.1 0.05 -eN 0.2 0.5 -O || exit 1
test_scrm 10 2 -r 10 100 -I 3 3 3 4 0.5 -eN 0.1 0.05 -eN 0.2 0.5 -l 10 || exit 1
echo ""
echo "Testing Splits & Merges"
test_scrm 5 3 -r 20 200 -I 2 3 2 0.4 -ej 1.1 2 1 -l 25 || exit 1
test_scrm 6 3 -r 20 200 -I 3 2 2 2 -ej 0.2 2 1 -ej 0.25 3 1 -l 25 || exit 1
test_scrm 20 2 -r 5 200 -I 2 10 10 1.5 -es 1.6 2 0.5 -eM 2.0 1 -l 25 || exit 1
test_scrm 20 2 -r 5 200 -I 2 10 10 1.5 -es 0.9 1 0.8 -es 1.6 2 0.5 -eM 2.0 1 -l 50 || exit 1
test_scrm 20 2 -r 5 200 -I 2 10 10 -eps 1.0 2 1 0.0 || exit 1
echo ""
echo "Testing Growth"
test_scrm 5 10 -r 20 200 -G 1.5 -l 25 || exit 1
test_scrm 5 6 -r 2 200 -G -2.5 -eG 1 0.0 -l 25 || exit 1
test_scrm 4 3 -r 2 200 -G -2.5 -eN 1 0.25 -eG 2 0.0 -eN 2.5 0.25 -l 25 || exit 1
echo ""
echo "Testing Variable Rates"
test_scrm 3 2 -r 2 100 -t 5 -st 10 10 -sr 20 5 -st 30 1 -sr 40 0 -st 50 20 -T || exit 1
echo ""
echo "Various Edge Cases"
test_scrm 6 1 -I 2 3 3 0.5 -r 1 100 -es 0 2 0.5 -ej 1 3 1 -t 1 || exit 1
test_scrm 10 1 -es 1.0 1 0.5 -ej 1.0 2 1 || exit 1
test_scrm 4 1 -t 1.0 -I 2 2 0 -ej 1.0 1 2 -eI 3.5 0 2 || exit 1
echo ""
|
