sortalgo.cc

Go to the documentation of this file.

#include "sortalgo.h"
#include <qinputdialog.h>
#include <set>
#include <vector>

void SortAlgo::start()
{

  if ( getMAl()->get_num_seq() == 0 ) {
    return;
  }
  bool ok;
  int res = QInputDialog::getInteger(
                                     "Compact reads", "Enter start row for selected reads:", 0, 1, numeric_limits<int>::max(), 1,
                                     &ok, 0 );
  
  size_t row(1);
  if ( ok ) {
    // user entered something and pressed OK
    row = res;
  } else {
    // user pressed Cancel
    return;
  }
  
  multiset<sort_struct, compBeg> set_begin;

  for( size_t i = 0; i < getMAl()->get_num_seq(); i++ ) {
    sort_struct s;
    s.begin = getMAl()->get_seq_begin_global(i);
    s.end = getMAl()->get_seq_end_global(i);
    s.ID = i;
    
    set_begin.insert(s);
    
  }
  
    multiset<sort_struct, compBeg>::iterator it = set_begin.begin();
    sort_struct tmp(*it);
    tmp.begin = tmp.end;
    vector<sort_struct> row_vec;
    row_vec.push_back(tmp);
    getMAl()->set_seq_row( (*it).ID, row );
    set_begin.erase( set_begin.begin() );
    
    while ( !set_begin.empty() ) {
      it = set_begin.lower_bound( row_vec[ row_vec.size() - 1 ] );
      
      if ( it == set_begin.end() ) {
        //New row

        row_vec.clear();
        ++row;
        it = set_begin.begin();
      }

      tmp = *it;
      tmp.begin = tmp.end;
      row_vec.push_back(tmp);
      getMAl()->set_seq_row( (*it).ID, row );
      set_begin.erase(it);
    } 
    
}

void SortAlgo::old_start()
{

  if ( getMAl()->get_num_seq() == 0 ) {
    return;
  }
  bool ok;
  int res = QInputDialog::getInteger(
                                     "Compact reads", "Enter start row for selected reads:", 0, 1, numeric_limits<int>::max(), 1,
                                     &ok, 0 );
  
  size_t row(1);
  if ( ok ) {
    // user entered something and pressed OK
    row = res;
  } else {
    // user pressed Cancel
    return;
  }
  
  multiset<sort_struct, compBeg> set_begin;

  for( size_t i = 0; i < getMAl()->get_num_seq(); i++ ) {
    sort_struct s;
    s.begin = getMAl()->get_seq_begin_global(i);
    s.end = getMAl()->get_seq_end_global(i);
    s.ID = i;
    
    set_begin.insert(s);
    
  }
  


  multiset<sort_struct, compBeg>::iterator it = set_begin.begin();
  multiset<sort_struct, compEnd> set_end;
  sort_struct tmp(*it);
  tmp.row = row;
  set_end.insert( tmp );
  getMAl()->set_seq_row( (*it).ID, row );
  ++it;
  ++row;

  int debugnum(0);
  
  for( ; it != set_begin.end(); ++it ) {
    
    if ( (*it).begin >= (*(set_end.begin())).end ) {
      getMAl()->set_seq_row( (*it).ID, (*(set_end.begin())).row );
      tmp = *(set_end.begin());
      set_end.erase( set_end.begin() );
      tmp.end = (*it).end;
      set_end.insert(tmp);
      debugnum++;
//       if (debugnum > 4 ) return;
    }
    else {
      getMAl()->set_seq_row( (*it).ID, row );
      tmp = *it;
      tmp.row = row;
      set_end.insert(tmp);
      row++;
    }
  }
}

---