[Seqan-dev] usage of setStepSize

[Johannes Merkle <merkle@in.tum.de>]

Hi,

I'm experiencing some odd behaviour when using the q-gram-index with a 
different stepsize.
The attached code (which is very similar to your tutorial, only with the 
added setStepSize command and the OpenAddressing spec) returns wrong 
results (0 and 7). Without the OpenAddressing spec it doesn't work at 
all (see attached error message).

Am i using the setStepSize command incorrect or is this a bug?

Thanks,
Johannes



Code:
typedef Index<DnaString, IndexQGram< UngappedShape<3>, OpenAddressing > 
> TIndex;
    TIndex index("CATGATTACATA");
    setStepSize(index,2);
    hash(indexShape(index), "CAT");
    for (unsigned i = 0; i < length(getOccurrences(index, 
indexShape(index))); ++i)
        std::cout << getOccurrences(index, indexShape(index))[i] << 
std::endl;

Error Message:
...SeqAn\core\include\seqan/sequence/string_base.h:227 Assertion failed 
: static_cast<TStringPos>(pos) < static_cast<TStringPos>(length(me)) 
was: 254 >= 65 (Trying to access an element behind the last one!)