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
function test_all_impls {
# program states taxa sites rates reps rsrc rescaling precision sse ctips rseed rfreq root derivs lnl_exp d1_exp d2_exp lscalers ecount ecomplex ievect smatrix
echo -n " testing resource=$R precision=SINGLE " 1>&2;
./parse_test.sh "synthetictest" "${1}" "${2}" "${3}" "${4}" "${5}" "$R" "${10}" "single" "no" "${6}" "${7}" "${11}" "${8}" "${9}" "${17}" "${18}" "${19}" "${12}" "${13}" "${14}" "${15}" "${16}" >> test_results.csv
echo -n " testing resource=$R precision=DOUBLE " 1>&2;
./parse_test.sh "synthetictest" "${1}" "${2}" "${3}" "${4}" "${5}" "$R" "${10}" "double" "no" "${6}" "${7}" "${11}" "${8}" "${9}" "${17}" "${18}" "${19}" "${12}" "${13}" "${14}" "${15}" "${16}" >> test_results.csv
}
if [ -z "${1}" ];
then
R=0
else
R=${1}
fi
if [ ! -f test_results.csv ]
then
echo "program,states,taxa,sites,rates,reps,rsrc,rescaling,precision,sse,ctips,rseed,rfreq,root,derivs,lscalers,ecount,ecomplex,ievect,smatrix,rsrc_name,impl_name,lnl,lnl_diff,d1,d1_diff,d2,d2_diff,best_run,time_real,time_user,time_sys,cpu,gcc_version,revision,date" >> test_results.csv
fi
set -v
# states taxa sites rates reps ctips rseed root derivs rescale rfreq lscalers ecount ecomplex ievect smatrix lnl_exp d1_exp d2_exp
test_all_impls "4" "14" "1240" "4" "2" "7" "1" "yes" "no" "manual" "2" "no" "1" "no" "no" "no" "-3528.89396" "0" "0"
test_all_impls "4" "14" "1240" "4" "2" "7" "1" "yes" "no" "manual" "1" "yes" "1" "no" "no" "no" "-3528.89396" "0" "0"
test_all_impls "4" "14" "1240" "4" "2" "7" "1" "yes" "no" "none" "2" "no" "1" "no" "no" "no" "-3528.89396" "0" "0"
test_all_impls "4" "17" "695" "4" "1" "7" "1" "yes" "no" "manual" "1" "yes" "1" "no" "no" "yes" "-294.51375" "0" "0"
test_all_impls "4" "7" "739" "5" "2" "7" "1" "yes" "no" "manual" "2" "no" "3" "no" "yes" "no" "-4415.64392" "0" "0"
test_all_impls "4" "7" "739" "5" "2" "7" "1" "yes" "no" "none" "2" "no" "3" "no" "yes" "no" "-4415.64392" "0" "0"
test_all_impls "4" "10" "587" "4" "2" "0" "1" "no" "no" "manual" "2" "no" "1" "no" "no" "no" "-4916.38623" "0" "0"
test_all_impls "4" "19" "1853" "4" "2" "0" "1" "yes" "no" "manual" "2" "no" "1" "yes" "no" "no" "-30002.00553" "0" "0"
test_all_impls "4" "9" "900" "2" "2" "9" "1" "no" "no" "manual" "2" "no" "1" "no" "no" "no" "-2726.08630" "0" "0"
test_all_impls "4" "5" "625" "2" "2" "3" "1" "no" "yes" "none" "2" "no" "1" "no" "no" "no" "-3661.78867" "-25.00759" "7.92562"
test_all_impls "4" "8" "1001" "4" "2" "8" "1" "no" "yes" "manual" "2" "no" "1" "no" "no" "no" "-7665.31520" "-111.65185" "91.72033"
test_all_impls "4" "21" "352" "3" "2" "0" "1" "no" "yes" "manual" "1" "no" "1" "no" "no" "no" "-788.68141" "-19.57312" "0.67422"
test_all_impls "5" "8" "456" "2" "2" "8" "1" "no" "no" "manual" "2" "no" "4" "no" "yes" "no" "-3523.14699" "0" "0"
test_all_impls "5" "8" "456" "2" "2" "8" "1" "no" "no" "none" "2" "no" "4" "no" "yes" "no" "-3523.14699" "0" "0"
test_all_impls "8" "13" "637" "2" "2" "13" "1" "yes" "no" "manual" "2" "no" "1" "no" "no" "no" "-14028.08481" "0" "0"
test_all_impls "8" "13" "637" "2" "2" "13" "1" "yes" "no" "none" "2" "no" "1" "no" "no" "no" "-14028.08481" "0" "0"
test_all_impls "11" "15" "854" "1" "2" "1" "1" "yes" "no" "manual" "2" "no" "2" "no" "no" "no" "-18672.38823" "0" "0"
test_all_impls "11" "15" "854" "1" "2" "1" "1" "yes" "no" "none" "2" "no" "2" "no" "no" "no" "-18672.38823" "0" "0"
test_all_impls "13" "5" "425" "1" "1" "5" "1" "no" "yes" "none" "2" "no" "1" "no" "no" "yes" "680.04980" "211.95093" "-15.52490"
test_all_impls "20" "11" "315" "4" "2" "8" "1" "no" "no" "manual" "2" "no" "2" "no" "no" "no" "-5541.25325" "0" "0"
test_all_impls "20" "11" "315" "4" "2" "8" "1" "no" "yes" "manual" "2" "no" "1" "no" "no" "no" "-5873.09996" "-58.63180" "90.59944"
test_all_impls "61" "6" "664" "1" "2" "6" "1" "no" "yes" "manual" "2" "no" "1" "no" "yes" "no" "-9872.86554" "-1.64439" "13.20485"
test_all_impls "61" "6" "664" "1" "2" "6" "1" "no" "yes" "none" "2" "no" "1" "no" "yes" "no" "-9872.86554" "-1.64439" "13.20485"
test_all_impls "64" "12" "399" "3" "2" "12" "1" "no" "yes" "manual" "2" "no" "1" "no" "no" "no" "-11662.60981" "-91.70365" "32.13057"
set +v
######################
# TODO: investigate issue with setmatrix + partial-tips
#test_all_impls "13" "5" "425" "1" "1" "1" "0" "no" "yes" "none" "2" "no" "1" "no" "no" "yes" "-16714.00237" "0" "0"
# TODO: investigate issues with eigencomplex + unrooted and eigencomplex + ratecount != 4
#test_all_impls "4" "10" "587" "4" "2" "0" "0" "no" "no" "manual" "2" "no" "1" "yes" "no" "no" "-9427.94740" "0" "0"
# program states taxa sites rates reps rsrc rescaling precision sse ctips rseed rfreq root derivs lnl_exp d1_exp d2_exp
#./parse_test.sh "fourtaxon" "4" "4" "1314" "4" "500" "0" "no" "single" "no" "0" "0" "1" "no" "no" "-3970.10422" "0" "0"
|
