popgenmethods / smcpp

SMC++ infers population history from whole-genome sequence data.
GNU General Public License v3.0
149 stars 32 forks source link

split --timepoints command is useless #243

Open xiaobaiymq opened 1 year ago

xiaobaiymq commented 1 year ago

hello, When I use the split --timepoints command, no matter how I set the start and end time points, the smc++ guessed separation time is always the farthest, and the content of the model.final.json file does not change,even if it is not possible. I don't know how to modify it.

It's my debug.txt 5975 smcpp.commands.command DEBUG ['/usr/local/bin/smc++', 'split', '--em-iterations', '20', '--timepoints', '1', '10000', '--cores', '1', './analysis5/WM-1/model.final.json', './analysis5/ST-2/model.final.json', 'split_smc_WM1_ST-2/ST-2-WM-1.1.smc.gz', 'split_smc_WM1_ST-2/ST-2-WM-1.2.smc.gz', 'split_smc_WM1_ST-2/ST-2-WM-1.3.smc.gz', 'split_smc_WM1_ST-2/ST-2-WM-1.4.smc.gz', 'split_smc_WM1_ST-2/ST-2-WM-1.5.smc.gz', 'split_smc_WM1_ST-2/ST-2-WM-1.6.smc.gz', 'split_smc_WM1_ST-2/ST-2-WM-1.7.smc.gz', 'split_smc_WM1_ST-2/ST-2-WM-1.8.smc.gz', 'split_smc_WM1_ST-2/ST-2-WM-1.9.smc.gz', 'split_smc_WM1_ST-2/WM-1-ST-2.1.smc.gz', 'split_smc_WM1_ST-2/WM-1-ST-2.2.smc.gz', 'split_smc_WM1_ST-2/WM-1-ST-2.3.smc.gz', 'split_smc_WM1_ST-2/WM-1-ST-2.4.smc.gz', 'split_smc_WM1_ST-2/WM-1-ST-2.5.smc.gz', 'split_smc_WM1_ST-2/WM-1-ST-2.6.smc.gz', 'split_smc_WM1_ST-2/WM-1-ST-2.7.smc.gz', 'split_smc_WM1_ST-2/WM-1-ST-2.8.smc.gz', 'split_smc_WM1_ST-2/WM-1-ST-2.9.smc.gz', '-o', 'split_WM1_ST-2_analysis/'] 5977 smcpp.commands.command DEBUG Namespace(Nmax=1000.0, Nmin=0.001, algorithm='L-BFGS-B', base='model', command='split', cores=1, data=['split_smc_WM1_ST-2/ST-2-WM-1.1.smc.gz', 'split_smc_WM1_ST-2/ST-2-WM-1.2.smc.gz', 'split_smc_WM1_ST-2/ST-2-WM-1.3.smc.gz', 'split_smc_WM1_ST-2/ST-2-WM-1.4.smc.gz', 'split_smc_WM1_ST-2/ST-2-WM-1.5.smc.gz', 'split_smc_WM1_ST-2/ST-2-WM-1.6.smc.gz', 'split_smc_WM1_ST-2/ST-2-WM-1.7.smc.gz', 'split_smc_WM1_ST-2/ST-2-WM-1.8.smc.gz', 'split_smc_WM1_ST-2/ST-2-WM-1.9.smc.gz', 'split_smc_WM1_ST-2/WM-1-ST-2.1.smc.gz', 'split_smc_WM1_ST-2/WM-1-ST-2.2.smc.gz', 'split_smc_WM1_ST-2/WM-1-ST-2.3.smc.gz', 'split_smc_WM1_ST-2/WM-1-ST-2.4.smc.gz', 'split_smc_WM1_ST-2/WM-1-ST-2.5.smc.gz', 'split_smc_WM1_ST-2/WM-1-ST-2.6.smc.gz', 'split_smc_WM1_ST-2/WM-1-ST-2.7.smc.gz', 'split_smc_WM1_ST-2/WM-1-ST-2.8.smc.gz', 'split_smc_WM1_ST-2/WM-1-ST-2.9.smc.gz'], emiterations=20, ftol=0.0001, lambda=None, length_cutoff=None, multi=False, no_initialize=False, nonseg_cutoff=None, outdir='split_WM1_ST-2_analysis/', polarization_error=0.5, pop1='./analysis5/WM-1/model.final.json', pop2='./analysis5/ST-2/model.final.json', regularization_penalty=6, seed=0, thinning=None, timepoints=[1.0, 10000.0], unfold=False, verbose=0, w=100, xtol=0.1) 5979 smcpp.analysis.base INFO theta: 0.000100 5981 smcpp.analysis.base INFO rho: 0.000100 5981 smcpp.analysis.base DEBUG Polarization error p=0.500000 5981 smcpp.data_filter DEBUG LoadData() 5981 smcpp.data_filter INFO Loading data... 32585 smcpp.data_filter INFO 0.62 Gb of data 32586 smcpp.data_filter DEBUG Contig(pid=('ST-2', 'WM-1'), data=array([[ 468, -1, 0, ..., -1, 0, 0], [ 1, 0, 4, ..., 0, 0, 4], [ 1, 0, 0, ..., 0, 0, 4], ..., [ 11, -1, 0, ..., -1, 0, 0], [ 2, 0, 0, ..., 0, 0, 4], [35861, -1, 0, ..., -1, 0, 0]], dtype=int32), n=array([26, 4]), a=array([2, 0]), fn='split_smc_WM1_ST-2/ST-2-WM-1.6.smc.gz') 32589 smcpp.data_filter DEBUG Contig(pid=('WM-1', 'ST-2'), data=array([[10665, -1, 0, ..., -1, 0, 0], [ 1, -1, 0, ..., 0, 0, 12], [ 1, 0, 0, ..., 0, 0, 28], ..., [ 1, 1, 1, ..., 0, 0, 28], [ 1, 0, 0, ..., 0, 0, 28], [ 726, -1, 0, ..., -1, 0, 0]], dtype=int32), n=array([ 2, 28]), a=array([2, 0]), fn='split_smc_WM1_ST-2/WM-1-ST-2.7.smc.gz') 32589 smcpp.data_filter DEBUG Contig(pid=('ST-2', 'WM-1'), data=array([[10665, -1, 0, ..., -1, 0, 0], [ 1, -1, 0, ..., 0, 0, 0], [ 1, 0, 0, ..., 0, 0, 4], ..., [ 1, 0, 0, ..., 0, 2, 4], [ 1, 0, 0, ..., 0, 0, 4], [ 726, -1, 0, ..., -1, 0, 0]], dtype=int32), n=array([26, 4]), a=array([2, 0]), fn='split_smc_WM1_ST-2/ST-2-WM-1.7.smc.gz') 32593 smcpp.data_filter DEBUG Contig(pid=('WM-1', 'ST-2'), data=array([[ 532, -1, 0, ..., -1, 0, 0], [ 1, -1, 0, ..., 0, 0, 28], [ 1, 0, 0, ..., 0, 0, 28], ..., [ 1, -1, 0, ..., 0, 0, 28], [ 1, 0, 0, ..., 0, 0, 28], [24364, -1, 0, ..., -1, 0, 0]], dtype=int32), n=array([ 2, 28]), a=array([2, 0]), fn='split_smc_WM1_ST-2/WM-1-ST-2.1.smc.gz') 32593 smcpp.data_filter DEBUG Contig(pid=('ST-2', 'WM-1'), data=array([[133017, -1, 0, ..., -1, 0, 0], [ 1, 1, 7, ..., 0, 0, 4], [ 1, 0, 0, ..., 0, 0, 4], ..., [ 1, 1, 5, ..., 0, 0, 4], [ 1, 0, 0, ..., 0, 0, 4], [ 9658, -1, 0, ..., -1, 0, 0]], dtype=int32), n=array([26, 4]), a=array([2, 0]), fn='split_smc_WM1_ST-2/ST-2-WM-1.5.smc.gz') 32593 smcpp.data_filter DEBUG Contig(pid=('ST-2', 'WM-1'), data=array([[23366, -1, 0, ..., -1, 0, 0], [ 1, 0, 0, ..., 0, 3, 4], [ 1, 0, 0, ..., 0, 0, 4], ..., [ 1, 0, 0, ..., 0, 0, 2], [ 1, 0, 0, ..., 0, 0, 4], [ 2992, -1, 0, ..., -1, 0, 0]], dtype=int32), n=array([26, 4]), a=array([2, 0]), fn='split_smc_WM1_ST-2/ST-2-WM-1.3.smc.gz') 32594 smcpp.data_filter DEBUG Contig(pid=('WM-1', 'ST-2'), data=array([[133017, -1, 0, ..., -1, 0, 0], [ 1, 0, 0, ..., 0, 8, 26], [ 1, 0, 0, ..., 0, 0, 28], ..., [ 1, 0, 0, ..., 0, 6, 24], [ 1, 0, 0, ..., 0, 0, 28], [ 9658, -1, 0, ..., -1, 0, 0]], dtype=int32), n=array([ 2, 28]), a=array([2, 0]), fn='split_smc_WM1_ST-2/WM-1-ST-2.5.smc.gz') 32594 smcpp.data_filter DEBUG Contig(pid=('ST-2', 'WM-1'), data=array([[ 532, -1, 0, ..., -1, 0, 0], [ 1, 0, 0, ..., 0, 0, 2], [ 1, 0, 0, ..., 0, 0, 4], ..., [ 1, 0, 0, ..., 0, 0, 2], [ 1, 0, 0, ..., 0, 0, 4], [24364, -1, 0, ..., -1, 0, 0]], dtype=int32), n=array([26, 4]), a=array([2, 0]), fn='split_smc_WM1_ST-2/ST-2-WM-1.1.smc.gz') 32594 smcpp.data_filter DEBUG Contig(pid=('ST-2', 'WM-1'), data=array([[ 1207, -1, 0, ..., -1, 0, 0], [ 1, 0, 0, ..., 0, 0, 0], [ 1, 0, 0, ..., 0, 0, 4], ..., [ 2, -1, 0, ..., -1, 0, 0], [ 2, 0, 0, ..., 0, 0, 4], [105981, -1, 0, ..., -1, 0, 0]], dtype=int32), n=array([26, 4]), a=array([2, 0]), fn='split_smc_WM1_ST-2/ST-2-WM-1.9.smc.gz') 32595 smcpp.data_filter DEBUG Contig(pid=('ST-2', 'WM-1'), data=array([[3213, -1, 0, ..., -1, 0, 0], [ 1, 0, 0, ..., 0, 0, 4], [ 1, 0, 0, ..., 0, 0, 4], ..., [ 76, -1, 0, ..., -1, 0, 0], [ 2, 0, 0, ..., 0, 0, 4], [8347, -1, 0, ..., -1, 0, 0]], dtype=int32), n=array([26, 4]), a=array([2, 0]), fn='split_smc_WM1_ST-2/ST-2-WM-1.2.smc.gz') 32595 smcpp.data_filter DEBUG Contig(pid=('WM-1', 'ST-2'), data=array([[ 1207, -1, 0, ..., -1, 0, 0], [ 1, -1, 0, ..., 0, 0, 26], [ 1, 0, 0, ..., 0, 0, 28], ..., [ 2, -1, 0, ..., -1, 0, 0], [ 2, 0, 0, ..., 0, 0, 28], [105981, -1, 0, ..., -1, 0, 0]], dtype=int32), n=array([ 2, 28]), a=array([2, 0]), fn='split_smc_WM1_ST-2/WM-1-ST-2.9.smc.gz') 32596 smcpp.data_filter DEBUG Contig(pid=('WM-1', 'ST-2'), data=array([[ 468, -1, 0, ..., -1, 0, 0], [ 1, 0, 0, ..., 0, 4, 18], [ 1, 0, 0, ..., 0, 0, 28], ..., [ 11, -1, 0, ..., -1, 0, 0], [ 2, 0, 0, ..., 0, 0, 28], [35861, -1, 0, ..., -1, 0, 0]], dtype=int32), n=array([ 2, 28]), a=array([2, 0]), fn='split_smc_WM1_ST-2/WM-1-ST-2.6.smc.gz') 32596 smcpp.data_filter DEBUG Contig(pid=('WM-1', 'ST-2'), data=array([[7611, -1, 0, ..., -1, 0, 0], [ 1, -1, 0, ..., 0, 0, 16], [ 1, 0, 0, ..., 0, 0, 28], ..., [ 1, 1, 0, ..., 0, 14, 28], [ 1, 0, 0, ..., 0, 0, 28], [ 189, -1, 0, ..., -1, 0, 0]], dtype=int32), n=array([ 2, 28]), a=array([2, 0]), fn='split_smc_WM1_ST-2/WM-1-ST-2.8.smc.gz') 32596 smcpp.data_filter DEBUG Contig(pid=('WM-1', 'ST-2'), data=array([[3109, -1, 0, ..., -1, 0, 0], [ 1, 1, 0, ..., 0, 7, 28], [ 1, 0, 0, ..., 0, 0, 28], ..., [ 1, 1, 2, ..., 0, 24, 28], [ 1, 0, 0, ..., 0, 0, 28], [2796, -1, 0, ..., -1, 0, 0]], dtype=int32), n=array([ 2, 28]), a=array([2, 0]), fn='split_smc_WM1_ST-2/WM-1-ST-2.4.smc.gz') 32597 smcpp.data_filter DEBUG Contig(pid=('ST-2', 'WM-1'), data=array([[3109, -1, 0, ..., -1, 0, 0], [ 1, 0, 7, ..., 0, 1, 4], [ 1, 0, 0, ..., 0, 0, 4], ..., [ 1, 2, 22, ..., 0, 3, 4], [ 1, 0, 0, ..., 0, 0, 4], [2796, -1, 0, ..., -1, 0, 0]], dtype=int32), n=array([26, 4]), a=array([2, 0]), fn='split_smc_WM1_ST-2/ST-2-WM-1.4.smc.gz') 32597 smcpp.data_filter DEBUG Contig(pid=('ST-2', 'WM-1'), data=array([[7611, -1, 0, ..., -1, 0, 0], [ 1, -1, 0, ..., 0, 0, 0], [ 1, 0, 0, ..., 0, 0, 4], ..., [ 1, 1, 13, ..., 0, 1, 4], [ 1, 0, 0, ..., 0, 0, 4], [ 189, -1, 0, ..., -1, 0, 0]], dtype=int32), n=array([26, 4]), a=array([2, 0]), fn='split_smc_WM1_ST-2/ST-2-WM-1.8.smc.gz') 32598 smcpp.data_filter DEBUG Contig(pid=('WM-1', 'ST-2'), data=array([[23366, -1, 0, ..., -1, 0, 0], [ 1, 2, 1, ..., 0, 0, 28], [ 1, 0, 0, ..., 0, 0, 28], ..., [ 1, -1, 0, ..., 0, 0, 26], [ 1, 0, 0, ..., 0, 0, 28], [ 2992, -1, 0, ..., -1, 0, 0]], dtype=int32), n=array([ 2, 28]), a=array([2, 0]), fn='split_smc_WM1_ST-2/WM-1-ST-2.3.smc.gz') 32598 smcpp.data_filter DEBUG Contig(pid=('WM-1', 'ST-2'), data=array([[3213, -1, 0, ..., -1, 0, 0], [ 1, 0, 0, ..., 0, 0, 26], [ 1, 0, 0, ..., 0, 0, 28], ..., [ 76, -1, 0, ..., -1, 0, 0], [ 2, 0, 0, ..., 0, 0, 28], [8347, -1, 0, ..., -1, 0, 0]], dtype=int32), n=array([ 2, 28]), a=array([2, 0]), fn='split_smc_WM1_ST-2/WM-1-ST-2.2.smc.gz') 32598 smcpp.data_filter INFO 2 populations 32598 smcpp.data_filter DEBUG RecodeNonseg(cutoff=None) 36421 smcpp.data_filter DEBUG Compress() 45411 smcpp.data_filter DEBUG BreakLongSpans(cutoff=100000) 47112 smcpp.estimation_tools DEBUG Long missing spans: [133017 173952] (base pairs) 47622 smcpp.estimation_tools DEBUG Long missing spans: [133859 102743 127176] (base pairs) 48074 smcpp.estimation_tools DEBUG Long missing spans: [133017 173952] (base pairs) 48883 smcpp.estimation_tools DEBUG Long missing spans: [105981] (base pairs) 49575 smcpp.estimation_tools DEBUG Long missing spans: [105981] (base pairs) 50194 smcpp.estimation_tools DEBUG Long missing spans: [1490394 259916 236463] (base pairs) 50524 smcpp.estimation_tools DEBUG Long missing spans: [106424 247157] (base pairs) 50877 smcpp.estimation_tools DEBUG Long missing spans: [106424 247157] (base pairs) 51224 smcpp.estimation_tools DEBUG Long missing spans: [1490394 259916 236463] (base pairs) 51594 smcpp.estimation_tools DEBUG Long missing spans: [133859 102743 127176] (base pairs) 52109 smcpp.data_filter DEBUG DropSmallContigs(cutoff=100000) 52356 smcpp.data_filter DEBUG Watterson() 60940 smcpp.data_filter DEBUG sites: 6526515/214775751 watterson:0.030388 60942 smcpp.data_filter DEBUG CountMutations(w=45454) 61118 smcpp.analysis.base DEBUG Creating inference manager... 61119 smcpp.analysis.base DEBUG Creating inference manager for ('ST-2', 'WM-1') 78286 smcpp.analysis.base DEBUG Creating inference manager for ('WM-1', 'ST-2') 94893 smcpp.optimize.plugins.progress_printer INFO Starting EM algorithm... 94894 smcpp.optimize.plugins.progress_printer INFO EM iteration 1 of 1... 94894 smcpp.analysis.base INFO Running E-step 251142 smcpp.analysis.base INFO E-step completed 251172 smcpp.optimize.plugins.loglikelihood_monitor INFO Loglik: -179958373.076126 251173 smcpp.optimize.plugins.parameter_optimizer INFO Updating split, bounds (0.000000, 22000.000000) 251844 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 251844 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.5 251844 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -89861897.9 251844 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 252401 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 252401 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 252401 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -90092992.4 252401 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 252404 smcpp.analysis.base DEBUG reg: 0.0 [] 252405 smcpp.analysis.base DEBUG Q: -179958373.3 [] 252405 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(11000.000000)=179958373.281117 252405 smcpp.optimize.plugins.parameter_optimizer DEBUG Old split: Q(11000.000000)=179958373.281117 253046 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 253046 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.5 253046 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -77455278.3 253046 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 253588 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 253588 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 253588 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -77686384.3 253588 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 253591 smcpp.analysis.base DEBUG reg: 0.0 [] 253592 smcpp.analysis.base DEBUG Q: -155145145.6 [] 253592 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(8403.252248)=155145145.603265 254251 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 254251 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.5 254251 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -102234672.6 254251 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 254813 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 254813 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 254813 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -102465756.9 254813 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 254816 smcpp.analysis.base DEBUG reg: 0.0 [] 254817 smcpp.analysis.base DEBUG Q: -204703912.5 [] 254817 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(13596.747752)=204703912.503458 255439 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 255439 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.5 255439 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -62081441.6 255439 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 255968 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 255968 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 255968 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -62312554.8 255968 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 255971 smcpp.analysis.base DEBUG reg: 0.0 [] 255971 smcpp.analysis.base DEBUG Q: -124397479.4 [] 255972 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(5193.495505)=124397479.360655 256582 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 256582 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.5 256582 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -52487670.2 256582 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 257091 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 257092 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 257092 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -52710093.2 257092 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 257095 smcpp.analysis.base DEBUG reg: 0.0 [] 257095 smcpp.analysis.base DEBUG Q: -105201246.4 [] 257095 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(3209.756743)=105201246.374038 257694 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 257694 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.5 257695 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -45445061.2 257695 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 258195 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 258196 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 258196 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -45697528.5 258196 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 258199 smcpp.analysis.base DEBUG reg: 0.0 [] 258199 smcpp.analysis.base DEBUG Q: -91146072.8 [] 258199 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(1983.738762)=91146072.763554 258781 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 258782 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.5 258782 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -41102286.6 258782 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 259265 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 259266 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 259266 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -41354509.6 259266 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 259269 smcpp.analysis.base DEBUG reg: 0.0 [] 259269 smcpp.analysis.base DEBUG Q: -82460279.2 [] 259269 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(1226.017980)=82460279.222820 259841 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 259841 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.5 259841 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -38435436.5 259841 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 260312 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 260312 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 260312 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -38687437.4 260312 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 260315 smcpp.analysis.base DEBUG reg: 0.0 [] 260315 smcpp.analysis.base DEBUG Q: -77126356.8 [] 260315 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(757.720782)=77126356.840759 260880 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 260880 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.5 260881 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -36798769.9 260881 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 261348 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 261348 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 261348 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -37050589.7 261348 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 261351 smcpp.analysis.base DEBUG reg: 0.0 [] 261351 smcpp.analysis.base DEBUG Q: -73852842.6 [] 261351 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(468.297198)=73852842.567145 261904 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 261904 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.5 261904 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -35794237.1 261904 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 262361 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 262361 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 262361 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -36045938.0 262361 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 262364 smcpp.analysis.base DEBUG reg: 0.0 [] 262365 smcpp.analysis.base DEBUG Q: -71843658.0 [] 262365 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(289.423585)=71843658.039864 262916 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 262916 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.5 262916 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -34784914.2 262916 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 263370 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 263370 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 263370 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -35033139.8 263370 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 263373 smcpp.analysis.base DEBUG reg: 0.0 [] 263373 smcpp.analysis.base DEBUG Q: -69821537.1 [] 263373 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(178.873613)=69821537.051811 263923 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 263924 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.4 263924 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -32246829.0 263924 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 264374 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 264374 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 264374 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -32493610.0 264374 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 264377 smcpp.analysis.base DEBUG reg: 0.0 [] 264377 smcpp.analysis.base DEBUG Q: -64743921.9 [] 264377 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(110.549972)=64743921.887973 264925 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 264925 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.0 264925 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -29979687.1 264925 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 265379 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 265380 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 265380 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -30220454.5 265380 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 265383 smcpp.analysis.base DEBUG reg: 0.0 [] 265383 smcpp.analysis.base DEBUG Q: -60203624.0 [] 265383 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(68.323640)=60203624.016530 265931 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 265932 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -6.2 265932 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -35152511.0 265932 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 266392 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 266392 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 266392 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -33187214.3 266392 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 266395 smcpp.analysis.base DEBUG reg: 0.0 [] 266395 smcpp.analysis.base DEBUG Q: -68343209.0 [] 266396 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(42.226332)=68343209.005841 266946 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 266947 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.1 266947 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -30716993.8 266947 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 267398 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 267399 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 267399 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -30889337.2 267399 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 267402 smcpp.analysis.base DEBUG reg: 0.0 [] 267402 smcpp.analysis.base DEBUG Q: -61609813.6 [] 267402 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(80.678978)=61609813.601536 267956 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 267956 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.0 267956 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -29775642.6 267956 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 268411 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 268411 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 268411 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -30017555.7 268411 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 268415 smcpp.analysis.base DEBUG reg: 0.0 [] 268415 smcpp.analysis.base DEBUG Q: -59796680.9 [] 268415 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(58.355356)=59796680.857300 268968 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 268968 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.0 268968 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -29758008.1 268968 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 269423 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 269423 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 269423 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -29999864.3 269423 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 269427 smcpp.analysis.base DEBUG reg: 0.0 [] 269427 smcpp.analysis.base DEBUG Q: -59761355.0 [] 269427 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(57.096497)=59761354.959181 269992 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 269992 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.1 269992 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -29598298.3 269992 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 270449 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 270449 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 270449 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -29827449.2 270449 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 270453 smcpp.analysis.base DEBUG reg: 0.0 [] 270453 smcpp.analysis.base DEBUG Q: -59429230.1 [] 270453 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(51.416599)=59429230.125341 271007 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 271007 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.2 271008 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -29606237.7 271008 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 271463 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 271463 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 271463 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -29779396.8 271463 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 271466 smcpp.analysis.base DEBUG reg: 0.0 [] 271466 smcpp.analysis.base DEBUG Q: -59389117.3 [] 271466 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(47.906229)=59389117.273042 272018 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 272018 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.2 272018 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -29588876.7 272018 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 272475 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 272475 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 272475 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -29779110.8 272475 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 272479 smcpp.analysis.base DEBUG reg: 0.0 [] 272479 smcpp.analysis.base DEBUG Q: -59371470.2 [] 272479 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(48.545325)=59371470.178223 273034 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 273034 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.2 273034 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -29580660.7 273034 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 273490 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 273490 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 273490 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -29784881.4 273490 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 273494 smcpp.analysis.base DEBUG reg: 0.0 [] 273494 smcpp.analysis.base DEBUG Q: -59369024.7 [] 273494 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(49.241201)=59369024.738190 274047 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 274048 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.2 274048 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -29582545.8 274048 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 274502 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 274502 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 274503 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -29782191.5 274503 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 274506 smcpp.analysis.base DEBUG reg: 0.0 [] 274506 smcpp.analysis.base DEBUG Q: -59368219.9 [] 274506 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(48.990600)=59368219.937724 275061 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 275062 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.2 275062 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -29582363.6 275062 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 275519 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 275519 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 275520 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -29782373.0 275520 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 275523 smcpp.analysis.base DEBUG reg: 0.0 [] 275523 smcpp.analysis.base DEBUG Q: -59368219.2 [] 275523 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(49.009592)=59368219.240006 276076 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 276077 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.2 276077 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -29582441.3 276077 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 276534 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 276534 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 276534 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -29782294.3 276535 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 276538 smcpp.analysis.base DEBUG reg: 0.0 [] 276538 smcpp.analysis.base DEBUG Q: -59368218.3 [] 276538 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(49.001406)=59368218.253824 277088 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 277088 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.2 277088 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -29582441.8 277088 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 277545 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 277546 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 277546 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -29782293.8 277546 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 277549 smcpp.analysis.base DEBUG reg: 0.0 [] 277549 smcpp.analysis.base DEBUG Q: -59368218.3 [] 277549 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(49.001358)=59368218.253782 278105 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 278105 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.2 278105 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -29582441.8 278105 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 278561 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 278561 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 278561 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -29782293.8 278561 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 278564 smcpp.analysis.base DEBUG reg: 0.0 [] 278564 smcpp.analysis.base DEBUG Q: -59368218.3 [] 278564 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(49.001352)=59368218.253841 279119 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 279119 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.2 279119 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -29582441.5 279119 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 279577 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 279577 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 279577 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -29782294.0 279577 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 279580 smcpp.analysis.base DEBUG reg: 0.0 [] 279580 smcpp.analysis.base DEBUG Q: -59368218.3 [] 279580 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(49.001380)=59368218.253793 280135 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 280136 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.2 280136 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -29582441.7 280136 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 280594 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 280594 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -398.1 280594 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -29782293.9 280594 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 280597 smcpp.analysis.base DEBUG reg: 0.0 [] 280597 smcpp.analysis.base DEBUG Q: -59368218.3 [] 280597 smcpp.optimize.plugins.parameter_optimizer DEBUG split f(49.001366)=59368218.253844 280598 smcpp.optimize.plugins.parameter_optimizer INFO New split: 49.0014 281185 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 281185 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -5.8 281185 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -31491627.9 281185 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 281647 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 281648 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -433.4 281648 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -31660691.5 281648 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 281652 smcpp.analysis.base DEBUG reg: 0.0 [] 281652 smcpp.analysis.base DEBUG Q: -63155838.0 [] 281652 smcpp.optimize.plugins.scale_optimizer DEBUG scale Q(-0.236068)=-63155838.037279 282219 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 282219 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -4.6 282219 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -28040607.0 282219 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 282672 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 282673 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -359.8 282673 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -28307649.1 282673 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 282676 smcpp.analysis.base DEBUG reg: 0.0 [] 282676 smcpp.analysis.base DEBUG Q: -56351699.7 [] 282676 smcpp.optimize.plugins.scale_optimizer DEBUG scale Q(0.236068)=-56351699.734827 283237 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 283237 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -4.0 283237 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -26867190.1 283237 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 283687 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 283687 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -319.2 283688 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -27231475.8 283688 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 283691 smcpp.analysis.base DEBUG reg: 0.0 [] 283691 smcpp.analysis.base DEBUG Q: -54102068.5 [] 283691 smcpp.optimize.plugins.scale_optimizer DEBUG scale Q(0.527864)=-54102068.510383 284247 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 284247 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -3.7 284247 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -26568168.6 284247 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 284696 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 284696 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -294.5 284696 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -27002885.8 284696 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 284699 smcpp.analysis.base DEBUG reg: 0.0 [] 284699 smcpp.analysis.base DEBUG Q: -53574431.9 [] 284699 smcpp.optimize.plugins.scale_optimizer DEBUG scale Q(0.821371)=-53574431.868902 285258 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 285259 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -3.7 285259 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -26561100.3 285259 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 285708 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 285708 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -297.9 285708 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -26985294.5 285708 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 285711 smcpp.analysis.base DEBUG reg: 0.0 [] 285712 smcpp.analysis.base DEBUG Q: -53549775.8 [] 285712 smcpp.optimize.plugins.scale_optimizer DEBUG scale Q(0.763607)=-53549775.803717 286273 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 286273 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -3.7 286273 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -26560389.1 286273 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 286719 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 286719 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -297.5 286720 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -26985808.5 286720 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 286723 smcpp.analysis.base DEBUG reg: 0.0 [] 286723 smcpp.analysis.base DEBUG Q: -53549578.2 [] 286723 smcpp.optimize.plugins.scale_optimizer DEBUG scale Q(0.769872)=-53549578.171092 287280 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 287280 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -3.7 287280 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -26560512.5 287280 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 287732 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 287732 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -297.6 287732 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -26985667.6 287732 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 287736 smcpp.analysis.base DEBUG reg: 0.0 [] 287736 smcpp.analysis.base DEBUG Q: -53549560.8 [] 287736 smcpp.optimize.plugins.scale_optimizer DEBUG scale Q(0.768512)=-53549560.831791 288306 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 288306 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -3.7 288306 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -26560515.0 288306 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 288755 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 288756 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -297.6 288756 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -26985665.1 288756 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 288759 smcpp.analysis.base DEBUG reg: 0.0 [] 288759 smcpp.analysis.base DEBUG Q: -53549560.8 [] 288759 smcpp.optimize.plugins.scale_optimizer DEBUG scale Q(0.768486)=-53549560.825472 289319 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 289320 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -3.7 289320 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -26560515.3 289320 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 289768 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 289768 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -297.6 289769 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -26985664.8 289769 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 289772 smcpp.analysis.base DEBUG reg: 0.0 [] 289772 smcpp.analysis.base DEBUG Q: -53549560.8 [] 289772 smcpp.optimize.plugins.scale_optimizer DEBUG scale Q(0.768483)=-53549560.825433 290333 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q1: 0.0 290334 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q2: -3.7 290334 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q3: -26560515.6 290334 smcpp._smcpp DEBUG im(('ST-2', 'WM-1')).q4: -1539.7 290783 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q1: 0.0 290783 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q2: -297.6 290783 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q3: -26985664.5 290783 smcpp._smcpp DEBUG im(('WM-1', 'ST-2')).q4: -1539.7 290787 smcpp.analysis.base DEBUG reg: 0.0 [] 290787 smcpp.analysis.base DEBUG Q: -53549560.8 [] 290787 smcpp.optimize.plugins.scale_optimizer DEBUG scale Q(0.768480)=-53549560.825748 290818 smcpp.optimize.plugins.progress_printer INFO Current model:

