Skip to content

Instantly share code, notes, and snippets.

@cbergman
Last active January 3, 2016 05:59
Show Gist options
  • Save cbergman/8419306 to your computer and use it in GitHub Desktop.
Save cbergman/8419306 to your computer and use it in GitHub Desktop.
#!/bin/bash
cd ../example/ecoli/output/local/
for l in 100; do
for o in 0; do
for s in 0.02; do
for i in 0; do
for d in 0; do
rm -rf ecoli.l$l.o$o.qA.s$s.i$i.d$d;
mkdir ecoli.l$l.o$o.qA.s$s.i$i.d$d;
cd ecoli.l$l.o$o.qA.s$s.i$i.d$d;
rm -fr ./*
date >> ../timestamp.txt
echo "l:$l" >> ../timestamp.txt
echo "q:$q" >> ../timestamp.txt
echo "s:$s" >> ../timestamp.txt
echo "i:$i" >> ../timestamp.txt
echo "d:$d" >> ../timestamp.txt
echo "o:$o" >> ../timestamp.txt
echo "t:$t" >> ../timestamp.txt
echo "f:$f" >> ../timestamp.txt
echo "b:$b" >> ../timestamp.txt
echo "x:$x" >> ../timestamp.txt
echo "m:$m" >> ../timestamp.txt
echo "p:$p" >> ../timestamp.txt
echo "--" >> ../timestamp.txt
cat ../../../input/ppd/NC_000913.fa.ppd | ../../../../../bin/mapper runall -l $l -q A -s $s -i $i -d $d -o $o -t 20 -f ref.fa -b 70 -x ../../../input/index/NC_000913.fna -p ../../../../../bin 1> map.txt
echo "sorting..."
cat map.txt | sort > mapsort.txt
echo "complete"
echo "====================================="
echo "analyzing..."
cat mapsort.txt | ../../../../../bin/reducer analyzer -l $l -t 20 -o $o 1> mapred.txt
echo "====================================="
echo "reducer is done"
echo "====================================="
done;
done;
done;
done;
done;
cd ../../../../../script/
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment