luntergroup / octopus

Bayesian haplotype-based mutation calling
MIT License
301 stars 37 forks source link

free() error with cancer calling model #150

Closed jbedo closed 3 years ago

jbedo commented 3 years ago

Describe the bug

Calling commences but after some progress octopus terminates with an free() error (see log). Seems to be pretty common, happens on ~10% of our cohort. The entire cohort called fine with 0.7.0.

Version 0.7.1

octopus version 0.7.1
Target: x86_64 Linux 3.10.0-957.el7.x86_64
SIMD extension: AVX2
Compiler: GNU 9.3.0
Boost: 1_69

Command Command line to install octopus:

$ nix build nixpkgs#octopus-caller

Command line to run octopus:

$ octopus -R ref.fa -I *.bam -o $out \\
  --bamout $evidence \\
  --threads=$NIX_BUILD_CORES \\
  --fast \\
  --max-genotypes 1000 \\
  -T chr22 \\
  -N $normal \\
  --annotations AF

Additional context Add any other context about the problem here, e.g.

dancooke commented 3 years ago

Thanks. I've already encountered this one myself and have started working on a fix.

jbedo commented 3 years ago

Thanks for your efforts! Let me know if you need me to test a patch.

dancooke commented 3 years ago

Please can you try cf636091cf9b022a6ed13bc0add4e76b9176f4bd?

jbedo commented 3 years ago

Testing now, I'm unfortunatly stuck in our computing queue waiting for an allocation. Hopefully won't be too long.

jbedo commented 3 years ago

The specific issue appears to be resolved, though I now seem to have a hanging issue where it gets stuck waiting for a futex. If I can get more details I'll open a separate issue for that.

dancooke commented 3 years ago

Ok - I'm no longer seeing this issue in my runs either so I'll close. Please do open another issue for the hanging problem.

lordkev commented 2 years ago

I'm hitting the hanging waiting for a futex also. Will see if I can find a way to replicate it and open a new issue.