Pop. 1: 5.17 5.17 5.17 5.16 5.17 5.18 5.20 5.20 5.10 4.84 4.35 3.65 2.82 2.00 1.47 1.95 6.09 9.41 9.52 3.85 5.03 5.14 6.52 7.82 8.78 9.53 9.98 176.72 176.68 176.55 176.34 176.04 175.69 175.30 174.93 174.60 174.35 174.24 174.29 174.55 175.03 175.75 176.70 177.83 179.11 180.41 181.58 182.44 182.76 182.21 180.48 177.26 172.25 165.22 156.06 144.89 131.99 117.76 102.81 87.82 73.37 60.02 48.24 38.22 29.96 23.25 17.89 13.67 10.50 8.23 6.71 5.67 4.88 4.22 3.80 4.19 6.77 18.08 67.93 281.67 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 256.23 45.09 18.15 15.29 21.53 40.54 81.53 140.06 169.46 160.64 156.81 205.26 370.55 754.53 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 Pop. 2: 0.23 0.24 0.26 0.32 0.38 0.39 0.30 0.22 0.42 0.32 -0.36 0.15 1.91 3.74 5.45 6.74 7.31 6.30 3.53 3.50 1.26 1.26 1.26 1.27 1.27 1.28 1.29 1.31 1.32 1.35 1.38 1.41 1.43 1.46 1.48 1.49 1.49 1.47 1.43 1.38 1.33 1.27 1.24 1.25 1.31 1.41 1.51 1.58 1.58 1.47 1.23 0.97 0.77 0.68 0.70 0.84 1.17 1.86 3.23 5.83 10.42 18.28 31.68 54.74 94.00 158.16 257.01 399.02 590.51 831.93 1000.00 1000.00 1000.00 1000.00 1000.00 753.14 420.84 197.87 80.94 30.01 17.10 17.81 27.31 49.77 86.84 117.02 102.56 68.85 48.71 49.17 72.72 129.29 239.28 453.54 868.04 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 1000.00 Split: 49.001

