Re: [Seqan-dev] index creation fails


Fixed in [10241]. Thanks for reporting.

David

Am 08.08.2011 um 18:02 schrieb Daniel James:

> Hi
> 
> I emailed a month ago, got no response. I can't seem to create a
> ticket on the trac. Thought I might try again here.
> 
> The following code snippet is crashing. I think it's a bug though
> can't be sure since I'm not an expert. The code might not look very
> minimal but I needed to generate enough data to reproduce it, thus the
> rather large example. The point at which the crash occurs is here:
> 
> indexRequire(index, QGramSA());
> 
> I'm on an up to date svn. I'm running OSX Snow Leopard. Please let me
> know if there's any other information I can provide or otherwise.
> 
> Thank you in advance,
> 
> Daniel
> 
> 
> 
> #include <algorithm>
> #include <seqan.h>
> 
> using namespace seqan;
> 
> // Generates random nucleotides.
> struct MyGenerator : std::unary_function<char, void>
> {
>   std::string syms;
>   MyGenerator (std::string syms = "ACGT") : syms(syms) { srand(time(NULL)); }
>  char operator()(void) { return syms[rand() % syms.size()]; }
> };
> 
> int main(int argc, char** argv)
> {
>  typedef StringSet<DnaString>                                  TMyStringSet;
>  typedef Index<TMyStringSet, IndexQGram<UngappedShape<8> > >   TMyIndex;
> 
>  StringSet<DnaString> myStringSet;
> 
>  for (unsigned i = 0; i < 100; ++i) {
>      DnaString input;
>      resize(input, 60);
>      generate_n(begin(input), 60, MyGenerator());
>      appendValue(myStringSet, input);
>  }
> 
>  std::cout << myStringSet[0] << std::endl;
>  TMyIndex index(myStringSet);
>  std::cout << "requiring QGramSA..." << std::endl;
> 
>  indexRequire(index, QGramSA());
> 
>  return 0;
> }
> 
> _______________________________________________
> seqan-dev mailing list
> seqan-dev@lists.fu-berlin.de
> https://lists.fu-berlin.de/listinfo/seqan-dev