File: bloom_filter.sh

package info (click to toggle)
jellyfish 2.2.10-2
  • links: PTS, VCS
  • area: main
  • in suites: buster
  • size: 3,320 kB
  • sloc: cpp: 35,615; sh: 963; ruby: 578; makefile: 372; python: 165; perl: 36
file content (32 lines) | stat: -rw-r--r-- 904 bytes parent folder | download
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
#! /bin/sh

cd tests
. ./compat.sh

# sort -k2,2 > ${pref}.md5sum <<EOF
# EOF

$JF count --bf-size 10M --bf-fp 0.001 -t $nCPUs -o ${pref}_10m.jf -s 1M -m 40 seq10m.fa
$JF histo ${pref}_10m.jf > ${pref}_10m.histo
if bc < /dev/null 2> /dev/null; then
    COLLISIONS=$(cut -d\  -f2 ${pref}_10m.histo | paste -sd+ - | bc)
    [ $((COLLISIONS > 10000)) = 0 ] || {
        echo >&2 "Too many collisions"
        false
    }
else
    echo > "bc missing: Skip collisions test"
fi

$JF count --bf-size 3M --bf-fp 0.001 -t $nCPUs -o ${pref}_3m.jf -s 1M -m 40 seq1m_0.fa seq1m_1.fa seq1m_0.fa seq1m_2.fa
$JF histo ${pref}_3m.jf > ${pref}_3m.histo

if bc < /dev/null 2> /dev/null; then
    COLLISIONS=$(cut -d\  -f2 ${pref}_3m.histo | paste -sd+ - | bc)
    [ $((COLLISIONS - 1000000 > 20000)) = 0 ] || {
        echo >&2 "Too many collisions"
        false
    }
else
    echo "bc missing: Skip collisions test"
fi