[khmer] keep kmers stranded?
C. Titus Brown
ctb at msu.edu
Wed Jul 31 05:16:58 PDT 2013
On Wed, Jul 31, 2013 at 11:54:22AM +0200, Karl Nordstr?m wrote:
> Hi again,
>
> I'm advancing by piecing together snippets from all the available scripts.
> One question though.
>
> Does K have an upper limit? Loading a counting hash with K=31 and
> everything is fine, but with K=41 the following error appears:
>
> python: ktable.cc:21: khmer::HashIntoType khmer::_hash(const char*,
> khmer::WordLength, khmer::HashIntoType&, khmer::HashIntoType&): Assertion
> `k <= sizeof(HashIntoType)*4' failed.
> Aborted
>
> Or is this due to only using one strand? The code looked like a limitation.
Hi Karl, at the moment you can only go up to k=32. Fixing this is rapidly
rising to the top of our stack -- it's not hard, just something we need to
do carefully.
cheers,
--titus
> On Mon, Jul 29, 2013 at 6:25 PM, C. Titus Brown <ctb at msu.edu> wrote:
>
> > On Mon, Jul 29, 2013 at 05:35:13PM +0200, Karl Nordstr?m wrote:
> > > thank you very much. My intent was to write a simple filter with the help
> > > of the data structures, not anything fancy. Learn how to make use of this
> > > resource.
> > >
> > > I will try and see what happens.
> >
> > Cool -- let us know how it works out ;)
> >
> > cheers,
> > --titus
> >
--
C. Titus Brown, ctb at msu.edu
More information about the khmer
mailing list