File: shulength.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 (29 lines) | stat: -rwxr-xr-x 1,082 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
#!/bin/sh
if test $# -ne 2
then
  echo "Usage: <file1> <file2>"
  exit 1
fi

WOTD="${TEACHDIR}/repertoire/exercises/effiziente_algorithmen_auf_sequenzen/suffix_trees/suffixtree_code/wotd.rb"
file1=$1
file2=$2

echo "# pairwise comparisons"
env -i ${GTDIR}/bin/gt suffixerator -db ${file1} -indexname idx1 -dna -suf -tis
env -i ${GTDIR}/bin/gt shulengthdist -ii idx1 -q ${file2} > result.pairwise
env -i ${GTDIR}/bin/gt suffixerator -db ${file2} -indexname idx2 -dna -suf -tis
env -i ${GTDIR}/bin/gt shulengthdist -ii idx2 -q ${file1} >> result.pairwise
env -i ${GTDIR}/bin/gt suffixerator -db ${file1} ${file2} -dna -suf -tis -lcp -indexname both 
echo "# multi comparisons"
env -i ${GTDIR}/bin/gt shulengthdist -ii both > result.multi
diff result.pairwise result.multi
if test $? -ne 0
then
  echo "failure: cmp -s result.pairwise result.multi"
  # echo "`tail -1 ${file1}`x`tail -1 ${file2}`y" | ${WOTD} | dot -Tpdf > dot.pdf
  # mkvtree -db ${file1} ${file2} -dna -suf -lcp -tis -indexname mkv-tmp
  # vstree2tex -suf -lcp -tis -s mkv-tmp > tmp.tex
  # latex tmp.tex
  exit 1
fi