File: cmplcpintervals.sh

package info (click to toggle)
genometools 1.6.1%2Bds-3
  • links: PTS, VCS
  • area: main
  • in suites: bullseye
  • size: 50,412 kB
  • sloc: ansic: 271,241; ruby: 30,339; python: 4,880; sh: 3,193; makefile: 1,194; perl: 219; pascal: 159; haskell: 37; sed: 5
file content (50 lines) | stat: -rwxr-xr-x 1,183 bytes parent folder | download | duplicates (9)
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
#!/bin/sh 

checkerror() {
if test $? -ne 0
then
  echo "failure: ${cmd}"
  exit 1
fi
}

if test $# -eq 1
then
  filenamelist="$1"
else
  filenamelist=`find testdata/ -name '[A-Z]*.fna'`
fi

for filename in ${filenamelist}
do
  echo "${filename}"
  cmd="env -i bin/gt suffixerator -tis -suf -lcp -dna -indexname sfx -db ${filename}"
  ${cmd}
  checkerror
  cmd="env -i RUBYLIB=gtruby LD_LIBRARY_PATH=lib scripts/lcpintervals.rb --itv sfx"
  ${cmd} > itvs.result1
  checkerror
  cmd="env -i bin/gt dev sfxmap -enumlcpitvs -esa sfx"
  ${cmd} > itvs.result2
  checkerror
  cmd="cmp -s itvs.result1 itvs.result2"
  ${cmd}
  checkerror
  cmd="env -i RUBYLIB=gtruby LD_LIBRARY_PATH=lib scripts/lcpintervals.rb --tree sfx"
  ${cmd} > itvtree.result1
  checkerror
  cmd="env -i bin/gt dev sfxmap -enumlcpitvtree -esa sfx"
  ${cmd} > itvtree.result2
  checkerror
  cmd="env -i bin/gt dev sfxmap -enumlcpitvtreeBU -esa sfx"
  ${cmd} > itvtree.result3
  checkerror
  cmd="cmp -s itvtree.result1 itvtree.result2"
  ${cmd}
  checkerror
  cmd="cmp -s itvtree.result2 itvtree.result3"
  ${cmd}
  checkerror
done
rm -f itvs.result1 itvs.result2
rm -f itvtree.result1 itvtree.result2 itvtree.result3