Skip to content

Bug if the bed files contains too few individuals #9

@olivierlabayle

Description

@olivierlabayle

I don't know if this repo is maintained but there is a bug if the bed file contains too few individuals at this line : https://github.com/sriramlab/SCOPE/blob/e09c6546c67915f87d9aaf2581fbb300ad9c47a3/src/genotype.cpp#L300C2-L300C19

Specifically,

segment_size_hori = floor(log(Nindv)/log(3)) - 2;

will be 0 (whenever Nindv<27) and result in

terminate called after throwing an instance of 'std::length_error'
  what():  vector::_M_fill_insert
Aborted

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions