[sfscode-users] convertSFS_CODE error

[Daniel R. <dr...@sa...>]

Hi,

I have stumbled across a memory corruption error when running convertSFS_CODE:

*** glibc detected *** /nfs/team143/software/sfs_code/bin/convertSFS_CODE: malloc(): memory corruption: 0x0000000000fda470 ***

The command used to generate the simulation output is:

sfs_code 1 100 -N 7000 -n 5000 -A -L 1 5000 -t 0.001 -r 0.001 -TW 0 P 0 1 140.0 0 1 -Td 0.265753 P 0 0.128575 -Td 0.342466 P 0 0.554541 -Tg 0.342466 P 0 44.822 -Tg 0.391465 P 0 284.7 -TW 0.391465 P 0 1 140.0 0 1 -TE 0.405479 --seed 22732920 --outfile small.test.o --popFreq small.test.freq

The command used in the attempt to generate the SFS is:

convertSFS_CODE small.test.o --SFS > small.test.sfs

I have created a repository https://github.com/dlrice/sfs_coder_small_test that contains the following files to recreate the error:


  *   Command used to create the simulation is located in small.test.sh.
  *   Output is located in small.test.o
  *   Frequency output is located in small.test.freq
  *   convertSFS_CODE command is located in to_sfs.sh
  *   VERSION contains the version of sfs_code / convertSFS_CODE

The converstSFS_CODE binary has worked in the past, it just seems to struggle with this particular simulation. Any advice would be appreciated.

Cheers,
Dan



-- 
 The Wellcome Trust Sanger Institute is operated by Genome Research 
 Limited, a charity registered in England with number 1021457 and a 
 company registered in England with number 2742969, whose registered 
 office is 215 Euston Road, London, NW1 2BE.