diff --git a/regression/run_fwbw.sh b/regression/run_fwbw.sh new file mode 100644 index 0000000..4f2f610 --- /dev/null +++ b/regression/run_fwbw.sh @@ -0,0 +1,34 @@ +#!/bin/sh -e +fake_pref() { + QDB="$1" + TDB="$2" + RES="$3" + # create link to data file which contains a list of all targets that should be aligned + ln -s "${TDB}.index" "${RES}" + # create new index repeatedly pointing to same entry + INDEX_SIZE="$(echo $(wc -c < "${TDB}.index"))" + awk -v size=$INDEX_SIZE '{ print $1"\t0\t"size; }' "${QDB}.index" > "${RES}.index" + # create dbtype (7) + awk 'BEGIN { printf("%c%c%c%c",7,0,0,0); exit; }' > "${RES}.dbtype" +} + +mkdir -p "${RESULTS}/fwbw" #makes sure index files do not get mixed + +QUERY="${DATADIR}/query.fasta" +QUERYDB="${RESULTS}/fwbw/query_fwbw" +"${MMSEQS}" createdb "${QUERY}" "${QUERYDB}" + +TARGET="${DATADIR}/targetannotation.fasta" +TARGETDB="${RESULTS}/fwbw/targetannotation_fwbw" +"${MMSEQS}" createdb "${TARGET}" "${TARGETDB}" + +"${MMSEQS}" prefilter "$QUERYDB" "${TARGETDB}" "$RESULTS/fwbw/pref" +awk 'BEGIN { printf("%c%c%c%c",5,0,0,0); exit; }' > "$RESULTS/fwbw/pref.dbtype" +"${MMSEQS}" fwbw "$QUERYDB" "$TARGETDB" "$RESULTS/fwbw/pref" "$RESULTS/fwbw/results_fwbw" -e 10000 +"${MMSEQS}" convertalis "$QUERYDB" "$TARGETDB" "$RESULTS/fwbw/results_fwbw" "$RESULTS/fwbw/results_fwbw.m8" +"${EVALUATE}" "$QUERY" "$TARGET" "$RESULTS/fwbw/results_fwbw.m8" "${RESULTS}/fwbw/evaluation_roc5.dat" 4000 1 | tee "${RESULTS}/fwbw/evaluation.log" +ACTUAL=$(grep "^ROC5 AUC:" "${RESULTS}/fwbw/evaluation.log" | cut -d" " -f3) +TARGET="0.207452" +awk -v actual="$ACTUAL" -v target="$TARGET" \ + 'BEGIN { print (actual == target) ? "GOOD" : "BAD"; print "Expected: ", target; print "Actual: ", actual; }' \ + > "${RESULTS}.report" diff --git a/run_regression.sh b/run_regression.sh index eadc524..fbc266b 100755 --- a/run_regression.sh +++ b/run_regression.sh @@ -95,6 +95,7 @@ run_test LINCLUST_UPDATE "run_cluster_update.sh" run_test EASYNUCLNUCLTAX_SEARCH "run_easy_nuclnucltax.sh" run_test EXTRACTORFS "run_extractorfs.sh" run_test RBH "run_rbh.sh" +run_test FWBW "run_fwbw.sh" case "$(uname -s)" in CYGWIN*|MINGW32*|MSYS*|MINGW*) ;;