0
|
1 #!/usr/bin/env bash
|
|
2 DIR=$1
|
|
3 DIR2=$2
|
|
4 EXT=$3
|
|
5 DIM=$4
|
|
6 KMER=$5
|
|
7
|
|
8
|
|
9
|
|
10 if [ $# != 5 ]; then
|
|
11 echo "***ERROR*** Use: $0 genomesDirectory1 genomesDirectory2 extension dim kmer"
|
|
12 exit -1
|
|
13 fi
|
|
14
|
|
15 indexnameA=$(basename "$DIR")
|
|
16 indexnameB=$(basename "$DIR2")
|
|
17
|
|
18 BINDIR="$( cd "$( dirname "${BASH_SOURCE[0]}" )" && pwd )"
|
|
19
|
|
20
|
|
21 array=()
|
|
22 x=0
|
|
23 array2=()
|
|
24
|
|
25 for elem in $(ls -d $DIR/*.$EXT | awk -F "/" '{print $NF}' | awk -F ".$EXT" '{print $1}')
|
|
26 do
|
|
27 array[$x]=$elem
|
|
28 x=`expr $x + 1`
|
|
29 #echo "X: $elem"
|
|
30 done
|
|
31
|
|
32 x=0
|
|
33
|
|
34 for elem in $(ls -d $DIR2/*.$EXT | awk -F "/" '{print $NF}' | awk -F ".$EXT" '{print $1}')
|
|
35 do
|
|
36 array2[$x]=$elem
|
|
37 x=`expr $x + 1`
|
|
38 #echo "X: $elem"
|
|
39 done
|
|
40
|
|
41 for ((i=0 ; i < ${#array[@]} ; i++))
|
|
42 do
|
|
43 for ((j=0 ; j < ${#array2[@]} ; j++))
|
|
44 do
|
|
45 seqX=${array[$i]}
|
|
46 seqY=${array2[$j]}
|
|
47 echo "----------${seqX}-${seqY}-----------"
|
|
48
|
|
49
|
|
50 #echo "$BINDIR/run_and_plot_chromeister.sh $DIR/${seqX}.$EXT $DIR/${seqY}.$EXT 30 10000"
|
|
51 if [[ ! -f ${seqX}.$EXT-${seqY}.$EXT.mat ]]; then
|
|
52 $BINDIR/run_and_plot_chromeister.sh $DIR/${seqX}.$EXT $DIR2/${seqY}.$EXT $KMER $DIM
|
|
53 fi
|
|
54
|
|
55 done
|
|
56 done
|
|
57
|
|
58
|
|
59 # generate index
|
|
60 if [[ ! -f index.csv.temp ]] && [ ! -f index-$indexnameA-$indexnameB.csv ]; then
|
|
61
|
|
62 echo "Launching... $BINDIR/index_chromeister_solo.sh . index-$indexnameA-$indexnameB.csv $DIR $DIR2"
|
|
63 $BINDIR/index_chromeister_solo.sh . index-$indexnameA-$indexnameB.csv $DIR $DIR2
|
|
64 fi
|
|
65
|
|
66
|
|
67
|