290980 smcpp.optimize.plugins.ascii_plotter INFO Plot of current model:

  1e+07 +-------------------------------------------------------------+   
        |+      +       +      +       +       +       +      +      +|   
        |+                             ####****        ***************|   
        |+                            ## *#   *       *Pop. 2 #######+|   
  1e+06 |-+                          #   * #  *       *             +-|   
        |*****************          ##   *  # *    ***               +|   
        |+                ***       #   *   #  *##**                 +|   
 100000 |-+                  **    #    *   #  #  *                 +-|   

N_e |+ # # # * +|
|+ # # +|
|+ # +|
10000 |-+ # ***** +-|
|+ # +|
|################### # +|
|+ + + ####+ + + + + +|
1000 +-------------------------------------------------------------+
1 10 100 1000 10000 100000 1e+06 1e+07 1e+08 Generations

It's my model.final.json { "alpha": 1, "hidden_states": { "ST-2": [ 0.0, 12.656215235611679, 26.250756614130886, 40.93397988813272, 56.89555231429243, 74.37947672199772, 93.70705846720283, 115.31345661013056, 139.80877053667945, 168.0865351892005, 201.5320320230655, 242.46601191119817, 295.23909049026827, 369.6185672122657, 496.7711225133337, Infinity ], "WM-1": [ 0.0, 22.39300894795084, 46.446225574347444, 72.42567867611761, 100.66695203536335, 131.60176694001194, 165.7986182818587, 204.0274456161297, 247.3677154937037, 297.4003852218705, 356.5764734857648, 429.00215216188246, 522.3750917676234, 653.976858707635, 878.9515652533883, Infinity ] }, "model": { "class": "SMCTwoPopulationModel", "model1": { "N0": 2272.727272727273, "class": "SMCModel", "knots": [ 0.00021999999999999998, 0.0004042964845327641, 0.00074298021547978, 0.001365383132708503, 0.0025091800027017288, 0.004611148427964639, 0.008473959541294934, 0.01557269114826553, 0.02861811039071931, 0.05259182465881085, 0.09664859011236712, 0.17761220553399135, 0.32639995594319754, 0.5998288851794722, 1.1023123163603075, 2.0257317925529112, 3.7227101924336217, 6.841266562432742, 12.572272822468664, 22.39300894795084, 72.42567867611761, 131.60176694001194, 204.0274456161297, 297.4003852218705, 429.00215216188246, 653.976858707635, 21999.999999999996 ], "pid": "WM-1", "spline_class": "CubicSpline", "y": [ 5.174560617875921, 5.171992322251927, 5.165342225012986, 5.160471248067957, 5.165227699564784, 5.183272945600572, 5.204915722770314, 5.1969949553190355, 5.099049013721095, 4.83584443948981, 4.3509291588789525, 3.649492518218441, 2.8184860346990055, 1.996920434721217, 1.4710016027910264, 1.9549862740174748, 6.086236735867133, 9.414356313350984, 9.516910834981413, 3.8494943339022667, 5.034423872673192, 5.1367463963431925, 6.516313528030586, 7.821115392722951, 8.776651763939654, 9.533031451419488, 9.978816725686231 ] }, "model2": { "N0": 2272.727272727273, "class": "SMCModel", "knots": [ 0.00021999999999999998, 0.0004042964845327641, 0.00074298021547978, 0.001365383132708503, 0.0025091800027017288, 0.004611148427964639, 0.008473959541294934, 0.01557269114826553, 0.02861811039071931, 0.05259182465881085, 0.09664859011236712, 0.17761220553399135, 0.32639995594319754, 0.5998288851794722, 1.1023123163603075, 2.0257317925529112, 3.7227101924336217, 6.841266562432742, 12.572272822468664, 12.656215235611679, 40.93397988813272, 74.37947672199772, 115.31345661013056, 168.0865351892005, 242.46601191119817, 369.6185672122657, 21999.999999999996 ], "pid": "ST-2", "spline_class": "CubicSpline", "y": [ 0.2334777864895975, 0.23687880918861937, 0.26030115463356185, 0.31529008122213, 0.3797100901323612, 0.39223061165963297, 0.29894529192081276, 0.2190830388926932, 0.4222918564056228, 0.32023464058056633, -0.36043009906355394, 0.14638909873746686, 1.9072959774249063, 3.7377864707874346, 5.448891926368258, 6.741617622305777, 7.305253415452289, 6.29889134901193, 3.532379931622749, 3.4957224861758593, 4.771948176950006, 3.833081019278895, 4.728818644662528, 5.989604149270964, 7.245898508355312, 8.419742106364035, 9.978816689650431 ] }, "split": 49.00135840123682 }, "rho": 0.005287922137082097, "theta": 0.0001