DReichLab / DReichLab/EIG

Smartpca memory error

Open
#87 2 comments 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

Dominant language
C
Stars
207
Forks
61
PR merge metrics
No merged PRs in 30d

Description

I am running smartpca with 81 million mutations, 2504 samples, .snp, .ind, and .eigenstratgeno input files, and fastmode enabled.
I use the following command to run it:

smartpca -p filepath/parameter_file

After 20 hours, smartpca fails and I get the following error message in my log file:

*** Error in `smartpca': free(): invalid size: 0x00002bbd4662b010 ***
======= Backtrace: =========
...
======= Memory map: ========
...
number of samples used: 2504 number of snps used: 81255192

I am not sure what could be causing this error. Any help is appreciated!

Contributor guide

No contributing guide indexed for this repository

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. Open a pull request that references the issue number.

Research direction

Start with the smartpca entry point and the parameter file used by the reported command, then review the log around the invalid free error. Reproduce with the reported input scale of 81,255,192 SNPs and 2,504 samples if the data is available; done means identifying the cause of the memory failure and confirming a run that completes without the allocator error.

Written by the indexing model from the issue text.

Assessment

Tech stack
c
Domain
bioinformatics
Issue type
Bug
Difficulty
4/5
Estimated time
3-5 days
Activity status
Stale
Clarity
Needs clarification
Newbie friendliness
30/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.