New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

aai-matrix is really slow #25

sbridel opened this Issue Apr 4, 2017 · 0 comments


None yet
1 participant

sbridel commented Apr 4, 2017

I can't re-open my issue about multithreading.
So I re post my message here

Unfortunately i can't read MiGA documentation. I see only the titles on the doc.

I was wondering how blast is called within aai.rb.
Then, I wonder how to add mpiBLAST ( in aai.rb to do the job using BLASt on several CPUs
Or, how to add parallel / xargs to to parallelize the loop. I don't know if it's better to launch the main loop in several jobs or launch the nested loop as several jobs. Then, aai.rb will be launched as several jobs, each jobs calling BLAST. This solution doesn't allow faster BLAST pairwise comparisons, but make aai.rb between some pairs of genomes in parallel at the same time so it may reduce the computing time.
But, I haven't done something like this before and I don't really know how to do that

for i in "${SEQS[@]}" ; do
  for j in "${SEQS[@]}" ; do
    echo -n " o $i vs $j: "
    AAI=$(aai.rb -1 "$i" -2 "$j" -S "$OUT.db" -t "$THR" \
      --no-save-rbm --auto --quiet)
    echo ${AAI:-Below detection}
    [[ "$i" == "$j" ]] && break

I want to add that I have launched AAI comparaison using aai-matrix. I have 35 genomes.
from 11h30 am to 5p.m, it have only done 2 species againts the 35 species (so 70 comparisons).
I wonder how much time will take the full analysis.

I'm pretty sure that parallel can perform well by parallelizing the loops. I tried to make it works but I have never used parallel before so I hoped for some help now.

This will be a great features for your tool.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment