NVIDIA-Genomics-Research / GenomeWorks

SDK for GPU accelerated genome assembly and analysis
https://clara-parabricks.github.io/GenomeWorks/
Apache License 2.0
286 stars 76 forks source link

[cudamapper] IndexGPU correctly skips too short reads and does not copy basepairs twice #550

Closed mimaric closed 4 years ago

mimaric commented 4 years ago

IndexGPU was segfaulting when passed reads that are shorter than one full window (window_size_ + kmer_size_ - 1). This was not exposed because parser already filters out such reads, but this adds an additional level of robustness.

IndexGPU was making an unnecessary copy of reads in fasta_reads which had a measurable performance cost.