Segmentation fault

[dreamy...@gmail.com]

Hi Laurent,

I tried to do simulation using fsc25, but met the "Segmentation fault" error. I did simulation following your manual, and here is the input file I used (1popDNA.par):

//Number of population samples (demes) 

1 

//Population effective sizes (number of genes) 

1000000 

//Sample sizes 

100

//Growth rates  : negative growth implies population expansion 

0 

//Number of migration matrices : 0 implies no migration between demes 

0 

//historical event: time, source, sink, migrants, new size, growth rate, migr. matrix  

0  historical event  

//Number of independent loci [chromosome]  

1 0 

//Per chromosome: Number of linkage blocks 

1 

//per Block: data type, num loci, rec. rate and mut rate + optional parameters 

DNA 50000 0.00000001 0.000000005 0.33 

And I used this command line:    ./fsc25 -i 1popDNA.par -n 10

to get 10 samples. However, the programme can only give 5 samples and then give the error which is as below:

        Genealogy # 5/10

Segmentation fault (core dumped)

What's more, even the 5 samples given by fsc25 looks strange. Here is part of one of the five output files (1popDNA_1_5.arp):

#Number of independent chromosomes: 1

#Total number of polymorphic sites: 1843

# 1843 polymorphic positions on chromosome 1

#16, 21, 37, 80, 124, 140, 177, 195, 228, 229,............

#Total number of recombination events: 17523

#Positions of recombination events: 

# Chromosome 1

# 9, 10, 24, 58, 66, 73, 85, 93, 94, 113,............

# Chromosome 2850

# -1

# Chromosome 35094385

# 8322104, 0, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,............

As you can see, the part of recombination information is very strange, for example, the chromosome numbers and the "-1"s in the positions of recombination events.

I have no idea why this happened, therefore I am here asking for your help. It would great if you could find the reason for this error.

Best regards,

Yang

[Laurent Excoffier]

Dear Yang,

I think you have spotted a bug that I discovered my self a few weeks ago, and that I have now corrected. I will post a new version of fsc2 addressing this issue by the end of the week.

best and sorry about this

laurent