Note
Click here to download the full example code
TRAM on a 1D double well¶
This example shows how to use the transition-based reweighting analysis method (TRAM) to estimate the free energies and Markov model of a simple double-well potential, sampled using umbrella sampling.
For more information see the TRAM
estimator and
its respective TRAM tutorial.
0%| | 0/1000 [00:00<?, ?it/s]
Initializing free energies using MBAR: 0%| | 0/1000 [00:00<?, ?it/s]
Initializing free energies using MBAR - [inc: 2.9e-01]: 0%| | 1/1000 [00:00<00:38, 25.98it/s]
Initializing free energies using MBAR - [inc: 7.1e-02]: 0%| | 2/1000 [00:00<00:27, 35.84it/s]
Initializing free energies using MBAR - [inc: 2.5e-02]: 0%| | 3/1000 [00:00<00:22, 44.21it/s]
Initializing free energies using MBAR - [inc: 1.6e-02]: 0%| | 4/1000 [00:00<00:19, 49.94it/s]
Initializing free energies using MBAR - [inc: 1.2e-02]: 0%| | 5/1000 [00:00<00:18, 53.85it/s]
Initializing free energies using MBAR - [inc: 1.2e-02]: 1%| | 6/1000 [00:00<00:17, 57.08it/s]
Initializing free energies using MBAR - [inc: 9.0e-03]: 1%| | 6/1000 [00:00<00:17, 57.08it/s]
Initializing free energies using MBAR - [inc: 6.7e-03]: 1%| | 7/1000 [00:00<00:17, 57.08it/s]
Initializing free energies using MBAR - [inc: 5.0e-03]: 1%| | 8/1000 [00:00<00:17, 57.08it/s]
Initializing free energies using MBAR - [inc: 3.7e-03]: 1%| | 9/1000 [00:00<00:17, 57.08it/s]
Initializing free energies using MBAR - [inc: 2.7e-03]: 1%|1 | 10/1000 [00:00<00:17, 57.08it/s]
Initializing free energies using MBAR - [inc: 2.0e-03]: 1%|1 | 11/1000 [00:00<00:17, 57.08it/s]
Initializing free energies using MBAR - [inc: 1.5e-03]: 1%|1 | 12/1000 [00:00<00:17, 57.08it/s]
Initializing free energies using MBAR - [inc: 1.1e-03]: 1%|1 | 13/1000 [00:00<00:17, 57.08it/s]
Initializing free energies using MBAR - [inc: 8.3e-04]: 1%|1 | 14/1000 [00:00<00:17, 57.08it/s]
Initializing free energies using MBAR - [inc: 6.2e-04]: 2%|1 | 15/1000 [00:00<00:17, 57.08it/s]
Initializing free energies using MBAR - [inc: 6.2e-04]: 2%|1 | 16/1000 [00:00<00:12, 78.67it/s]
Initializing free energies using MBAR - [inc: 4.6e-04]: 2%|1 | 16/1000 [00:00<00:12, 78.67it/s]
Initializing free energies using MBAR - [inc: 3.4e-04]: 2%|1 | 17/1000 [00:00<00:12, 78.67it/s]
Initializing free energies using MBAR - [inc: 2.5e-04]: 2%|1 | 18/1000 [00:00<00:12, 78.67it/s]
Initializing free energies using MBAR - [inc: 1.9e-04]: 2%|1 | 19/1000 [00:00<00:12, 78.67it/s]
Initializing free energies using MBAR - [inc: 1.4e-04]: 2%|2 | 20/1000 [00:00<00:12, 78.67it/s]
Initializing free energies using MBAR - [inc: 1.0e-04]: 2%|2 | 21/1000 [00:00<00:12, 78.67it/s]
Initializing free energies using MBAR - [inc: 7.7e-05]: 2%|2 | 22/1000 [00:00<00:12, 78.67it/s]
Initializing free energies using MBAR - [inc: 5.7e-05]: 2%|2 | 23/1000 [00:00<00:12, 78.67it/s]
Initializing free energies using MBAR - [inc: 4.2e-05]: 2%|2 | 24/1000 [00:00<00:12, 78.67it/s]
Initializing free energies using MBAR - [inc: 3.1e-05]: 2%|2 | 25/1000 [00:00<00:12, 78.67it/s]
Initializing free energies using MBAR - [inc: 2.3e-05]: 3%|2 | 26/1000 [00:00<00:12, 78.67it/s]
Initializing free energies using MBAR - [inc: 1.7e-05]: 3%|2 | 27/1000 [00:00<00:12, 78.67it/s]
Initializing free energies using MBAR - [inc: 1.3e-05]: 3%|2 | 28/1000 [00:00<00:12, 78.67it/s]
Initializing free energies using MBAR - [inc: 1.3e-05]: 3%|2 | 29/1000 [00:00<00:09, 98.09it/s]
Initializing free energies using MBAR - [inc: 9.6e-06]: 3%|2 | 29/1000 [00:00<00:09, 98.09it/s]
Initializing free energies using MBAR - [inc: 7.1e-06]: 3%|3 | 30/1000 [00:00<00:09, 98.09it/s]
Initializing free energies using MBAR - [inc: 5.3e-06]: 3%|3 | 31/1000 [00:00<00:09, 98.09it/s]
Initializing free energies using MBAR - [inc: 3.9e-06]: 3%|3 | 32/1000 [00:00<00:09, 98.09it/s]
Initializing free energies using MBAR - [inc: 2.9e-06]: 3%|3 | 33/1000 [00:00<00:09, 98.09it/s]
Initializing free energies using MBAR - [inc: 2.2e-06]: 3%|3 | 34/1000 [00:00<00:09, 98.09it/s]
Initializing free energies using MBAR - [inc: 1.6e-06]: 4%|3 | 35/1000 [00:00<00:09, 98.09it/s]
Initializing free energies using MBAR - [inc: 1.2e-06]: 4%|3 | 36/1000 [00:00<00:09, 98.09it/s]
Initializing free energies using MBAR - [inc: 8.8e-07]: 4%|3 | 37/1000 [00:00<00:09, 98.09it/s]
Initializing free energies using MBAR - [inc: 6.6e-07]: 4%|3 | 38/1000 [00:00<00:09, 98.09it/s]
Initializing free energies using MBAR - [inc: 4.9e-07]: 4%|3 | 39/1000 [00:00<00:09, 98.09it/s]
Initializing free energies using MBAR - [inc: 4.9e-07]: 4%|4 | 40/1000 [00:00<00:09, 101.00it/s]
Initializing free energies using MBAR - [inc: 3.6e-07]: 4%|4 | 40/1000 [00:00<00:09, 101.00it/s]
Initializing free energies using MBAR - [inc: 2.7e-07]: 4%|4 | 41/1000 [00:00<00:09, 101.00it/s]
Initializing free energies using MBAR - [inc: 2.0e-07]: 4%|4 | 42/1000 [00:00<00:09, 101.00it/s]
Initializing free energies using MBAR - [inc: 1.5e-07]: 4%|4 | 43/1000 [00:00<00:09, 101.00it/s]
Initializing free energies using MBAR - [inc: 1.1e-07]: 4%|4 | 44/1000 [00:00<00:09, 101.00it/s]
Initializing free energies using MBAR - [inc: 8.2e-08]: 4%|4 | 45/1000 [00:00<00:09, 101.00it/s]
Initializing free energies using MBAR - [inc: 6.1e-08]: 5%|4 | 46/1000 [00:00<00:09, 101.00it/s]
Initializing free energies using MBAR - [inc: 4.5e-08]: 5%|4 | 47/1000 [00:00<00:09, 101.00it/s]
Initializing free energies using MBAR - [inc: 3.3e-08]: 5%|4 | 48/1000 [00:00<00:09, 101.00it/s]
Initializing free energies using MBAR - [inc: 2.5e-08]: 5%|4 | 49/1000 [00:00<00:09, 101.00it/s]
Initializing free energies using MBAR - [inc: 1.8e-08]: 5%|5 | 50/1000 [00:00<00:09, 101.00it/s]
Initializing free energies using MBAR - [inc: 1.8e-08]: 5%|5 | 51/1000 [00:00<00:09, 99.33it/s]
Initializing free energies using MBAR - [inc: 1.4e-08]: 5%|5 | 51/1000 [00:00<00:09, 99.33it/s]
Initializing free energies using MBAR - [inc: 1.0e-08]: 5%|5 | 52/1000 [00:00<00:09, 99.33it/s]
Initializing free energies using MBAR - [inc: 7.5e-09]: 5%|5 | 53/1000 [00:00<00:09, 99.33it/s]
Initializing free energies using MBAR - [inc: 7.5e-09]: 100%|##########| 53/53 [00:00<00:00, 94.48it/s]
0%| | 0/1000 [00:00<?, ?it/s]
Running TRAM estimate: 0%| | 0/1000 [00:00<?, ?it/s]
Running TRAM estimate - [inc: 9.9e-01]: 0%| | 1/1000 [00:00<00:15, 66.10it/s]
Running TRAM estimate - [inc: 3.9e-01]: 0%| | 2/1000 [00:00<00:11, 86.39it/s]
Running TRAM estimate - [inc: 2.6e-01]: 0%| | 3/1000 [00:00<00:10, 96.45it/s]
Running TRAM estimate - [inc: 1.9e-01]: 0%| | 4/1000 [00:00<00:09, 102.28it/s]
Running TRAM estimate - [inc: 1.4e-01]: 0%| | 5/1000 [00:00<00:09, 106.10it/s]
Running TRAM estimate - [inc: 1.1e-01]: 1%| | 6/1000 [00:00<00:09, 108.81it/s]
Running TRAM estimate - [inc: 8.4e-02]: 1%| | 7/1000 [00:00<00:08, 110.74it/s]
Running TRAM estimate - [inc: 6.6e-02]: 1%| | 8/1000 [00:00<00:08, 112.20it/s]
Running TRAM estimate - [inc: 5.2e-02]: 1%| | 9/1000 [00:00<00:10, 95.10it/s]
Running TRAM estimate - [inc: 5.2e-02]: 1%|1 | 10/1000 [00:00<00:10, 97.58it/s]
Running TRAM estimate - [inc: 4.1e-02]: 1%|1 | 10/1000 [00:00<00:10, 97.58it/s]
Running TRAM estimate - [inc: 3.2e-02]: 1%|1 | 11/1000 [00:00<00:10, 97.58it/s]
Running TRAM estimate - [inc: 2.6e-02]: 1%|1 | 12/1000 [00:00<00:10, 97.58it/s]
Running TRAM estimate - [inc: 2.0e-02]: 1%|1 | 13/1000 [00:00<00:10, 97.58it/s]
Running TRAM estimate - [inc: 1.7e-02]: 1%|1 | 14/1000 [00:00<00:10, 97.58it/s]
Running TRAM estimate - [inc: 1.4e-02]: 2%|1 | 15/1000 [00:00<00:10, 97.58it/s]
Running TRAM estimate - [inc: 1.1e-02]: 2%|1 | 16/1000 [00:00<00:10, 97.58it/s]
Running TRAM estimate - [inc: 9.4e-03]: 2%|1 | 17/1000 [00:00<00:10, 97.58it/s]
Running TRAM estimate - [inc: 7.9e-03]: 2%|1 | 18/1000 [00:00<00:10, 97.58it/s]
Running TRAM estimate - [inc: 6.7e-03]: 2%|1 | 19/1000 [00:00<00:10, 97.58it/s]
Running TRAM estimate - [inc: 5.8e-03]: 2%|2 | 20/1000 [00:00<00:10, 97.58it/s]
Running TRAM estimate - [inc: 5.8e-03]: 2%|2 | 21/1000 [00:00<00:09, 102.07it/s]
Running TRAM estimate - [inc: 5.0e-03]: 2%|2 | 21/1000 [00:00<00:09, 102.07it/s]
Running TRAM estimate - [inc: 4.4e-03]: 2%|2 | 22/1000 [00:00<00:09, 102.07it/s]
Running TRAM estimate - [inc: 3.9e-03]: 2%|2 | 23/1000 [00:00<00:09, 102.07it/s]
Running TRAM estimate - [inc: 3.4e-03]: 2%|2 | 24/1000 [00:00<00:09, 102.07it/s]
Running TRAM estimate - [inc: 3.1e-03]: 2%|2 | 25/1000 [00:00<00:09, 102.07it/s]
Running TRAM estimate - [inc: 2.8e-03]: 3%|2 | 26/1000 [00:00<00:09, 102.07it/s]
Running TRAM estimate - [inc: 2.6e-03]: 3%|2 | 27/1000 [00:00<00:09, 102.07it/s]
Running TRAM estimate - [inc: 2.4e-03]: 3%|2 | 28/1000 [00:00<00:09, 102.07it/s]
Running TRAM estimate - [inc: 2.3e-03]: 3%|2 | 29/1000 [00:00<00:09, 102.07it/s]
Running TRAM estimate - [inc: 2.2e-03]: 3%|3 | 30/1000 [00:00<00:09, 102.07it/s]
Running TRAM estimate - [inc: 2.1e-03]: 3%|3 | 31/1000 [00:00<00:09, 102.07it/s]
Running TRAM estimate - [inc: 2.0e-03]: 3%|3 | 32/1000 [00:00<00:09, 102.07it/s]
Running TRAM estimate - [inc: 1.9e-03]: 3%|3 | 33/1000 [00:00<00:09, 102.07it/s]
Running TRAM estimate - [inc: 1.9e-03]: 3%|3 | 34/1000 [00:00<00:08, 111.68it/s]
Running TRAM estimate - [inc: 1.9e-03]: 3%|3 | 34/1000 [00:00<00:08, 111.68it/s]
Running TRAM estimate - [inc: 1.8e-03]: 4%|3 | 35/1000 [00:00<00:08, 111.68it/s]
Running TRAM estimate - [inc: 1.8e-03]: 4%|3 | 36/1000 [00:00<00:08, 111.68it/s]
Running TRAM estimate - [inc: 1.8e-03]: 4%|3 | 37/1000 [00:00<00:08, 111.68it/s]
Running TRAM estimate - [inc: 1.8e-03]: 4%|3 | 38/1000 [00:00<00:08, 111.68it/s]
Running TRAM estimate - [inc: 1.7e-03]: 4%|3 | 39/1000 [00:00<00:08, 111.68it/s]
Running TRAM estimate - [inc: 1.7e-03]: 4%|4 | 40/1000 [00:00<00:08, 111.68it/s]
Running TRAM estimate - [inc: 1.7e-03]: 4%|4 | 41/1000 [00:00<00:08, 111.68it/s]
Running TRAM estimate - [inc: 1.7e-03]: 4%|4 | 42/1000 [00:00<00:08, 111.68it/s]
Running TRAM estimate - [inc: 1.7e-03]: 4%|4 | 43/1000 [00:00<00:08, 111.68it/s]
Running TRAM estimate - [inc: 1.7e-03]: 4%|4 | 44/1000 [00:00<00:08, 111.68it/s]
Running TRAM estimate - [inc: 1.7e-03]: 4%|4 | 45/1000 [00:00<00:08, 111.68it/s]
Running TRAM estimate - [inc: 1.7e-03]: 5%|4 | 46/1000 [00:00<00:08, 109.37it/s]
Running TRAM estimate - [inc: 1.7e-03]: 5%|4 | 46/1000 [00:00<00:08, 109.37it/s]
Running TRAM estimate - [inc: 1.7e-03]: 5%|4 | 47/1000 [00:00<00:08, 109.37it/s]
Running TRAM estimate - [inc: 1.7e-03]: 5%|4 | 48/1000 [00:00<00:08, 109.37it/s]
Running TRAM estimate - [inc: 1.7e-03]: 5%|4 | 49/1000 [00:00<00:08, 109.37it/s]
Running TRAM estimate - [inc: 1.7e-03]: 5%|5 | 50/1000 [00:00<00:08, 109.37it/s]
Running TRAM estimate - [inc: 1.7e-03]: 5%|5 | 51/1000 [00:00<00:08, 109.37it/s]
Running TRAM estimate - [inc: 1.7e-03]: 5%|5 | 52/1000 [00:00<00:08, 109.37it/s]
Running TRAM estimate - [inc: 1.7e-03]: 5%|5 | 53/1000 [00:00<00:08, 109.37it/s]
Running TRAM estimate - [inc: 1.7e-03]: 5%|5 | 54/1000 [00:00<00:08, 109.37it/s]
Running TRAM estimate - [inc: 1.7e-03]: 6%|5 | 55/1000 [00:00<00:08, 109.37it/s]
Running TRAM estimate - [inc: 1.7e-03]: 6%|5 | 56/1000 [00:00<00:08, 109.37it/s]
Running TRAM estimate - [inc: 1.7e-03]: 6%|5 | 57/1000 [00:00<00:08, 109.37it/s]
Running TRAM estimate - [inc: 1.7e-03]: 6%|5 | 58/1000 [00:00<00:08, 109.37it/s]
Running TRAM estimate - [inc: 1.7e-03]: 6%|5 | 59/1000 [00:00<00:08, 115.34it/s]
Running TRAM estimate - [inc: 1.7e-03]: 6%|5 | 59/1000 [00:00<00:08, 115.34it/s]
Running TRAM estimate - [inc: 1.6e-03]: 6%|6 | 60/1000 [00:00<00:08, 115.34it/s]
Running TRAM estimate - [inc: 1.6e-03]: 6%|6 | 61/1000 [00:00<00:08, 115.34it/s]
Running TRAM estimate - [inc: 1.6e-03]: 6%|6 | 62/1000 [00:00<00:08, 115.34it/s]
Running TRAM estimate - [inc: 1.6e-03]: 6%|6 | 63/1000 [00:00<00:08, 115.34it/s]
Running TRAM estimate - [inc: 1.6e-03]: 6%|6 | 64/1000 [00:00<00:08, 115.34it/s]
Running TRAM estimate - [inc: 1.6e-03]: 6%|6 | 65/1000 [00:00<00:08, 115.34it/s]
Running TRAM estimate - [inc: 1.6e-03]: 7%|6 | 66/1000 [00:00<00:08, 115.34it/s]
Running TRAM estimate - [inc: 1.6e-03]: 7%|6 | 67/1000 [00:00<00:08, 115.34it/s]
Running TRAM estimate - [inc: 1.6e-03]: 7%|6 | 68/1000 [00:00<00:08, 115.34it/s]
Running TRAM estimate - [inc: 1.6e-03]: 7%|6 | 69/1000 [00:00<00:08, 115.34it/s]
Running TRAM estimate - [inc: 1.6e-03]: 7%|7 | 70/1000 [00:00<00:08, 115.34it/s]
Running TRAM estimate - [inc: 1.6e-03]: 7%|7 | 71/1000 [00:00<00:08, 106.51it/s]
Running TRAM estimate - [inc: 1.6e-03]: 7%|7 | 71/1000 [00:00<00:08, 106.51it/s]
Running TRAM estimate - [inc: 1.6e-03]: 7%|7 | 72/1000 [00:00<00:08, 106.51it/s]
Running TRAM estimate - [inc: 1.6e-03]: 7%|7 | 73/1000 [00:00<00:08, 106.51it/s]
Running TRAM estimate - [inc: 1.6e-03]: 7%|7 | 74/1000 [00:00<00:08, 106.51it/s]
Running TRAM estimate - [inc: 1.6e-03]: 8%|7 | 75/1000 [00:00<00:08, 106.51it/s]
Running TRAM estimate - [inc: 1.6e-03]: 8%|7 | 76/1000 [00:00<00:08, 106.51it/s]
Running TRAM estimate - [inc: 1.6e-03]: 8%|7 | 77/1000 [00:00<00:08, 106.51it/s]
Running TRAM estimate - [inc: 1.6e-03]: 8%|7 | 78/1000 [00:00<00:08, 106.51it/s]
Running TRAM estimate - [inc: 1.6e-03]: 8%|7 | 79/1000 [00:00<00:08, 106.51it/s]
Running TRAM estimate - [inc: 1.6e-03]: 8%|8 | 80/1000 [00:00<00:08, 106.51it/s]
Running TRAM estimate - [inc: 1.6e-03]: 8%|8 | 81/1000 [00:00<00:08, 106.51it/s]
Running TRAM estimate - [inc: 1.6e-03]: 8%|8 | 82/1000 [00:00<00:08, 106.51it/s]
Running TRAM estimate - [inc: 1.6e-03]: 8%|8 | 83/1000 [00:00<00:08, 106.51it/s]
Running TRAM estimate - [inc: 1.6e-03]: 8%|8 | 84/1000 [00:00<00:08, 112.25it/s]
Running TRAM estimate - [inc: 1.6e-03]: 8%|8 | 84/1000 [00:00<00:08, 112.25it/s]
Running TRAM estimate - [inc: 1.6e-03]: 8%|8 | 85/1000 [00:00<00:08, 112.25it/s]
Running TRAM estimate - [inc: 1.6e-03]: 9%|8 | 86/1000 [00:00<00:08, 112.25it/s]
Running TRAM estimate - [inc: 1.6e-03]: 9%|8 | 87/1000 [00:00<00:08, 112.25it/s]
Running TRAM estimate - [inc: 1.6e-03]: 9%|8 | 88/1000 [00:00<00:08, 112.25it/s]
Running TRAM estimate - [inc: 1.6e-03]: 9%|8 | 89/1000 [00:00<00:08, 112.25it/s]
Running TRAM estimate - [inc: 1.6e-03]: 9%|9 | 90/1000 [00:00<00:08, 112.25it/s]
Running TRAM estimate - [inc: 1.6e-03]: 9%|9 | 91/1000 [00:00<00:08, 112.25it/s]
Running TRAM estimate - [inc: 1.6e-03]: 9%|9 | 92/1000 [00:00<00:08, 112.25it/s]
Running TRAM estimate - [inc: 1.6e-03]: 9%|9 | 93/1000 [00:00<00:08, 112.25it/s]
Running TRAM estimate - [inc: 1.6e-03]: 9%|9 | 94/1000 [00:00<00:08, 112.25it/s]
Running TRAM estimate - [inc: 1.6e-03]: 10%|9 | 95/1000 [00:00<00:08, 112.25it/s]
Running TRAM estimate - [inc: 1.6e-03]: 10%|9 | 96/1000 [00:00<00:08, 112.25it/s]
Running TRAM estimate - [inc: 1.6e-03]: 10%|9 | 97/1000 [00:00<00:07, 116.73it/s]
Running TRAM estimate - [inc: 1.6e-03]: 10%|9 | 97/1000 [00:00<00:07, 116.73it/s]
Running TRAM estimate - [inc: 1.6e-03]: 10%|9 | 98/1000 [00:00<00:07, 116.73it/s]
Running TRAM estimate - [inc: 1.6e-03]: 10%|9 | 99/1000 [00:00<00:07, 116.73it/s]
Running TRAM estimate - [inc: 1.6e-03]: 10%|# | 100/1000 [00:00<00:07, 116.73it/s]
Running TRAM estimate - [inc: 1.6e-03]: 10%|# | 101/1000 [00:00<00:07, 116.73it/s]
Running TRAM estimate - [inc: 1.6e-03]: 10%|# | 102/1000 [00:00<00:07, 116.73it/s]
Running TRAM estimate - [inc: 1.6e-03]: 10%|# | 103/1000 [00:00<00:07, 116.73it/s]
Running TRAM estimate - [inc: 1.6e-03]: 10%|# | 104/1000 [00:00<00:07, 116.73it/s]
Running TRAM estimate - [inc: 1.6e-03]: 10%|# | 105/1000 [00:00<00:07, 116.73it/s]
Running TRAM estimate - [inc: 1.6e-03]: 11%|# | 106/1000 [00:00<00:07, 116.73it/s]
Running TRAM estimate - [inc: 1.6e-03]: 11%|# | 107/1000 [00:00<00:07, 116.73it/s]
Running TRAM estimate - [inc: 1.6e-03]: 11%|# | 108/1000 [00:00<00:07, 116.73it/s]
Running TRAM estimate - [inc: 1.6e-03]: 11%|# | 109/1000 [00:00<00:07, 116.73it/s]
Running TRAM estimate - [inc: 1.6e-03]: 11%|#1 | 110/1000 [00:00<00:07, 118.95it/s]
Running TRAM estimate - [inc: 1.6e-03]: 11%|#1 | 110/1000 [00:00<00:07, 118.95it/s]
Running TRAM estimate - [inc: 1.6e-03]: 11%|#1 | 111/1000 [00:00<00:07, 118.95it/s]
Running TRAM estimate - [inc: 1.6e-03]: 11%|#1 | 112/1000 [00:00<00:07, 118.95it/s]
Running TRAM estimate - [inc: 1.6e-03]: 11%|#1 | 113/1000 [00:01<00:07, 118.95it/s]
Running TRAM estimate - [inc: 1.6e-03]: 11%|#1 | 114/1000 [00:01<00:07, 118.95it/s]
Running TRAM estimate - [inc: 1.6e-03]: 12%|#1 | 115/1000 [00:01<00:07, 118.95it/s]
Running TRAM estimate - [inc: 1.6e-03]: 12%|#1 | 116/1000 [00:01<00:07, 118.95it/s]
Running TRAM estimate - [inc: 1.6e-03]: 12%|#1 | 117/1000 [00:01<00:07, 118.95it/s]
Running TRAM estimate - [inc: 1.6e-03]: 12%|#1 | 118/1000 [00:01<00:07, 118.95it/s]
Running TRAM estimate - [inc: 1.6e-03]: 12%|#1 | 119/1000 [00:01<00:07, 118.95it/s]
Running TRAM estimate - [inc: 1.6e-03]: 12%|#2 | 120/1000 [00:01<00:07, 118.95it/s]
Running TRAM estimate - [inc: 1.6e-03]: 12%|#2 | 121/1000 [00:01<00:07, 118.95it/s]
Running TRAM estimate - [inc: 1.6e-03]: 12%|#2 | 122/1000 [00:01<00:07, 115.10it/s]
Running TRAM estimate - [inc: 1.6e-03]: 12%|#2 | 122/1000 [00:01<00:07, 115.10it/s]
Running TRAM estimate - [inc: 1.6e-03]: 12%|#2 | 123/1000 [00:01<00:07, 115.10it/s]
Running TRAM estimate - [inc: 1.6e-03]: 12%|#2 | 124/1000 [00:01<00:07, 115.10it/s]
Running TRAM estimate - [inc: 1.6e-03]: 12%|#2 | 125/1000 [00:01<00:07, 115.10it/s]
Running TRAM estimate - [inc: 1.6e-03]: 13%|#2 | 126/1000 [00:01<00:07, 115.10it/s]
Running TRAM estimate - [inc: 1.6e-03]: 13%|#2 | 127/1000 [00:01<00:07, 115.10it/s]
Running TRAM estimate - [inc: 1.6e-03]: 13%|#2 | 128/1000 [00:01<00:07, 115.10it/s]
Running TRAM estimate - [inc: 1.6e-03]: 13%|#2 | 129/1000 [00:01<00:07, 115.10it/s]
Running TRAM estimate - [inc: 1.6e-03]: 13%|#3 | 130/1000 [00:01<00:07, 115.10it/s]
Running TRAM estimate - [inc: 1.6e-03]: 13%|#3 | 131/1000 [00:01<00:07, 115.10it/s]
Running TRAM estimate - [inc: 1.6e-03]: 13%|#3 | 132/1000 [00:01<00:07, 115.10it/s]
Running TRAM estimate - [inc: 1.6e-03]: 13%|#3 | 133/1000 [00:01<00:07, 115.10it/s]
Running TRAM estimate - [inc: 1.6e-03]: 13%|#3 | 134/1000 [00:01<00:07, 115.10it/s]
Running TRAM estimate - [inc: 1.6e-03]: 14%|#3 | 135/1000 [00:01<00:07, 119.24it/s]
Running TRAM estimate - [inc: 1.6e-03]: 14%|#3 | 135/1000 [00:01<00:07, 119.24it/s]
Running TRAM estimate - [inc: 1.6e-03]: 14%|#3 | 136/1000 [00:01<00:07, 119.24it/s]
Running TRAM estimate - [inc: 1.6e-03]: 14%|#3 | 137/1000 [00:01<00:07, 119.24it/s]
Running TRAM estimate - [inc: 1.6e-03]: 14%|#3 | 138/1000 [00:01<00:07, 119.24it/s]
Running TRAM estimate - [inc: 1.6e-03]: 14%|#3 | 139/1000 [00:01<00:07, 119.24it/s]
Running TRAM estimate - [inc: 1.6e-03]: 14%|#4 | 140/1000 [00:01<00:07, 119.24it/s]
Running TRAM estimate - [inc: 1.6e-03]: 14%|#4 | 141/1000 [00:01<00:07, 119.24it/s]
Running TRAM estimate - [inc: 1.6e-03]: 14%|#4 | 142/1000 [00:01<00:07, 119.24it/s]
Running TRAM estimate - [inc: 1.6e-03]: 14%|#4 | 143/1000 [00:01<00:07, 119.24it/s]
Running TRAM estimate - [inc: 1.6e-03]: 14%|#4 | 144/1000 [00:01<00:07, 119.24it/s]
Running TRAM estimate - [inc: 1.6e-03]: 14%|#4 | 145/1000 [00:01<00:07, 119.24it/s]
Running TRAM estimate - [inc: 1.6e-03]: 15%|#4 | 146/1000 [00:01<00:07, 119.24it/s]
Running TRAM estimate - [inc: 1.6e-03]: 15%|#4 | 147/1000 [00:01<00:07, 119.24it/s]
Running TRAM estimate - [inc: 1.6e-03]: 15%|#4 | 148/1000 [00:01<00:07, 120.81it/s]
Running TRAM estimate - [inc: 1.6e-03]: 15%|#4 | 148/1000 [00:01<00:07, 120.81it/s]
Running TRAM estimate - [inc: 1.6e-03]: 15%|#4 | 149/1000 [00:01<00:07, 120.81it/s]
Running TRAM estimate - [inc: 1.5e-03]: 15%|#5 | 150/1000 [00:01<00:07, 120.81it/s]
Running TRAM estimate - [inc: 1.5e-03]: 15%|#5 | 151/1000 [00:01<00:07, 120.81it/s]
Running TRAM estimate - [inc: 1.5e-03]: 15%|#5 | 152/1000 [00:01<00:07, 120.81it/s]
Running TRAM estimate - [inc: 1.5e-03]: 15%|#5 | 153/1000 [00:01<00:07, 120.81it/s]
Running TRAM estimate - [inc: 1.5e-03]: 15%|#5 | 154/1000 [00:01<00:07, 120.81it/s]
Running TRAM estimate - [inc: 1.5e-03]: 16%|#5 | 155/1000 [00:01<00:06, 120.81it/s]
Running TRAM estimate - [inc: 1.5e-03]: 16%|#5 | 156/1000 [00:01<00:06, 120.81it/s]
Running TRAM estimate - [inc: 1.5e-03]: 16%|#5 | 157/1000 [00:01<00:06, 120.81it/s]
Running TRAM estimate - [inc: 1.5e-03]: 16%|#5 | 158/1000 [00:01<00:06, 120.81it/s]
Running TRAM estimate - [inc: 1.5e-03]: 16%|#5 | 159/1000 [00:01<00:06, 120.81it/s]
Running TRAM estimate - [inc: 1.5e-03]: 16%|#6 | 160/1000 [00:01<00:06, 120.81it/s]
Running TRAM estimate - [inc: 1.5e-03]: 16%|#6 | 161/1000 [00:01<00:06, 122.89it/s]
Running TRAM estimate - [inc: 1.5e-03]: 16%|#6 | 161/1000 [00:01<00:06, 122.89it/s]
Running TRAM estimate - [inc: 1.5e-03]: 16%|#6 | 162/1000 [00:01<00:06, 122.89it/s]
Running TRAM estimate - [inc: 1.5e-03]: 16%|#6 | 163/1000 [00:01<00:06, 122.89it/s]
Running TRAM estimate - [inc: 1.5e-03]: 16%|#6 | 164/1000 [00:01<00:06, 122.89it/s]
Running TRAM estimate - [inc: 1.5e-03]: 16%|#6 | 165/1000 [00:01<00:06, 122.89it/s]
Running TRAM estimate - [inc: 1.5e-03]: 17%|#6 | 166/1000 [00:01<00:06, 122.89it/s]
Running TRAM estimate - [inc: 1.5e-03]: 17%|#6 | 167/1000 [00:01<00:06, 122.89it/s]
Running TRAM estimate - [inc: 1.5e-03]: 17%|#6 | 168/1000 [00:01<00:06, 122.89it/s]
Running TRAM estimate - [inc: 1.5e-03]: 17%|#6 | 169/1000 [00:01<00:06, 122.89it/s]
Running TRAM estimate - [inc: 1.5e-03]: 17%|#7 | 170/1000 [00:01<00:06, 122.89it/s]
Running TRAM estimate - [inc: 1.5e-03]: 17%|#7 | 171/1000 [00:01<00:06, 122.89it/s]
Running TRAM estimate - [inc: 1.5e-03]: 17%|#7 | 172/1000 [00:01<00:06, 122.89it/s]
Running TRAM estimate - [inc: 1.5e-03]: 17%|#7 | 173/1000 [00:01<00:06, 122.89it/s]
Running TRAM estimate - [inc: 1.5e-03]: 17%|#7 | 174/1000 [00:01<00:06, 124.72it/s]
Running TRAM estimate - [inc: 1.5e-03]: 17%|#7 | 174/1000 [00:01<00:06, 124.72it/s]
Running TRAM estimate - [inc: 1.5e-03]: 18%|#7 | 175/1000 [00:01<00:06, 124.72it/s]
Running TRAM estimate - [inc: 1.5e-03]: 18%|#7 | 176/1000 [00:01<00:06, 124.72it/s]
Running TRAM estimate - [inc: 1.5e-03]: 18%|#7 | 177/1000 [00:01<00:06, 124.72it/s]
Running TRAM estimate - [inc: 1.5e-03]: 18%|#7 | 178/1000 [00:01<00:06, 124.72it/s]
Running TRAM estimate - [inc: 1.5e-03]: 18%|#7 | 179/1000 [00:01<00:06, 124.72it/s]
Running TRAM estimate - [inc: 1.5e-03]: 18%|#8 | 180/1000 [00:01<00:06, 124.72it/s]
Running TRAM estimate - [inc: 1.5e-03]: 18%|#8 | 181/1000 [00:01<00:06, 124.72it/s]
Running TRAM estimate - [inc: 1.5e-03]: 18%|#8 | 182/1000 [00:01<00:06, 124.72it/s]
Running TRAM estimate - [inc: 1.5e-03]: 18%|#8 | 183/1000 [00:01<00:06, 124.72it/s]
Running TRAM estimate - [inc: 1.5e-03]: 18%|#8 | 184/1000 [00:01<00:06, 124.72it/s]
Running TRAM estimate - [inc: 1.5e-03]: 18%|#8 | 185/1000 [00:01<00:06, 124.72it/s]
Running TRAM estimate - [inc: 1.5e-03]: 19%|#8 | 186/1000 [00:01<00:06, 124.72it/s]
Running TRAM estimate - [inc: 1.5e-03]: 19%|#8 | 187/1000 [00:01<00:06, 125.67it/s]
Running TRAM estimate - [inc: 1.5e-03]: 19%|#8 | 187/1000 [00:01<00:06, 125.67it/s]
Running TRAM estimate - [inc: 1.5e-03]: 19%|#8 | 188/1000 [00:01<00:06, 125.67it/s]
Running TRAM estimate - [inc: 1.5e-03]: 19%|#8 | 189/1000 [00:01<00:06, 125.67it/s]
Running TRAM estimate - [inc: 1.5e-03]: 19%|#9 | 190/1000 [00:01<00:06, 125.67it/s]
Running TRAM estimate - [inc: 1.5e-03]: 19%|#9 | 191/1000 [00:01<00:06, 125.67it/s]
Running TRAM estimate - [inc: 1.5e-03]: 19%|#9 | 192/1000 [00:01<00:06, 125.67it/s]
Running TRAM estimate - [inc: 1.5e-03]: 19%|#9 | 193/1000 [00:01<00:06, 125.67it/s]
Running TRAM estimate - [inc: 1.5e-03]: 19%|#9 | 194/1000 [00:01<00:06, 125.67it/s]
Running TRAM estimate - [inc: 1.5e-03]: 20%|#9 | 195/1000 [00:01<00:06, 125.67it/s]
Running TRAM estimate - [inc: 1.5e-03]: 20%|#9 | 196/1000 [00:01<00:06, 125.67it/s]
Running TRAM estimate - [inc: 1.5e-03]: 20%|#9 | 197/1000 [00:01<00:06, 125.67it/s]
Running TRAM estimate - [inc: 1.5e-03]: 20%|#9 | 198/1000 [00:01<00:06, 125.67it/s]
Running TRAM estimate - [inc: 1.5e-03]: 20%|#9 | 199/1000 [00:01<00:06, 125.67it/s]
Running TRAM estimate - [inc: 1.5e-03]: 20%|## | 200/1000 [00:01<00:06, 119.08it/s]
Running TRAM estimate - [inc: 1.5e-03]: 20%|## | 200/1000 [00:01<00:06, 119.08it/s]
Running TRAM estimate - [inc: 1.5e-03]: 20%|## | 201/1000 [00:01<00:06, 119.08it/s]
Running TRAM estimate - [inc: 1.5e-03]: 20%|## | 202/1000 [00:01<00:06, 119.08it/s]
Running TRAM estimate - [inc: 1.5e-03]: 20%|## | 203/1000 [00:01<00:06, 119.08it/s]
Running TRAM estimate - [inc: 1.5e-03]: 20%|## | 204/1000 [00:01<00:06, 119.08it/s]
Running TRAM estimate - [inc: 1.5e-03]: 20%|## | 205/1000 [00:01<00:06, 119.08it/s]
Running TRAM estimate - [inc: 1.5e-03]: 21%|## | 206/1000 [00:01<00:06, 119.08it/s]
Running TRAM estimate - [inc: 1.5e-03]: 21%|## | 207/1000 [00:01<00:06, 119.08it/s]
Running TRAM estimate - [inc: 1.5e-03]: 21%|## | 208/1000 [00:01<00:06, 119.08it/s]
Running TRAM estimate - [inc: 1.5e-03]: 21%|## | 209/1000 [00:01<00:06, 119.08it/s]
Running TRAM estimate - [inc: 1.5e-03]: 21%|##1 | 210/1000 [00:01<00:06, 119.08it/s]
Running TRAM estimate - [inc: 1.5e-03]: 21%|##1 | 211/1000 [00:01<00:06, 119.08it/s]
Running TRAM estimate - [inc: 1.5e-03]: 21%|##1 | 212/1000 [00:01<00:08, 96.29it/s]
Running TRAM estimate - [inc: 1.5e-03]: 21%|##1 | 212/1000 [00:01<00:08, 96.29it/s]
Running TRAM estimate - [inc: 1.5e-03]: 21%|##1 | 213/1000 [00:01<00:08, 96.29it/s]
Running TRAM estimate - [inc: 1.5e-03]: 21%|##1 | 214/1000 [00:01<00:08, 96.29it/s]
Running TRAM estimate - [inc: 1.5e-03]: 22%|##1 | 215/1000 [00:01<00:08, 96.29it/s]
Running TRAM estimate - [inc: 1.5e-03]: 22%|##1 | 216/1000 [00:01<00:08, 96.29it/s]
Running TRAM estimate - [inc: 1.5e-03]: 22%|##1 | 217/1000 [00:01<00:08, 96.29it/s]
Running TRAM estimate - [inc: 1.5e-03]: 22%|##1 | 218/1000 [00:01<00:08, 96.29it/s]
Running TRAM estimate - [inc: 1.5e-03]: 22%|##1 | 219/1000 [00:01<00:08, 96.29it/s]
Running TRAM estimate - [inc: 1.5e-03]: 22%|##2 | 220/1000 [00:01<00:08, 96.29it/s]
Running TRAM estimate - [inc: 1.5e-03]: 22%|##2 | 221/1000 [00:02<00:08, 96.29it/s]
Running TRAM estimate - [inc: 1.5e-03]: 22%|##2 | 222/1000 [00:02<00:08, 96.29it/s]
Running TRAM estimate - [inc: 1.5e-03]: 22%|##2 | 223/1000 [00:02<00:08, 92.43it/s]
Running TRAM estimate - [inc: 1.5e-03]: 22%|##2 | 223/1000 [00:02<00:08, 92.43it/s]
Running TRAM estimate - [inc: 1.5e-03]: 22%|##2 | 224/1000 [00:02<00:08, 92.43it/s]
Running TRAM estimate - [inc: 1.5e-03]: 22%|##2 | 225/1000 [00:02<00:08, 92.43it/s]
Running TRAM estimate - [inc: 1.5e-03]: 23%|##2 | 226/1000 [00:02<00:08, 92.43it/s]
Running TRAM estimate - [inc: 1.5e-03]: 23%|##2 | 227/1000 [00:02<00:08, 92.43it/s]
Running TRAM estimate - [inc: 1.5e-03]: 23%|##2 | 228/1000 [00:02<00:08, 92.43it/s]
Running TRAM estimate - [inc: 1.4e-03]: 23%|##2 | 229/1000 [00:02<00:08, 92.43it/s]
Running TRAM estimate - [inc: 1.4e-03]: 23%|##3 | 230/1000 [00:02<00:08, 92.43it/s]
Running TRAM estimate - [inc: 1.4e-03]: 23%|##3 | 231/1000 [00:02<00:08, 92.43it/s]
Running TRAM estimate - [inc: 1.4e-03]: 23%|##3 | 232/1000 [00:02<00:08, 92.43it/s]
Running TRAM estimate - [inc: 1.4e-03]: 23%|##3 | 233/1000 [00:02<00:08, 87.98it/s]
Running TRAM estimate - [inc: 1.4e-03]: 23%|##3 | 233/1000 [00:02<00:08, 87.98it/s]
Running TRAM estimate - [inc: 1.4e-03]: 23%|##3 | 234/1000 [00:02<00:08, 87.98it/s]
Running TRAM estimate - [inc: 1.4e-03]: 24%|##3 | 235/1000 [00:02<00:08, 87.98it/s]
Running TRAM estimate - [inc: 1.4e-03]: 24%|##3 | 236/1000 [00:02<00:08, 87.98it/s]
Running TRAM estimate - [inc: 1.4e-03]: 24%|##3 | 237/1000 [00:02<00:08, 87.98it/s]
Running TRAM estimate - [inc: 1.4e-03]: 24%|##3 | 238/1000 [00:02<00:08, 87.98it/s]
Running TRAM estimate - [inc: 1.4e-03]: 24%|##3 | 239/1000 [00:02<00:08, 87.98it/s]
Running TRAM estimate - [inc: 1.4e-03]: 24%|##4 | 240/1000 [00:02<00:08, 87.98it/s]
Running TRAM estimate - [inc: 1.4e-03]: 24%|##4 | 241/1000 [00:02<00:08, 87.98it/s]
Running TRAM estimate - [inc: 1.4e-03]: 24%|##4 | 242/1000 [00:02<00:08, 87.98it/s]
Running TRAM estimate - [inc: 1.4e-03]: 24%|##4 | 243/1000 [00:02<00:08, 87.98it/s]
Running TRAM estimate - [inc: 1.4e-03]: 24%|##4 | 244/1000 [00:02<00:08, 87.98it/s]
Running TRAM estimate - [inc: 1.4e-03]: 24%|##4 | 245/1000 [00:02<00:08, 87.98it/s]
Running TRAM estimate - [inc: 1.4e-03]: 25%|##4 | 246/1000 [00:02<00:07, 97.73it/s]
Running TRAM estimate - [inc: 1.4e-03]: 25%|##4 | 246/1000 [00:02<00:07, 97.73it/s]
Running TRAM estimate - [inc: 1.4e-03]: 25%|##4 | 247/1000 [00:02<00:07, 97.73it/s]
Running TRAM estimate - [inc: 1.4e-03]: 25%|##4 | 248/1000 [00:02<00:07, 97.73it/s]
Running TRAM estimate - [inc: 1.4e-03]: 25%|##4 | 249/1000 [00:02<00:07, 97.73it/s]
Running TRAM estimate - [inc: 1.4e-03]: 25%|##5 | 250/1000 [00:02<00:07, 97.73it/s]
Running TRAM estimate - [inc: 1.4e-03]: 25%|##5 | 251/1000 [00:02<00:07, 97.73it/s]
Running TRAM estimate - [inc: 1.4e-03]: 25%|##5 | 252/1000 [00:02<00:07, 97.73it/s]
Running TRAM estimate - [inc: 1.4e-03]: 25%|##5 | 253/1000 [00:02<00:07, 97.73it/s]
Running TRAM estimate - [inc: 1.4e-03]: 25%|##5 | 254/1000 [00:02<00:07, 97.73it/s]
Running TRAM estimate - [inc: 1.4e-03]: 26%|##5 | 255/1000 [00:02<00:07, 97.73it/s]
Running TRAM estimate - [inc: 1.4e-03]: 26%|##5 | 256/1000 [00:02<00:07, 97.73it/s]
Running TRAM estimate - [inc: 1.4e-03]: 26%|##5 | 257/1000 [00:02<00:07, 97.73it/s]
Running TRAM estimate - [inc: 1.4e-03]: 26%|##5 | 258/1000 [00:02<00:07, 97.73it/s]
Running TRAM estimate - [inc: 1.4e-03]: 26%|##5 | 259/1000 [00:02<00:07, 105.50it/s]
Running TRAM estimate - [inc: 1.4e-03]: 26%|##5 | 259/1000 [00:02<00:07, 105.50it/s]
Running TRAM estimate - [inc: 1.4e-03]: 26%|##6 | 260/1000 [00:02<00:07, 105.50it/s]
Running TRAM estimate - [inc: 1.4e-03]: 26%|##6 | 261/1000 [00:02<00:07, 105.50it/s]
Running TRAM estimate - [inc: 1.4e-03]: 26%|##6 | 262/1000 [00:02<00:06, 105.50it/s]
Running TRAM estimate - [inc: 1.4e-03]: 26%|##6 | 263/1000 [00:02<00:06, 105.50it/s]
Running TRAM estimate - [inc: 1.4e-03]: 26%|##6 | 264/1000 [00:02<00:06, 105.50it/s]
Running TRAM estimate - [inc: 1.4e-03]: 26%|##6 | 265/1000 [00:02<00:06, 105.50it/s]
Running TRAM estimate - [inc: 1.4e-03]: 27%|##6 | 266/1000 [00:02<00:06, 105.50it/s]
Running TRAM estimate - [inc: 1.4e-03]: 27%|##6 | 267/1000 [00:02<00:06, 105.50it/s]
Running TRAM estimate - [inc: 1.4e-03]: 27%|##6 | 268/1000 [00:02<00:06, 105.50it/s]
Running TRAM estimate - [inc: 1.4e-03]: 27%|##6 | 269/1000 [00:02<00:06, 105.50it/s]
Running TRAM estimate - [inc: 1.4e-03]: 27%|##7 | 270/1000 [00:02<00:08, 81.46it/s]
Running TRAM estimate - [inc: 1.4e-03]: 27%|##7 | 270/1000 [00:02<00:08, 81.46it/s]
Running TRAM estimate - [inc: 1.4e-03]: 27%|##7 | 271/1000 [00:02<00:08, 81.46it/s]
Running TRAM estimate - [inc: 1.4e-03]: 27%|##7 | 272/1000 [00:02<00:08, 81.46it/s]
Running TRAM estimate - [inc: 1.4e-03]: 27%|##7 | 273/1000 [00:02<00:08, 81.46it/s]
Running TRAM estimate - [inc: 1.4e-03]: 27%|##7 | 274/1000 [00:02<00:08, 81.46it/s]
Running TRAM estimate - [inc: 1.4e-03]: 28%|##7 | 275/1000 [00:02<00:08, 81.46it/s]
Running TRAM estimate - [inc: 1.4e-03]: 28%|##7 | 276/1000 [00:02<00:08, 81.46it/s]
Running TRAM estimate - [inc: 1.4e-03]: 28%|##7 | 277/1000 [00:02<00:08, 81.46it/s]
Running TRAM estimate - [inc: 1.4e-03]: 28%|##7 | 278/1000 [00:02<00:08, 81.46it/s]
Running TRAM estimate - [inc: 1.4e-03]: 28%|##7 | 279/1000 [00:02<00:08, 81.46it/s]
Running TRAM estimate - [inc: 1.4e-03]: 28%|##8 | 280/1000 [00:02<00:12, 57.58it/s]
Running TRAM estimate - [inc: 1.4e-03]: 28%|##8 | 280/1000 [00:02<00:12, 57.58it/s]
Running TRAM estimate - [inc: 1.4e-03]: 28%|##8 | 281/1000 [00:02<00:12, 57.58it/s]
Running TRAM estimate - [inc: 1.4e-03]: 28%|##8 | 282/1000 [00:02<00:12, 57.58it/s]
Running TRAM estimate - [inc: 1.4e-03]: 28%|##8 | 283/1000 [00:02<00:12, 57.58it/s]
Running TRAM estimate - [inc: 1.4e-03]: 28%|##8 | 284/1000 [00:02<00:12, 57.58it/s]
Running TRAM estimate - [inc: 1.4e-03]: 28%|##8 | 285/1000 [00:02<00:12, 57.58it/s]
Running TRAM estimate - [inc: 1.4e-03]: 29%|##8 | 286/1000 [00:03<00:12, 57.58it/s]
Running TRAM estimate - [inc: 1.4e-03]: 29%|##8 | 287/1000 [00:03<00:12, 57.58it/s]
Running TRAM estimate - [inc: 1.4e-03]: 29%|##8 | 288/1000 [00:03<00:12, 56.90it/s]
Running TRAM estimate - [inc: 1.4e-03]: 29%|##8 | 288/1000 [00:03<00:12, 56.90it/s]
Running TRAM estimate - [inc: 1.4e-03]: 29%|##8 | 289/1000 [00:03<00:12, 56.90it/s]
Running TRAM estimate - [inc: 1.4e-03]: 29%|##9 | 290/1000 [00:03<00:12, 56.90it/s]
Running TRAM estimate - [inc: 1.4e-03]: 29%|##9 | 291/1000 [00:03<00:12, 56.90it/s]
Running TRAM estimate - [inc: 1.4e-03]: 29%|##9 | 292/1000 [00:03<00:12, 56.90it/s]
Running TRAM estimate - [inc: 1.4e-03]: 29%|##9 | 293/1000 [00:03<00:12, 56.90it/s]
Running TRAM estimate - [inc: 1.4e-03]: 29%|##9 | 294/1000 [00:03<00:12, 56.90it/s]
Running TRAM estimate - [inc: 1.4e-03]: 30%|##9 | 295/1000 [00:03<00:14, 48.65it/s]
Running TRAM estimate - [inc: 1.4e-03]: 30%|##9 | 295/1000 [00:03<00:14, 48.65it/s]
Running TRAM estimate - [inc: 1.4e-03]: 30%|##9 | 296/1000 [00:03<00:14, 48.65it/s]
Running TRAM estimate - [inc: 1.4e-03]: 30%|##9 | 297/1000 [00:03<00:14, 48.65it/s]
Running TRAM estimate - [inc: 1.4e-03]: 30%|##9 | 298/1000 [00:03<00:14, 48.65it/s]
Running TRAM estimate - [inc: 1.4e-03]: 30%|##9 | 299/1000 [00:03<00:14, 48.65it/s]
Running TRAM estimate - [inc: 1.3e-03]: 30%|### | 300/1000 [00:03<00:14, 48.65it/s]
Running TRAM estimate - [inc: 1.3e-03]: 30%|### | 301/1000 [00:03<00:14, 48.44it/s]
Running TRAM estimate - [inc: 1.3e-03]: 30%|### | 301/1000 [00:03<00:14, 48.44it/s]
Running TRAM estimate - [inc: 1.3e-03]: 30%|### | 302/1000 [00:03<00:14, 48.44it/s]
Running TRAM estimate - [inc: 1.3e-03]: 30%|### | 303/1000 [00:03<00:14, 48.44it/s]
Running TRAM estimate - [inc: 1.3e-03]: 30%|### | 304/1000 [00:03<00:14, 48.44it/s]
Running TRAM estimate - [inc: 1.3e-03]: 30%|### | 305/1000 [00:03<00:14, 48.44it/s]
Running TRAM estimate - [inc: 1.3e-03]: 31%|### | 306/1000 [00:03<00:14, 48.44it/s]
Running TRAM estimate - [inc: 1.3e-03]: 31%|### | 307/1000 [00:03<00:14, 48.24it/s]
Running TRAM estimate - [inc: 1.3e-03]: 31%|### | 307/1000 [00:03<00:14, 48.24it/s]
Running TRAM estimate - [inc: 1.3e-03]: 31%|### | 308/1000 [00:03<00:14, 48.24it/s]
Running TRAM estimate - [inc: 1.3e-03]: 31%|### | 309/1000 [00:03<00:14, 48.24it/s]
Running TRAM estimate - [inc: 1.3e-03]: 31%|###1 | 310/1000 [00:03<00:14, 48.24it/s]
Running TRAM estimate - [inc: 1.3e-03]: 31%|###1 | 311/1000 [00:03<00:14, 48.24it/s]
Running TRAM estimate - [inc: 1.3e-03]: 31%|###1 | 312/1000 [00:03<00:14, 48.24it/s]
Running TRAM estimate - [inc: 1.3e-03]: 31%|###1 | 313/1000 [00:03<00:14, 48.24it/s]
Running TRAM estimate - [inc: 1.3e-03]: 31%|###1 | 314/1000 [00:03<00:14, 48.24it/s]
Running TRAM estimate - [inc: 1.3e-03]: 32%|###1 | 315/1000 [00:03<00:14, 48.24it/s]
Running TRAM estimate - [inc: 1.3e-03]: 32%|###1 | 316/1000 [00:03<00:14, 48.24it/s]
Running TRAM estimate - [inc: 1.3e-03]: 32%|###1 | 317/1000 [00:03<00:14, 48.24it/s]
Running TRAM estimate - [inc: 1.3e-03]: 32%|###1 | 318/1000 [00:03<00:11, 60.69it/s]
Running TRAM estimate - [inc: 1.3e-03]: 32%|###1 | 318/1000 [00:03<00:11, 60.69it/s]
Running TRAM estimate - [inc: 1.3e-03]: 32%|###1 | 319/1000 [00:03<00:11, 60.69it/s]
Running TRAM estimate - [inc: 1.3e-03]: 32%|###2 | 320/1000 [00:03<00:11, 60.69it/s]
Running TRAM estimate - [inc: 1.3e-03]: 32%|###2 | 321/1000 [00:03<00:11, 60.69it/s]
Running TRAM estimate - [inc: 1.3e-03]: 32%|###2 | 322/1000 [00:03<00:11, 60.69it/s]
Running TRAM estimate - [inc: 1.3e-03]: 32%|###2 | 323/1000 [00:03<00:11, 60.69it/s]
Running TRAM estimate - [inc: 1.3e-03]: 32%|###2 | 324/1000 [00:03<00:11, 60.69it/s]
Running TRAM estimate - [inc: 1.3e-03]: 32%|###2 | 325/1000 [00:03<00:11, 60.69it/s]
Running TRAM estimate - [inc: 1.3e-03]: 33%|###2 | 326/1000 [00:03<00:11, 60.69it/s]
Running TRAM estimate - [inc: 1.3e-03]: 33%|###2 | 327/1000 [00:03<00:11, 60.69it/s]
Running TRAM estimate - [inc: 1.3e-03]: 33%|###2 | 328/1000 [00:03<00:11, 60.69it/s]
Running TRAM estimate - [inc: 1.3e-03]: 33%|###2 | 329/1000 [00:03<00:11, 60.69it/s]
Running TRAM estimate - [inc: 1.3e-03]: 33%|###3 | 330/1000 [00:03<00:11, 60.69it/s]
Running TRAM estimate - [inc: 1.3e-03]: 33%|###3 | 331/1000 [00:03<00:08, 76.28it/s]
Running TRAM estimate - [inc: 1.3e-03]: 33%|###3 | 331/1000 [00:03<00:08, 76.28it/s]
Running TRAM estimate - [inc: 1.3e-03]: 33%|###3 | 332/1000 [00:03<00:08, 76.28it/s]
Running TRAM estimate - [inc: 1.3e-03]: 33%|###3 | 333/1000 [00:03<00:08, 76.28it/s]
Running TRAM estimate - [inc: 1.3e-03]: 33%|###3 | 334/1000 [00:03<00:08, 76.28it/s]
Running TRAM estimate - [inc: 1.3e-03]: 34%|###3 | 335/1000 [00:03<00:08, 76.28it/s]
Running TRAM estimate - [inc: 1.3e-03]: 34%|###3 | 336/1000 [00:03<00:08, 76.28it/s]
Running TRAM estimate - [inc: 1.3e-03]: 34%|###3 | 337/1000 [00:03<00:08, 76.28it/s]
Running TRAM estimate - [inc: 1.3e-03]: 34%|###3 | 338/1000 [00:03<00:08, 76.28it/s]
Running TRAM estimate - [inc: 1.3e-03]: 34%|###3 | 339/1000 [00:03<00:08, 76.28it/s]
Running TRAM estimate - [inc: 1.3e-03]: 34%|###4 | 340/1000 [00:03<00:08, 76.28it/s]
Running TRAM estimate - [inc: 1.3e-03]: 34%|###4 | 341/1000 [00:03<00:08, 76.28it/s]
Running TRAM estimate - [inc: 1.3e-03]: 34%|###4 | 342/1000 [00:03<00:08, 76.28it/s]
Running TRAM estimate - [inc: 1.3e-03]: 34%|###4 | 343/1000 [00:03<00:08, 76.28it/s]
Running TRAM estimate - [inc: 1.3e-03]: 34%|###4 | 344/1000 [00:03<00:07, 88.44it/s]
Running TRAM estimate - [inc: 1.3e-03]: 34%|###4 | 344/1000 [00:03<00:07, 88.44it/s]
Running TRAM estimate - [inc: 1.3e-03]: 34%|###4 | 345/1000 [00:03<00:07, 88.44it/s]
Running TRAM estimate - [inc: 1.3e-03]: 35%|###4 | 346/1000 [00:03<00:07, 88.44it/s]
Running TRAM estimate - [inc: 1.3e-03]: 35%|###4 | 347/1000 [00:03<00:07, 88.44it/s]
Running TRAM estimate - [inc: 1.3e-03]: 35%|###4 | 348/1000 [00:03<00:07, 88.44it/s]
Running TRAM estimate - [inc: 1.3e-03]: 35%|###4 | 349/1000 [00:03<00:07, 88.44it/s]
Running TRAM estimate - [inc: 1.3e-03]: 35%|###5 | 350/1000 [00:03<00:07, 88.44it/s]
Running TRAM estimate - [inc: 1.3e-03]: 35%|###5 | 351/1000 [00:03<00:07, 88.44it/s]
Running TRAM estimate - [inc: 1.3e-03]: 35%|###5 | 352/1000 [00:03<00:07, 88.44it/s]
Running TRAM estimate - [inc: 1.3e-03]: 35%|###5 | 353/1000 [00:03<00:07, 88.44it/s]
Running TRAM estimate - [inc: 1.3e-03]: 35%|###5 | 354/1000 [00:03<00:07, 88.44it/s]
Running TRAM estimate - [inc: 1.3e-03]: 36%|###5 | 355/1000 [00:03<00:07, 88.44it/s]
Running TRAM estimate - [inc: 1.3e-03]: 36%|###5 | 356/1000 [00:03<00:06, 95.72it/s]
Running TRAM estimate - [inc: 1.3e-03]: 36%|###5 | 356/1000 [00:03<00:06, 95.72it/s]
Running TRAM estimate - [inc: 1.3e-03]: 36%|###5 | 357/1000 [00:03<00:06, 95.72it/s]
Running TRAM estimate - [inc: 1.3e-03]: 36%|###5 | 358/1000 [00:03<00:06, 95.72it/s]
Running TRAM estimate - [inc: 1.3e-03]: 36%|###5 | 359/1000 [00:03<00:06, 95.72it/s]
Running TRAM estimate - [inc: 1.3e-03]: 36%|###6 | 360/1000 [00:03<00:06, 95.72it/s]
Running TRAM estimate - [inc: 1.3e-03]: 36%|###6 | 361/1000 [00:03<00:06, 95.72it/s]
Running TRAM estimate - [inc: 1.3e-03]: 36%|###6 | 362/1000 [00:03<00:06, 95.72it/s]
Running TRAM estimate - [inc: 1.3e-03]: 36%|###6 | 363/1000 [00:03<00:06, 95.72it/s]
Running TRAM estimate - [inc: 1.3e-03]: 36%|###6 | 364/1000 [00:03<00:06, 95.72it/s]
Running TRAM estimate - [inc: 1.3e-03]: 36%|###6 | 365/1000 [00:03<00:06, 95.72it/s]
Running TRAM estimate - [inc: 1.3e-03]: 37%|###6 | 366/1000 [00:04<00:06, 95.72it/s]
Running TRAM estimate - [inc: 1.2e-03]: 37%|###6 | 367/1000 [00:04<00:06, 95.72it/s]
Running TRAM estimate - [inc: 1.2e-03]: 37%|###6 | 368/1000 [00:04<00:06, 95.72it/s]
Running TRAM estimate - [inc: 1.2e-03]: 37%|###6 | 369/1000 [00:04<00:06, 103.70it/s]
Running TRAM estimate - [inc: 1.2e-03]: 37%|###6 | 369/1000 [00:04<00:06, 103.70it/s]
Running TRAM estimate - [inc: 1.2e-03]: 37%|###7 | 370/1000 [00:04<00:06, 103.70it/s]
Running TRAM estimate - [inc: 1.2e-03]: 37%|###7 | 371/1000 [00:04<00:06, 103.70it/s]
Running TRAM estimate - [inc: 1.2e-03]: 37%|###7 | 372/1000 [00:04<00:06, 103.70it/s]
Running TRAM estimate - [inc: 1.2e-03]: 37%|###7 | 373/1000 [00:04<00:06, 103.70it/s]
Running TRAM estimate - [inc: 1.2e-03]: 37%|###7 | 374/1000 [00:04<00:06, 103.70it/s]
Running TRAM estimate - [inc: 1.2e-03]: 38%|###7 | 375/1000 [00:04<00:06, 103.70it/s]
Running TRAM estimate - [inc: 1.2e-03]: 38%|###7 | 376/1000 [00:04<00:06, 103.70it/s]
Running TRAM estimate - [inc: 1.2e-03]: 38%|###7 | 377/1000 [00:04<00:06, 103.70it/s]
Running TRAM estimate - [inc: 1.2e-03]: 38%|###7 | 378/1000 [00:04<00:05, 103.70it/s]
Running TRAM estimate - [inc: 1.2e-03]: 38%|###7 | 379/1000 [00:04<00:05, 103.70it/s]
Running TRAM estimate - [inc: 1.2e-03]: 38%|###8 | 380/1000 [00:04<00:05, 103.70it/s]
Running TRAM estimate - [inc: 1.2e-03]: 38%|###8 | 381/1000 [00:04<00:05, 103.70it/s]
Running TRAM estimate - [inc: 1.2e-03]: 38%|###8 | 382/1000 [00:04<00:05, 109.96it/s]
Running TRAM estimate - [inc: 1.2e-03]: 38%|###8 | 382/1000 [00:04<00:05, 109.96it/s]
Running TRAM estimate - [inc: 1.2e-03]: 38%|###8 | 383/1000 [00:04<00:05, 109.96it/s]
Running TRAM estimate - [inc: 1.2e-03]: 38%|###8 | 384/1000 [00:04<00:05, 109.96it/s]
Running TRAM estimate - [inc: 1.2e-03]: 38%|###8 | 385/1000 [00:04<00:05, 109.96it/s]
Running TRAM estimate - [inc: 1.2e-03]: 39%|###8 | 386/1000 [00:04<00:05, 109.96it/s]
Running TRAM estimate - [inc: 1.2e-03]: 39%|###8 | 387/1000 [00:04<00:05, 109.96it/s]
Running TRAM estimate - [inc: 1.2e-03]: 39%|###8 | 388/1000 [00:04<00:05, 109.96it/s]
Running TRAM estimate - [inc: 1.2e-03]: 39%|###8 | 389/1000 [00:04<00:05, 109.96it/s]
Running TRAM estimate - [inc: 1.2e-03]: 39%|###9 | 390/1000 [00:04<00:05, 109.96it/s]
Running TRAM estimate - [inc: 1.2e-03]: 39%|###9 | 391/1000 [00:04<00:05, 109.96it/s]
Running TRAM estimate - [inc: 1.2e-03]: 39%|###9 | 392/1000 [00:04<00:05, 109.96it/s]
Running TRAM estimate - [inc: 1.2e-03]: 39%|###9 | 393/1000 [00:04<00:05, 109.96it/s]
Running TRAM estimate - [inc: 1.2e-03]: 39%|###9 | 394/1000 [00:04<00:05, 109.96it/s]
Running TRAM estimate - [inc: 1.2e-03]: 40%|###9 | 395/1000 [00:04<00:05, 113.55it/s]
Running TRAM estimate - [inc: 1.2e-03]: 40%|###9 | 395/1000 [00:04<00:05, 113.55it/s]
Running TRAM estimate - [inc: 1.2e-03]: 40%|###9 | 396/1000 [00:04<00:05, 113.55it/s]
Running TRAM estimate - [inc: 1.2e-03]: 40%|###9 | 397/1000 [00:04<00:05, 113.55it/s]
Running TRAM estimate - [inc: 1.2e-03]: 40%|###9 | 398/1000 [00:04<00:05, 113.55it/s]
Running TRAM estimate - [inc: 1.2e-03]: 40%|###9 | 399/1000 [00:04<00:05, 113.55it/s]
Running TRAM estimate - [inc: 1.2e-03]: 40%|#### | 400/1000 [00:04<00:05, 113.55it/s]
Running TRAM estimate - [inc: 1.2e-03]: 40%|#### | 401/1000 [00:04<00:05, 113.55it/s]
Running TRAM estimate - [inc: 1.2e-03]: 40%|#### | 402/1000 [00:04<00:05, 113.55it/s]
Running TRAM estimate - [inc: 1.2e-03]: 40%|#### | 403/1000 [00:04<00:05, 113.55it/s]
Running TRAM estimate - [inc: 1.2e-03]: 40%|#### | 404/1000 [00:04<00:05, 113.55it/s]
Running TRAM estimate - [inc: 1.2e-03]: 40%|#### | 405/1000 [00:04<00:05, 113.55it/s]
Running TRAM estimate - [inc: 1.2e-03]: 41%|#### | 406/1000 [00:04<00:05, 113.55it/s]
Running TRAM estimate - [inc: 1.2e-03]: 41%|#### | 407/1000 [00:04<00:05, 113.55it/s]
Running TRAM estimate - [inc: 1.2e-03]: 41%|#### | 408/1000 [00:04<00:05, 115.72it/s]
Running TRAM estimate - [inc: 1.2e-03]: 41%|#### | 408/1000 [00:04<00:05, 115.72it/s]
Running TRAM estimate - [inc: 1.2e-03]: 41%|#### | 409/1000 [00:04<00:05, 115.72it/s]
Running TRAM estimate - [inc: 1.2e-03]: 41%|####1 | 410/1000 [00:04<00:05, 115.72it/s]
Running TRAM estimate - [inc: 1.2e-03]: 41%|####1 | 411/1000 [00:04<00:05, 115.72it/s]
Running TRAM estimate - [inc: 1.2e-03]: 41%|####1 | 412/1000 [00:04<00:05, 115.72it/s]
Running TRAM estimate - [inc: 1.2e-03]: 41%|####1 | 413/1000 [00:04<00:05, 115.72it/s]
Running TRAM estimate - [inc: 1.2e-03]: 41%|####1 | 414/1000 [00:04<00:05, 115.72it/s]
Running TRAM estimate - [inc: 1.2e-03]: 42%|####1 | 415/1000 [00:04<00:05, 115.72it/s]
Running TRAM estimate - [inc: 1.2e-03]: 42%|####1 | 416/1000 [00:04<00:05, 115.72it/s]
Running TRAM estimate - [inc: 1.2e-03]: 42%|####1 | 417/1000 [00:04<00:05, 115.72it/s]
Running TRAM estimate - [inc: 1.2e-03]: 42%|####1 | 418/1000 [00:04<00:05, 115.72it/s]
Running TRAM estimate - [inc: 1.2e-03]: 42%|####1 | 419/1000 [00:04<00:05, 115.72it/s]
Running TRAM estimate - [inc: 1.2e-03]: 42%|####2 | 420/1000 [00:04<00:05, 105.67it/s]
Running TRAM estimate - [inc: 1.2e-03]: 42%|####2 | 420/1000 [00:04<00:05, 105.67it/s]
Running TRAM estimate - [inc: 1.2e-03]: 42%|####2 | 421/1000 [00:04<00:05, 105.67it/s]
Running TRAM estimate - [inc: 1.2e-03]: 42%|####2 | 422/1000 [00:04<00:05, 105.67it/s]
Running TRAM estimate - [inc: 1.2e-03]: 42%|####2 | 423/1000 [00:04<00:05, 105.67it/s]
Running TRAM estimate - [inc: 1.2e-03]: 42%|####2 | 424/1000 [00:04<00:05, 105.67it/s]
Running TRAM estimate - [inc: 1.2e-03]: 42%|####2 | 425/1000 [00:04<00:05, 105.67it/s]
Running TRAM estimate - [inc: 1.2e-03]: 43%|####2 | 426/1000 [00:04<00:05, 105.67it/s]
Running TRAM estimate - [inc: 1.2e-03]: 43%|####2 | 427/1000 [00:04<00:05, 105.67it/s]
Running TRAM estimate - [inc: 1.2e-03]: 43%|####2 | 428/1000 [00:04<00:05, 105.67it/s]
Running TRAM estimate - [inc: 1.2e-03]: 43%|####2 | 429/1000 [00:04<00:05, 105.67it/s]
Running TRAM estimate - [inc: 1.2e-03]: 43%|####3 | 430/1000 [00:04<00:05, 105.67it/s]
Running TRAM estimate - [inc: 1.2e-03]: 43%|####3 | 431/1000 [00:04<00:05, 105.67it/s]
Running TRAM estimate - [inc: 1.1e-03]: 43%|####3 | 432/1000 [00:04<00:05, 105.67it/s]
Running TRAM estimate - [inc: 1.1e-03]: 43%|####3 | 433/1000 [00:04<00:05, 111.32it/s]
Running TRAM estimate - [inc: 1.1e-03]: 43%|####3 | 433/1000 [00:04<00:05, 111.32it/s]
Running TRAM estimate - [inc: 1.1e-03]: 43%|####3 | 434/1000 [00:04<00:05, 111.32it/s]
Running TRAM estimate - [inc: 1.1e-03]: 44%|####3 | 435/1000 [00:04<00:05, 111.32it/s]
Running TRAM estimate - [inc: 1.1e-03]: 44%|####3 | 436/1000 [00:04<00:05, 111.32it/s]
Running TRAM estimate - [inc: 1.1e-03]: 44%|####3 | 437/1000 [00:04<00:05, 111.32it/s]
Running TRAM estimate - [inc: 1.1e-03]: 44%|####3 | 438/1000 [00:04<00:05, 111.32it/s]
Running TRAM estimate - [inc: 1.1e-03]: 44%|####3 | 439/1000 [00:04<00:05, 111.32it/s]
Running TRAM estimate - [inc: 1.1e-03]: 44%|####4 | 440/1000 [00:04<00:05, 111.32it/s]
Running TRAM estimate - [inc: 1.1e-03]: 44%|####4 | 441/1000 [00:04<00:05, 111.32it/s]
Running TRAM estimate - [inc: 1.1e-03]: 44%|####4 | 442/1000 [00:04<00:05, 111.32it/s]
Running TRAM estimate - [inc: 1.1e-03]: 44%|####4 | 443/1000 [00:04<00:05, 111.32it/s]
Running TRAM estimate - [inc: 1.1e-03]: 44%|####4 | 444/1000 [00:04<00:04, 111.32it/s]
Running TRAM estimate - [inc: 1.1e-03]: 44%|####4 | 445/1000 [00:04<00:04, 111.32it/s]
Running TRAM estimate - [inc: 1.1e-03]: 45%|####4 | 446/1000 [00:04<00:04, 115.42it/s]
Running TRAM estimate - [inc: 1.1e-03]: 45%|####4 | 446/1000 [00:04<00:04, 115.42it/s]
Running TRAM estimate - [inc: 1.1e-03]: 45%|####4 | 447/1000 [00:04<00:04, 115.42it/s]
Running TRAM estimate - [inc: 1.1e-03]: 45%|####4 | 448/1000 [00:04<00:04, 115.42it/s]
Running TRAM estimate - [inc: 1.1e-03]: 45%|####4 | 449/1000 [00:04<00:04, 115.42it/s]
Running TRAM estimate - [inc: 1.1e-03]: 45%|####5 | 450/1000 [00:04<00:04, 115.42it/s]
Running TRAM estimate - [inc: 1.1e-03]: 45%|####5 | 451/1000 [00:04<00:04, 115.42it/s]
Running TRAM estimate - [inc: 1.1e-03]: 45%|####5 | 452/1000 [00:04<00:04, 115.42it/s]
Running TRAM estimate - [inc: 1.1e-03]: 45%|####5 | 453/1000 [00:04<00:04, 115.42it/s]
Running TRAM estimate - [inc: 1.1e-03]: 45%|####5 | 454/1000 [00:04<00:04, 115.42it/s]
Running TRAM estimate - [inc: 1.1e-03]: 46%|####5 | 455/1000 [00:04<00:04, 115.42it/s]
Running TRAM estimate - [inc: 1.1e-03]: 46%|####5 | 456/1000 [00:04<00:04, 115.42it/s]
Running TRAM estimate - [inc: 1.1e-03]: 46%|####5 | 457/1000 [00:04<00:04, 115.42it/s]
Running TRAM estimate - [inc: 1.1e-03]: 46%|####5 | 458/1000 [00:04<00:04, 115.42it/s]
Running TRAM estimate - [inc: 1.1e-03]: 46%|####5 | 459/1000 [00:04<00:04, 118.08it/s]
Running TRAM estimate - [inc: 1.1e-03]: 46%|####5 | 459/1000 [00:04<00:04, 118.08it/s]
Running TRAM estimate - [inc: 1.1e-03]: 46%|####6 | 460/1000 [00:04<00:04, 118.08it/s]
Running TRAM estimate - [inc: 1.1e-03]: 46%|####6 | 461/1000 [00:04<00:04, 118.08it/s]
Running TRAM estimate - [inc: 1.1e-03]: 46%|####6 | 462/1000 [00:04<00:04, 118.08it/s]
Running TRAM estimate - [inc: 1.1e-03]: 46%|####6 | 463/1000 [00:04<00:04, 118.08it/s]
Running TRAM estimate - [inc: 1.1e-03]: 46%|####6 | 464/1000 [00:04<00:04, 118.08it/s]
Running TRAM estimate - [inc: 1.1e-03]: 46%|####6 | 465/1000 [00:04<00:04, 118.08it/s]
Running TRAM estimate - [inc: 1.1e-03]: 47%|####6 | 466/1000 [00:04<00:04, 118.08it/s]
Running TRAM estimate - [inc: 1.1e-03]: 47%|####6 | 467/1000 [00:04<00:04, 118.08it/s]
Running TRAM estimate - [inc: 1.1e-03]: 47%|####6 | 468/1000 [00:04<00:04, 118.08it/s]
Running TRAM estimate - [inc: 1.1e-03]: 47%|####6 | 469/1000 [00:04<00:04, 118.08it/s]
Running TRAM estimate - [inc: 1.1e-03]: 47%|####6 | 470/1000 [00:04<00:04, 118.08it/s]
Running TRAM estimate - [inc: 1.1e-03]: 47%|####7 | 471/1000 [00:04<00:04, 118.08it/s]
Running TRAM estimate - [inc: 1.1e-03]: 47%|####7 | 472/1000 [00:04<00:04, 121.04it/s]
Running TRAM estimate - [inc: 1.1e-03]: 47%|####7 | 472/1000 [00:04<00:04, 121.04it/s]
Running TRAM estimate - [inc: 1.1e-03]: 47%|####7 | 473/1000 [00:04<00:04, 121.04it/s]
Running TRAM estimate - [inc: 1.1e-03]: 47%|####7 | 474/1000 [00:04<00:04, 121.04it/s]
Running TRAM estimate - [inc: 1.1e-03]: 48%|####7 | 475/1000 [00:04<00:04, 121.04it/s]
Running TRAM estimate - [inc: 1.1e-03]: 48%|####7 | 476/1000 [00:04<00:04, 121.04it/s]
Running TRAM estimate - [inc: 1.1e-03]: 48%|####7 | 477/1000 [00:04<00:04, 121.04it/s]
Running TRAM estimate - [inc: 1.1e-03]: 48%|####7 | 478/1000 [00:04<00:04, 121.04it/s]
Running TRAM estimate - [inc: 1.1e-03]: 48%|####7 | 479/1000 [00:04<00:04, 121.04it/s]
Running TRAM estimate - [inc: 1.1e-03]: 48%|####8 | 480/1000 [00:04<00:04, 121.04it/s]
Running TRAM estimate - [inc: 1.1e-03]: 48%|####8 | 481/1000 [00:04<00:04, 121.04it/s]
Running TRAM estimate - [inc: 1.1e-03]: 48%|####8 | 482/1000 [00:04<00:04, 121.04it/s]
Running TRAM estimate - [inc: 1.1e-03]: 48%|####8 | 483/1000 [00:05<00:04, 121.04it/s]
Running TRAM estimate - [inc: 1.1e-03]: 48%|####8 | 484/1000 [00:05<00:04, 121.04it/s]
Running TRAM estimate - [inc: 1.1e-03]: 48%|####8 | 485/1000 [00:05<00:04, 107.94it/s]
Running TRAM estimate - [inc: 1.1e-03]: 48%|####8 | 485/1000 [00:05<00:04, 107.94it/s]
Running TRAM estimate - [inc: 1.1e-03]: 49%|####8 | 486/1000 [00:05<00:04, 107.94it/s]
Running TRAM estimate - [inc: 1.1e-03]: 49%|####8 | 487/1000 [00:05<00:04, 107.94it/s]
Running TRAM estimate - [inc: 1.1e-03]: 49%|####8 | 488/1000 [00:05<00:04, 107.94it/s]
Running TRAM estimate - [inc: 1.1e-03]: 49%|####8 | 489/1000 [00:05<00:04, 107.94it/s]
Running TRAM estimate - [inc: 1.1e-03]: 49%|####9 | 490/1000 [00:05<00:04, 107.94it/s]
Running TRAM estimate - [inc: 1.1e-03]: 49%|####9 | 491/1000 [00:05<00:04, 107.94it/s]
Running TRAM estimate - [inc: 1.1e-03]: 49%|####9 | 492/1000 [00:05<00:04, 107.94it/s]
Running TRAM estimate - [inc: 1.1e-03]: 49%|####9 | 493/1000 [00:05<00:04, 107.94it/s]
Running TRAM estimate - [inc: 1.1e-03]: 49%|####9 | 494/1000 [00:05<00:04, 107.94it/s]
Running TRAM estimate - [inc: 1.1e-03]: 50%|####9 | 495/1000 [00:05<00:04, 107.94it/s]
Running TRAM estimate - [inc: 1.0e-03]: 50%|####9 | 496/1000 [00:05<00:04, 107.94it/s]
Running TRAM estimate - [inc: 1.0e-03]: 50%|####9 | 497/1000 [00:05<00:05, 94.94it/s]
Running TRAM estimate - [inc: 1.0e-03]: 50%|####9 | 497/1000 [00:05<00:05, 94.94it/s]
Running TRAM estimate - [inc: 1.0e-03]: 50%|####9 | 498/1000 [00:05<00:05, 94.94it/s]
Running TRAM estimate - [inc: 1.0e-03]: 50%|####9 | 499/1000 [00:05<00:05, 94.94it/s]
Running TRAM estimate - [inc: 1.0e-03]: 50%|##### | 500/1000 [00:05<00:05, 94.94it/s]
Running TRAM estimate - [inc: 1.0e-03]: 50%|##### | 501/1000 [00:05<00:05, 94.94it/s]
Running TRAM estimate - [inc: 1.0e-03]: 50%|##### | 502/1000 [00:05<00:05, 94.94it/s]
Running TRAM estimate - [inc: 1.0e-03]: 50%|##### | 503/1000 [00:05<00:05, 94.94it/s]
Running TRAM estimate - [inc: 1.0e-03]: 50%|##### | 504/1000 [00:05<00:05, 94.94it/s]
Running TRAM estimate - [inc: 1.0e-03]: 50%|##### | 505/1000 [00:05<00:05, 94.94it/s]
Running TRAM estimate - [inc: 1.0e-03]: 51%|##### | 506/1000 [00:05<00:05, 94.94it/s]
Running TRAM estimate - [inc: 1.0e-03]: 51%|##### | 507/1000 [00:05<00:05, 94.94it/s]
Running TRAM estimate - [inc: 1.0e-03]: 51%|##### | 508/1000 [00:05<00:05, 95.84it/s]
Running TRAM estimate - [inc: 1.0e-03]: 51%|##### | 508/1000 [00:05<00:05, 95.84it/s]
Running TRAM estimate - [inc: 1.0e-03]: 51%|##### | 509/1000 [00:05<00:05, 95.84it/s]
Running TRAM estimate - [inc: 1.0e-03]: 51%|#####1 | 510/1000 [00:05<00:05, 95.84it/s]
Running TRAM estimate - [inc: 1.0e-03]: 51%|#####1 | 511/1000 [00:05<00:05, 95.84it/s]
Running TRAM estimate - [inc: 1.0e-03]: 51%|#####1 | 512/1000 [00:05<00:05, 95.84it/s]
Running TRAM estimate - [inc: 1.0e-03]: 51%|#####1 | 513/1000 [00:05<00:05, 95.84it/s]
Running TRAM estimate - [inc: 1.0e-03]: 51%|#####1 | 514/1000 [00:05<00:05, 95.84it/s]
Running TRAM estimate - [inc: 1.0e-03]: 52%|#####1 | 515/1000 [00:05<00:05, 95.84it/s]
Running TRAM estimate - [inc: 1.0e-03]: 52%|#####1 | 516/1000 [00:05<00:05, 95.84it/s]
Running TRAM estimate - [inc: 1.0e-03]: 52%|#####1 | 517/1000 [00:05<00:05, 95.84it/s]
Running TRAM estimate - [inc: 1.0e-03]: 52%|#####1 | 518/1000 [00:05<00:05, 95.84it/s]
Running TRAM estimate - [inc: 1.0e-03]: 52%|#####1 | 519/1000 [00:05<00:05, 95.84it/s]
Running TRAM estimate - [inc: 1.0e-03]: 52%|#####2 | 520/1000 [00:05<00:05, 95.84it/s]
Running TRAM estimate - [inc: 1.0e-03]: 52%|#####2 | 521/1000 [00:05<00:04, 103.58it/s]
Running TRAM estimate - [inc: 1.0e-03]: 52%|#####2 | 521/1000 [00:05<00:04, 103.58it/s]
Running TRAM estimate - [inc: 1.0e-03]: 52%|#####2 | 522/1000 [00:05<00:04, 103.58it/s]
Running TRAM estimate - [inc: 1.0e-03]: 52%|#####2 | 523/1000 [00:05<00:04, 103.58it/s]
Running TRAM estimate - [inc: 1.0e-03]: 52%|#####2 | 524/1000 [00:05<00:04, 103.58it/s]
Running TRAM estimate - [inc: 1.0e-03]: 52%|#####2 | 525/1000 [00:05<00:04, 103.58it/s]
Running TRAM estimate - [inc: 1.0e-03]: 53%|#####2 | 526/1000 [00:05<00:04, 103.58it/s]
Running TRAM estimate - [inc: 1.0e-03]: 53%|#####2 | 527/1000 [00:05<00:04, 103.58it/s]
Running TRAM estimate - [inc: 1.0e-03]: 53%|#####2 | 528/1000 [00:05<00:04, 103.58it/s]
Running TRAM estimate - [inc: 1.0e-03]: 53%|#####2 | 529/1000 [00:05<00:04, 103.58it/s]
Running TRAM estimate - [inc: 1.0e-03]: 100%|##########| 529/529 [00:05<00:00, 96.01it/s]
12 import numpy as np
13 import matplotlib.pyplot as plt
14
15 from deeptime.data import tmatrix_metropolis1d
16 from deeptime.markov.msm import MarkovStateModel, TRAM
17 from deeptime.clustering import ClusterModel
18
19 xs = np.linspace(-1.5, 1.5, num=100)
20 n_samples = 10000
21 bias_centers = [-1, -0.5, 0.0, 0.5, 1]
22
23
24 def harmonic(x0, x):
25 return 2 * (x - x0) ** 4
26
27
28 def plot_contour_with_colourbar(data, vmin=None, vmax=None):
29 if vmin is None:
30 vmin = np.min(data)
31 if vmax is None:
32 vmax = np.max(data)
33 fig, (ax1) = plt.subplots(1, figsize=(3.5, 3))
34 im = ax1.contourf(data, vmin=vmin, vmax=vmax, levels=50, cmap='jet')
35 plt.colorbar(im)
36 plt.show()
37
38
39 def get_bias_functions():
40 bias_functions = []
41 for i, bias_center in enumerate(bias_centers):
42 bias = lambda x, x0=bias_center: harmonic(x0, x)
43 bias_functions.append(bias)
44 return bias_functions
45
46
47 def sample_trajectories(bias_functions):
48 trajs = np.zeros((len(bias_centers), n_samples), dtype=np.int32)
49
50 for i, bias in enumerate(bias_functions):
51 biased_energies = (xs - 1) ** 4 * (xs + 1) ** 4 - 0.1 * xs + bias(xs)
52
53 biased_energies /= np.max(biased_energies)
54 transition_matrix = tmatrix_metropolis1d(biased_energies)
55
56 msm = MarkovStateModel(transition_matrix)
57 trajs[i] = msm.simulate(n_steps=n_samples)
58 return trajs
59
60
61 if __name__ == "__main__":
62 bias_functions = get_bias_functions()
63 trajectories = sample_trajectories(bias_functions)
64
65 # move from trajectory over 100 bins back to the space of the xs: (-1.5, 1.5)
66 trajectories = trajectories / 100 * 3 - 1.5
67
68 bias_matrices = np.zeros((len(bias_centers), n_samples, len(bias_centers)))
69 for i, traj in enumerate(trajectories):
70 for j, bias_function in enumerate(bias_functions):
71 bias_matrices[i, :, j] = bias_function(traj)
72
73 # discretize the trajectories into two Markov states (centered around the two wells)
74 clustering = ClusterModel(cluster_centers=np.asarray([-0.75, 0.75]), metric='euclidean')
75
76 dtrajs = clustering.transform(trajectories.flatten()).reshape((len(bias_matrices), n_samples))
77
78 from tqdm import tqdm
79 tram = TRAM(lagtime=1, maxiter=1000, maxerr=1e-3, progress=tqdm, init_strategy="MBAR")
80
81 # For every simulation frame seen in trajectory i and time step t, btrajs[i][t,k] is the
82 # bias energy of that frame evaluated in the k'th thermodynamic state (i.e. at the k'th
83 # Umbrella/Hamiltonian/temperature).
84 model = tram.fit_fetch((dtrajs, bias_matrices))
85
86 plot_contour_with_colourbar(model.biased_conf_energies)
Total running time of the script: ( 0 minutes 6.523 seconds)
Estimated memory usage: 8 MB