1#! /bin/sh
2
3mkdir -p tests-data; cd tests-data
4. ../compat.sh
5
6sort -k2,2 > ${pref}.md5sum <<EOF
772f1913b3503114c7df7a4dcc68ce867 ${pref}_m40_s16m.histo
872f1913b3503114c7df7a4dcc68ce867 ${pref}_automerge_m40_s1m.histo
972f1913b3503114c7df7a4dcc68ce867 ${pref}_m40_s1m_merged.histo
1072f1913b3503114c7df7a4dcc68ce867 ${pref}_m40_s1m_text.histo
114199aa97e646281b9c36a03564f082ee ${pref}_m9_min.histo
1210b5ca10bcf85183f82837ea10638a8d ${pref}_m9_max.histo
130ff4b7a2f3f67fd26011a41f820bdf36 ${pref}_m9_jaccard
14EOF
15
16# Compare counts without merging, with explicit merging and implicit merging
17FILES="seq1m_0.fa seq1m_1.fa seq1m_0.fa seq1m_2.fa seq1m_2.fa"
18echo $FILES | xargs $JF count -t $nCPUs -o ${pref}_m40_s16m.jf -s 4M -C -m 40
19$JF histo ${pref}_m40_s16m.jf > ${pref}_m40_s16m.histo
20
21ls | grep "^${pref}_m40_s1m[0-9].*" | xargs rm -f
22echo $FILES | xargs $JF count -t $nCPUs -o ${pref}_m40_s1m -s 1M --disk --no-merge -C -m 40
23$JF merge -o ${pref}_m40_s1m_merged.jf ${pref}_m40_s1m[0-9]*
24ls | grep "^${pref}_m40_s1m[0-9].*" | xargs rm -f
25
26echo $FILES | xargs $JF count -t $nCPUs -o ${pref}_automerge_m40_s1m.jf -s 1M --disk -C -m 40
27
28echo $FILES | xargs $JF count -t $nCPUs -o ${pref}_m40_s1m_text.jf -s 1M --text --disk -C -m 40
29
30$JF histo ${pref}_automerge_m40_s1m.jf > ${pref}_automerge_m40_s1m.histo
31$JF histo ${pref}_m40_s1m_merged.jf > ${pref}_m40_s1m_merged.histo
32$JF histo ${pref}_m40_s1m_text.jf > ${pref}_m40_s1m_text.histo
33
34# Compute weighted intersection (min), weighted uniont (max) and Jaccard
35$JF count -t $nCPUs -o ${pref}_m9_2.jf -C -m 9 -s 4M seq1m_2.fa
36$JF count -t $nCPUs -o ${pref}_m9_3.jf -C -m 9 -s 4M seq1m_3.fa
37$JF merge -m -o ${pref}_m9_min.jf ${pref}_m9_2.jf ${pref}_m9_3.jf
38$JF merge -M -o ${pref}_m9_max.jf ${pref}_m9_2.jf ${pref}_m9_3.jf
39$JF merge -j -o ${pref}_m9_jaccard ${pref}_m9_2.jf ${pref}_m9_3.jf
40$JF histo ${pref}_m9_min.jf > ${pref}_m9_min.histo
41$JF histo ${pref}_m9_max.jf > ${pref}_m9_max.histo
42
43
44check ${pref}.md5sum
45