From e722744689bbda4b4cc5f78cb3ba35ee8a1554c0 Mon Sep 17 00:00:00 2001 From: jaseg Date: Mon, 17 Feb 2020 11:17:03 +0000 Subject: MLE decoder: Make distance fun handle skipped bits --- lab-windows/dsss_experiments.ipynb | 663 ++++--------------------------------- 1 file changed, 71 insertions(+), 592 deletions(-) diff --git a/lab-windows/dsss_experiments.ipynb b/lab-windows/dsss_experiments.ipynb index d11e02e..dbc06e4 100644 --- a/lab-windows/dsss_experiments.ipynb +++ b/lab-windows/dsss_experiments.ipynb @@ -289,7 +289,7 @@ }, { "cell_type": "code", - "execution_count": 20, + "execution_count": 76, "metadata": {}, "outputs": [ { @@ -302,10 +302,18 @@ "(31,) (31,)\n" ] }, + { + "name": "stderr", + "output_type": "stream", + "text": [ + ":27: RuntimeWarning: More than 20 figures have been opened. Figures created through the pyplot interface (`matplotlib.pyplot.figure`) are retained until explicitly closed and may consume too much memory. (To control this warning, see the rcParam `figure.max_open_warning`).\n", + " fig, ((ax1, ax3, ax5), (ax2, ax4, ax6)) = plt.subplots(2, 3, figsize=(16, 9))\n" + ] + }, { "data": { "application/vnd.jupyter.widget-view+json": { - "model_id": "97fba8da8bd74c1aa46c4fabb9e68088", + "model_id": "c24adc84c7294295a47a58aa7d5914f9", "version_major": 2, "version_minor": 0 }, @@ -322,7 +330,7 @@ "(0.001, 0.010294564)" ] }, - "execution_count": 20, + "execution_count": 76, "metadata": {}, "output_type": "execute_result" } @@ -498,13 +506,21 @@ }, { "cell_type": "code", - "execution_count": 24, + "execution_count": 77, "metadata": {}, "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + ":1: RuntimeWarning: More than 20 figures have been opened. Figures created through the pyplot interface (`matplotlib.pyplot.figure`) are retained until explicitly closed and may consume too much memory. (To control this warning, see the rcParam `figure.max_open_warning`).\n", + " fig, ax = plt.subplots()\n" + ] + }, { "data": { "application/vnd.jupyter.widget-view+json": { - "model_id": "61142286f07245fea6b3d080bd6d963b", + "model_id": "298c99e458364426829979eeaf01af6e", "version_major": 2, "version_minor": 0 }, @@ -515,609 +531,72 @@ "metadata": {}, "output_type": "display_data" }, + { + "data": { + "text/plain": [ + "[]" + ] + }, + "execution_count": 77, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "fig, ax = plt.subplots()\n", + "nonlinear_distance = lambda x: 100**(2*np.abs(0.5-x%1)) / (np.abs(x)+7)**2\n", + "x = np.linspace(-1.5, 5.5, 10000)\n", + "ax.plot(x, nonlinear_distance(x))" + ] + }, + { + "cell_type": "code", + "execution_count": 78, + "metadata": {}, + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + ":9: RuntimeWarning: More than 20 figures have been opened. Figures created through the pyplot interface (`matplotlib.pyplot.figure`) are retained until explicitly closed and may consume too much memory. (To control this warning, see the rcParam `figure.max_open_warning`).\n", + " fig, (ax1, ax2, ax3) = plt.subplots(3, figsize=(12, 12))\n" + ] + }, { "data": { "application/vnd.jupyter.widget-view+json": { - "model_id": "591cc5914cfd4812bad96502207bc01f", + "model_id": "126c195977aa4639a8026e0a857a8ab8", "version_major": 2, "version_minor": 0 }, "text/plain": [ - "interactive(children=(FloatSlider(value=10.0, description='w', max=30.0, min=-10.0), Output()), _dom_classes=(…" + "Canvas(toolbar=Toolbar(toolitems=[('Home', 'Reset original view', 'home', 'home'), ('Back', 'Back to previous …" ] }, "metadata": {}, "output_type": "display_data" }, { - "name": "stdout", - "output_type": "stream", - "text": [ - " candidates for 32.5, -0.30796614799228145:\n", - " 7.6945 309.50 -0.51\n", - " 4.9788 419.50 -0.41\n", - " 2.4515 211.00 0.38\n", - " 1.8240 158.50 0.35\n", - " candidates for 158.5, 0.34903632894813835:\n", - " 5.6264 419.50 -0.41\n", - " 2.5562 618.50 0.54\n", - " 2.1462 309.50 -0.51\n", - " 1.3863 211.00 0.38\n", - " candidates for 211.0, 0.3800947813243649:\n", - " 3.9276 618.50 0.54\n", - " 3.0726 419.50 -0.41\n", - " 1.6974 309.50 -0.51\n", - " candidates for 309.5, -0.5053314278026041:\n", - " 29.3614 618.50 0.54\n", - " 1.7265 419.50 -0.41\n", - " candidates for 419.5, -0.40531218065962255:\n", - " 4.0528 823.00 -0.44\n", - " 2.9151 618.50 0.54\n", - "chain candidates:\n", - " [(29.36136019270347, [('309.50', '-0.51'), ('618.50', '0.54')])]\n", - " [(7.69452617901758, [('32.50', '-0.31'), ('309.50', '-0.51')])]\n", - " [(5.626384903889135, [('158.50', '0.35'), ('419.50', '-0.41')])]\n", - " [(4.978775493314884, [('32.50', '-0.31'), ('419.50', '-0.41')])]\n", - " [(4.052840585124125, [('419.50', '-0.41'), ('823.00', '-0.44')])]\n", - " [(3.92759395893727, [('211.00', '0.38'), ('618.50', '0.54')])]\n", - " [(3.0726112472804843, [('211.00', '0.38'), ('419.50', '-0.41')])]\n", - " [(2.9150792009679662, [('419.50', '-0.41'), ('618.50', '0.54')])]\n", - " [(2.5561653875086985, [('158.50', '0.35'), ('618.50', '0.54')])]\n", - " [(2.451531489602112, [('32.50', '-0.31'), ('211.00', '0.38')])]\n", - " [(2.1461568381616574, [('158.50', '0.35'), ('309.50', '-0.51')])]\n", - " [(1.726465614139335, [('309.50', '-0.51'), ('419.50', '-0.41')])]\n", - " [(1.6973609734949815, [('211.00', '0.38'), ('309.50', '-0.51')])]\n", - " candidates for 618.5, 0.540902203588826:\n", - " 35.8652 929.50 -0.62\n", - " 4.0145 1022.50 0.41\n", - " 3.0106 1055.00 0.53\n", - " 2.9911 823.00 -0.44\n", - " candidates for 309.5, -0.5053314278026041:\n", - " 29.3614 618.50 0.54\n", - " 1.7265 419.50 -0.41\n", - " candidates for 419.5, -0.40531218065962255:\n", - " 4.0528 823.00 -0.44\n", - " 2.9151 618.50 0.54\n", - " candidates for 419.5, -0.40531218065962255:\n", - " 4.0528 823.00 -0.44\n", - " 2.9151 618.50 0.54\n", - " candidates for 823.0, -0.4410105115656168:\n", - " 30.6933 1153.50 -0.43\n", - " 3.9002 1055.00 0.53\n", - " 3.5615 1240.50 -0.56\n", - " 2.8605 1022.50 0.41\n", - " 1.8085 929.50 -0.62\n", - " candidates for 618.5, 0.540902203588826:\n", - " 35.8652 929.50 -0.62\n", - " 4.0145 1022.50 0.41\n", - " 3.0106 1055.00 0.53\n", - " 2.9911 823.00 -0.44\n", - " candidates for 419.5, -0.40531218065962255:\n", - " 4.0528 823.00 -0.44\n", - " 2.9151 618.50 0.54\n", - " candidates for 618.5, 0.540902203588826:\n", - " 35.8652 929.50 -0.62\n", - " 4.0145 1022.50 0.41\n", - " 3.0106 1055.00 0.53\n", - " 2.9911 823.00 -0.44\n", - " candidates for 618.5, 0.540902203588826:\n", - " 35.8652 929.50 -0.62\n", - " 4.0145 1022.50 0.41\n", - " 3.0106 1055.00 0.53\n", - " 2.9911 823.00 -0.44\n", - " candidates for 211.0, 0.3800947813243649:\n", - " 3.9276 618.50 0.54\n", - " 3.0726 419.50 -0.41\n", - " 1.6974 309.50 -0.51\n", - "chain candidates:\n", - " [(35.86521493289496, [('419.50', '-0.41'), ('618.50', '0.54'), ('929.50', '-0.62')])]\n", - " [(35.86521493289496, [('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62')])]\n", - " [(35.86521493289496, [('211.00', '0.38'), ('618.50', '0.54'), ('929.50', '-0.62')])]\n", - " [(35.86521493289496, [('158.50', '0.35'), ('618.50', '0.54'), ('929.50', '-0.62')])]\n", - " [(30.693332126158438, [('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43')])]\n", - " [(29.36136019270347, [('32.50', '-0.31'), ('309.50', '-0.51'), ('618.50', '0.54')])]\n", - " [(4.052840585124125, [('211.00', '0.38'), ('419.50', '-0.41'), ('823.00', '-0.44')])]\n", - " [(4.052840585124125, [('158.50', '0.35'), ('419.50', '-0.41'), ('823.00', '-0.44')])]\n", - " [(4.052840585124125, [('32.50', '-0.31'), ('419.50', '-0.41'), ('823.00', '-0.44')])]\n", - " [(4.014470627816998, [('419.50', '-0.41'), ('618.50', '0.54'), ('1022.50', '0.41')])]\n", - " [(4.014470627816998, [('309.50', '-0.51'), ('618.50', '0.54'), ('1022.50', '0.41')])]\n", - " [(4.014470627816998, [('211.00', '0.38'), ('618.50', '0.54'), ('1022.50', '0.41')])]\n", - " [(4.014470627816998, [('158.50', '0.35'), ('618.50', '0.54'), ('1022.50', '0.41')])]\n", - " [(3.92759395893727, [('32.50', '-0.31'), ('211.00', '0.38'), ('618.50', '0.54')])]\n", - " [(3.900204331202999, [('419.50', '-0.41'), ('823.00', '-0.44'), ('1055.00', '0.53')])]\n", - " [(3.5614986168976417, [('419.50', '-0.41'), ('823.00', '-0.44'), ('1240.50', '-0.56')])]\n", - " [(3.0726112472804843, [('32.50', '-0.31'), ('211.00', '0.38'), ('419.50', '-0.41')])]\n", - " [(3.0106218107191913, [('419.50', '-0.41'), ('618.50', '0.54'), ('1055.00', '0.53')])]\n", - " [(3.0106218107191913, [('309.50', '-0.51'), ('618.50', '0.54'), ('1055.00', '0.53')])]\n", - " [(3.0106218107191913, [('211.00', '0.38'), ('618.50', '0.54'), ('1055.00', '0.53')])]\n", - " [(3.0106218107191913, [('158.50', '0.35'), ('618.50', '0.54'), ('1055.00', '0.53')])]\n", - " [(2.9150792009679662, [('211.00', '0.38'), ('419.50', '-0.41'), ('618.50', '0.54')])]\n", - " [(2.9150792009679662, [('158.50', '0.35'), ('419.50', '-0.41'), ('618.50', '0.54')])]\n", - " [(2.9150792009679662, [('32.50', '-0.31'), ('419.50', '-0.41'), ('618.50', '0.54')])]\n", - " [(1.726465614139335, [('32.50', '-0.31'), ('309.50', '-0.51'), ('419.50', '-0.41')])]\n", - " [(1.6973609734949815, [('32.50', '-0.31'), ('211.00', '0.38'), ('309.50', '-0.51')])]\n", - " candidates for 929.5, -0.6193187546786915:\n", - " 35.8350 1240.50 -0.56\n", - " 3.5505 1153.50 -0.43\n", - " 3.0349 1366.00 0.58\n", - " 1.9091 1055.00 0.53\n", - " 1.6146 1022.50 0.41\n", - " candidates for 929.5, -0.6193187546786915:\n", - " 35.8350 1240.50 -0.56\n", - " 3.5505 1153.50 -0.43\n", - " 3.0349 1366.00 0.58\n", - " 1.9091 1055.00 0.53\n", - " 1.6146 1022.50 0.41\n", - " candidates for 929.5, -0.6193187546786915:\n", - " 35.8350 1240.50 -0.56\n", - " 3.5505 1153.50 -0.43\n", - " 3.0349 1366.00 0.58\n", - " 1.9091 1055.00 0.53\n", - " 1.6146 1022.50 0.41\n", - " candidates for 929.5, -0.6193187546786915:\n", - " 35.8350 1240.50 -0.56\n", - " 3.5505 1153.50 -0.43\n", - " 3.0349 1366.00 0.58\n", - " 1.9091 1055.00 0.53\n", - " 1.6146 1022.50 0.41\n", - " candidates for 1153.5, -0.4342832999357559:\n", - " 6.9736 1426.00 -0.47\n", - " 4.5113 1549.50 0.60\n", - " 3.2649 1366.00 0.58\n", - " 1.6528 1240.50 -0.56\n", - " candidates for 618.5, 0.540902203588826:\n", - " 35.8652 929.50 -0.62\n", - " 4.0145 1022.50 0.41\n", - " 3.0106 1055.00 0.53\n", - " 2.9911 823.00 -0.44\n", - " candidates for 823.0, -0.4410105115656168:\n", - " 30.6933 1153.50 -0.43\n", - " 3.9002 1055.00 0.53\n", - " 3.5615 1240.50 -0.56\n", - " 2.8605 1022.50 0.41\n", - " 1.8085 929.50 -0.62\n", - " candidates for 823.0, -0.4410105115656168:\n", - " 30.6933 1153.50 -0.43\n", - " 3.9002 1055.00 0.53\n", - " 3.5615 1240.50 -0.56\n", - " 2.8605 1022.50 0.41\n", - " 1.8085 929.50 -0.62\n", - " candidates for 823.0, -0.4410105115656168:\n", - " 30.6933 1153.50 -0.43\n", - " 3.9002 1055.00 0.53\n", - " 3.5615 1240.50 -0.56\n", - " 2.8605 1022.50 0.41\n", - " 1.8085 929.50 -0.62\n", - " candidates for 1022.5, 0.4098936365863757:\n", - " 13.9051 1366.00 0.58\n", - " 4.0691 1426.00 -0.47\n", - " 3.4167 1240.50 -0.56\n", - " 1.9103 1153.50 -0.43\n", - " 1.3769 1055.00 0.53\n", - "chain candidates:\n", - " [(35.86521493289496, [('32.50', '-0.31'), ('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62')])]\n", - " [(35.835002890401974, [('419.50', '-0.41'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56')])]\n", - " [(35.835002890401974, [('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56')])]\n", - " [(35.835002890401974, [('211.00', '0.38'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56')])]\n", - " [(35.835002890401974, [('158.50', '0.35'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56')])]\n", - " [(30.693332126158438, [('211.00', '0.38'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43')])]\n", - " [(30.693332126158438, [('158.50', '0.35'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43')])]\n", - " [(30.693332126158438, [('32.50', '-0.31'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43')])]\n", - " [(13.905133395544146, [('419.50', '-0.41'), ('618.50', '0.54'), ('1022.50', '0.41'), ('1366.00', '0.58')])]\n", - " [(6.9736079278832195, [('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47')])]\n", - " [(4.511327725367169, [('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1549.50', '0.60')])]\n", - " [(4.069101151961379, [('419.50', '-0.41'), ('618.50', '0.54'), ('1022.50', '0.41'), ('1426.00', '-0.47')])]\n", - " [(4.014470627816998, [('32.50', '-0.31'), ('309.50', '-0.51'), ('618.50', '0.54'), ('1022.50', '0.41')])]\n", - " [(3.900204331202999, [('211.00', '0.38'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1055.00', '0.53')])]\n", - " [(3.900204331202999, [('158.50', '0.35'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1055.00', '0.53')])]\n", - " [(3.900204331202999, [('32.50', '-0.31'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1055.00', '0.53')])]\n", - " [(3.5614986168976417, [('211.00', '0.38'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1240.50', '-0.56')])]\n", - " [(3.5614986168976417, [('158.50', '0.35'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1240.50', '-0.56')])]\n", - " [(3.5614986168976417, [('32.50', '-0.31'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1240.50', '-0.56')])]\n", - " [(3.550474983301211, [('419.50', '-0.41'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1153.50', '-0.43')])]\n", - " [(3.550474983301211, [('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1153.50', '-0.43')])]\n", - " [(3.550474983301211, [('211.00', '0.38'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1153.50', '-0.43')])]\n", - " [(3.550474983301211, [('158.50', '0.35'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1153.50', '-0.43')])]\n", - " [(3.4167022368071445, [('419.50', '-0.41'), ('618.50', '0.54'), ('1022.50', '0.41'), ('1240.50', '-0.56')])]\n", - " [(3.2648563049949213, [('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1366.00', '0.58')])]\n", - " [(3.0348932348282385, [('419.50', '-0.41'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1366.00', '0.58')])]\n", - " [(3.0348932348282385, [('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1366.00', '0.58')])]\n", - " [(3.0348932348282385, [('211.00', '0.38'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1366.00', '0.58')])]\n", - " [(3.0348932348282385, [('158.50', '0.35'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1366.00', '0.58')])]\n", - " [(3.0106218107191913, [('32.50', '-0.31'), ('309.50', '-0.51'), ('618.50', '0.54'), ('1055.00', '0.53')])]\n", - " candidates for 929.5, -0.6193187546786915:\n", - " 35.8350 1240.50 -0.56\n", - " 3.5505 1153.50 -0.43\n", - " 3.0349 1366.00 0.58\n", - " 1.9091 1055.00 0.53\n", - " 1.6146 1022.50 0.41\n", - " candidates for 1240.5, -0.5588946696927191:\n", - " 29.3917 1549.50 0.60\n", - " 2.6159 1426.00 -0.47\n", - " 2.2893 1716.00 0.46\n", - " 1.9334 1366.00 0.58\n", - " candidates for 1240.5, -0.5588946696927191:\n", - " 29.3917 1549.50 0.60\n", - " 2.6159 1426.00 -0.47\n", - " 2.2893 1716.00 0.46\n", - " 1.9334 1366.00 0.58\n", - " candidates for 1240.5, -0.5588946696927191:\n", - " 29.3917 1549.50 0.60\n", - " 2.6159 1426.00 -0.47\n", - " 2.2893 1716.00 0.46\n", - " 1.9334 1366.00 0.58\n", - " candidates for 1240.5, -0.5588946696927191:\n", - " 29.3917 1549.50 0.60\n", - " 2.6159 1426.00 -0.47\n", - " 2.2893 1716.00 0.46\n", - " 1.9334 1366.00 0.58\n", - " candidates for 1153.5, -0.4342832999357559:\n", - " 6.9736 1426.00 -0.47\n", - " 4.5113 1549.50 0.60\n", - " 3.2649 1366.00 0.58\n", - " 1.6528 1240.50 -0.56\n", - " candidates for 1153.5, -0.4342832999357559:\n", - " 6.9736 1426.00 -0.47\n", - " 4.5113 1549.50 0.60\n", - " 3.2649 1366.00 0.58\n", - " 1.6528 1240.50 -0.56\n", - " candidates for 1153.5, -0.4342832999357559:\n", - " 6.9736 1426.00 -0.47\n", - " 4.5113 1549.50 0.60\n", - " 3.2649 1366.00 0.58\n", - " 1.6528 1240.50 -0.56\n", - " candidates for 1366.0, 0.5762242378968192:\n", - " 10.8981 1716.00 0.46\n", - " 2.6451 1549.50 0.60\n", - " 1.4675 1426.00 -0.47\n", - " candidates for 1426.0, -0.47353164524012403:\n", - " 10.8981 1716.00 0.46\n", - " 3.0884 1860.50 0.59\n", - " 1.9293 1549.50 0.60\n", - "chain candidates:\n", - " [(35.835002890401974, [('32.50', '-0.31'), ('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56')])]\n", - " [(29.391710500486752, [('419.50', '-0.41'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60')])]\n", - " [(29.391710500486752, [('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60')])]\n", - " [(29.391710500486752, [('211.00', '0.38'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60')])]\n", - " [(29.391710500486752, [('158.50', '0.35'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60')])]\n", - " [(10.8980671118677, [('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46')])]\n", - " [(10.8980671118677, [('419.50', '-0.41'), ('618.50', '0.54'), ('1022.50', '0.41'), ('1366.00', '0.58'), ('1716.00', '0.46')])]\n", - " [(6.9736079278832195, [('211.00', '0.38'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47')])]\n", - " [(6.9736079278832195, [('158.50', '0.35'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47')])]\n", - " [(6.9736079278832195, [('32.50', '-0.31'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47')])]\n", - " [(4.511327725367169, [('211.00', '0.38'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1549.50', '0.60')])]\n", - " [(4.511327725367169, [('158.50', '0.35'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1549.50', '0.60')])]\n", - " [(4.511327725367169, [('32.50', '-0.31'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1549.50', '0.60')])]\n", - " [(3.550474983301211, [('32.50', '-0.31'), ('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1153.50', '-0.43')])]\n", - " [(3.2648563049949213, [('211.00', '0.38'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1366.00', '0.58')])]\n", - " [(3.2648563049949213, [('158.50', '0.35'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1366.00', '0.58')])]\n", - " [(3.2648563049949213, [('32.50', '-0.31'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1366.00', '0.58')])]\n", - " [(3.088380538661272, [('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1860.50', '0.59')])]\n", - " [(3.0348932348282385, [('32.50', '-0.31'), ('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1366.00', '0.58')])]\n", - " [(2.645123753900041, [('419.50', '-0.41'), ('618.50', '0.54'), ('1022.50', '0.41'), ('1366.00', '0.58'), ('1549.50', '0.60')])]\n", - " [(2.615947978753891, [('419.50', '-0.41'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1426.00', '-0.47')])]\n", - " [(2.615947978753891, [('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1426.00', '-0.47')])]\n", - " [(2.615947978753891, [('211.00', '0.38'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1426.00', '-0.47')])]\n", - " [(2.615947978753891, [('158.50', '0.35'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1426.00', '-0.47')])]\n", - " [(2.289278258705857, [('419.50', '-0.41'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1716.00', '0.46')])]\n", - " [(2.289278258705857, [('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1716.00', '0.46')])]\n", - " [(2.289278258705857, [('211.00', '0.38'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1716.00', '0.46')])]\n", - " [(2.289278258705857, [('158.50', '0.35'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1716.00', '0.46')])]\n", - " [(1.9293001373130654, [('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1549.50', '0.60')])]\n", - " [(1.4675350533892928, [('419.50', '-0.41'), ('618.50', '0.54'), ('1022.50', '0.41'), ('1366.00', '0.58'), ('1426.00', '-0.47')])]\n", - " candidates for 1240.5, -0.5588946696927191:\n", - " 29.3917 1549.50 0.60\n", - " 2.6159 1426.00 -0.47\n", - " 2.2893 1716.00 0.46\n", - " 1.9334 1366.00 0.58\n", - " candidates for 1549.5, 0.6016028191553928:\n", - " 35.8492 1860.50 0.59\n", - " 3.3339 1972.50 0.45\n", - " 2.3161 1716.00 0.46\n", - " candidates for 1549.5, 0.6016028191553928:\n", - " 35.8492 1860.50 0.59\n", - " 3.3339 1972.50 0.45\n", - " 2.3161 1716.00 0.46\n", - " candidates for 1549.5, 0.6016028191553928:\n", - " 35.8492 1860.50 0.59\n", - " 3.3339 1972.50 0.45\n", - " 2.3161 1716.00 0.46\n", - " candidates for 1549.5, 0.6016028191553928:\n", - " 35.8492 1860.50 0.59\n", - " 3.3339 1972.50 0.45\n", - " 2.3161 1716.00 0.46\n", - " candidates for 1716.0, 0.462800890402067:\n", - " 6.1740 2089.50 0.39\n", - " 5.2664 1972.50 0.45\n", - " 2.8998 2156.50 -0.49\n", - " 2.1170 1860.50 0.59\n", - " candidates for 1716.0, 0.462800890402067:\n", - " 6.1740 2089.50 0.39\n", - " 5.2664 1972.50 0.45\n", - " 2.8998 2156.50 -0.49\n", - " 2.1170 1860.50 0.59\n", - " candidates for 1426.0, -0.47353164524012403:\n", - " 10.8981 1716.00 0.46\n", - " 3.0884 1860.50 0.59\n", - " 1.9293 1549.50 0.60\n", - " candidates for 1426.0, -0.47353164524012403:\n", - " 10.8981 1716.00 0.46\n", - " 3.0884 1860.50 0.59\n", - " 1.9293 1549.50 0.60\n", - " candidates for 1426.0, -0.47353164524012403:\n", - " 10.8981 1716.00 0.46\n", - " 3.0884 1860.50 0.59\n", - " 1.9293 1549.50 0.60\n", - "chain candidates:\n", - " [(35.84917626888937, [('419.50', '-0.41'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59')])]\n", - " [(35.84917626888937, [('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59')])]\n", - " [(35.84917626888937, [('211.00', '0.38'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59')])]\n", - " [(35.84917626888937, [('158.50', '0.35'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59')])]\n", - " [(29.391710500486752, [('32.50', '-0.31'), ('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60')])]\n", - " [(10.8980671118677, [('211.00', '0.38'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46')])]\n", - " [(10.8980671118677, [('158.50', '0.35'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46')])]\n", - " [(10.8980671118677, [('32.50', '-0.31'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46')])]\n", - " [(6.173953539976651, [('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('2089.50', '0.39')])]\n", - " [(6.173953539976651, [('419.50', '-0.41'), ('618.50', '0.54'), ('1022.50', '0.41'), ('1366.00', '0.58'), ('1716.00', '0.46'), ('2089.50', '0.39')])]\n", - " [(5.266444098570901, [('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('1972.50', '0.45')])]\n", - " [(5.266444098570901, [('419.50', '-0.41'), ('618.50', '0.54'), ('1022.50', '0.41'), ('1366.00', '0.58'), ('1716.00', '0.46'), ('1972.50', '0.45')])]\n", - " [(3.333870136335599, [('419.50', '-0.41'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1972.50', '0.45')])]\n", - " [(3.333870136335599, [('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1972.50', '0.45')])]\n", - " [(3.333870136335599, [('211.00', '0.38'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1972.50', '0.45')])]\n", - " [(3.333870136335599, [('158.50', '0.35'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1972.50', '0.45')])]\n", - " [(3.088380538661272, [('211.00', '0.38'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1860.50', '0.59')])]\n", - " [(3.088380538661272, [('158.50', '0.35'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1860.50', '0.59')])]\n", - " [(3.088380538661272, [('32.50', '-0.31'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1860.50', '0.59')])]\n", - " [(2.8998303227280857, [('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('2156.50', '-0.49')])]\n", - " [(2.8998303227280857, [('419.50', '-0.41'), ('618.50', '0.54'), ('1022.50', '0.41'), ('1366.00', '0.58'), ('1716.00', '0.46'), ('2156.50', '-0.49')])]\n", - " [(2.615947978753891, [('32.50', '-0.31'), ('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1426.00', '-0.47')])]\n", - " [(2.3160909989469616, [('419.50', '-0.41'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1716.00', '0.46')])]\n", - " [(2.3160909989469616, [('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1716.00', '0.46')])]\n", - " [(2.3160909989469616, [('211.00', '0.38'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1716.00', '0.46')])]\n", - " [(2.3160909989469616, [('158.50', '0.35'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1716.00', '0.46')])]\n", - " [(2.289278258705857, [('32.50', '-0.31'), ('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1716.00', '0.46')])]\n", - " [(1.9293001373130654, [('211.00', '0.38'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1549.50', '0.60')])]\n", - " [(1.9293001373130654, [('158.50', '0.35'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1549.50', '0.60')])]\n", - " [(1.9293001373130654, [('32.50', '-0.31'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1549.50', '0.60')])]\n", - " candidates for 1860.5, 0.587241426667523:\n", - " 13.5776 2156.50 -0.49\n", - " 4.6432 2252.00 0.34\n", - " 3.7091 2089.50 0.39\n", - " 1.7655 1972.50 0.45\n", - " candidates for 1860.5, 0.587241426667523:\n", - " 13.5776 2156.50 -0.49\n", - " 4.6432 2252.00 0.34\n", - " 3.7091 2089.50 0.39\n", - " 1.7655 1972.50 0.45\n", - " candidates for 1860.5, 0.587241426667523:\n", - " 13.5776 2156.50 -0.49\n", - " 4.6432 2252.00 0.34\n", - " 3.7091 2089.50 0.39\n", - " 1.7655 1972.50 0.45\n", - " candidates for 1860.5, 0.587241426667523:\n", - " 13.5776 2156.50 -0.49\n", - " 4.6432 2252.00 0.34\n", - " 3.7091 2089.50 0.39\n", - " 1.7655 1972.50 0.45\n", - " candidates for 1549.5, 0.6016028191553928:\n", - " 35.8492 1860.50 0.59\n", - " 3.3339 1972.50 0.45\n", - " 2.3161 1716.00 0.46\n", - " candidates for 1716.0, 0.462800890402067:\n", - " 6.1740 2089.50 0.39\n", - " 5.2664 1972.50 0.45\n", - " 2.8998 2156.50 -0.49\n", - " 2.1170 1860.50 0.59\n", - " candidates for 1716.0, 0.462800890402067:\n", - " 6.1740 2089.50 0.39\n", - " 5.2664 1972.50 0.45\n", - " 2.8998 2156.50 -0.49\n", - " 2.1170 1860.50 0.59\n", - " candidates for 1716.0, 0.462800890402067:\n", - " 6.1740 2089.50 0.39\n", - " 5.2664 1972.50 0.45\n", - " 2.8998 2156.50 -0.49\n", - " 2.1170 1860.50 0.59\n", - " candidates for 2089.5, 0.3852902575233996:\n", - " 8.3442 2370.50 -0.28\n", - " 4.7167 2480.50 -0.42\n", - " 2.1995 2252.00 0.34\n", - " 1.5091 2156.50 -0.49\n", - " candidates for 2089.5, 0.3852902575233996:\n", - " 8.3442 2370.50 -0.28\n", - " 4.7167 2480.50 -0.42\n", - " 2.1995 2252.00 0.34\n", - " 1.5091 2156.50 -0.49\n", - "chain candidates:\n", - " [(35.84917626888937, [('32.50', '-0.31'), ('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59')])]\n", - " [(13.577561996310388, [('419.50', '-0.41'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2156.50', '-0.49')])]\n", - " [(13.577561996310388, [('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2156.50', '-0.49')])]\n", - " [(13.577561996310388, [('211.00', '0.38'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2156.50', '-0.49')])]\n", - " [(13.577561996310388, [('158.50', '0.35'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2156.50', '-0.49')])]\n", - " [(8.34422970205904, [('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('2089.50', '0.39'), ('2370.50', '-0.28')])]\n", - " [(8.34422970205904, [('419.50', '-0.41'), ('618.50', '0.54'), ('1022.50', '0.41'), ('1366.00', '0.58'), ('1716.00', '0.46'), ('2089.50', '0.39'), ('2370.50', '-0.28')])]\n", - " [(6.173953539976651, [('211.00', '0.38'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('2089.50', '0.39')])]\n", - " [(6.173953539976651, [('158.50', '0.35'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('2089.50', '0.39')])]\n", - " [(6.173953539976651, [('32.50', '-0.31'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('2089.50', '0.39')])]\n", - " [(5.266444098570901, [('211.00', '0.38'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('1972.50', '0.45')])]\n", - " [(5.266444098570901, [('158.50', '0.35'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('1972.50', '0.45')])]\n", - " [(5.266444098570901, [('32.50', '-0.31'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('1972.50', '0.45')])]\n", - " [(4.716677612654047, [('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('2089.50', '0.39'), ('2480.50', '-0.42')])]\n", - " [(4.716677612654047, [('419.50', '-0.41'), ('618.50', '0.54'), ('1022.50', '0.41'), ('1366.00', '0.58'), ('1716.00', '0.46'), ('2089.50', '0.39'), ('2480.50', '-0.42')])]\n", - " [(4.643240934821138, [('419.50', '-0.41'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2252.00', '0.34')])]\n", - " [(4.643240934821138, [('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2252.00', '0.34')])]\n", - " [(4.643240934821138, [('211.00', '0.38'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2252.00', '0.34')])]\n", - " [(4.643240934821138, [('158.50', '0.35'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2252.00', '0.34')])]\n", - " [(3.7091286452452157, [('419.50', '-0.41'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2089.50', '0.39')])]\n", - " [(3.7091286452452157, [('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2089.50', '0.39')])]\n", - " [(3.7091286452452157, [('211.00', '0.38'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2089.50', '0.39')])]\n", - " [(3.7091286452452157, [('158.50', '0.35'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2089.50', '0.39')])]\n", - " [(3.333870136335599, [('32.50', '-0.31'), ('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1972.50', '0.45')])]\n", - " [(2.8998303227280857, [('211.00', '0.38'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('2156.50', '-0.49')])]\n", - " [(2.8998303227280857, [('158.50', '0.35'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('2156.50', '-0.49')])]\n", - " [(2.8998303227280857, [('32.50', '-0.31'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('2156.50', '-0.49')])]\n", - " [(2.3160909989469616, [('32.50', '-0.31'), ('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1716.00', '0.46')])]\n", - " [(2.1994624910426936, [('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('2089.50', '0.39'), ('2252.00', '0.34')])]\n", - " [(2.1994624910426936, [('419.50', '-0.41'), ('618.50', '0.54'), ('1022.50', '0.41'), ('1366.00', '0.58'), ('1716.00', '0.46'), ('2089.50', '0.39'), ('2252.00', '0.34')])]\n", - " candidates for 1860.5, 0.587241426667523:\n", - " 13.5776 2156.50 -0.49\n", - " 4.6432 2252.00 0.34\n", - " 3.7091 2089.50 0.39\n", - " 1.7655 1972.50 0.45\n", - " candidates for 2156.5, -0.4884573259540951:\n", - " 80.2096 2480.50 -0.42\n", - " 3.1580 2370.50 -0.28\n", - " 2.6215 2605.50 0.28\n", - " 1.5931 2252.00 0.34\n", - " candidates for 2156.5, -0.4884573259540951:\n", - " 80.2096 2480.50 -0.42\n", - " 3.1580 2370.50 -0.28\n", - " 2.6215 2605.50 0.28\n", - " 1.5931 2252.00 0.34\n", - " candidates for 2156.5, -0.4884573259540951:\n", - " 80.2096 2480.50 -0.42\n", - " 3.1580 2370.50 -0.28\n", - " 2.6215 2605.50 0.28\n", - " 1.5931 2252.00 0.34\n", - " candidates for 2156.5, -0.4884573259540951:\n", - " 80.2096 2480.50 -0.42\n", - " 3.1580 2370.50 -0.28\n", - " 2.6215 2605.50 0.28\n", - " 1.5931 2252.00 0.34\n", - " candidates for 2370.5, -0.2782029938616636:\n", - " 8.4892 2728.50 -0.14\n", - " 8.1897 2651.00 -0.18\n", - " 4.5779 2761.50 0.14\n", - " 3.9056 2605.50 0.28\n", - " 1.7334 2480.50 -0.42\n", - " candidates for 2370.5, -0.2782029938616636:\n", - " 8.4892 2728.50 -0.14\n", - " 8.1897 2651.00 -0.18\n", - " 4.5779 2761.50 0.14\n", - " 3.9056 2605.50 0.28\n", - " 1.7334 2480.50 -0.42\n", - " candidates for 2089.5, 0.3852902575233996:\n", - " 8.3442 2370.50 -0.28\n", - " 4.7167 2480.50 -0.42\n", - " 2.1995 2252.00 0.34\n", - " 1.5091 2156.50 -0.49\n", - " candidates for 2089.5, 0.3852902575233996:\n", - " 8.3442 2370.50 -0.28\n", - " 4.7167 2480.50 -0.42\n", - " 2.1995 2252.00 0.34\n", - " 1.5091 2156.50 -0.49\n", - " candidates for 2089.5, 0.3852902575233996:\n", - " 8.3442 2370.50 -0.28\n", - " 4.7167 2480.50 -0.42\n", - " 2.1995 2252.00 0.34\n", - " 1.5091 2156.50 -0.49\n", - "chain candidates:\n", - " [(80.2096353591332, [('419.50', '-0.41'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2156.50', '-0.49'), ('2480.50', '-0.42')])]\n", - " [(80.2096353591332, [('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2156.50', '-0.49'), ('2480.50', '-0.42')])]\n", - " [(80.2096353591332, [('211.00', '0.38'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2156.50', '-0.49'), ('2480.50', '-0.42')])]\n", - " [(80.2096353591332, [('158.50', '0.35'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2156.50', '-0.49'), ('2480.50', '-0.42')])]\n", - " [(13.577561996310388, [('32.50', '-0.31'), ('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2156.50', '-0.49')])]\n", - " [(8.34422970205904, [('211.00', '0.38'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('2089.50', '0.39'), ('2370.50', '-0.28')])]\n", - " [(8.34422970205904, [('158.50', '0.35'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('2089.50', '0.39'), ('2370.50', '-0.28')])]\n", - " [(8.34422970205904, [('32.50', '-0.31'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('2089.50', '0.39'), ('2370.50', '-0.28')])]\n", - " [(4.716677612654047, [('211.00', '0.38'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('2089.50', '0.39'), ('2480.50', '-0.42')])]\n", - " [(4.716677612654047, [('158.50', '0.35'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('2089.50', '0.39'), ('2480.50', '-0.42')])]\n", - " [(4.716677612654047, [('32.50', '-0.31'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('2089.50', '0.39'), ('2480.50', '-0.42')])]\n", - " [(4.643240934821138, [('32.50', '-0.31'), ('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2252.00', '0.34')])]\n", - " [(3.7091286452452157, [('32.50', '-0.31'), ('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2089.50', '0.39')])]\n", - " [(3.1579694214591334, [('419.50', '-0.41'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2156.50', '-0.49'), ('2370.50', '-0.28')])]\n", - " [(3.1579694214591334, [('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2156.50', '-0.49'), ('2370.50', '-0.28')])]\n", - " [(3.1579694214591334, [('211.00', '0.38'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2156.50', '-0.49'), ('2370.50', '-0.28')])]\n", - " [(3.1579694214591334, [('158.50', '0.35'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2156.50', '-0.49'), ('2370.50', '-0.28')])]\n", - " [(2.6214860493499117, [('419.50', '-0.41'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2156.50', '-0.49'), ('2605.50', '0.28')])]\n", - " [(2.6214860493499117, [('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2156.50', '-0.49'), ('2605.50', '0.28')])]\n", - " [(2.6214860493499117, [('211.00', '0.38'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2156.50', '-0.49'), ('2605.50', '0.28')])]\n", - " [(2.6214860493499117, [('158.50', '0.35'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2156.50', '-0.49'), ('2605.50', '0.28')])]\n", - " [(2.1994624910426936, [('211.00', '0.38'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('2089.50', '0.39'), ('2252.00', '0.34')])]\n", - " [(2.1994624910426936, [('158.50', '0.35'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('2089.50', '0.39'), ('2252.00', '0.34')])]\n", - " [(2.1994624910426936, [('32.50', '-0.31'), ('419.50', '-0.41'), ('823.00', '-0.44'), ('1153.50', '-0.43'), ('1426.00', '-0.47'), ('1716.00', '0.46'), ('2089.50', '0.39'), ('2252.00', '0.34')])]\n", - " candidates for 2480.5, -0.41927071826584167:\n", - " 8.2759 2761.50 0.14\n", - " 4.5126 2728.50 -0.14\n", - " 2.2289 2651.00 -0.18\n", - " 1.7819 2605.50 0.28\n", - " candidates for 2480.5, -0.41927071826584167:\n", - " 8.2759 2761.50 0.14\n", - " 4.5126 2728.50 -0.14\n", - " 2.2289 2651.00 -0.18\n", - " 1.7819 2605.50 0.28\n", - " candidates for 2480.5, -0.41927071826584167:\n", - " 8.2759 2761.50 0.14\n", - " 4.5126 2728.50 -0.14\n", - " 2.2289 2651.00 -0.18\n", - " 1.7819 2605.50 0.28\n", - " candidates for 2480.5, -0.41927071826584167:\n", - " 8.2759 2761.50 0.14\n", - " 4.5126 2728.50 -0.14\n", - " 2.2289 2651.00 -0.18\n", - " 1.7819 2605.50 0.28\n", - " candidates for 2156.5, -0.4884573259540951:\n", - " 80.2096 2480.50 -0.42\n", - " 3.1580 2370.50 -0.28\n", - " 2.6215 2605.50 0.28\n", - " 1.5931 2252.00 0.34\n", - " candidates for 2370.5, -0.2782029938616636:\n", - " 8.4892 2728.50 -0.14\n", - " 8.1897 2651.00 -0.18\n", - " 4.5779 2761.50 0.14\n", - " 3.9056 2605.50 0.28\n", - " 1.7334 2480.50 -0.42\n", - " candidates for 2370.5, -0.2782029938616636:\n", - " 8.4892 2728.50 -0.14\n", - " 8.1897 2651.00 -0.18\n", - " 4.5779 2761.50 0.14\n", - " 3.9056 2605.50 0.28\n", - " 1.7334 2480.50 -0.42\n", - " candidates for 2370.5, -0.2782029938616636:\n", - " 8.4892 2728.50 -0.14\n", - " 8.1897 2651.00 -0.18\n", - " 4.5779 2761.50 0.14\n", - " 3.9056 2605.50 0.28\n", - " 1.7334 2480.50 -0.42\n", - " candidates for 2480.5, -0.41927071826584167:\n", - " 8.2759 2761.50 0.14\n", - " 4.5126 2728.50 -0.14\n", - " 2.2289 2651.00 -0.18\n", - " 1.7819 2605.50 0.28\n", - " candidates for 2480.5, -0.41927071826584167:\n", - " 8.2759 2761.50 0.14\n", - " 4.5126 2728.50 -0.14\n", - " 2.2289 2651.00 -0.18\n", - " 1.7819 2605.50 0.28\n", - "chain candidates:\n", - " [(80.2096353591332, [('32.50', '-0.31'), ('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2156.50', '-0.49'), ('2480.50', '-0.42')])]\n", - " [(3.1579694214591334, [('32.50', '-0.31'), ('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2156.50', '-0.49'), ('2370.50', '-0.28')])]\n", - " [(2.6214860493499117, [('32.50', '-0.31'), ('309.50', '-0.51'), ('618.50', '0.54'), ('929.50', '-0.62'), ('1240.50', '-0.56'), ('1549.50', '0.60'), ('1860.50', '0.59'), ('2156.50', '-0.49'), ('2605.50', '0.28')])]\n", - " candidates for 2480.5, -0.41927071826584167:\n", - " 8.2759 2761.50 0.14\n", - " 4.5126 2728.50 -0.14\n", - " 2.2289 2651.00 -0.18\n", - " 1.7819 2605.50 0.28\n", - " candidates for 2370.5, -0.2782029938616636:\n", - " 8.4892 2728.50 -0.14\n", - " 8.1897 2651.00 -0.18\n", - " 4.5779 2761.50 0.14\n", - " 3.9056 2605.50 0.28\n", - " 1.7334 2480.50 -0.42\n", - " candidates for 2605.5, 0.2817317886223028:\n", - " 2.0220 2761.50 0.14\n", - " 1.6925 2728.50 -0.14\n", - " 1.2542 2651.00 -0.18\n", - "chain candidates:\n" - ] + "data": { + "application/vnd.jupyter.widget-view+json": { + "model_id": "5a40fd3c41814b5f99e9f452c8923db4", + "version_major": 2, + "version_minor": 0 + }, + "text/plain": [ + "interactive(children=(FloatSlider(value=10.0, description='w', max=30.0, min=-10.0), Output()), _dom_classes=(…" + ] + }, + "metadata": {}, + "output_type": "display_data" }, { "data": { "text/plain": [ - "[]" + "[]" ] }, - "execution_count": 24, + "execution_count": 78, "metadata": {}, "output_type": "execute_result" } @@ -1190,8 +669,8 @@ " chain_candidates = []\n", " for chain in candidates:\n", " pos, ampl = chain[-1]\n", - " score_fun = lambda pos, npos, npol: abs(npol)/2 + 1/(abs((npos-pos)/bit_period - 1))\n", - " next_candidates = sorted([ (score_fun(pos, npos, npol), npos, npol) for npos, npol in peak_groups if pos < npos < pos + bit_period*1.5 ], reverse=True)\n", + " score_fun = lambda pos, npos, npol: abs(npol)/2 + nonlinear_distance((npos-pos)/bit_period)\n", + " next_candidates = sorted([ (score_fun(pos, npos, npol), npos, npol) for npos, npol in peak_groups if pos < npos < pos + bit_period*3.5 ], reverse=True)\n", " \n", " print(f' candidates for {pos}, {ampl}:')\n", " for score, npos, npol in next_candidates:\n", @@ -1210,7 +689,7 @@ " print(' ', [(score, [(f'{pos:.2f}', f'{pol:.2f}') for pos, pol in chain])])\n", " candidates = [ chain for _score, chain in sorted(chain_candidates, reverse=True)[:10] ]\n", "\n", - "res = sorted(mle_decode(peak_group), reverse=True)\n", + "res = sorted(mle_decode(peak_group, print=lambda *args, **kwargs: None), reverse=True)\n", "#for i, (score, chain) in enumerate(res):\n", "# print(f'Chain {i}@{score:.4f}: {chain}')\n", "(_score, chain), *_ = res\n", -- cgit