statgen / demuxlet

Genetic multiplexing of barcoded single cell RNA-seq
Apache License 2.0
117 stars 25 forks source link

Core dumped after/during bcf_filtered_reader #60

Open cramirezs opened 4 years ago

cramirezs commented 4 years ago

Dear Developers,

Thank you for your software, it has been very helpful.

I've been having an issue related to a core dumped while using the SAMFilteredReader function. After backtracing the core file I got the following messages:

(gdb) bt
#0  BCFFilteredReader::parse_genotypes (this=this@entry=0x7ffebc17fd70, hdr=hdr@entry=0x248ab00, v=v@entry=0xdb6ef80) at bcf_filtered_reader.cpp:237
#1  0x0000000000425b98 in BCFFilteredReader::passed_vfilter (this=this@entry=0x7ffebc17fd70, hdr=0x248ab00, hdr@entry=0x0, v=0xdb6ef80, v@entry=0x0)
    at bcf_filtered_reader.cpp:549
#2  0x00000000004270cc in BCFFilteredReader::read (this=this@entry=0x7ffebc17fd70) at bcf_filtered_reader.cpp:753
#3  0x0000000000407d2d in main (argc=<optimized out>, argv=<optimized out>) at cmd_cram_demuxlet.cpp:210
(gdb) list +
232     for(int32_t j=0; j < 2; ++j) {
233       tmp_gt = gts[sm_icols[i]*2+j];
234       if ( ( j < ploidies[i] ) && ( bcf_gt_allele(tmp_gt) >= 0 ) ) {
235     notice("%d %d %x %x", i, j, tmp_gt, bcf_gt_allele(tmp_gt));
236     ++an;
237     ++acs[bcf_gt_allele(tmp_gt)];
238       }
239     }
240   }
241   return true;

I have tried increasing the memory up to 400gb, but no luck after this. Any input on this would be greatly appreciated.

Thank you, Ciro