[sfscode-users] free(): invalid next size (fast): 0x00000000010052a0 Error
Brought to you by:
luricchio,
ryan_hernandez
From: Vasili <vas...@gm...> - 2019-10-16 13:22:55
|
Dear SFS_code users, I was trying to simulate a locus under negative selection within a neutrality evolving region in a way, similar to the one from Torres et al., 2018. Ultimately I want to go for a 2.2 Mb region, but to start with I'm using a 100x shorter region. My command is as follows (I'm using the sfs_code version from the aforementioned paper as I want to draw values for the selection coefficient from two negative gamma distributions): /mmg/vasili85/SFS_code/SFS_CODE_torres_et_al_2018/bin/sfs_code 1 1 \ -r 6.0443e-05 -N 18449 -s 111 -n 100 \ -TE 0.5994242 \ -Td 0 2.10709 -Td 0.437017 0.152957396219 -Td 0.546498 0.573964845871 -Tg 0.546498 60.0453856768 \ -t 0.000861 \ -L 11 2000 2000 2000 2000 2000 2000 2000 2000 2000 2000 2000 \ -v L A 2000 -v L 5 409 \ -W 0 -W L 5 2 -0.3394 0.184 0.00040244 0.0415 0.00515625 \ -a N \ -o 2000test.txt when running it I get an output that ends with a mutation in the 7th locus 7,A,1959,-14141,22117,GTG,G,0,X,X,0.0,156,0.0,0.2,0.3,0.4,0.8,0.9,0.10,0.12,0.13,0.14,0.16,0.18,0.19,0.20,0.21,0.22,0.24,0.25,0.26,0.27,0.28,0.29,0.31,0.32,0.33,0.34,0.35,0.37,0.38,0.39,0.40,0.42,0.43,0.44,0.47,0.49,0.51,0.53,0.54,0.55,0.56,0.57,0.58,0.59,0.60,0.61,0.62,0.63,0.64,0.66,0.68,0.71,0.72,0.73,0.74,0.75,0.76,0.77,0.78,0.79,0.80,0.81,0.82,0.83,0.84,0.85,0.86,0.87,0.88,0.89,0.90,0.91,0.93,0.94,0.95,0.97,0.100,0.101,0.102,0.103,0.104,0.106,0.107,0.108,0.109,0.111,0.113,0.114,0.115,0.116,0.117,0.118,0.119,0.120,0.121,0.122,0.123,0.124,0.125,0.127,0.128,0.130,0.131,0.132,0.133,0.134,0.135,0.136,0.137,0.139,0.140,0.141,0.144,0.145,0.146,0.147,0.149,0.150,0.151,0.152,0.153,0.154,0.155,0.156,0.159,0.160,0.161,0.162,0.163,0.164,0.166,0.167,0.168,0.169,0.172,0.173,0.174,0.175,0.176,0.177,0.178,0.179,0.182,0.183,0.184,0.185,0.186,0.187,0.190,0.191,0.192,0.193,0.195,0.196,0.197,0.199; so for some reason, it stops at the 8th locus, which isn't special in any way Below is the error I'm getting: *** Error in `/mmg/vasili85/SFS_code/SFS_CODE_torres_et_al_2018/bin/sfs_code': free(): invalid next size (fast): 0x00000000010052a0 *** ======= Backtrace: ========= /lib64/libc.so.6(+0x81489)[0x2ab6717bf489] /mmg/vasili85/SFS_code/SFS_CODE_torres_et_al_2018/bin/sfs_code[0x4025d8] /lib64/libc.so.6(__libc_start_main+0xf5)[0x2ab6717603d5] /mmg/vasili85/SFS_code/SFS_CODE_torres_et_al_2018/bin/sfs_code[0x40348f] ======= Memory map: ======== 00400000-0047a000 r-xp 00000000 00:2b 68985308 /mmg/vasili85/SFS_code/SFS_CODE_torres_et_al_2018/bin/sfs_code 00679000-0067a000 r--p 00079000 00:2b 68985308 /mmg/vasili85/SFS_code/SFS_CODE_torres_et_al_2018/bin/sfs_code 0067a000-0067b000 rw-p 0007a000 00:2b 68985308 /mmg/vasili85/SFS_code/SFS_CODE_torres_et_al_2018/bin/sfs_code 01004000-03531000 rw-p 00000000 00:00 0 [heap] 2ab671218000-2ab67123a000 r-xp 00000000 fd:01 201895039 /usr/lib64/ld-2.17.so 2ab67123a000-2ab67123e000 rw-p 00000000 00:00 0 2ab671249000-2ab67124c000 rw-p 00000000 00:00 0 2ab671371000-2ab671395000 rw-p 00000000 00:00 0 2ab6713b1000-2ab6713e2000 rw-p 00000000 00:00 0 2ab6713f8000-2ab67141f000 rw-p 00000000 00:00 0 2ab671439000-2ab67143a000 r--p 00021000 fd:01 201895039 /usr/lib64/ld-2.17.so 2ab67143a000-2ab67143b000 rw-p 00022000 fd:01 201895039 /usr/lib64/ld-2.17.so 2ab67143b000-2ab67143c000 rw-p 00000000 00:00 0 2ab67143c000-2ab67153d000 r-xp 00000000 fd:01 201340513 /usr/lib64/libm-2.17.so 2ab67153d000-2ab67173c000 ---p 00101000 fd:01 201340513 /usr/lib64/libm-2.17.so 2ab67173c000-2ab67173d000 r--p 00100000 fd:01 201340513 /usr/lib64/libm-2.17.so 2ab67173d000-2ab67173e000 rw-p 00101000 fd:01 201340513 /usr/lib64/libm-2.17.so 2ab67173e000-2ab671900000 r-xp 00000000 fd:01 201328540 /usr/lib64/libc-2.17.so 2ab671900000-2ab671b00000 ---p 001c2000 fd:01 201328540 /usr/lib64/libc-2.17.so 2ab671b00000-2ab671b04000 r--p 001c2000 fd:01 201328540 /usr/lib64/libc-2.17.so 2ab671b04000-2ab671b06000 rw-p 001c6000 fd:01 201328540 /usr/lib64/libc-2.17.so 2ab671b06000-2ab671b2d000 rw-p 00000000 00:00 0 2ab671b4f000-2ab671b93000 rw-p 00000000 00:00 0 2ab671ba9000-2ab671bed000 rw-p 00000000 00:00 0 2ab671bed000-2ab671c31000 rw-p 00000000 00:00 0 2ab671c38000-2ab671c7c000 rw-p 00000000 00:00 0 2ab671c89000-2ab671cba000 rw-p 00000000 00:00 0 2ab671cd2000-2ab671cf6000 rw-p 00000000 00:00 0 2ab671cf6000-2ab671d19000 rw-p 00000000 00:00 0 2ab671d37000-2ab671d7b000 rw-p 00000000 00:00 0 2ab671dd1000-2ab671e15000 rw-p 00000000 00:00 0 2ab671e15000-2ab671e2a000 r-xp 00000000 fd:01 201794459 /usr/lib64/libgcc_s-4.8.5-20150702.so.1 2ab671e2a000-2ab672029000 ---p 00015000 fd:01 201794459 /usr/lib64/libgcc_s-4.8.5-20150702.so.1 2ab672029000-2ab67202a000 r--p 00014000 fd:01 201794459 /usr/lib64/libgcc_s-4.8.5-20150702.so.1 2ab67202a000-2ab67202b000 rw-p 00015000 fd:01 201794459 /usr/lib64/libgcc_s-4.8.5-20150702.so.1 2ab674000000-2ab674021000 rw-p 00000000 00:00 0 2ab674021000-2ab678000000 ---p 00000000 00:00 0 7fff17b6a000-7fff17b8b000 rw-p 00000000 00:00 0 [stack] 7fff17bd1000-7fff17bd3000 r-xp 00000000 00:00 0 [vdso] ffffffffff600000-ffffffffff601000 r-xp 00000000 00:00 0 [vsyscall] /var/spool/slurm/slurmd/job08874/slurm_script: line 22: 4934 Aborted (core dumped) /mmg/vasili85/SFS_code/SFS_CODE_torres_et_al_2018/bin/sfs_code 1 1 -r 6.0443e-05 -N 18449 -s 111 -n 100 -TE 0.5994242 -Td 0 2.10709 -Td 0.437017 0.152957396219 -Td 0.546498 0.573964845871 -Tg 0.546498 60.0453856768 -t 0.000861 -L 11 2000 2000 2000 2000 2000 2000 2000 2000 2000 2000 2000 -v L A 2000 -v L 5 409 -W 0 -W L 5 2 -0.3394 0.184 0.00040244 0.0415 0.00515625 -a N -o 2000test.txt I've tried different random seeds and allocating up to 80g of RAM (some desperate moves). Most likely I'm doing something obviously wrong, but I can't see whether the problem lies. And yes, that's my first experience with sfs_code and forward simulators in general. Would greatly appreciate any hint. Kind regards, Vasili <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail> Без вирусов. www.avast.ru <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail> <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2> |