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