kr-colab / discoal

discoal is a coalescent simulation program capable of simulating models with recombination, selective sweeps, and demographic changes including population splits, admixture events, and ancient samples
GNU General Public License v3.0
48 stars 10 forks source link

trajectory too bigly #30

Open pengyan19 opened 1 year ago

pengyan19 commented 1 year ago

Hello, When i used the discoal, i used the question as follows: i also change the discoal.h, I can increase the define on line 119 to a larger number.But i still have the same question.

trajectory too bigly. step= 50000000000 freq = 1.000000. killing myself gently

this is our code: /public/home/pengyan/biosoft/discoal/discoal 18 1000 1100000 -Pt 11.057476220207999 49.14433875648 -Pr 6.450194461787999 40.236927356868 -ws 0 -Pa 739.957373322 991.263651054 -Pu 0 0.05 -Pf 0 0.2 -x 0.045454545454545456-en 0.0005123594455420342 0 1.0 -en 0.0006049945447130395 0 0.9999759043822755 -en 0.0007143781617403659 0 1.0001175973010572 -en 0.0008435384455698648 0 1.0029344631843478 -en 0.0009960510380363108 0 1.0133665133649477 -en 0.0011761380592139644 0 1.0387192901658302 -en 0.0013887849934739569 0 1.0858933719555524 -en 0.001639878705137193 0 1.1629591568885242 -en 0.001936370411402724 0 1.2882524614416981 -en 0.002286468113513557 0 1.48797233867264 -en 0.0026998638280081627 0 1.7695179130949956 -en 0.0031880019010682437 0 2.12247787175784 -en 0.003764395824957696 0 2.5250306571199475 -en 0.004445002343599473 0 2.970739432288111 -en 0.005248663201508407 0 3.4574194687113566 -en 0.006197626743972406 0 3.9814785202741185 -en 0.007318163839748038 0 4.537918146940054 -en 0.008641295161738328 0 5.120402501254934 -en 0.010203649945325195 0 5.721394463107419 -en 0.012048480032755063 0 6.332354421233643 -en 0.014226857239327694 0 6.943992688081553 -en 0.016799087224164594 0 7.546562939038391 -en 0.019836378938834744 0 8.130181295108308 -en 0.02342281603437426 0 8.685154228449564 -en 0.0276576845462343 0 9.202298112052004 -en 0.03265822149796196 0 9.673234260516272 -en 0.03856286050734812 0 10.090645320147127 -en 0.045535064150699 0 10.448481816121573 -en 0.05376784916296355 0 10.742111116853538 -en 0.06348912991905076 0 10.968404899059198 -en 0.07496802794592912 0 11.125764867130398 -en 0.08852232216070617 0 11.214089951407326 -en 0.10452724630306809 0 11.234691036293558 -en 0.12342587668851686 0 11.190161556883046 -en 0.1457413982752102 0 11.084213702335639 -en 0.17209158853612302 0 10.921490681927809 -en 0.20320591926660578 0 10.707365544085508 -en 0.23994575225178202 0 10.44773637115422 -en 0.28332818356385014 0 10.148826605437552 -en 0.3345542016489557 0 9.816997852441057 -en 0.39504193499643186 0 9.458580861298259 -en 0.46646591092727446 0 9.079728710610912 -en 0.5508034127806999 0 8.686294632113915 -en 0.650389219175433 0 8.283735376081774 -en 0.7679802387977698 0 7.877039800458334 -en 0.906831831933108 0 7.47068130992789 -en 1.0707879315572497 0 7.068591977532621 -en 1.2643874575662444 0 6.674155707278724 -en 1.4929899715442896 0 6.290217442342519 -en 1.7629240467779466 0 5.919105364754637 -en 2.0816624724269146 0 5.562663082578437 -en 2.458029123724286 0 5.234511473696049 -en 2.9024432420182307 0 5.506560850643605 -en 3.427207874522306 0 6.3459278355979185 -en 4.046850476232654 0 7.22699384954018 -en 4.77852507847275 0 7.935169380554136 -en 5.642487179059596 0 8.304416097339134 -en 6.662654488820973 0 8.293892528792151 -en 7.867269065790819 0 8.007797749589406 -en 9.289679158948962 0 7.625175385724137 -en 10.969262416833702 0 7.280676891769579 -en 12.952516002768839 0 7.00133644204606 -en 15.294343809822708 0 6.706180402692751 -en 18.059576425066336 0 6.360982768384097 -en 21.32476585449663 0 5.99743151473244 -en 25.1803047942492 0 5.657142438742077 -en 29.73292901647457 0 5.395809736221858 -en 35.10867224657684 0 5.262542308833054 -en 41.45635519662705 0 5.292219047780075 -en 48.95170556431167 0 5.485757659385585 -en 57.80222276775352 0 5.8495725717114215 -en 68.2529222874877 0 6.4037220052107 -en 80.59311871070311 0 7.182799474054486 -en 95.16443496530583 0 8.238346651283726 -en 112.37025971199664 0 9.642808531127395 -en 132.68691474376973 0 11.495193845657576 -en 156.67684122873123 0 13.928695452697328 -en 185.00417043135357 0 17.120564780273064 -en 218.45310902477905 0 21.304499979541323 -en 257.949648998146 0 26.785658733040314 -en 304.58720279813053 0 33.958082369517115 -en 359.65687279960446 0 43.323739641756035 -en 424.68319412941645 0 55.51147712305491 -en 501.46633927861114 0 71.29282274080903 -en 592.1319535785387 0 91.58980409930733 -en 699.190001339254 0 117.46780571666169 -en 825.6042509899493 0 150.1042891583789 -en 974.8743229457496 0 190.7225217155683 -en 1151.1325728858374 0 240.479247750621 -en 1359.2584901053735 0 300.2977569564728 -en 1605.013780953549 0 370.64448651510696 -en 1895.2018727837517 0 451.25921707748773 -en 2237.856261592261 0 540.866179335767 -en 2642.4629066695265 0 636.9141207474848 -en 3120.2228358028515 0 735.4130877903922 -en 3684.3622360576646 0 830.9470871025904 -en 4350.498603166347 0 916.936026141531 -en 5137.073090297949 0 986.1897121260888 -en 6065.861027736246 0 1031.7388324716105 -en 7162.574752388678 0 1047.84990137343 |gzip>HN.soft_0.msOut.gz /public/home/pengyan/biosoft/discoal/discoal 18 1000 1100000 -Pt 11.057476220207999 49.14433875648 -Pr 6.450194461787999 40.236927356868 -ws 0 -Pa 739.957373322 991.263651054 -Pu 0 0.05 -x 0.045454545454545456-en 0.0005123594455420342 0 1.0 -en 0.0006049945447130395 0 0.9999759043822755 -en 0.0007143781617403659 0 1.0001175973010572 -en 0.0008435384455698648 0 1.0029344631843478 -en 0.0009960510380363108 0 1.0133665133649477 -en 0.0011761380592139644 0 1.0387192901658302 -en 0.0013887849934739569 0 1.0858933719555524 -en 0.001639878705137193 0 1.1629591568885242 -en 0.001936370411402724 0 1.2882524614416981 -en 0.002286468113513557 0 1.48797233867264 -en 0.0026998638280081627 0 1.7695179130949956 -en 0.0031880019010682437 0 2.12247787175784 -en 0.003764395824957696 0 2.5250306571199475 -en 0.004445002343599473 0 2.970739432288111 -en 0.005248663201508407 0 3.4574194687113566 -en 0.006197626743972406 0 3.9814785202741185 -en 0.007318163839748038 0 4.537918146940054 -en 0.008641295161738328 0 5.120402501254934 -en 0.010203649945325195 0 5.721394463107419 -en 0.012048480032755063 0 6.332354421233643 -en 0.014226857239327694 0 6.943992688081553 -en 0.016799087224164594 0 7.546562939038391 -en 0.019836378938834744 0 8.130181295108308 -en 0.02342281603437426 0 8.685154228449564 -en 0.0276576845462343 0 9.202298112052004 -en 0.03265822149796196 0 9.673234260516272 -en 0.03856286050734812 0 10.090645320147127 -en 0.045535064150699 0 10.448481816121573 -en 0.05376784916296355 0 10.742111116853538 -en 0.06348912991905076 0 10.968404899059198 -en 0.07496802794592912 0 11.125764867130398 -en 0.08852232216070617 0 11.214089951407326 -en 0.10452724630306809 0 11.234691036293558 -en 0.12342587668851686 0 11.190161556883046 -en 0.1457413982752102 0 11.084213702335639 -en 0.17209158853612302 0 10.921490681927809 -en 0.20320591926660578 0 10.707365544085508 -en 0.23994575225178202 0 10.44773637115422 -en 0.28332818356385014 0 10.148826605437552 -en 0.3345542016489557 0 9.816997852441057 -en 0.39504193499643186 0 9.458580861298259 -en 0.46646591092727446 0 9.079728710610912 -en 0.5508034127806999 0 8.686294632113915 -en 0.650389219175433 0 8.283735376081774 -en 0.7679802387977698 0 7.877039800458334 -en 0.906831831933108 0 7.47068130992789 -en 1.0707879315572497 0 7.068591977532621 -en 1.2643874575662444 0 6.674155707278724 -en 1.4929899715442896 0 6.290217442342519 -en 1.7629240467779466 0 5.919105364754637 -en 2.0816624724269146 0 5.562663082578437 -en 2.458029123724286 0 5.234511473696049 -en 2.9024432420182307 0 5.506560850643605 -en 3.427207874522306 0 6.3459278355979185 -en 4.046850476232654 0 7.22699384954018 -en 4.77852507847275 0 7.935169380554136 -en 5.642487179059596 0 8.304416097339134 -en 6.662654488820973 0 8.293892528792151 -en 7.867269065790819 0 8.007797749589406 -en 9.289679158948962 0 7.625175385724137 -en 10.969262416833702 0 7.280676891769579 -en 12.952516002768839 0 7.00133644204606 -en 15.294343809822708 0 6.706180402692751 -en 18.059576425066336 0 6.360982768384097 -en 21.32476585449663 0 5.99743151473244 -en 25.1803047942492 0 5.657142438742077 -en 29.73292901647457 0 5.395809736221858 -en 35.10867224657684 0 5.262542308833054 -en 41.45635519662705 0 5.292219047780075 -en 48.95170556431167 0 5.485757659385585 -en 57.80222276775352 0 5.8495725717114215 -en 68.2529222874877 0 6.4037220052107 -en 80.59311871070311 0 7.182799474054486 -en 95.16443496530583 0 8.238346651283726 -en 112.37025971199664 0 9.642808531127395 -en 132.68691474376973 0 11.495193845657576 -en 156.67684122873123 0 13.928695452697328 -en 185.00417043135357 0 17.120564780273064 -en 218.45310902477905 0 21.304499979541323 -en 257.949648998146 0 26.785658733040314 -en 304.58720279813053 0 33.958082369517115 -en 359.65687279960446 0 43.323739641756035 -en 424.68319412941645 0 55.51147712305491 -en 501.46633927861114 0 71.29282274080903 -en 592.1319535785387 0 91.58980409930733 -en 699.190001339254 0 117.46780571666169 -en 825.6042509899493 0 150.1042891583789 -en 974.8743229457496 0 190.7225217155683 -en 1151.1325728858374 0 240.479247750621 -en 1359.2584901053735 0 300.2977569564728 -en 1605.013780953549 0 370.64448651510696 -en 1895.2018727837517 0 451.25921707748773 -en 2237.856261592261 0 540.866179335767 -en 2642.4629066695265 0 636.9141207474848 -en 3120.2228358028515 0 735.4130877903922 -en 3684.3622360576646 0 830.9470871025904 -en 4350.498603166347 0 916.936026141531 -en 5137.073090297949 0 986.1897121260888 -en 6065.861027736246 0 1031.7388324716105 -en 7162.574752388678 0 1047.84990137343 |gzip>HN.hard_0.msOut.gz