Backward toy: Bayesian inversion with Active-MCMC (documentation notebook)¶
This notebook-style tutorial demonstrates how to use gp_active_mcmc on a small toy
inverse problem.
We show two inference modes supported by the library:
MCMC-guided active learning (single posterior) The sampler uses only the coarse model (
ActiveMCMCModel.coarse). The model decides internally when to fall back to the high-fidelity (HF) forward model based on a simple uncertainty trigger.DA-MCMC guided active learning with an adaptive subchain (two posteriors, recommended) The sampler uses a coarse posterior and a fine posterior (delayed acceptance). An
AdaptiveSubchainmonitors LF-HF discrepancy during HF corrections and adjusts how often the fine level is called. BecausetinyDAuses a fixed subsampling rate per call, the run is performed in chunks.
Along the way we:
- define a Gaussian prior and generate a synthetic observation,
- train a POD-GP surrogate on an initial design,
- wrap surrogate + HF model into an
ActiveMCMCModel, - run short chains and plot lightweight diagnostics.
The budgets are intentionally small so the tutorial executes quickly in documentation builds.
In [1]:
Copied!
from __future__ import annotations
import copy
import numpy as np
import tinyDA as tda
from scipy.stats import multivariate_normal
# For documentation builds, import the toy helpers from the tutorial package.
# If you place this file at `docs/tutorials/backward_toy_notebook.py`, also place
# `toy.py` in the same directory so this import works.
from docs.tutorials.toy_tutorial import make_forward_model, make_observation, make_timeline
# Diagnostics return (fig, ax) and never force `plt.show()` unless `show=True`.
from gp_active_mcmc.diagnostics import (
plot_chain_2d,
plot_cumulative_hf_fraction,
plot_prediction_at_theta,
plot_subchain_length_history,
)
from gp_active_mcmc.inference import (
ActiveGPLogLike,
ActiveMCMCModel,
AdaptiveMetropolisShared,
AdaptiveSubchain,
AdaptiveSubchainControl,
AdaptiveSubchainState,
ChunkedMCMCConfig,
sample_active_chain,
sample_adaptive_active_chain,
)
from gp_active_mcmc.surrogates import POD, MultiOutputGP, PODGPSurrogate
from gp_active_mcmc.utils.rng import set_seed
from __future__ import annotations
import copy
import numpy as np
import tinyDA as tda
from scipy.stats import multivariate_normal
# For documentation builds, import the toy helpers from the tutorial package.
# If you place this file at `docs/tutorials/backward_toy_notebook.py`, also place
# `toy.py` in the same directory so this import works.
from docs.tutorials.toy_tutorial import make_forward_model, make_observation, make_timeline
# Diagnostics return (fig, ax) and never force `plt.show()` unless `show=True`.
from gp_active_mcmc.diagnostics import (
plot_chain_2d,
plot_cumulative_hf_fraction,
plot_prediction_at_theta,
plot_subchain_length_history,
)
from gp_active_mcmc.inference import (
ActiveGPLogLike,
ActiveMCMCModel,
AdaptiveMetropolisShared,
AdaptiveSubchain,
AdaptiveSubchainControl,
AdaptiveSubchainState,
ChunkedMCMCConfig,
sample_active_chain,
sample_adaptive_active_chain,
)
from gp_active_mcmc.surrogates import POD, MultiOutputGP, PODGPSurrogate
from gp_active_mcmc.utils.rng import set_seed
Configuration¶
tdefines the toy trajectory grid (the observation space).- The prior is a multivariate Normal over three parameters.
- The surrogate is trained on
n_initHF snapshots drawn from the prior. gamma_thresholdcontrols whenActiveMCMCModel.coarsetriggers HF: it compares the mean surrogate variance togamma_threshold**2.
In [2]:
Copied!
rng = set_seed(2)
# Time grid for the toy trajectory
t = make_timeline(T=500, t_end=0.05)
# High-fidelity forward model: theta -> y(t)
hf_forward = make_forward_model(t)
# Prior over theta = [A, f, tau]
prior_mean = np.array([0.8, 150.0, 0.01])
prior_cov = np.diag([0.5**2, 40.0**2, 0.01**2])
prior = multivariate_normal(mean=prior_mean, cov=prior_cov)
# Observation noise (standard deviation)
sigma_obs = 0.1
# Surrogate configuration (small budgets for docs)
n_init = 50
pod_rank = 20
gp_kernel = "matern52"
gp_ard = True
noise_variance_gp = 1e-6
# Active coupling: trigger HF when average predictive std exceeds gamma_threshold
gamma_threshold = 0.10
# MCMC budget (expressed as coarse evaluations)
n_coarse_evals = 1000
burn_in = 100
# Chunking for the adaptive workflow
chunk_size = 250
rng = set_seed(2)
# Time grid for the toy trajectory
t = make_timeline(T=500, t_end=0.05)
# High-fidelity forward model: theta -> y(t)
hf_forward = make_forward_model(t)
# Prior over theta = [A, f, tau]
prior_mean = np.array([0.8, 150.0, 0.01])
prior_cov = np.diag([0.5**2, 40.0**2, 0.01**2])
prior = multivariate_normal(mean=prior_mean, cov=prior_cov)
# Observation noise (standard deviation)
sigma_obs = 0.1
# Surrogate configuration (small budgets for docs)
n_init = 50
pod_rank = 20
gp_kernel = "matern52"
gp_ard = True
noise_variance_gp = 1e-6
# Active coupling: trigger HF when average predictive std exceeds gamma_threshold
gamma_threshold = 0.10
# MCMC budget (expressed as coarse evaluations)
n_coarse_evals = 1000
burn_in = 100
# Chunking for the adaptive workflow
chunk_size = 250
Synthetic observation¶
We generate a synthetic "truth" parameter theta_true, evaluate the HF model, and
corrupt with Gaussian noise. This is what the likelihood will condition on.
In [3]:
Copied!
theta_true = prior.rvs(random_state=rng)
y_obs = make_observation(rng, theta_true, t, sigma_obs)
theta_true = prior.rvs(random_state=rng)
y_obs = make_observation(rng, theta_true, t, sigma_obs)
Initial surrogate training set¶
We build an initial design by sampling parameters from the prior and evaluating the HF model.
This yields training pairs (theta_i, y_i(t)) for POD and GP fitting.
In [4]:
Copied!
theta_train = prior.rvs(size=n_init, random_state=rng)
y_train = np.asarray([hf_forward(th) for th in np.atleast_2d(theta_train)], dtype=float)
theta_train = prior.rvs(size=n_init, random_state=rng)
y_train = np.asarray([hf_forward(th) for th in np.atleast_2d(theta_train)], dtype=float)
Fit a POD-GP surrogate¶
- Fit POD on snapshot matrix
y_trainto obtain a reduced basis. - Project snapshots to POD coefficients.
- Fit a multi-output GP mapping
theta -> coefficients.
In [5]:
Copied!
pod = POD(rank=pod_rank).fit(y_train)
coeffs = pod.transform(y_train) # (n_init, pod_rank)
gp = MultiOutputGP(
theta_train,
coeffs,
kernel=gp_kernel,
ard=gp_ard,
noise_variance=noise_variance_gp,
update_every=100,
n_retrain_max=0,
)
# We create two surrogate *instances* so the two runs below do not share state.
# This keeps the tutorial results easier to interpret.
lf_surrogate_single = PODGPSurrogate(pod=copy.deepcopy(pod), gp=copy.deepcopy(gp))
lf_surrogate_adapt = PODGPSurrogate(pod=copy.deepcopy(pod), gp=copy.deepcopy(gp))
pod = POD(rank=pod_rank).fit(y_train)
coeffs = pod.transform(y_train) # (n_init, pod_rank)
gp = MultiOutputGP(
theta_train,
coeffs,
kernel=gp_kernel,
ard=gp_ard,
noise_variance=noise_variance_gp,
update_every=100,
n_retrain_max=0,
)
# We create two surrogate *instances* so the two runs below do not share state.
# This keeps the tutorial results easier to interpret.
lf_surrogate_single = PODGPSurrogate(pod=copy.deepcopy(pod), gp=copy.deepcopy(gp))
lf_surrogate_adapt = PODGPSurrogate(pod=copy.deepcopy(pod), gp=copy.deepcopy(gp))
Wrap LF + HF in an ActiveMCMCModel¶
ActiveMCMCModel is the heart of the inference API. It couples:
- an LF surrogate (
predict,update) - an HF forward model (
theta -> y)
It exposes two callables intended for MCMC:
coarse(theta): LF-first; triggers HF if uncertainty is too large.fine(theta): always HF; updates the surrogate.
The choice of posteriors determines the inference mode:
- single posterior (
Posterior(..., model.coarse)) → MCMC-guided active learning - two posteriors (
[Posterior(..., model.coarse), Posterior(..., model.fine)]) → DA-MCMC guided AL
In the adaptive subchain workflow, DA-MCMC is mandatory because subchains describe how often we attempt fine (HF) corrections relative to coarse steps.
In [6]:
Copied!
model_single = ActiveMCMCModel(
lf_model=lf_surrogate_single,
hf_model=hf_forward,
gamma_threshold=gamma_threshold,
)
adaptive_policy = AdaptiveSubchain(
state=AdaptiveSubchainState(subchain_length=20),
control=AdaptiveSubchainControl(
update_every=5,
target_error=0.05,
min_subchain=1,
max_subchain=500,
),
)
model_adapt = ActiveMCMCModel(
lf_model=lf_surrogate_adapt,
hf_model=hf_forward,
gamma_threshold=gamma_threshold,
adaptive=adaptive_policy,
)
model_single = ActiveMCMCModel(
lf_model=lf_surrogate_single,
hf_model=hf_forward,
gamma_threshold=gamma_threshold,
)
adaptive_policy = AdaptiveSubchain(
state=AdaptiveSubchainState(subchain_length=20),
control=AdaptiveSubchainControl(
update_every=5,
target_error=0.05,
min_subchain=1,
max_subchain=500,
),
)
model_adapt = ActiveMCMCModel(
lf_model=lf_surrogate_adapt,
hf_model=hf_forward,
gamma_threshold=gamma_threshold,
adaptive=adaptive_policy,
)
Likelihood and posterior objects¶
ActiveGPLogLike inflates the observation covariance when the model output carries a
predictive variance (the LF surrogate case).
For the fine level we use tinyDA.AdaptiveGaussianLogLike, which expects a plain mean vector.
In [7]:
Copied!
cov = (sigma_obs**2) * np.eye(len(y_obs))
loglike_coarse = ActiveGPLogLike(data=y_obs, covariance=cov)
loglike_fine = tda.AdaptiveGaussianLogLike(data=y_obs, covariance=cov)
# Single-level posterior (MCMC-guided active learning)
posterior_single = tda.Posterior(prior, loglike_coarse, model_single.coarse)
# Two-level posterior list (DA-MCMC)
posterior_adapt = [
tda.Posterior(prior, loglike_coarse, model_adapt.coarse),
tda.Posterior(prior, loglike_fine, model_adapt.fine),
]
cov = (sigma_obs**2) * np.eye(len(y_obs))
loglike_coarse = ActiveGPLogLike(data=y_obs, covariance=cov)
loglike_fine = tda.AdaptiveGaussianLogLike(data=y_obs, covariance=cov)
# Single-level posterior (MCMC-guided active learning)
posterior_single = tda.Posterior(prior, loglike_coarse, model_single.coarse)
# Two-level posterior list (DA-MCMC)
posterior_adapt = [
tda.Posterior(prior, loglike_coarse, model_adapt.coarse),
tda.Posterior(prior, loglike_fine, model_adapt.fine),
]
Proposal¶
We use an adaptive Metropolis proposal. The AdaptiveMetropolisShared wrapper ensures
proposal adaptation is not accidentally reset by deepcopies performed in chunked sampling.
In [8]:
Copied!
theta0 = prior_mean.copy()
proposal = AdaptiveMetropolisShared(
C0=0.1 * prior_cov,
period=100,
share_across_deepcopy=True,
adaptive=True,
sd=1,
)
theta0 = prior_mean.copy()
proposal = AdaptiveMetropolisShared(
C0=0.1 * prior_cov,
period=100,
share_across_deepcopy=True,
adaptive=True,
sd=1,
)
Sanity check: surrogate prediction before sampling¶
This plot is only meant to provide intuition: the surrogate has been trained on a small initial design and may be inaccurate away from the training region.
In [9]:
Copied!
plot_prediction_at_theta(
model=lf_surrogate_single,
theta=theta_true,
t=t,
y_obs=y_obs,
y_true=hf_forward(theta_true),
title="Surrogate prediction (before sampling)",
show=True,
)
plot_prediction_at_theta(
model=lf_surrogate_single,
theta=theta_true,
t=t,
y_obs=y_obs,
y_true=hf_forward(theta_true),
title="Surrogate prediction (before sampling)",
show=True,
)
Out[9]:
(<Figure size 800x400 with 1 Axes>,
<Axes: title={'center': 'Surrogate prediction (before sampling)'}, xlabel='t', ylabel='y'>)
Part 1 — MCMC-guided active learning (single posterior)¶
Here we pass a single posterior using model.coarse. The chain is driven by the coarse model.
HF calls occur only when the uncertainty trigger fires inside model.coarse.
In [10]:
Copied!
result_single = sample_active_chain(
model=model_single,
posterior=posterior_single,
proposal=proposal,
subsampling_rate=1,
iterations=n_coarse_evals,
initial_parameters=theta0,
chain_key="chain_0",
)
chain_single = result_single.chain
chain_single.summary(theta_true=theta_true, burn_in=burn_in)
result_single = sample_active_chain(
model=model_single,
posterior=posterior_single,
proposal=proposal,
subsampling_rate=1,
iterations=n_coarse_evals,
initial_parameters=theta0,
chain_key="chain_0",
)
chain_single = result_single.chain
chain_single.summary(theta_true=theta_true, burn_in=burn_in)
/opt/hostedtoolcache/Python/3.10.19/x64/lib/python3.10/site-packages/tinyDA/sampler.py:114: UserWarning: subsampling_rate has been deprecated in favour of subchain_length.
Sampling chain 1/1
0%| | 0/1000 [00:00<?, ?it/s]
Running chain, α = 1.00: 0%| | 0/1000 [00:00<?, ?it/s]
Running chain, α = 0.50: 0%| | 0/1000 [00:00<?, ?it/s]
Running chain, α = 0.67: 0%| | 0/1000 [00:00<?, ?it/s]
Running chain, α = 0.50: 0%| | 0/1000 [00:00<?, ?it/s]
Running chain, α = 0.40: 0%| | 0/1000 [00:00<?, ?it/s]
Running chain, α = 0.40: 0%| | 5/1000 [00:00<00:21, 46.55it/s]
Running chain, α = 0.33: 0%| | 5/1000 [00:00<00:21, 46.55it/s]
Running chain, α = 0.29: 0%| | 5/1000 [00:00<00:21, 46.55it/s]
Running chain, α = 0.25: 0%| | 5/1000 [00:00<00:21, 46.55it/s]
Running chain, α = 0.33: 0%| | 5/1000 [00:00<00:21, 46.55it/s]
Running chain, α = 0.30: 0%| | 5/1000 [00:00<00:21, 46.55it/s]
Running chain, α = 0.30: 1%| | 10/1000 [00:00<00:21, 46.63it/s]
Running chain, α = 0.27: 1%| | 10/1000 [00:00<00:21, 46.63it/s]
Running chain, α = 0.25: 1%| | 10/1000 [00:00<00:21, 46.63it/s]
Running chain, α = 0.23: 1%| | 10/1000 [00:00<00:21, 46.63it/s]
Running chain, α = 0.21: 1%| | 10/1000 [00:00<00:21, 46.63it/s]
Running chain, α = 0.20: 1%| | 10/1000 [00:00<00:21, 46.63it/s]
Running chain, α = 0.20: 2%|▏ | 15/1000 [00:00<00:21, 46.48it/s]
Running chain, α = 0.19: 2%|▏ | 15/1000 [00:00<00:21, 46.48it/s]
Running chain, α = 0.18: 2%|▏ | 15/1000 [00:00<00:21, 46.48it/s]
Running chain, α = 0.17: 2%|▏ | 15/1000 [00:00<00:21, 46.48it/s]
Running chain, α = 0.16: 2%|▏ | 15/1000 [00:00<00:21, 46.48it/s]
Running chain, α = 0.15: 2%|▏ | 15/1000 [00:00<00:21, 46.48it/s]
Running chain, α = 0.15: 2%|▏ | 20/1000 [00:00<00:21, 46.13it/s]
Running chain, α = 0.14: 2%|▏ | 20/1000 [00:00<00:21, 46.13it/s]
Running chain, α = 0.18: 2%|▏ | 20/1000 [00:00<00:21, 46.13it/s]
Running chain, α = 0.17: 2%|▏ | 20/1000 [00:00<00:21, 46.13it/s]
Running chain, α = 0.17: 2%|▏ | 20/1000 [00:00<00:21, 46.13it/s]
Running chain, α = 0.16: 2%|▏ | 20/1000 [00:00<00:21, 46.13it/s]
Running chain, α = 0.16: 2%|▎ | 25/1000 [00:00<00:29, 33.58it/s]
Running chain, α = 0.15: 2%|▎ | 25/1000 [00:00<00:29, 33.58it/s]
Running chain, α = 0.19: 2%|▎ | 25/1000 [00:00<00:29, 33.58it/s]
Running chain, α = 0.18: 2%|▎ | 25/1000 [00:00<00:29, 33.58it/s]
Running chain, α = 0.17: 2%|▎ | 25/1000 [00:00<00:29, 33.58it/s]
Running chain, α = 0.17: 3%|▎ | 29/1000 [00:00<00:28, 33.50it/s]
Running chain, α = 0.17: 3%|▎ | 29/1000 [00:00<00:28, 33.50it/s]
Running chain, α = 0.16: 3%|▎ | 29/1000 [00:00<00:28, 33.50it/s]
Running chain, α = 0.16: 3%|▎ | 29/1000 [00:00<00:28, 33.50it/s]
Running chain, α = 0.18: 3%|▎ | 29/1000 [00:00<00:28, 33.50it/s]
Running chain, α = 0.18: 3%|▎ | 29/1000 [00:00<00:28, 33.50it/s]
Running chain, α = 0.18: 3%|▎ | 34/1000 [00:00<00:26, 36.56it/s]
Running chain, α = 0.17: 3%|▎ | 34/1000 [00:00<00:26, 36.56it/s]
Running chain, α = 0.17: 3%|▎ | 34/1000 [00:00<00:26, 36.56it/s]
Running chain, α = 0.16: 3%|▎ | 34/1000 [00:00<00:26, 36.56it/s]
Running chain, α = 0.16: 3%|▎ | 34/1000 [00:00<00:26, 36.56it/s]
Running chain, α = 0.15: 3%|▎ | 34/1000 [00:00<00:26, 36.56it/s]
Running chain, α = 0.15: 4%|▍ | 39/1000 [00:00<00:24, 38.73it/s]
Running chain, α = 0.15: 4%|▍ | 39/1000 [00:01<00:24, 38.73it/s]
Running chain, α = 0.15: 4%|▍ | 39/1000 [00:01<00:24, 38.73it/s]
Running chain, α = 0.14: 4%|▍ | 39/1000 [00:01<00:24, 38.73it/s]
Running chain, α = 0.14: 4%|▍ | 39/1000 [00:01<00:24, 38.73it/s]
Running chain, α = 0.14: 4%|▍ | 39/1000 [00:01<00:24, 38.73it/s]
Running chain, α = 0.14: 4%|▍ | 44/1000 [00:01<00:23, 40.61it/s]
Running chain, α = 0.13: 4%|▍ | 44/1000 [00:01<00:23, 40.61it/s]
Running chain, α = 0.13: 4%|▍ | 44/1000 [00:01<00:23, 40.61it/s]
Running chain, α = 0.13: 4%|▍ | 44/1000 [00:01<00:23, 40.61it/s]
Running chain, α = 0.12: 4%|▍ | 44/1000 [00:01<00:23, 40.61it/s]
Running chain, α = 0.14: 4%|▍ | 44/1000 [00:01<00:23, 40.61it/s]
Running chain, α = 0.14: 5%|▍ | 49/1000 [00:01<00:22, 42.22it/s]
Running chain, α = 0.14: 5%|▍ | 49/1000 [00:01<00:22, 42.22it/s]
Running chain, α = 0.14: 5%|▍ | 49/1000 [00:01<00:22, 42.22it/s]
Running chain, α = 0.13: 5%|▍ | 49/1000 [00:01<00:22, 42.22it/s]
Running chain, α = 0.13: 5%|▍ | 49/1000 [00:01<00:22, 42.22it/s]
Running chain, α = 0.13: 5%|▍ | 49/1000 [00:01<00:22, 42.22it/s]
Running chain, α = 0.13: 5%|▌ | 54/1000 [00:01<00:21, 43.42it/s]
Running chain, α = 0.13: 5%|▌ | 54/1000 [00:01<00:21, 43.42it/s]
Running chain, α = 0.12: 5%|▌ | 54/1000 [00:01<00:21, 43.42it/s]
Running chain, α = 0.12: 5%|▌ | 54/1000 [00:01<00:21, 43.42it/s]
Running chain, α = 0.12: 5%|▌ | 54/1000 [00:01<00:21, 43.42it/s]
Running chain, α = 0.12: 5%|▌ | 54/1000 [00:01<00:21, 43.42it/s]
Running chain, α = 0.12: 6%|▌ | 59/1000 [00:01<00:21, 44.34it/s]
Running chain, α = 0.12: 6%|▌ | 59/1000 [00:01<00:21, 44.34it/s]
Running chain, α = 0.11: 6%|▌ | 59/1000 [00:01<00:21, 44.34it/s]
Running chain, α = 0.11: 6%|▌ | 59/1000 [00:01<00:21, 44.34it/s]
Running chain, α = 0.11: 6%|▌ | 59/1000 [00:01<00:21, 44.34it/s]
Running chain, α = 0.11: 6%|▌ | 59/1000 [00:01<00:21, 44.34it/s]
Running chain, α = 0.11: 6%|▋ | 64/1000 [00:01<00:20, 44.94it/s]
Running chain, α = 0.11: 6%|▋ | 64/1000 [00:01<00:20, 44.94it/s]
Running chain, α = 0.11: 6%|▋ | 64/1000 [00:01<00:20, 44.94it/s]
Running chain, α = 0.10: 6%|▋ | 64/1000 [00:01<00:20, 44.94it/s]
Running chain, α = 0.10: 6%|▋ | 64/1000 [00:01<00:20, 44.94it/s]
Running chain, α = 0.10: 6%|▋ | 64/1000 [00:01<00:20, 44.94it/s]
Running chain, α = 0.10: 7%|▋ | 69/1000 [00:01<00:20, 45.38it/s]
Running chain, α = 0.10: 7%|▋ | 69/1000 [00:01<00:20, 45.38it/s]
Running chain, α = 0.10: 7%|▋ | 69/1000 [00:01<00:20, 45.38it/s]
Running chain, α = 0.10: 7%|▋ | 69/1000 [00:01<00:20, 45.38it/s]
Running chain, α = 0.10: 7%|▋ | 69/1000 [00:01<00:20, 45.38it/s]
Running chain, α = 0.09: 7%|▋ | 69/1000 [00:01<00:20, 45.38it/s]
Running chain, α = 0.09: 7%|▋ | 74/1000 [00:01<00:20, 45.63it/s]
Running chain, α = 0.09: 7%|▋ | 74/1000 [00:01<00:20, 45.63it/s]
Running chain, α = 0.09: 7%|▋ | 74/1000 [00:01<00:20, 45.63it/s]
Running chain, α = 0.09: 7%|▋ | 74/1000 [00:01<00:20, 45.63it/s]
Running chain, α = 0.09: 7%|▋ | 74/1000 [00:01<00:20, 45.63it/s]
Running chain, α = 0.09: 7%|▋ | 74/1000 [00:01<00:20, 45.63it/s]
Running chain, α = 0.09: 8%|▊ | 79/1000 [00:01<00:20, 45.87it/s]
Running chain, α = 0.09: 8%|▊ | 79/1000 [00:01<00:20, 45.87it/s]
Running chain, α = 0.09: 8%|▊ | 79/1000 [00:01<00:20, 45.87it/s]
Running chain, α = 0.09: 8%|▊ | 79/1000 [00:01<00:20, 45.87it/s]
Running chain, α = 0.08: 8%|▊ | 79/1000 [00:01<00:20, 45.87it/s]
Running chain, α = 0.08: 8%|▊ | 79/1000 [00:01<00:20, 45.87it/s]
Running chain, α = 0.08: 8%|▊ | 84/1000 [00:01<00:20, 44.94it/s]
Running chain, α = 0.08: 8%|▊ | 84/1000 [00:01<00:20, 44.94it/s]
Running chain, α = 0.08: 8%|▊ | 84/1000 [00:02<00:20, 44.94it/s]
Running chain, α = 0.08: 8%|▊ | 84/1000 [00:02<00:20, 44.94it/s]
Running chain, α = 0.08: 8%|▊ | 84/1000 [00:02<00:20, 44.94it/s]
Running chain, α = 0.08: 8%|▊ | 84/1000 [00:02<00:20, 44.94it/s]
Running chain, α = 0.08: 9%|▉ | 89/1000 [00:02<00:20, 45.51it/s]
Running chain, α = 0.08: 9%|▉ | 89/1000 [00:02<00:20, 45.51it/s]
Running chain, α = 0.08: 9%|▉ | 89/1000 [00:02<00:20, 45.51it/s]
Running chain, α = 0.08: 9%|▉ | 89/1000 [00:02<00:20, 45.51it/s]
Running chain, α = 0.08: 9%|▉ | 89/1000 [00:02<00:20, 45.51it/s]
Running chain, α = 0.07: 9%|▉ | 89/1000 [00:02<00:20, 45.51it/s]
Running chain, α = 0.07: 9%|▉ | 94/1000 [00:02<00:19, 45.59it/s]
Running chain, α = 0.07: 9%|▉ | 94/1000 [00:02<00:19, 45.59it/s]
Running chain, α = 0.07: 9%|▉ | 94/1000 [00:02<00:19, 45.59it/s]
Running chain, α = 0.07: 9%|▉ | 94/1000 [00:02<00:19, 45.59it/s]
Running chain, α = 0.07: 9%|▉ | 94/1000 [00:02<00:19, 45.59it/s]
Running chain, α = 0.07: 9%|▉ | 94/1000 [00:02<00:19, 45.59it/s]
Running chain, α = 0.07: 10%|▉ | 99/1000 [00:02<00:19, 45.82it/s]
Running chain, α = 0.07: 10%|▉ | 99/1000 [00:02<00:19, 45.82it/s]
Running chain, α = 0.06: 10%|▉ | 99/1000 [00:02<00:19, 45.82it/s]
Running chain, α = 0.07: 10%|▉ | 99/1000 [00:02<00:19, 45.82it/s]
Running chain, α = 0.06: 10%|▉ | 99/1000 [00:02<00:19, 45.82it/s]
Running chain, α = 0.06: 10%|▉ | 99/1000 [00:02<00:19, 45.82it/s]
Running chain, α = 0.06: 10%|█ | 104/1000 [00:02<00:19, 45.82it/s]
Running chain, α = 0.07: 10%|█ | 104/1000 [00:02<00:19, 45.82it/s]
Running chain, α = 0.08: 10%|█ | 104/1000 [00:02<00:19, 45.82it/s]
Running chain, α = 0.08: 10%|█ | 104/1000 [00:02<00:19, 45.82it/s]
Running chain, α = 0.08: 10%|█ | 104/1000 [00:02<00:19, 45.82it/s]
Running chain, α = 0.08: 10%|█ | 104/1000 [00:02<00:19, 45.82it/s]
Running chain, α = 0.08: 11%|█ | 109/1000 [00:02<00:19, 45.60it/s]
Running chain, α = 0.08: 11%|█ | 109/1000 [00:02<00:19, 45.60it/s]
Running chain, α = 0.08: 11%|█ | 109/1000 [00:02<00:19, 45.60it/s]
Running chain, α = 0.08: 11%|█ | 109/1000 [00:02<00:19, 45.60it/s]
Running chain, α = 0.09: 11%|█ | 109/1000 [00:02<00:19, 45.60it/s]
Running chain, α = 0.09: 11%|█ | 109/1000 [00:02<00:19, 45.60it/s]
Running chain, α = 0.09: 11%|█▏ | 114/1000 [00:02<00:19, 45.83it/s]
Running chain, α = 0.10: 11%|█▏ | 114/1000 [00:02<00:19, 45.83it/s]
Running chain, α = 0.10: 11%|█▏ | 114/1000 [00:02<00:19, 45.83it/s]
Running chain, α = 0.10: 11%|█▏ | 114/1000 [00:02<00:19, 45.83it/s]
Running chain, α = 0.10: 11%|█▏ | 114/1000 [00:02<00:19, 45.83it/s]
Running chain, α = 0.10: 11%|█▏ | 114/1000 [00:02<00:19, 45.83it/s]
Running chain, α = 0.10: 12%|█▏ | 119/1000 [00:02<00:19, 45.93it/s]
Running chain, α = 0.11: 12%|█▏ | 119/1000 [00:02<00:19, 45.93it/s]
Running chain, α = 0.11: 12%|█▏ | 119/1000 [00:02<00:19, 45.93it/s]
Running chain, α = 0.10: 12%|█▏ | 119/1000 [00:02<00:19, 45.93it/s]
Running chain, α = 0.10: 12%|█▏ | 119/1000 [00:02<00:19, 45.93it/s]
Running chain, α = 0.10: 12%|█▏ | 119/1000 [00:02<00:19, 45.93it/s]
Running chain, α = 0.10: 12%|█▏ | 124/1000 [00:02<00:19, 46.01it/s]
Running chain, α = 0.10: 12%|█▏ | 124/1000 [00:02<00:19, 46.01it/s]
Running chain, α = 0.11: 12%|█▏ | 124/1000 [00:02<00:19, 46.01it/s]
Running chain, α = 0.10: 12%|█▏ | 124/1000 [00:02<00:19, 46.01it/s]
Running chain, α = 0.10: 12%|█▏ | 124/1000 [00:02<00:19, 46.01it/s]
Running chain, α = 0.10: 12%|█▏ | 124/1000 [00:02<00:19, 46.01it/s]
Running chain, α = 0.10: 13%|█▎ | 129/1000 [00:02<00:19, 45.39it/s]
Running chain, α = 0.10: 13%|█▎ | 129/1000 [00:02<00:19, 45.39it/s]
Running chain, α = 0.11: 13%|█▎ | 129/1000 [00:02<00:19, 45.39it/s]
Running chain, α = 0.11: 13%|█▎ | 129/1000 [00:03<00:19, 45.39it/s]
Running chain, α = 0.11: 13%|█▎ | 129/1000 [00:03<00:19, 45.39it/s]
Running chain, α = 0.11: 13%|█▎ | 129/1000 [00:03<00:19, 45.39it/s]
Running chain, α = 0.11: 13%|█▎ | 134/1000 [00:03<00:18, 45.75it/s]
Running chain, α = 0.11: 13%|█▎ | 134/1000 [00:03<00:18, 45.75it/s]
Running chain, α = 0.11: 13%|█▎ | 134/1000 [00:03<00:18, 45.75it/s]
Running chain, α = 0.11: 13%|█▎ | 134/1000 [00:03<00:18, 45.75it/s]
Running chain, α = 0.11: 13%|█▎ | 134/1000 [00:03<00:18, 45.75it/s]
Running chain, α = 0.11: 13%|█▎ | 134/1000 [00:03<00:18, 45.75it/s]
Running chain, α = 0.11: 14%|█▍ | 139/1000 [00:03<00:18, 45.60it/s]
Running chain, α = 0.11: 14%|█▍ | 139/1000 [00:03<00:18, 45.60it/s]
Running chain, α = 0.11: 14%|█▍ | 139/1000 [00:03<00:18, 45.60it/s]
Running chain, α = 0.11: 14%|█▍ | 139/1000 [00:03<00:18, 45.60it/s]
Running chain, α = 0.11: 14%|█▍ | 139/1000 [00:03<00:18, 45.60it/s]
Running chain, α = 0.11: 14%|█▍ | 139/1000 [00:03<00:18, 45.60it/s]
Running chain, α = 0.11: 14%|█▍ | 144/1000 [00:03<00:18, 45.41it/s]
Running chain, α = 0.12: 14%|█▍ | 144/1000 [00:03<00:18, 45.41it/s]
Running chain, α = 0.12: 14%|█▍ | 144/1000 [00:03<00:18, 45.41it/s]
Running chain, α = 0.12: 14%|█▍ | 144/1000 [00:03<00:18, 45.41it/s]
Running chain, α = 0.12: 14%|█▍ | 144/1000 [00:03<00:18, 45.41it/s]
Running chain, α = 0.11: 14%|█▍ | 144/1000 [00:03<00:18, 45.41it/s]
Running chain, α = 0.11: 15%|█▍ | 149/1000 [00:03<00:18, 45.69it/s]
Running chain, α = 0.11: 15%|█▍ | 149/1000 [00:03<00:18, 45.69it/s]
Running chain, α = 0.12: 15%|█▍ | 149/1000 [00:03<00:18, 45.69it/s]
Running chain, α = 0.12: 15%|█▍ | 149/1000 [00:03<00:18, 45.69it/s]
Running chain, α = 0.12: 15%|█▍ | 149/1000 [00:03<00:18, 45.69it/s]
Running chain, α = 0.12: 15%|█▍ | 149/1000 [00:03<00:18, 45.69it/s]
Running chain, α = 0.12: 15%|█▌ | 154/1000 [00:03<00:18, 45.88it/s]
Running chain, α = 0.12: 15%|█▌ | 154/1000 [00:03<00:18, 45.88it/s]
Running chain, α = 0.12: 15%|█▌ | 154/1000 [00:03<00:18, 45.88it/s]
Running chain, α = 0.12: 15%|█▌ | 154/1000 [00:03<00:18, 45.88it/s]
Running chain, α = 0.12: 15%|█▌ | 154/1000 [00:03<00:18, 45.88it/s]
Running chain, α = 0.12: 15%|█▌ | 154/1000 [00:03<00:18, 45.88it/s]
Running chain, α = 0.12: 16%|█▌ | 159/1000 [00:03<00:18, 46.04it/s]
Running chain, α = 0.12: 16%|█▌ | 159/1000 [00:03<00:18, 46.04it/s]
Running chain, α = 0.12: 16%|█▌ | 159/1000 [00:03<00:18, 46.04it/s]
Running chain, α = 0.12: 16%|█▌ | 159/1000 [00:03<00:18, 46.04it/s]
Running chain, α = 0.13: 16%|█▌ | 159/1000 [00:03<00:18, 46.04it/s]
Running chain, α = 0.13: 16%|█▌ | 159/1000 [00:03<00:18, 46.04it/s]
Running chain, α = 0.13: 16%|█▋ | 164/1000 [00:03<00:18, 45.67it/s]
Running chain, α = 0.13: 16%|█▋ | 164/1000 [00:03<00:18, 45.67it/s]
Running chain, α = 0.13: 16%|█▋ | 164/1000 [00:03<00:18, 45.67it/s]
Running chain, α = 0.13: 16%|█▋ | 164/1000 [00:03<00:18, 45.67it/s]
Running chain, α = 0.13: 16%|█▋ | 164/1000 [00:03<00:18, 45.67it/s]
Running chain, α = 0.14: 16%|█▋ | 164/1000 [00:03<00:18, 45.67it/s]
Running chain, α = 0.14: 17%|█▋ | 169/1000 [00:03<00:18, 45.78it/s]
Running chain, α = 0.14: 17%|█▋ | 169/1000 [00:03<00:18, 45.78it/s]
Running chain, α = 0.14: 17%|█▋ | 169/1000 [00:03<00:18, 45.78it/s]
Running chain, α = 0.15: 17%|█▋ | 169/1000 [00:03<00:18, 45.78it/s]
Running chain, α = 0.15: 17%|█▋ | 169/1000 [00:03<00:18, 45.78it/s]
Running chain, α = 0.16: 17%|█▋ | 169/1000 [00:03<00:18, 45.78it/s]
Running chain, α = 0.16: 17%|█▋ | 174/1000 [00:03<00:17, 45.89it/s]
Running chain, α = 0.17: 17%|█▋ | 174/1000 [00:03<00:17, 45.89it/s]
Running chain, α = 0.17: 17%|█▋ | 174/1000 [00:03<00:17, 45.89it/s]
Running chain, α = 0.17: 17%|█▋ | 174/1000 [00:03<00:17, 45.89it/s]
Running chain, α = 0.17: 17%|█▋ | 174/1000 [00:04<00:17, 45.89it/s]
Running chain, α = 0.17: 17%|█▋ | 174/1000 [00:04<00:17, 45.89it/s]
Running chain, α = 0.17: 18%|█▊ | 179/1000 [00:04<00:17, 46.04it/s]
Running chain, α = 0.17: 18%|█▊ | 179/1000 [00:04<00:17, 46.04it/s]
Running chain, α = 0.17: 18%|█▊ | 179/1000 [00:04<00:17, 46.04it/s]
Running chain, α = 0.17: 18%|█▊ | 179/1000 [00:04<00:17, 46.04it/s]
Running chain, α = 0.17: 18%|█▊ | 179/1000 [00:04<00:17, 46.04it/s]
Running chain, α = 0.17: 18%|█▊ | 179/1000 [00:04<00:17, 46.04it/s]
Running chain, α = 0.17: 18%|█▊ | 184/1000 [00:04<00:17, 46.06it/s]
Running chain, α = 0.17: 18%|█▊ | 184/1000 [00:04<00:17, 46.06it/s]
Running chain, α = 0.17: 18%|█▊ | 184/1000 [00:04<00:17, 46.06it/s]
Running chain, α = 0.17: 18%|█▊ | 184/1000 [00:04<00:17, 46.06it/s]
Running chain, α = 0.17: 18%|█▊ | 184/1000 [00:04<00:17, 46.06it/s]
Running chain, α = 0.17: 18%|█▊ | 184/1000 [00:04<00:17, 46.06it/s]
Running chain, α = 0.17: 19%|█▉ | 189/1000 [00:04<00:17, 46.03it/s]
Running chain, α = 0.18: 19%|█▉ | 189/1000 [00:04<00:17, 46.03it/s]
Running chain, α = 0.19: 19%|█▉ | 189/1000 [00:04<00:17, 46.03it/s]
Running chain, α = 0.19: 19%|█▉ | 189/1000 [00:04<00:17, 46.03it/s]
Running chain, α = 0.19: 19%|█▉ | 189/1000 [00:04<00:17, 46.03it/s]
Running chain, α = 0.20: 19%|█▉ | 189/1000 [00:04<00:17, 46.03it/s]
Running chain, α = 0.20: 19%|█▉ | 194/1000 [00:04<00:17, 45.23it/s]
Running chain, α = 0.20: 19%|█▉ | 194/1000 [00:04<00:17, 45.23it/s]
Running chain, α = 0.20: 19%|█▉ | 194/1000 [00:04<00:17, 45.23it/s]
Running chain, α = 0.21: 19%|█▉ | 194/1000 [00:04<00:17, 45.23it/s]
Running chain, α = 0.21: 19%|█▉ | 194/1000 [00:04<00:17, 45.23it/s]
Running chain, α = 0.22: 19%|█▉ | 194/1000 [00:04<00:17, 45.23it/s]
Running chain, α = 0.22: 20%|█▉ | 199/1000 [00:04<00:17, 45.39it/s]
Running chain, α = 0.23: 20%|█▉ | 199/1000 [00:04<00:17, 45.39it/s]
Running chain, α = 0.23: 20%|█▉ | 199/1000 [00:04<00:17, 45.39it/s]
Running chain, α = 0.22: 20%|█▉ | 199/1000 [00:04<00:17, 45.39it/s]
Running chain, α = 0.23: 20%|█▉ | 199/1000 [00:04<00:17, 45.39it/s]
Running chain, α = 0.23: 20%|█▉ | 199/1000 [00:04<00:17, 45.39it/s]
Running chain, α = 0.23: 20%|██ | 204/1000 [00:04<00:18, 44.15it/s]
Running chain, α = 0.22: 20%|██ | 204/1000 [00:04<00:18, 44.15it/s]
Running chain, α = 0.21: 20%|██ | 204/1000 [00:04<00:18, 44.15it/s]
Running chain, α = 0.22: 20%|██ | 204/1000 [00:04<00:18, 44.15it/s]
Running chain, α = 0.22: 20%|██ | 204/1000 [00:04<00:18, 44.15it/s]
Running chain, α = 0.21: 20%|██ | 204/1000 [00:04<00:18, 44.15it/s]
Running chain, α = 0.21: 21%|██ | 209/1000 [00:04<00:19, 40.97it/s]
Running chain, α = 0.22: 21%|██ | 209/1000 [00:04<00:19, 40.97it/s]
Running chain, α = 0.23: 21%|██ | 209/1000 [00:04<00:19, 40.97it/s]
Running chain, α = 0.24: 21%|██ | 209/1000 [00:04<00:19, 40.97it/s]
Running chain, α = 0.24: 21%|██ | 209/1000 [00:04<00:19, 40.97it/s]
Running chain, α = 0.25: 21%|██ | 209/1000 [00:04<00:19, 40.97it/s]
Running chain, α = 0.25: 21%|██▏ | 214/1000 [00:04<00:23, 33.90it/s]
Running chain, α = 0.24: 21%|██▏ | 214/1000 [00:04<00:23, 33.90it/s]
Running chain, α = 0.24: 21%|██▏ | 214/1000 [00:05<00:23, 33.90it/s]
Running chain, α = 0.24: 21%|██▏ | 214/1000 [00:05<00:23, 33.90it/s]
Running chain, α = 0.24: 21%|██▏ | 214/1000 [00:05<00:23, 33.90it/s]
Running chain, α = 0.24: 22%|██▏ | 218/1000 [00:05<00:23, 32.80it/s]
Running chain, α = 0.24: 22%|██▏ | 218/1000 [00:05<00:23, 32.80it/s]
Running chain, α = 0.23: 22%|██▏ | 218/1000 [00:05<00:23, 32.80it/s]
Running chain, α = 0.23: 22%|██▏ | 218/1000 [00:05<00:23, 32.80it/s]
Running chain, α = 0.23: 22%|██▏ | 218/1000 [00:05<00:23, 32.80it/s]
Running chain, α = 0.23: 22%|██▏ | 218/1000 [00:05<00:23, 32.80it/s]
Running chain, α = 0.23: 22%|██▏ | 223/1000 [00:05<00:21, 36.02it/s]
Running chain, α = 0.23: 22%|██▏ | 223/1000 [00:05<00:21, 36.02it/s]
Running chain, α = 0.23: 22%|██▏ | 223/1000 [00:05<00:21, 36.02it/s]
Running chain, α = 0.22: 22%|██▏ | 223/1000 [00:05<00:21, 36.02it/s]
Running chain, α = 0.22: 22%|██▏ | 223/1000 [00:05<00:21, 36.02it/s]
Running chain, α = 0.22: 22%|██▏ | 223/1000 [00:05<00:21, 36.02it/s]
Running chain, α = 0.22: 23%|██▎ | 228/1000 [00:05<00:19, 38.66it/s]
Running chain, α = 0.22: 23%|██▎ | 228/1000 [00:05<00:19, 38.66it/s]
Running chain, α = 0.22: 23%|██▎ | 228/1000 [00:05<00:19, 38.66it/s]
Running chain, α = 0.21: 23%|██▎ | 228/1000 [00:05<00:19, 38.66it/s]
Running chain, α = 0.21: 23%|██▎ | 228/1000 [00:05<00:19, 38.66it/s]
Running chain, α = 0.20: 23%|██▎ | 228/1000 [00:05<00:19, 38.66it/s]
Running chain, α = 0.20: 23%|██▎ | 233/1000 [00:05<00:18, 40.65it/s]
Running chain, α = 0.20: 23%|██▎ | 233/1000 [00:05<00:18, 40.65it/s]
Running chain, α = 0.20: 23%|██▎ | 233/1000 [00:05<00:18, 40.65it/s]
Running chain, α = 0.20: 23%|██▎ | 233/1000 [00:05<00:18, 40.65it/s]
Running chain, α = 0.20: 23%|██▎ | 233/1000 [00:05<00:18, 40.65it/s]
Running chain, α = 0.20: 23%|██▎ | 233/1000 [00:05<00:18, 40.65it/s]
Running chain, α = 0.20: 24%|██▍ | 238/1000 [00:05<00:18, 42.06it/s]
Running chain, α = 0.20: 24%|██▍ | 238/1000 [00:05<00:18, 42.06it/s]
Running chain, α = 0.20: 24%|██▍ | 238/1000 [00:05<00:18, 42.06it/s]
Running chain, α = 0.20: 24%|██▍ | 238/1000 [00:05<00:18, 42.06it/s]
Running chain, α = 0.20: 24%|██▍ | 238/1000 [00:05<00:18, 42.06it/s]
Running chain, α = 0.20: 24%|██▍ | 238/1000 [00:05<00:18, 42.06it/s]
Running chain, α = 0.20: 24%|██▍ | 243/1000 [00:05<00:18, 41.88it/s]
Running chain, α = 0.20: 24%|██▍ | 243/1000 [00:05<00:18, 41.88it/s]
Running chain, α = 0.19: 24%|██▍ | 243/1000 [00:05<00:18, 41.88it/s]
Running chain, α = 0.19: 24%|██▍ | 243/1000 [00:05<00:18, 41.88it/s]
Running chain, α = 0.19: 24%|██▍ | 243/1000 [00:05<00:18, 41.88it/s]
Running chain, α = 0.19: 24%|██▍ | 243/1000 [00:05<00:18, 41.88it/s]
Running chain, α = 0.19: 25%|██▍ | 248/1000 [00:05<00:17, 43.00it/s]
Running chain, α = 0.19: 25%|██▍ | 248/1000 [00:05<00:17, 43.00it/s]
Running chain, α = 0.19: 25%|██▍ | 248/1000 [00:05<00:17, 43.00it/s]
Running chain, α = 0.18: 25%|██▍ | 248/1000 [00:05<00:17, 43.00it/s]
Running chain, α = 0.18: 25%|██▍ | 248/1000 [00:05<00:17, 43.00it/s]
Running chain, α = 0.18: 25%|██▍ | 248/1000 [00:05<00:17, 43.00it/s]
Running chain, α = 0.18: 25%|██▌ | 253/1000 [00:05<00:16, 43.97it/s]
Running chain, α = 0.18: 25%|██▌ | 253/1000 [00:05<00:16, 43.97it/s]
Running chain, α = 0.18: 25%|██▌ | 253/1000 [00:05<00:16, 43.97it/s]
Running chain, α = 0.18: 25%|██▌ | 253/1000 [00:05<00:16, 43.97it/s]
Running chain, α = 0.18: 25%|██▌ | 253/1000 [00:05<00:16, 43.97it/s]
Running chain, α = 0.18: 25%|██▌ | 253/1000 [00:05<00:16, 43.97it/s]
Running chain, α = 0.18: 26%|██▌ | 258/1000 [00:05<00:16, 44.61it/s]
Running chain, α = 0.18: 26%|██▌ | 258/1000 [00:05<00:16, 44.61it/s]
Running chain, α = 0.18: 26%|██▌ | 258/1000 [00:05<00:16, 44.61it/s]
Running chain, α = 0.18: 26%|██▌ | 258/1000 [00:06<00:16, 44.61it/s]
Running chain, α = 0.18: 26%|██▌ | 258/1000 [00:06<00:16, 44.61it/s]
Running chain, α = 0.17: 26%|██▌ | 258/1000 [00:06<00:16, 44.61it/s]
Running chain, α = 0.17: 26%|██▋ | 263/1000 [00:06<00:16, 45.11it/s]
Running chain, α = 0.17: 26%|██▋ | 263/1000 [00:06<00:16, 45.11it/s]
Running chain, α = 0.17: 26%|██▋ | 263/1000 [00:06<00:16, 45.11it/s]
Running chain, α = 0.17: 26%|██▋ | 263/1000 [00:06<00:16, 45.11it/s]
Running chain, α = 0.17: 26%|██▋ | 263/1000 [00:06<00:16, 45.11it/s]
Running chain, α = 0.17: 26%|██▋ | 263/1000 [00:06<00:16, 45.11it/s]
Running chain, α = 0.17: 27%|██▋ | 268/1000 [00:06<00:16, 45.26it/s]
Running chain, α = 0.16: 27%|██▋ | 268/1000 [00:06<00:16, 45.26it/s]
Running chain, α = 0.16: 27%|██▋ | 268/1000 [00:06<00:16, 45.26it/s]
Running chain, α = 0.16: 27%|██▋ | 268/1000 [00:06<00:16, 45.26it/s]
Running chain, α = 0.15: 27%|██▋ | 268/1000 [00:06<00:16, 45.26it/s]
Running chain, α = 0.15: 27%|██▋ | 268/1000 [00:06<00:16, 45.26it/s]
Running chain, α = 0.15: 27%|██▋ | 273/1000 [00:06<00:16, 45.41it/s]
Running chain, α = 0.14: 27%|██▋ | 273/1000 [00:06<00:16, 45.41it/s]
Running chain, α = 0.13: 27%|██▋ | 273/1000 [00:06<00:16, 45.41it/s]
Running chain, α = 0.13: 27%|██▋ | 273/1000 [00:06<00:16, 45.41it/s]
Running chain, α = 0.13: 27%|██▋ | 273/1000 [00:06<00:16, 45.41it/s]
Running chain, α = 0.13: 27%|██▋ | 273/1000 [00:06<00:16, 45.41it/s]
Running chain, α = 0.13: 28%|██▊ | 278/1000 [00:06<00:15, 45.39it/s]
Running chain, α = 0.13: 28%|██▊ | 278/1000 [00:06<00:15, 45.39it/s]
Running chain, α = 0.13: 28%|██▊ | 278/1000 [00:06<00:15, 45.39it/s]
Running chain, α = 0.13: 28%|██▊ | 278/1000 [00:06<00:15, 45.39it/s]
Running chain, α = 0.13: 28%|██▊ | 278/1000 [00:06<00:15, 45.39it/s]
Running chain, α = 0.13: 28%|██▊ | 278/1000 [00:06<00:15, 45.39it/s]
Running chain, α = 0.13: 28%|██▊ | 283/1000 [00:06<00:15, 45.03it/s]
Running chain, α = 0.13: 28%|██▊ | 283/1000 [00:06<00:15, 45.03it/s]
Running chain, α = 0.13: 28%|██▊ | 283/1000 [00:06<00:15, 45.03it/s]
Running chain, α = 0.13: 28%|██▊ | 283/1000 [00:06<00:15, 45.03it/s]
Running chain, α = 0.13: 28%|██▊ | 283/1000 [00:06<00:15, 45.03it/s]
Running chain, α = 0.13: 28%|██▊ | 283/1000 [00:06<00:15, 45.03it/s]
Running chain, α = 0.13: 29%|██▉ | 288/1000 [00:06<00:15, 45.44it/s]
Running chain, α = 0.13: 29%|██▉ | 288/1000 [00:06<00:15, 45.44it/s]
Running chain, α = 0.12: 29%|██▉ | 288/1000 [00:06<00:15, 45.44it/s]
Running chain, α = 0.11: 29%|██▉ | 288/1000 [00:06<00:15, 45.44it/s]
Running chain, α = 0.11: 29%|██▉ | 288/1000 [00:06<00:15, 45.44it/s]
Running chain, α = 0.11: 29%|██▉ | 288/1000 [00:06<00:15, 45.44it/s]
Running chain, α = 0.11: 29%|██▉ | 293/1000 [00:06<00:15, 45.77it/s]
Running chain, α = 0.10: 29%|██▉ | 293/1000 [00:06<00:15, 45.77it/s]
Running chain, α = 0.10: 29%|██▉ | 293/1000 [00:06<00:15, 45.77it/s]
Running chain, α = 0.10: 29%|██▉ | 293/1000 [00:06<00:15, 45.77it/s]
Running chain, α = 0.09: 29%|██▉ | 293/1000 [00:06<00:15, 45.77it/s]
Running chain, α = 0.09: 29%|██▉ | 293/1000 [00:06<00:15, 45.77it/s]
Running chain, α = 0.09: 30%|██▉ | 298/1000 [00:06<00:15, 45.15it/s]
Running chain, α = 0.08: 30%|██▉ | 298/1000 [00:06<00:15, 45.15it/s]
Running chain, α = 0.07: 30%|██▉ | 298/1000 [00:06<00:15, 45.15it/s]
Running chain, α = 0.07: 30%|██▉ | 298/1000 [00:06<00:15, 45.15it/s]
Running chain, α = 0.07: 30%|██▉ | 298/1000 [00:06<00:15, 45.15it/s]
Running chain, α = 0.06: 30%|██▉ | 298/1000 [00:06<00:15, 45.15it/s]
Running chain, α = 0.06: 30%|███ | 303/1000 [00:06<00:15, 45.36it/s]
Running chain, α = 0.06: 30%|███ | 303/1000 [00:06<00:15, 45.36it/s]
Running chain, α = 0.06: 30%|███ | 303/1000 [00:06<00:15, 45.36it/s]
Running chain, α = 0.06: 30%|███ | 303/1000 [00:07<00:15, 45.36it/s]
Running chain, α = 0.05: 30%|███ | 303/1000 [00:07<00:15, 45.36it/s]
Running chain, α = 0.05: 30%|███ | 303/1000 [00:07<00:15, 45.36it/s]
Running chain, α = 0.05: 31%|███ | 308/1000 [00:07<00:15, 45.60it/s]
Running chain, α = 0.05: 31%|███ | 308/1000 [00:07<00:15, 45.60it/s]
Running chain, α = 0.04: 31%|███ | 308/1000 [00:07<00:15, 45.60it/s]
Running chain, α = 0.03: 31%|███ | 308/1000 [00:07<00:15, 45.60it/s]
Running chain, α = 0.02: 31%|███ | 308/1000 [00:07<00:15, 45.60it/s]
Running chain, α = 0.01: 31%|███ | 308/1000 [00:07<00:15, 45.60it/s]
Running chain, α = 0.01: 31%|███▏ | 313/1000 [00:07<00:14, 45.84it/s]
Running chain, α = 0.00: 31%|███▏ | 313/1000 [00:07<00:14, 45.84it/s]
Running chain, α = 0.00: 31%|███▏ | 313/1000 [00:07<00:14, 45.84it/s]
Running chain, α = 0.00: 31%|███▏ | 313/1000 [00:07<00:14, 45.84it/s]
Running chain, α = 0.00: 31%|███▏ | 313/1000 [00:07<00:14, 45.84it/s]
Running chain, α = 0.00: 31%|███▏ | 313/1000 [00:07<00:14, 45.84it/s]
Running chain, α = 0.00: 32%|███▏ | 318/1000 [00:07<00:14, 45.94it/s]
Running chain, α = 0.00: 32%|███▏ | 318/1000 [00:07<00:14, 45.94it/s]
Running chain, α = 0.00: 32%|███▏ | 318/1000 [00:07<00:14, 45.94it/s]
Running chain, α = 0.00: 32%|███▏ | 318/1000 [00:07<00:14, 45.94it/s]
Running chain, α = 0.00: 32%|███▏ | 318/1000 [00:07<00:14, 45.94it/s]
Running chain, α = 0.00: 32%|███▏ | 318/1000 [00:07<00:14, 45.94it/s]
Running chain, α = 0.00: 32%|███▏ | 323/1000 [00:07<00:14, 45.15it/s]
Running chain, α = 0.00: 32%|███▏ | 323/1000 [00:07<00:14, 45.15it/s]
Running chain, α = 0.00: 32%|███▏ | 323/1000 [00:07<00:14, 45.15it/s]
Running chain, α = 0.00: 32%|███▏ | 323/1000 [00:07<00:14, 45.15it/s]
Running chain, α = 0.00: 32%|███▏ | 323/1000 [00:07<00:14, 45.15it/s]
Running chain, α = 0.00: 32%|███▏ | 323/1000 [00:07<00:14, 45.15it/s]
Running chain, α = 0.00: 33%|███▎ | 328/1000 [00:07<00:14, 45.38it/s]
Running chain, α = 0.00: 33%|███▎ | 328/1000 [00:07<00:14, 45.38it/s]
Running chain, α = 0.00: 33%|███▎ | 328/1000 [00:07<00:14, 45.38it/s]
Running chain, α = 0.00: 33%|███▎ | 328/1000 [00:07<00:14, 45.38it/s]
Running chain, α = 0.00: 33%|███▎ | 328/1000 [00:07<00:14, 45.38it/s]
Running chain, α = 0.00: 33%|███▎ | 328/1000 [00:07<00:14, 45.38it/s]
Running chain, α = 0.00: 33%|███▎ | 333/1000 [00:07<00:14, 45.49it/s]
Running chain, α = 0.00: 33%|███▎ | 333/1000 [00:07<00:14, 45.49it/s]
Running chain, α = 0.00: 33%|███▎ | 333/1000 [00:07<00:14, 45.49it/s]
Running chain, α = 0.00: 33%|███▎ | 333/1000 [00:07<00:14, 45.49it/s]
Running chain, α = 0.00: 33%|███▎ | 333/1000 [00:07<00:14, 45.49it/s]
Running chain, α = 0.00: 33%|███▎ | 333/1000 [00:07<00:14, 45.49it/s]
Running chain, α = 0.00: 34%|███▍ | 338/1000 [00:07<00:14, 45.71it/s]
Running chain, α = 0.00: 34%|███▍ | 338/1000 [00:07<00:14, 45.71it/s]
Running chain, α = 0.00: 34%|███▍ | 338/1000 [00:07<00:14, 45.71it/s]
Running chain, α = 0.00: 34%|███▍ | 338/1000 [00:07<00:14, 45.71it/s]
Running chain, α = 0.00: 34%|███▍ | 338/1000 [00:07<00:14, 45.71it/s]
Running chain, α = 0.00: 34%|███▍ | 338/1000 [00:07<00:14, 45.71it/s]
Running chain, α = 0.00: 34%|███▍ | 343/1000 [00:07<00:14, 45.85it/s]
Running chain, α = 0.00: 34%|███▍ | 343/1000 [00:07<00:14, 45.85it/s]
Running chain, α = 0.00: 34%|███▍ | 343/1000 [00:07<00:14, 45.85it/s]
Running chain, α = 0.00: 34%|███▍ | 343/1000 [00:07<00:14, 45.85it/s]
Running chain, α = 0.00: 34%|███▍ | 343/1000 [00:07<00:14, 45.85it/s]
Running chain, α = 0.00: 34%|███▍ | 343/1000 [00:07<00:14, 45.85it/s]
Running chain, α = 0.00: 35%|███▍ | 348/1000 [00:07<00:14, 45.90it/s]
Running chain, α = 0.00: 35%|███▍ | 348/1000 [00:07<00:14, 45.90it/s]
Running chain, α = 0.00: 35%|███▍ | 348/1000 [00:07<00:14, 45.90it/s]
Running chain, α = 0.00: 35%|███▍ | 348/1000 [00:07<00:14, 45.90it/s]
Running chain, α = 0.00: 35%|███▍ | 348/1000 [00:08<00:14, 45.90it/s]
Running chain, α = 0.00: 35%|███▍ | 348/1000 [00:08<00:14, 45.90it/s]
Running chain, α = 0.00: 35%|███▌ | 353/1000 [00:08<00:14, 46.02it/s]
Running chain, α = 0.00: 35%|███▌ | 353/1000 [00:08<00:14, 46.02it/s]
Running chain, α = 0.00: 35%|███▌ | 353/1000 [00:08<00:14, 46.02it/s]
Running chain, α = 0.00: 35%|███▌ | 353/1000 [00:08<00:14, 46.02it/s]
Running chain, α = 0.00: 35%|███▌ | 353/1000 [00:08<00:14, 46.02it/s]
Running chain, α = 0.00: 35%|███▌ | 353/1000 [00:08<00:14, 46.02it/s]
Running chain, α = 0.00: 36%|███▌ | 358/1000 [00:08<00:13, 46.17it/s]
Running chain, α = 0.00: 36%|███▌ | 358/1000 [00:08<00:13, 46.17it/s]
Running chain, α = 0.00: 36%|███▌ | 358/1000 [00:08<00:13, 46.17it/s]
Running chain, α = 0.00: 36%|███▌ | 358/1000 [00:08<00:13, 46.17it/s]
Running chain, α = 0.00: 36%|███▌ | 358/1000 [00:08<00:13, 46.17it/s]
Running chain, α = 0.00: 36%|███▌ | 358/1000 [00:08<00:13, 46.17it/s]
Running chain, α = 0.00: 36%|███▋ | 363/1000 [00:08<00:13, 46.34it/s]
Running chain, α = 0.00: 36%|███▋ | 363/1000 [00:08<00:13, 46.34it/s]
Running chain, α = 0.00: 36%|███▋ | 363/1000 [00:08<00:13, 46.34it/s]
Running chain, α = 0.00: 36%|███▋ | 363/1000 [00:08<00:13, 46.34it/s]
Running chain, α = 0.00: 36%|███▋ | 363/1000 [00:08<00:13, 46.34it/s]
Running chain, α = 0.00: 36%|███▋ | 363/1000 [00:08<00:13, 46.34it/s]
Running chain, α = 0.00: 37%|███▋ | 368/1000 [00:08<00:13, 46.28it/s]
Running chain, α = 0.00: 37%|███▋ | 368/1000 [00:08<00:13, 46.28it/s]
Running chain, α = 0.00: 37%|███▋ | 368/1000 [00:08<00:13, 46.28it/s]
Running chain, α = 0.00: 37%|███▋ | 368/1000 [00:08<00:13, 46.28it/s]
Running chain, α = 0.00: 37%|███▋ | 368/1000 [00:08<00:13, 46.28it/s]
Running chain, α = 0.00: 37%|███▋ | 368/1000 [00:08<00:13, 46.28it/s]
Running chain, α = 0.00: 37%|███▋ | 373/1000 [00:08<00:13, 45.84it/s]
Running chain, α = 0.00: 37%|███▋ | 373/1000 [00:08<00:13, 45.84it/s]
Running chain, α = 0.00: 37%|███▋ | 373/1000 [00:08<00:13, 45.84it/s]
Running chain, α = 0.00: 37%|███▋ | 373/1000 [00:08<00:13, 45.84it/s]
Running chain, α = 0.00: 37%|███▋ | 373/1000 [00:08<00:13, 45.84it/s]
Running chain, α = 0.00: 37%|███▋ | 373/1000 [00:08<00:13, 45.84it/s]
Running chain, α = 0.00: 38%|███▊ | 378/1000 [00:08<00:13, 46.08it/s]
Running chain, α = 0.00: 38%|███▊ | 378/1000 [00:08<00:13, 46.08it/s]
Running chain, α = 0.00: 38%|███▊ | 378/1000 [00:08<00:13, 46.08it/s]
Running chain, α = 0.00: 38%|███▊ | 378/1000 [00:08<00:13, 46.08it/s]
Running chain, α = 0.00: 38%|███▊ | 378/1000 [00:08<00:13, 46.08it/s]
Running chain, α = 0.00: 38%|███▊ | 378/1000 [00:08<00:13, 46.08it/s]
Running chain, α = 0.00: 38%|███▊ | 383/1000 [00:08<00:13, 46.27it/s]
Running chain, α = 0.00: 38%|███▊ | 383/1000 [00:08<00:13, 46.27it/s]
Running chain, α = 0.00: 38%|███▊ | 383/1000 [00:08<00:13, 46.27it/s]
Running chain, α = 0.00: 38%|███▊ | 383/1000 [00:08<00:13, 46.27it/s]
Running chain, α = 0.00: 38%|███▊ | 383/1000 [00:08<00:13, 46.27it/s]
Running chain, α = 0.00: 38%|███▊ | 383/1000 [00:08<00:13, 46.27it/s]
Running chain, α = 0.00: 39%|███▉ | 388/1000 [00:08<00:13, 46.44it/s]
Running chain, α = 0.00: 39%|███▉ | 388/1000 [00:08<00:13, 46.44it/s]
Running chain, α = 0.00: 39%|███▉ | 388/1000 [00:08<00:13, 46.44it/s]
Running chain, α = 0.00: 39%|███▉ | 388/1000 [00:08<00:13, 46.44it/s]
Running chain, α = 0.00: 39%|███▉ | 388/1000 [00:08<00:13, 46.44it/s]
Running chain, α = 0.00: 39%|███▉ | 388/1000 [00:08<00:13, 46.44it/s]
Running chain, α = 0.00: 39%|███▉ | 393/1000 [00:08<00:13, 46.48it/s]
Running chain, α = 0.00: 39%|███▉ | 393/1000 [00:08<00:13, 46.48it/s]
Running chain, α = 0.00: 39%|███▉ | 393/1000 [00:08<00:13, 46.48it/s]
Running chain, α = 0.00: 39%|███▉ | 393/1000 [00:08<00:13, 46.48it/s]
Running chain, α = 0.00: 39%|███▉ | 393/1000 [00:08<00:13, 46.48it/s]
Running chain, α = 0.00: 39%|███▉ | 393/1000 [00:09<00:13, 46.48it/s]
Running chain, α = 0.00: 40%|███▉ | 398/1000 [00:09<00:12, 46.48it/s]
Running chain, α = 0.00: 40%|███▉ | 398/1000 [00:09<00:12, 46.48it/s]
Running chain, α = 0.00: 40%|███▉ | 398/1000 [00:09<00:12, 46.48it/s]
Running chain, α = 0.00: 40%|███▉ | 398/1000 [00:09<00:12, 46.48it/s]
Running chain, α = 0.00: 40%|███▉ | 398/1000 [00:09<00:12, 46.48it/s]
Running chain, α = 0.00: 40%|███▉ | 398/1000 [00:09<00:12, 46.48it/s]
Running chain, α = 0.00: 40%|████ | 403/1000 [00:09<00:12, 46.50it/s]
Running chain, α = 0.00: 40%|████ | 403/1000 [00:09<00:12, 46.50it/s]
Running chain, α = 0.00: 40%|████ | 403/1000 [00:09<00:12, 46.50it/s]
Running chain, α = 0.00: 40%|████ | 403/1000 [00:09<00:12, 46.50it/s]
Running chain, α = 0.00: 40%|████ | 403/1000 [00:09<00:12, 46.50it/s]
Running chain, α = 0.00: 40%|████ | 403/1000 [00:09<00:12, 46.50it/s]
Running chain, α = 0.00: 41%|████ | 408/1000 [00:09<00:12, 46.54it/s]
Running chain, α = 0.00: 41%|████ | 408/1000 [00:09<00:12, 46.54it/s]
Running chain, α = 0.00: 41%|████ | 408/1000 [00:09<00:12, 46.54it/s]
Running chain, α = 0.00: 41%|████ | 408/1000 [00:09<00:12, 46.54it/s]
Running chain, α = 0.00: 41%|████ | 408/1000 [00:09<00:12, 46.54it/s]
Running chain, α = 0.00: 41%|████ | 408/1000 [00:09<00:12, 46.54it/s]
Running chain, α = 0.00: 41%|████▏ | 413/1000 [00:09<00:12, 46.59it/s]
Running chain, α = 0.00: 41%|████▏ | 413/1000 [00:09<00:12, 46.59it/s]
Running chain, α = 0.00: 41%|████▏ | 413/1000 [00:09<00:12, 46.59it/s]
Running chain, α = 0.00: 41%|████▏ | 413/1000 [00:09<00:12, 46.59it/s]
Running chain, α = 0.00: 41%|████▏ | 413/1000 [00:09<00:12, 46.59it/s]
Running chain, α = 0.00: 41%|████▏ | 413/1000 [00:09<00:12, 46.59it/s]
Running chain, α = 0.00: 42%|████▏ | 418/1000 [00:09<00:12, 45.63it/s]
Running chain, α = 0.00: 42%|████▏ | 418/1000 [00:09<00:12, 45.63it/s]
Running chain, α = 0.00: 42%|████▏ | 418/1000 [00:09<00:12, 45.63it/s]
Running chain, α = 0.00: 42%|████▏ | 418/1000 [00:09<00:12, 45.63it/s]
Running chain, α = 0.00: 42%|████▏ | 418/1000 [00:09<00:12, 45.63it/s]
Running chain, α = 0.00: 42%|████▏ | 418/1000 [00:09<00:12, 45.63it/s]
Running chain, α = 0.00: 42%|████▏ | 423/1000 [00:09<00:12, 45.99it/s]
Running chain, α = 0.00: 42%|████▏ | 423/1000 [00:09<00:12, 45.99it/s]
Running chain, α = 0.00: 42%|████▏ | 423/1000 [00:09<00:12, 45.99it/s]
Running chain, α = 0.00: 42%|████▏ | 423/1000 [00:09<00:12, 45.99it/s]
Running chain, α = 0.00: 42%|████▏ | 423/1000 [00:09<00:12, 45.99it/s]
Running chain, α = 0.00: 42%|████▏ | 423/1000 [00:09<00:12, 45.99it/s]
Running chain, α = 0.00: 43%|████▎ | 428/1000 [00:09<00:12, 46.28it/s]
Running chain, α = 0.00: 43%|████▎ | 428/1000 [00:09<00:12, 46.28it/s]
Running chain, α = 0.00: 43%|████▎ | 428/1000 [00:09<00:12, 46.28it/s]
Running chain, α = 0.00: 43%|████▎ | 428/1000 [00:09<00:12, 46.28it/s]
Running chain, α = 0.00: 43%|████▎ | 428/1000 [00:09<00:12, 46.28it/s]
Running chain, α = 0.00: 43%|████▎ | 428/1000 [00:09<00:12, 46.28it/s]
Running chain, α = 0.00: 43%|████▎ | 433/1000 [00:09<00:12, 46.42it/s]
Running chain, α = 0.00: 43%|████▎ | 433/1000 [00:09<00:12, 46.42it/s]
Running chain, α = 0.00: 43%|████▎ | 433/1000 [00:09<00:12, 46.42it/s]
Running chain, α = 0.00: 43%|████▎ | 433/1000 [00:09<00:12, 46.42it/s]
Running chain, α = 0.00: 43%|████▎ | 433/1000 [00:09<00:12, 46.42it/s]
Running chain, α = 0.00: 43%|████▎ | 433/1000 [00:09<00:12, 46.42it/s]
Running chain, α = 0.00: 44%|████▍ | 438/1000 [00:09<00:12, 46.46it/s]
Running chain, α = 0.00: 44%|████▍ | 438/1000 [00:09<00:12, 46.46it/s]
Running chain, α = 0.00: 44%|████▍ | 438/1000 [00:09<00:12, 46.46it/s]
Running chain, α = 0.00: 44%|████▍ | 438/1000 [00:09<00:12, 46.46it/s]
Running chain, α = 0.00: 44%|████▍ | 438/1000 [00:09<00:12, 46.46it/s]
Running chain, α = 0.00: 44%|████▍ | 438/1000 [00:09<00:12, 46.46it/s]
Running chain, α = 0.00: 44%|████▍ | 443/1000 [00:09<00:11, 46.44it/s]
Running chain, α = 0.00: 44%|████▍ | 443/1000 [00:09<00:11, 46.44it/s]
Running chain, α = 0.00: 44%|████▍ | 443/1000 [00:10<00:11, 46.44it/s]
Running chain, α = 0.00: 44%|████▍ | 443/1000 [00:10<00:11, 46.44it/s]
Running chain, α = 0.00: 44%|████▍ | 443/1000 [00:10<00:11, 46.44it/s]
Running chain, α = 0.00: 44%|████▍ | 443/1000 [00:10<00:11, 46.44it/s]
Running chain, α = 0.00: 45%|████▍ | 448/1000 [00:10<00:12, 45.35it/s]
Running chain, α = 0.00: 45%|████▍ | 448/1000 [00:10<00:12, 45.35it/s]
Running chain, α = 0.00: 45%|████▍ | 448/1000 [00:10<00:12, 45.35it/s]
Running chain, α = 0.00: 45%|████▍ | 448/1000 [00:10<00:12, 45.35it/s]
Running chain, α = 0.00: 45%|████▍ | 448/1000 [00:10<00:12, 45.35it/s]
Running chain, α = 0.00: 45%|████▍ | 448/1000 [00:10<00:12, 45.35it/s]
Running chain, α = 0.00: 45%|████▌ | 453/1000 [00:10<00:12, 45.09it/s]
Running chain, α = 0.00: 45%|████▌ | 453/1000 [00:10<00:12, 45.09it/s]
Running chain, α = 0.00: 45%|████▌ | 453/1000 [00:10<00:12, 45.09it/s]
Running chain, α = 0.00: 45%|████▌ | 453/1000 [00:10<00:12, 45.09it/s]
Running chain, α = 0.00: 45%|████▌ | 453/1000 [00:10<00:12, 45.09it/s]
Running chain, α = 0.00: 45%|████▌ | 453/1000 [00:10<00:12, 45.09it/s]
Running chain, α = 0.00: 46%|████▌ | 458/1000 [00:10<00:13, 41.66it/s]
Running chain, α = 0.00: 46%|████▌ | 458/1000 [00:10<00:13, 41.66it/s]
Running chain, α = 0.00: 46%|████▌ | 458/1000 [00:10<00:13, 41.66it/s]
Running chain, α = 0.00: 46%|████▌ | 458/1000 [00:10<00:13, 41.66it/s]
Running chain, α = 0.00: 46%|████▌ | 458/1000 [00:10<00:13, 41.66it/s]
Running chain, α = 0.00: 46%|████▌ | 458/1000 [00:10<00:13, 41.66it/s]
Running chain, α = 0.00: 46%|████▋ | 463/1000 [00:10<00:12, 41.75it/s]
Running chain, α = 0.00: 46%|████▋ | 463/1000 [00:10<00:12, 41.75it/s]
Running chain, α = 0.00: 46%|████▋ | 463/1000 [00:10<00:12, 41.75it/s]
Running chain, α = 0.00: 46%|████▋ | 463/1000 [00:10<00:12, 41.75it/s]
Running chain, α = 0.00: 46%|████▋ | 463/1000 [00:10<00:12, 41.75it/s]
Running chain, α = 0.00: 46%|████▋ | 463/1000 [00:10<00:12, 41.75it/s]
Running chain, α = 0.00: 47%|████▋ | 468/1000 [00:10<00:12, 41.91it/s]
Running chain, α = 0.00: 47%|████▋ | 468/1000 [00:10<00:12, 41.91it/s]
Running chain, α = 0.00: 47%|████▋ | 468/1000 [00:10<00:12, 41.91it/s]
Running chain, α = 0.00: 47%|████▋ | 468/1000 [00:10<00:12, 41.91it/s]
Running chain, α = 0.00: 47%|████▋ | 468/1000 [00:10<00:12, 41.91it/s]
Running chain, α = 0.00: 47%|████▋ | 468/1000 [00:10<00:12, 41.91it/s]
Running chain, α = 0.00: 47%|████▋ | 473/1000 [00:10<00:12, 41.89it/s]
Running chain, α = 0.00: 47%|████▋ | 473/1000 [00:10<00:12, 41.89it/s]
Running chain, α = 0.00: 47%|████▋ | 473/1000 [00:10<00:12, 41.89it/s]
Running chain, α = 0.00: 47%|████▋ | 473/1000 [00:10<00:12, 41.89it/s]
Running chain, α = 0.00: 47%|████▋ | 473/1000 [00:10<00:12, 41.89it/s]
Running chain, α = 0.00: 47%|████▋ | 473/1000 [00:10<00:12, 41.89it/s]
Running chain, α = 0.00: 48%|████▊ | 478/1000 [00:10<00:12, 42.14it/s]
Running chain, α = 0.00: 48%|████▊ | 478/1000 [00:10<00:12, 42.14it/s]
Running chain, α = 0.00: 48%|████▊ | 478/1000 [00:10<00:12, 42.14it/s]
Running chain, α = 0.00: 48%|████▊ | 478/1000 [00:10<00:12, 42.14it/s]
Running chain, α = 0.00: 48%|████▊ | 478/1000 [00:10<00:12, 42.14it/s]
Running chain, α = 0.00: 48%|████▊ | 478/1000 [00:10<00:12, 42.14it/s]
Running chain, α = 0.00: 48%|████▊ | 483/1000 [00:10<00:12, 42.36it/s]
Running chain, α = 0.00: 48%|████▊ | 483/1000 [00:10<00:12, 42.36it/s]
Running chain, α = 0.00: 48%|████▊ | 483/1000 [00:10<00:12, 42.36it/s]
Running chain, α = 0.00: 48%|████▊ | 483/1000 [00:11<00:12, 42.36it/s]
Running chain, α = 0.00: 48%|████▊ | 483/1000 [00:11<00:12, 42.36it/s]
Running chain, α = 0.00: 48%|████▊ | 483/1000 [00:11<00:12, 42.36it/s]
Running chain, α = 0.00: 49%|████▉ | 488/1000 [00:11<00:11, 42.95it/s]
Running chain, α = 0.00: 49%|████▉ | 488/1000 [00:11<00:11, 42.95it/s]
Running chain, α = 0.00: 49%|████▉ | 488/1000 [00:11<00:11, 42.95it/s]
Running chain, α = 0.00: 49%|████▉ | 488/1000 [00:11<00:11, 42.95it/s]
Running chain, α = 0.00: 49%|████▉ | 488/1000 [00:11<00:11, 42.95it/s]
Running chain, α = 0.00: 49%|████▉ | 488/1000 [00:11<00:11, 42.95it/s]
Running chain, α = 0.00: 49%|████▉ | 493/1000 [00:11<00:11, 43.68it/s]
Running chain, α = 0.00: 49%|████▉ | 493/1000 [00:11<00:11, 43.68it/s]
Running chain, α = 0.00: 49%|████▉ | 493/1000 [00:11<00:11, 43.68it/s]
Running chain, α = 0.00: 49%|████▉ | 493/1000 [00:11<00:11, 43.68it/s]
Running chain, α = 0.00: 49%|████▉ | 493/1000 [00:11<00:11, 43.68it/s]
Running chain, α = 0.00: 49%|████▉ | 493/1000 [00:11<00:11, 43.68it/s]
Running chain, α = 0.00: 50%|████▉ | 498/1000 [00:11<00:11, 44.04it/s]
Running chain, α = 0.00: 50%|████▉ | 498/1000 [00:11<00:11, 44.04it/s]
Running chain, α = 0.00: 50%|████▉ | 498/1000 [00:11<00:11, 44.04it/s]
Running chain, α = 0.00: 50%|████▉ | 498/1000 [00:11<00:11, 44.04it/s]
Running chain, α = 0.00: 50%|████▉ | 498/1000 [00:11<00:11, 44.04it/s]
Running chain, α = 0.00: 50%|████▉ | 498/1000 [00:11<00:11, 44.04it/s]
Running chain, α = 0.00: 50%|█████ | 503/1000 [00:11<00:11, 44.30it/s]
Running chain, α = 0.00: 50%|█████ | 503/1000 [00:11<00:11, 44.30it/s]
Running chain, α = 0.00: 50%|█████ | 503/1000 [00:11<00:11, 44.30it/s]
Running chain, α = 0.00: 50%|█████ | 503/1000 [00:11<00:11, 44.30it/s]
Running chain, α = 0.00: 50%|█████ | 503/1000 [00:11<00:11, 44.30it/s]
Running chain, α = 0.00: 50%|█████ | 503/1000 [00:11<00:11, 44.30it/s]
Running chain, α = 0.00: 51%|█████ | 508/1000 [00:11<00:11, 44.56it/s]
Running chain, α = 0.00: 51%|█████ | 508/1000 [00:11<00:11, 44.56it/s]
Running chain, α = 0.00: 51%|█████ | 508/1000 [00:11<00:11, 44.56it/s]
Running chain, α = 0.00: 51%|█████ | 508/1000 [00:11<00:11, 44.56it/s]
Running chain, α = 0.00: 51%|█████ | 508/1000 [00:11<00:11, 44.56it/s]
Running chain, α = 0.00: 51%|█████ | 508/1000 [00:11<00:11, 44.56it/s]
Running chain, α = 0.00: 51%|█████▏ | 513/1000 [00:11<00:10, 45.11it/s]
Running chain, α = 0.00: 51%|█████▏ | 513/1000 [00:11<00:10, 45.11it/s]
Running chain, α = 0.00: 51%|█████▏ | 513/1000 [00:11<00:10, 45.11it/s]
Running chain, α = 0.00: 51%|█████▏ | 513/1000 [00:11<00:10, 45.11it/s]
Running chain, α = 0.00: 51%|█████▏ | 513/1000 [00:11<00:10, 45.11it/s]
Running chain, α = 0.00: 51%|█████▏ | 513/1000 [00:11<00:10, 45.11it/s]
Running chain, α = 0.00: 52%|█████▏ | 518/1000 [00:11<00:10, 45.52it/s]
Running chain, α = 0.00: 52%|█████▏ | 518/1000 [00:11<00:10, 45.52it/s]
Running chain, α = 0.00: 52%|█████▏ | 518/1000 [00:11<00:10, 45.52it/s]
Running chain, α = 0.00: 52%|█████▏ | 518/1000 [00:11<00:10, 45.52it/s]
Running chain, α = 0.00: 52%|█████▏ | 518/1000 [00:11<00:10, 45.52it/s]
Running chain, α = 0.00: 52%|█████▏ | 518/1000 [00:11<00:10, 45.52it/s]
Running chain, α = 0.00: 52%|█████▏ | 523/1000 [00:11<00:10, 45.82it/s]
Running chain, α = 0.00: 52%|█████▏ | 523/1000 [00:11<00:10, 45.82it/s]
Running chain, α = 0.00: 52%|█████▏ | 523/1000 [00:11<00:10, 45.82it/s]
Running chain, α = 0.00: 52%|█████▏ | 523/1000 [00:11<00:10, 45.82it/s]
Running chain, α = 0.00: 52%|█████▏ | 523/1000 [00:11<00:10, 45.82it/s]
Running chain, α = 0.00: 52%|█████▏ | 523/1000 [00:11<00:10, 45.82it/s]
Running chain, α = 0.00: 53%|█████▎ | 528/1000 [00:11<00:10, 46.03it/s]
Running chain, α = 0.00: 53%|█████▎ | 528/1000 [00:11<00:10, 46.03it/s]
Running chain, α = 0.00: 53%|█████▎ | 528/1000 [00:11<00:10, 46.03it/s]
Running chain, α = 0.00: 53%|█████▎ | 528/1000 [00:11<00:10, 46.03it/s]
Running chain, α = 0.00: 53%|█████▎ | 528/1000 [00:12<00:10, 46.03it/s]
Running chain, α = 0.00: 53%|█████▎ | 528/1000 [00:12<00:10, 46.03it/s]
Running chain, α = 0.00: 53%|█████▎ | 533/1000 [00:12<00:10, 46.24it/s]
Running chain, α = 0.00: 53%|█████▎ | 533/1000 [00:12<00:10, 46.24it/s]
Running chain, α = 0.00: 53%|█████▎ | 533/1000 [00:12<00:10, 46.24it/s]
Running chain, α = 0.00: 53%|█████▎ | 533/1000 [00:12<00:10, 46.24it/s]
Running chain, α = 0.00: 53%|█████▎ | 533/1000 [00:12<00:10, 46.24it/s]
Running chain, α = 0.00: 53%|█████▎ | 533/1000 [00:12<00:10, 46.24it/s]
Running chain, α = 0.00: 54%|█████▍ | 538/1000 [00:12<00:09, 46.26it/s]
Running chain, α = 0.00: 54%|█████▍ | 538/1000 [00:12<00:09, 46.26it/s]
Running chain, α = 0.00: 54%|█████▍ | 538/1000 [00:12<00:09, 46.26it/s]
Running chain, α = 0.00: 54%|█████▍ | 538/1000 [00:12<00:09, 46.26it/s]
Running chain, α = 0.00: 54%|█████▍ | 538/1000 [00:12<00:09, 46.26it/s]
Running chain, α = 0.00: 54%|█████▍ | 538/1000 [00:12<00:09, 46.26it/s]
Running chain, α = 0.00: 54%|█████▍ | 543/1000 [00:12<00:09, 46.20it/s]
Running chain, α = 0.00: 54%|█████▍ | 543/1000 [00:12<00:09, 46.20it/s]
Running chain, α = 0.00: 54%|█████▍ | 543/1000 [00:12<00:09, 46.20it/s]
Running chain, α = 0.00: 54%|█████▍ | 543/1000 [00:12<00:09, 46.20it/s]
Running chain, α = 0.00: 54%|█████▍ | 543/1000 [00:12<00:09, 46.20it/s]
Running chain, α = 0.00: 54%|█████▍ | 543/1000 [00:12<00:09, 46.20it/s]
Running chain, α = 0.00: 55%|█████▍ | 548/1000 [00:12<00:09, 46.24it/s]
Running chain, α = 0.00: 55%|█████▍ | 548/1000 [00:12<00:09, 46.24it/s]
Running chain, α = 0.00: 55%|█████▍ | 548/1000 [00:12<00:09, 46.24it/s]
Running chain, α = 0.00: 55%|█████▍ | 548/1000 [00:12<00:09, 46.24it/s]
Running chain, α = 0.00: 55%|█████▍ | 548/1000 [00:12<00:09, 46.24it/s]
Running chain, α = 0.00: 55%|█████▍ | 548/1000 [00:12<00:09, 46.24it/s]
Running chain, α = 0.00: 55%|█████▌ | 553/1000 [00:12<00:09, 46.12it/s]
Running chain, α = 0.00: 55%|█████▌ | 553/1000 [00:12<00:09, 46.12it/s]
Running chain, α = 0.00: 55%|█████▌ | 553/1000 [00:12<00:09, 46.12it/s]
Running chain, α = 0.00: 55%|█████▌ | 553/1000 [00:12<00:09, 46.12it/s]
Running chain, α = 0.00: 55%|█████▌ | 553/1000 [00:12<00:09, 46.12it/s]
Running chain, α = 0.00: 55%|█████▌ | 553/1000 [00:12<00:09, 46.12it/s]
Running chain, α = 0.00: 56%|█████▌ | 558/1000 [00:12<00:09, 46.22it/s]
Running chain, α = 0.00: 56%|█████▌ | 558/1000 [00:12<00:09, 46.22it/s]
Running chain, α = 0.00: 56%|█████▌ | 558/1000 [00:12<00:09, 46.22it/s]
Running chain, α = 0.00: 56%|█████▌ | 558/1000 [00:12<00:09, 46.22it/s]
Running chain, α = 0.00: 56%|█████▌ | 558/1000 [00:12<00:09, 46.22it/s]
Running chain, α = 0.00: 56%|█████▌ | 558/1000 [00:12<00:09, 46.22it/s]
Running chain, α = 0.00: 56%|█████▋ | 563/1000 [00:12<00:09, 46.28it/s]
Running chain, α = 0.00: 56%|█████▋ | 563/1000 [00:12<00:09, 46.28it/s]
Running chain, α = 0.00: 56%|█████▋ | 563/1000 [00:12<00:09, 46.28it/s]
Running chain, α = 0.00: 56%|█████▋ | 563/1000 [00:12<00:09, 46.28it/s]
Running chain, α = 0.00: 56%|█████▋ | 563/1000 [00:12<00:09, 46.28it/s]
Running chain, α = 0.00: 56%|█████▋ | 563/1000 [00:12<00:09, 46.28it/s]
Running chain, α = 0.00: 57%|█████▋ | 568/1000 [00:12<00:09, 46.40it/s]
Running chain, α = 0.00: 57%|█████▋ | 568/1000 [00:12<00:09, 46.40it/s]
Running chain, α = 0.00: 57%|█████▋ | 568/1000 [00:12<00:09, 46.40it/s]
Running chain, α = 0.00: 57%|█████▋ | 568/1000 [00:12<00:09, 46.40it/s]
Running chain, α = 0.00: 57%|█████▋ | 568/1000 [00:12<00:09, 46.40it/s]
Running chain, α = 0.00: 57%|█████▋ | 568/1000 [00:12<00:09, 46.40it/s]
Running chain, α = 0.00: 57%|█████▋ | 573/1000 [00:12<00:09, 46.40it/s]
Running chain, α = 0.00: 57%|█████▋ | 573/1000 [00:12<00:09, 46.40it/s]
Running chain, α = 0.00: 57%|█████▋ | 573/1000 [00:12<00:09, 46.40it/s]
Running chain, α = 0.00: 57%|█████▋ | 573/1000 [00:12<00:09, 46.40it/s]
Running chain, α = 0.00: 57%|█████▋ | 573/1000 [00:12<00:09, 46.40it/s]
Running chain, α = 0.00: 57%|█████▋ | 573/1000 [00:13<00:09, 46.40it/s]
Running chain, α = 0.00: 58%|█████▊ | 578/1000 [00:13<00:09, 46.36it/s]
Running chain, α = 0.00: 58%|█████▊ | 578/1000 [00:13<00:09, 46.36it/s]
Running chain, α = 0.00: 58%|█████▊ | 578/1000 [00:13<00:09, 46.36it/s]
Running chain, α = 0.00: 58%|█████▊ | 578/1000 [00:13<00:09, 46.36it/s]
Running chain, α = 0.00: 58%|█████▊ | 578/1000 [00:13<00:09, 46.36it/s]
Running chain, α = 0.00: 58%|█████▊ | 578/1000 [00:13<00:09, 46.36it/s]
Running chain, α = 0.00: 58%|█████▊ | 583/1000 [00:13<00:08, 46.39it/s]
Running chain, α = 0.00: 58%|█████▊ | 583/1000 [00:13<00:08, 46.39it/s]
Running chain, α = 0.00: 58%|█████▊ | 583/1000 [00:13<00:08, 46.39it/s]
Running chain, α = 0.00: 58%|█████▊ | 583/1000 [00:13<00:08, 46.39it/s]
Running chain, α = 0.00: 58%|█████▊ | 583/1000 [00:13<00:08, 46.39it/s]
Running chain, α = 0.00: 58%|█████▊ | 583/1000 [00:13<00:08, 46.39it/s]
Running chain, α = 0.00: 59%|█████▉ | 588/1000 [00:13<00:08, 46.39it/s]
Running chain, α = 0.00: 59%|█████▉ | 588/1000 [00:13<00:08, 46.39it/s]
Running chain, α = 0.00: 59%|█████▉ | 588/1000 [00:13<00:08, 46.39it/s]
Running chain, α = 0.00: 59%|█████▉ | 588/1000 [00:13<00:08, 46.39it/s]
Running chain, α = 0.00: 59%|█████▉ | 588/1000 [00:13<00:08, 46.39it/s]
Running chain, α = 0.00: 59%|█████▉ | 588/1000 [00:13<00:08, 46.39it/s]
Running chain, α = 0.00: 59%|█████▉ | 593/1000 [00:13<00:08, 46.35it/s]
Running chain, α = 0.00: 59%|█████▉ | 593/1000 [00:13<00:08, 46.35it/s]
Running chain, α = 0.00: 59%|█████▉ | 593/1000 [00:13<00:08, 46.35it/s]
Running chain, α = 0.00: 59%|█████▉ | 593/1000 [00:13<00:08, 46.35it/s]
Running chain, α = 0.00: 59%|█████▉ | 593/1000 [00:13<00:08, 46.35it/s]
Running chain, α = 0.00: 59%|█████▉ | 593/1000 [00:13<00:08, 46.35it/s]
Running chain, α = 0.00: 60%|█████▉ | 598/1000 [00:13<00:08, 46.34it/s]
Running chain, α = 0.00: 60%|█████▉ | 598/1000 [00:13<00:08, 46.34it/s]
Running chain, α = 0.00: 60%|█████▉ | 598/1000 [00:13<00:08, 46.34it/s]
Running chain, α = 0.00: 60%|█████▉ | 598/1000 [00:13<00:08, 46.34it/s]
Running chain, α = 0.00: 60%|█████▉ | 598/1000 [00:13<00:08, 46.34it/s]
Running chain, α = 0.00: 60%|█████▉ | 598/1000 [00:13<00:08, 46.34it/s]
Running chain, α = 0.00: 60%|██████ | 603/1000 [00:13<00:08, 46.40it/s]
Running chain, α = 0.00: 60%|██████ | 603/1000 [00:13<00:08, 46.40it/s]
Running chain, α = 0.00: 60%|██████ | 603/1000 [00:13<00:08, 46.40it/s]
Running chain, α = 0.00: 60%|██████ | 603/1000 [00:13<00:08, 46.40it/s]
Running chain, α = 0.00: 60%|██████ | 603/1000 [00:13<00:08, 46.40it/s]
Running chain, α = 0.00: 60%|██████ | 603/1000 [00:13<00:08, 46.40it/s]
Running chain, α = 0.00: 61%|██████ | 608/1000 [00:13<00:08, 46.44it/s]
Running chain, α = 0.00: 61%|██████ | 608/1000 [00:13<00:08, 46.44it/s]
Running chain, α = 0.00: 61%|██████ | 608/1000 [00:13<00:08, 46.44it/s]
Running chain, α = 0.00: 61%|██████ | 608/1000 [00:13<00:08, 46.44it/s]
Running chain, α = 0.00: 61%|██████ | 608/1000 [00:13<00:08, 46.44it/s]
Running chain, α = 0.00: 61%|██████ | 608/1000 [00:13<00:08, 46.44it/s]
Running chain, α = 0.00: 61%|██████▏ | 613/1000 [00:13<00:08, 46.39it/s]
Running chain, α = 0.00: 61%|██████▏ | 613/1000 [00:13<00:08, 46.39it/s]
Running chain, α = 0.00: 61%|██████▏ | 613/1000 [00:13<00:08, 46.39it/s]
Running chain, α = 0.00: 61%|██████▏ | 613/1000 [00:13<00:08, 46.39it/s]
Running chain, α = 0.00: 61%|██████▏ | 613/1000 [00:13<00:08, 46.39it/s]
Running chain, α = 0.00: 61%|██████▏ | 613/1000 [00:13<00:08, 46.39it/s]
Running chain, α = 0.00: 62%|██████▏ | 618/1000 [00:13<00:08, 46.50it/s]
Running chain, α = 0.00: 62%|██████▏ | 618/1000 [00:13<00:08, 46.50it/s]
Running chain, α = 0.00: 62%|██████▏ | 618/1000 [00:13<00:08, 46.50it/s]
Running chain, α = 0.00: 62%|██████▏ | 618/1000 [00:13<00:08, 46.50it/s]
Running chain, α = 0.00: 62%|██████▏ | 618/1000 [00:13<00:08, 46.50it/s]
Running chain, α = 0.00: 62%|██████▏ | 618/1000 [00:13<00:08, 46.50it/s]
Running chain, α = 0.00: 62%|██████▏ | 623/1000 [00:13<00:08, 46.45it/s]
Running chain, α = 0.00: 62%|██████▏ | 623/1000 [00:13<00:08, 46.45it/s]
Running chain, α = 0.00: 62%|██████▏ | 623/1000 [00:14<00:08, 46.45it/s]
Running chain, α = 0.00: 62%|██████▏ | 623/1000 [00:14<00:08, 46.45it/s]
Running chain, α = 0.00: 62%|██████▏ | 623/1000 [00:14<00:08, 46.45it/s]
Running chain, α = 0.00: 62%|██████▏ | 623/1000 [00:14<00:08, 46.45it/s]
Running chain, α = 0.00: 63%|██████▎ | 628/1000 [00:14<00:07, 46.54it/s]
Running chain, α = 0.00: 63%|██████▎ | 628/1000 [00:14<00:07, 46.54it/s]
Running chain, α = 0.00: 63%|██████▎ | 628/1000 [00:14<00:07, 46.54it/s]
Running chain, α = 0.00: 63%|██████▎ | 628/1000 [00:14<00:07, 46.54it/s]
Running chain, α = 0.00: 63%|██████▎ | 628/1000 [00:14<00:07, 46.54it/s]
Running chain, α = 0.00: 63%|██████▎ | 628/1000 [00:14<00:07, 46.54it/s]
Running chain, α = 0.00: 63%|██████▎ | 633/1000 [00:14<00:07, 46.45it/s]
Running chain, α = 0.00: 63%|██████▎ | 633/1000 [00:14<00:07, 46.45it/s]
Running chain, α = 0.00: 63%|██████▎ | 633/1000 [00:14<00:07, 46.45it/s]
Running chain, α = 0.00: 63%|██████▎ | 633/1000 [00:14<00:07, 46.45it/s]
Running chain, α = 0.00: 63%|██████▎ | 633/1000 [00:14<00:07, 46.45it/s]
Running chain, α = 0.00: 63%|██████▎ | 633/1000 [00:14<00:07, 46.45it/s]
Running chain, α = 0.00: 64%|██████▍ | 638/1000 [00:14<00:07, 46.43it/s]
Running chain, α = 0.00: 64%|██████▍ | 638/1000 [00:14<00:07, 46.43it/s]
Running chain, α = 0.00: 64%|██████▍ | 638/1000 [00:14<00:07, 46.43it/s]
Running chain, α = 0.00: 64%|██████▍ | 638/1000 [00:14<00:07, 46.43it/s]
Running chain, α = 0.00: 64%|██████▍ | 638/1000 [00:14<00:07, 46.43it/s]
Running chain, α = 0.00: 64%|██████▍ | 638/1000 [00:14<00:07, 46.43it/s]
Running chain, α = 0.00: 64%|██████▍ | 643/1000 [00:14<00:07, 46.30it/s]
Running chain, α = 0.00: 64%|██████▍ | 643/1000 [00:14<00:07, 46.30it/s]
Running chain, α = 0.00: 64%|██████▍ | 643/1000 [00:14<00:07, 46.30it/s]
Running chain, α = 0.00: 64%|██████▍ | 643/1000 [00:14<00:07, 46.30it/s]
Running chain, α = 0.00: 64%|██████▍ | 643/1000 [00:14<00:07, 46.30it/s]
Running chain, α = 0.00: 64%|██████▍ | 643/1000 [00:14<00:07, 46.30it/s]
Running chain, α = 0.00: 65%|██████▍ | 648/1000 [00:14<00:07, 45.97it/s]
Running chain, α = 0.00: 65%|██████▍ | 648/1000 [00:14<00:07, 45.97it/s]
Running chain, α = 0.00: 65%|██████▍ | 648/1000 [00:14<00:07, 45.97it/s]
Running chain, α = 0.00: 65%|██████▍ | 648/1000 [00:14<00:07, 45.97it/s]
Running chain, α = 0.00: 65%|██████▍ | 648/1000 [00:14<00:07, 45.97it/s]
Running chain, α = 0.00: 65%|██████▍ | 648/1000 [00:14<00:07, 45.97it/s]
Running chain, α = 0.00: 65%|██████▌ | 653/1000 [00:14<00:07, 46.14it/s]
Running chain, α = 0.00: 65%|██████▌ | 653/1000 [00:14<00:07, 46.14it/s]
Running chain, α = 0.00: 65%|██████▌ | 653/1000 [00:14<00:07, 46.14it/s]
Running chain, α = 0.00: 65%|██████▌ | 653/1000 [00:14<00:07, 46.14it/s]
Running chain, α = 0.00: 65%|██████▌ | 653/1000 [00:14<00:07, 46.14it/s]
Running chain, α = 0.00: 65%|██████▌ | 653/1000 [00:14<00:07, 46.14it/s]
Running chain, α = 0.00: 66%|██████▌ | 658/1000 [00:14<00:07, 46.23it/s]
Running chain, α = 0.00: 66%|██████▌ | 658/1000 [00:14<00:07, 46.23it/s]
Running chain, α = 0.00: 66%|██████▌ | 658/1000 [00:14<00:07, 46.23it/s]
Running chain, α = 0.00: 66%|██████▌ | 658/1000 [00:14<00:07, 46.23it/s]
Running chain, α = 0.00: 66%|██████▌ | 658/1000 [00:14<00:07, 46.23it/s]
Running chain, α = 0.00: 66%|██████▌ | 658/1000 [00:14<00:07, 46.23it/s]
Running chain, α = 0.00: 66%|██████▋ | 663/1000 [00:14<00:07, 46.38it/s]
Running chain, α = 0.00: 66%|██████▋ | 663/1000 [00:14<00:07, 46.38it/s]
Running chain, α = 0.00: 66%|██████▋ | 663/1000 [00:14<00:07, 46.38it/s]
Running chain, α = 0.00: 66%|██████▋ | 663/1000 [00:14<00:07, 46.38it/s]
Running chain, α = 0.00: 66%|██████▋ | 663/1000 [00:14<00:07, 46.38it/s]
Running chain, α = 0.00: 66%|██████▋ | 663/1000 [00:14<00:07, 46.38it/s]
Running chain, α = 0.00: 67%|██████▋ | 668/1000 [00:14<00:07, 46.45it/s]
Running chain, α = 0.00: 67%|██████▋ | 668/1000 [00:14<00:07, 46.45it/s]
Running chain, α = 0.00: 67%|██████▋ | 668/1000 [00:14<00:07, 46.45it/s]
Running chain, α = 0.00: 67%|██████▋ | 668/1000 [00:15<00:07, 46.45it/s]
Running chain, α = 0.00: 67%|██████▋ | 668/1000 [00:15<00:07, 46.45it/s]
Running chain, α = 0.00: 67%|██████▋ | 668/1000 [00:15<00:07, 46.45it/s]
Running chain, α = 0.00: 67%|██████▋ | 673/1000 [00:15<00:07, 46.40it/s]
Running chain, α = 0.00: 67%|██████▋ | 673/1000 [00:15<00:07, 46.40it/s]
Running chain, α = 0.00: 67%|██████▋ | 673/1000 [00:15<00:07, 46.40it/s]
Running chain, α = 0.00: 67%|██████▋ | 673/1000 [00:15<00:07, 46.40it/s]
Running chain, α = 0.00: 67%|██████▋ | 673/1000 [00:15<00:07, 46.40it/s]
Running chain, α = 0.00: 67%|██████▋ | 673/1000 [00:15<00:07, 46.40it/s]
Running chain, α = 0.00: 68%|██████▊ | 678/1000 [00:15<00:06, 46.32it/s]
Running chain, α = 0.00: 68%|██████▊ | 678/1000 [00:15<00:06, 46.32it/s]
Running chain, α = 0.00: 68%|██████▊ | 678/1000 [00:15<00:06, 46.32it/s]
Running chain, α = 0.00: 68%|██████▊ | 678/1000 [00:15<00:06, 46.32it/s]
Running chain, α = 0.00: 68%|██████▊ | 678/1000 [00:15<00:06, 46.32it/s]
Running chain, α = 0.00: 68%|██████▊ | 678/1000 [00:15<00:06, 46.32it/s]
Running chain, α = 0.00: 68%|██████▊ | 683/1000 [00:15<00:06, 46.33it/s]
Running chain, α = 0.00: 68%|██████▊ | 683/1000 [00:15<00:06, 46.33it/s]
Running chain, α = 0.00: 68%|██████▊ | 683/1000 [00:15<00:06, 46.33it/s]
Running chain, α = 0.00: 68%|██████▊ | 683/1000 [00:15<00:06, 46.33it/s]
Running chain, α = 0.00: 68%|██████▊ | 683/1000 [00:15<00:06, 46.33it/s]
Running chain, α = 0.00: 68%|██████▊ | 683/1000 [00:15<00:06, 46.33it/s]
Running chain, α = 0.00: 69%|██████▉ | 688/1000 [00:15<00:08, 38.17it/s]
Running chain, α = 0.00: 69%|██████▉ | 688/1000 [00:15<00:08, 38.17it/s]
Running chain, α = 0.00: 69%|██████▉ | 688/1000 [00:15<00:08, 38.17it/s]
Running chain, α = 0.00: 69%|██████▉ | 688/1000 [00:15<00:08, 38.17it/s]
Running chain, α = 0.00: 69%|██████▉ | 688/1000 [00:15<00:08, 38.17it/s]
Running chain, α = 0.00: 69%|██████▉ | 688/1000 [00:15<00:08, 38.17it/s]
Running chain, α = 0.00: 69%|██████▉ | 693/1000 [00:15<00:08, 36.82it/s]
Running chain, α = 0.00: 69%|██████▉ | 693/1000 [00:15<00:08, 36.82it/s]
Running chain, α = 0.00: 69%|██████▉ | 693/1000 [00:15<00:08, 36.82it/s]
Running chain, α = 0.00: 69%|██████▉ | 693/1000 [00:15<00:08, 36.82it/s]
Running chain, α = 0.00: 69%|██████▉ | 693/1000 [00:15<00:08, 36.82it/s]
Running chain, α = 0.00: 69%|██████▉ | 693/1000 [00:15<00:08, 36.82it/s]
Running chain, α = 0.00: 70%|██████▉ | 698/1000 [00:15<00:07, 38.35it/s]
Running chain, α = 0.00: 70%|██████▉ | 698/1000 [00:15<00:07, 38.35it/s]
Running chain, α = 0.00: 70%|██████▉ | 698/1000 [00:15<00:07, 38.35it/s]
Running chain, α = 0.00: 70%|██████▉ | 698/1000 [00:15<00:07, 38.35it/s]
Running chain, α = 0.00: 70%|██████▉ | 698/1000 [00:15<00:07, 38.35it/s]
Running chain, α = 0.00: 70%|██████▉ | 698/1000 [00:15<00:07, 38.35it/s]
Running chain, α = 0.00: 70%|███████ | 703/1000 [00:15<00:07, 40.12it/s]
Running chain, α = 0.00: 70%|███████ | 703/1000 [00:15<00:07, 40.12it/s]
Running chain, α = 0.00: 70%|███████ | 703/1000 [00:15<00:07, 40.12it/s]
Running chain, α = 0.00: 70%|███████ | 703/1000 [00:15<00:07, 40.12it/s]
Running chain, α = 0.00: 70%|███████ | 703/1000 [00:15<00:07, 40.12it/s]
Running chain, α = 0.00: 70%|███████ | 703/1000 [00:15<00:07, 40.12it/s]
Running chain, α = 0.00: 71%|███████ | 708/1000 [00:15<00:06, 41.97it/s]
Running chain, α = 0.00: 71%|███████ | 708/1000 [00:15<00:06, 41.97it/s]
Running chain, α = 0.00: 71%|███████ | 708/1000 [00:15<00:06, 41.97it/s]
Running chain, α = 0.00: 71%|███████ | 708/1000 [00:15<00:06, 41.97it/s]
Running chain, α = 0.00: 71%|███████ | 708/1000 [00:16<00:06, 41.97it/s]
Running chain, α = 0.00: 71%|███████ | 708/1000 [00:16<00:06, 41.97it/s]
Running chain, α = 0.00: 71%|███████▏ | 713/1000 [00:16<00:06, 43.33it/s]
Running chain, α = 0.00: 71%|███████▏ | 713/1000 [00:16<00:06, 43.33it/s]
Running chain, α = 0.00: 71%|███████▏ | 713/1000 [00:16<00:06, 43.33it/s]
Running chain, α = 0.00: 71%|███████▏ | 713/1000 [00:16<00:06, 43.33it/s]
Running chain, α = 0.00: 71%|███████▏ | 713/1000 [00:16<00:06, 43.33it/s]
Running chain, α = 0.00: 71%|███████▏ | 713/1000 [00:16<00:06, 43.33it/s]
Running chain, α = 0.00: 72%|███████▏ | 718/1000 [00:16<00:06, 44.34it/s]
Running chain, α = 0.00: 72%|███████▏ | 718/1000 [00:16<00:06, 44.34it/s]
Running chain, α = 0.00: 72%|███████▏ | 718/1000 [00:16<00:06, 44.34it/s]
Running chain, α = 0.00: 72%|███████▏ | 718/1000 [00:16<00:06, 44.34it/s]
Running chain, α = 0.00: 72%|███████▏ | 718/1000 [00:16<00:06, 44.34it/s]
Running chain, α = 0.00: 72%|███████▏ | 718/1000 [00:16<00:06, 44.34it/s]
Running chain, α = 0.00: 72%|███████▏ | 723/1000 [00:16<00:06, 44.97it/s]
Running chain, α = 0.00: 72%|███████▏ | 723/1000 [00:16<00:06, 44.97it/s]
Running chain, α = 0.00: 72%|███████▏ | 723/1000 [00:16<00:06, 44.97it/s]
Running chain, α = 0.00: 72%|███████▏ | 723/1000 [00:16<00:06, 44.97it/s]
Running chain, α = 0.00: 72%|███████▏ | 723/1000 [00:16<00:06, 44.97it/s]
Running chain, α = 0.00: 72%|███████▏ | 723/1000 [00:16<00:06, 44.97it/s]
Running chain, α = 0.00: 73%|███████▎ | 728/1000 [00:16<00:05, 45.35it/s]
Running chain, α = 0.00: 73%|███████▎ | 728/1000 [00:16<00:05, 45.35it/s]
Running chain, α = 0.00: 73%|███████▎ | 728/1000 [00:16<00:05, 45.35it/s]
Running chain, α = 0.00: 73%|███████▎ | 728/1000 [00:16<00:05, 45.35it/s]
Running chain, α = 0.00: 73%|███████▎ | 728/1000 [00:16<00:05, 45.35it/s]
Running chain, α = 0.00: 73%|███████▎ | 728/1000 [00:16<00:05, 45.35it/s]
Running chain, α = 0.00: 73%|███████▎ | 733/1000 [00:16<00:05, 45.44it/s]
Running chain, α = 0.00: 73%|███████▎ | 733/1000 [00:16<00:05, 45.44it/s]
Running chain, α = 0.00: 73%|███████▎ | 733/1000 [00:16<00:05, 45.44it/s]
Running chain, α = 0.00: 73%|███████▎ | 733/1000 [00:16<00:05, 45.44it/s]
Running chain, α = 0.00: 73%|███████▎ | 733/1000 [00:16<00:05, 45.44it/s]
Running chain, α = 0.00: 73%|███████▎ | 733/1000 [00:16<00:05, 45.44it/s]
Running chain, α = 0.00: 74%|███████▍ | 738/1000 [00:16<00:05, 45.67it/s]
Running chain, α = 0.00: 74%|███████▍ | 738/1000 [00:16<00:05, 45.67it/s]
Running chain, α = 0.00: 74%|███████▍ | 738/1000 [00:16<00:05, 45.67it/s]
Running chain, α = 0.00: 74%|███████▍ | 738/1000 [00:16<00:05, 45.67it/s]
Running chain, α = 0.00: 74%|███████▍ | 738/1000 [00:16<00:05, 45.67it/s]
Running chain, α = 0.00: 74%|███████▍ | 738/1000 [00:16<00:05, 45.67it/s]
Running chain, α = 0.00: 74%|███████▍ | 743/1000 [00:16<00:05, 45.98it/s]
Running chain, α = 0.00: 74%|███████▍ | 743/1000 [00:16<00:05, 45.98it/s]
Running chain, α = 0.00: 74%|███████▍ | 743/1000 [00:16<00:05, 45.98it/s]
Running chain, α = 0.00: 74%|███████▍ | 743/1000 [00:16<00:05, 45.98it/s]
Running chain, α = 0.00: 74%|███████▍ | 743/1000 [00:16<00:05, 45.98it/s]
Running chain, α = 0.00: 74%|███████▍ | 743/1000 [00:16<00:05, 45.98it/s]
Running chain, α = 0.00: 75%|███████▍ | 748/1000 [00:16<00:05, 46.17it/s]
Running chain, α = 0.00: 75%|███████▍ | 748/1000 [00:16<00:05, 46.17it/s]
Running chain, α = 0.00: 75%|███████▍ | 748/1000 [00:16<00:05, 46.17it/s]
Running chain, α = 0.00: 75%|███████▍ | 748/1000 [00:16<00:05, 46.17it/s]
Running chain, α = 0.00: 75%|███████▍ | 748/1000 [00:16<00:05, 46.17it/s]
Running chain, α = 0.00: 75%|███████▍ | 748/1000 [00:16<00:05, 46.17it/s]
Running chain, α = 0.00: 75%|███████▌ | 753/1000 [00:16<00:05, 46.25it/s]
Running chain, α = 0.00: 75%|███████▌ | 753/1000 [00:16<00:05, 46.25it/s]
Running chain, α = 0.00: 75%|███████▌ | 753/1000 [00:16<00:05, 46.25it/s]
Running chain, α = 0.00: 75%|███████▌ | 753/1000 [00:16<00:05, 46.25it/s]
Running chain, α = 0.00: 75%|███████▌ | 753/1000 [00:16<00:05, 46.25it/s]
Running chain, α = 0.00: 75%|███████▌ | 753/1000 [00:17<00:05, 46.25it/s]
Running chain, α = 0.00: 76%|███████▌ | 758/1000 [00:17<00:05, 46.39it/s]
Running chain, α = 0.00: 76%|███████▌ | 758/1000 [00:17<00:05, 46.39it/s]
Running chain, α = 0.00: 76%|███████▌ | 758/1000 [00:17<00:05, 46.39it/s]
Running chain, α = 0.00: 76%|███████▌ | 758/1000 [00:17<00:05, 46.39it/s]
Running chain, α = 0.00: 76%|███████▌ | 758/1000 [00:17<00:05, 46.39it/s]
Running chain, α = 0.00: 76%|███████▌ | 758/1000 [00:17<00:05, 46.39it/s]
Running chain, α = 0.00: 76%|███████▋ | 763/1000 [00:17<00:05, 46.50it/s]
Running chain, α = 0.00: 76%|███████▋ | 763/1000 [00:17<00:05, 46.50it/s]
Running chain, α = 0.00: 76%|███████▋ | 763/1000 [00:17<00:05, 46.50it/s]
Running chain, α = 0.00: 76%|███████▋ | 763/1000 [00:17<00:05, 46.50it/s]
Running chain, α = 0.00: 76%|███████▋ | 763/1000 [00:17<00:05, 46.50it/s]
Running chain, α = 0.00: 76%|███████▋ | 763/1000 [00:17<00:05, 46.50it/s]
Running chain, α = 0.00: 77%|███████▋ | 768/1000 [00:17<00:04, 46.48it/s]
Running chain, α = 0.00: 77%|███████▋ | 768/1000 [00:17<00:04, 46.48it/s]
Running chain, α = 0.00: 77%|███████▋ | 768/1000 [00:17<00:04, 46.48it/s]
Running chain, α = 0.00: 77%|███████▋ | 768/1000 [00:17<00:04, 46.48it/s]
Running chain, α = 0.00: 77%|███████▋ | 768/1000 [00:17<00:04, 46.48it/s]
Running chain, α = 0.00: 77%|███████▋ | 768/1000 [00:17<00:04, 46.48it/s]
Running chain, α = 0.00: 77%|███████▋ | 773/1000 [00:17<00:04, 46.49it/s]
Running chain, α = 0.00: 77%|███████▋ | 773/1000 [00:17<00:04, 46.49it/s]
Running chain, α = 0.00: 77%|███████▋ | 773/1000 [00:17<00:04, 46.49it/s]
Running chain, α = 0.00: 77%|███████▋ | 773/1000 [00:17<00:04, 46.49it/s]
Running chain, α = 0.00: 77%|███████▋ | 773/1000 [00:17<00:04, 46.49it/s]
Running chain, α = 0.00: 77%|███████▋ | 773/1000 [00:17<00:04, 46.49it/s]
Running chain, α = 0.00: 78%|███████▊ | 778/1000 [00:17<00:04, 46.53it/s]
Running chain, α = 0.00: 78%|███████▊ | 778/1000 [00:17<00:04, 46.53it/s]
Running chain, α = 0.00: 78%|███████▊ | 778/1000 [00:17<00:04, 46.53it/s]
Running chain, α = 0.00: 78%|███████▊ | 778/1000 [00:17<00:04, 46.53it/s]
Running chain, α = 0.00: 78%|███████▊ | 778/1000 [00:17<00:04, 46.53it/s]
Running chain, α = 0.00: 78%|███████▊ | 778/1000 [00:17<00:04, 46.53it/s]
Running chain, α = 0.00: 78%|███████▊ | 783/1000 [00:17<00:04, 46.46it/s]
Running chain, α = 0.00: 78%|███████▊ | 783/1000 [00:17<00:04, 46.46it/s]
Running chain, α = 0.00: 78%|███████▊ | 783/1000 [00:17<00:04, 46.46it/s]
Running chain, α = 0.00: 78%|███████▊ | 783/1000 [00:17<00:04, 46.46it/s]
Running chain, α = 0.00: 78%|███████▊ | 783/1000 [00:17<00:04, 46.46it/s]
Running chain, α = 0.00: 78%|███████▊ | 783/1000 [00:17<00:04, 46.46it/s]
Running chain, α = 0.00: 79%|███████▉ | 788/1000 [00:17<00:04, 46.38it/s]
Running chain, α = 0.00: 79%|███████▉ | 788/1000 [00:17<00:04, 46.38it/s]
Running chain, α = 0.00: 79%|███████▉ | 788/1000 [00:17<00:04, 46.38it/s]
Running chain, α = 0.00: 79%|███████▉ | 788/1000 [00:17<00:04, 46.38it/s]
Running chain, α = 0.00: 79%|███████▉ | 788/1000 [00:17<00:04, 46.38it/s]
Running chain, α = 0.00: 79%|███████▉ | 788/1000 [00:17<00:04, 46.38it/s]
Running chain, α = 0.00: 79%|███████▉ | 793/1000 [00:17<00:04, 46.31it/s]
Running chain, α = 0.00: 79%|███████▉ | 793/1000 [00:17<00:04, 46.31it/s]
Running chain, α = 0.00: 79%|███████▉ | 793/1000 [00:17<00:04, 46.31it/s]
Running chain, α = 0.00: 79%|███████▉ | 793/1000 [00:17<00:04, 46.31it/s]
Running chain, α = 0.00: 79%|███████▉ | 793/1000 [00:17<00:04, 46.31it/s]
Running chain, α = 0.00: 79%|███████▉ | 793/1000 [00:17<00:04, 46.31it/s]
Running chain, α = 0.00: 80%|███████▉ | 798/1000 [00:17<00:04, 46.41it/s]
Running chain, α = 0.00: 80%|███████▉ | 798/1000 [00:17<00:04, 46.41it/s]
Running chain, α = 0.00: 80%|███████▉ | 798/1000 [00:17<00:04, 46.41it/s]
Running chain, α = 0.00: 80%|███████▉ | 798/1000 [00:17<00:04, 46.41it/s]
Running chain, α = 0.00: 80%|███████▉ | 798/1000 [00:17<00:04, 46.41it/s]
Running chain, α = 0.00: 80%|███████▉ | 798/1000 [00:17<00:04, 46.41it/s]
Running chain, α = 0.00: 80%|████████ | 803/1000 [00:17<00:04, 46.48it/s]
Running chain, α = 0.00: 80%|████████ | 803/1000 [00:18<00:04, 46.48it/s]
Running chain, α = 0.00: 80%|████████ | 803/1000 [00:18<00:04, 46.48it/s]
Running chain, α = 0.00: 80%|████████ | 803/1000 [00:18<00:04, 46.48it/s]
Running chain, α = 0.00: 80%|████████ | 803/1000 [00:18<00:04, 46.48it/s]
Running chain, α = 0.00: 80%|████████ | 803/1000 [00:18<00:04, 46.48it/s]
Running chain, α = 0.00: 81%|████████ | 808/1000 [00:18<00:04, 46.55it/s]
Running chain, α = 0.00: 81%|████████ | 808/1000 [00:18<00:04, 46.55it/s]
Running chain, α = 0.00: 81%|████████ | 808/1000 [00:18<00:04, 46.55it/s]
Running chain, α = 0.00: 81%|████████ | 808/1000 [00:18<00:04, 46.55it/s]
Running chain, α = 0.00: 81%|████████ | 808/1000 [00:18<00:04, 46.55it/s]
Running chain, α = 0.00: 81%|████████ | 808/1000 [00:18<00:04, 46.55it/s]
Running chain, α = 0.00: 81%|████████▏ | 813/1000 [00:18<00:04, 46.55it/s]
Running chain, α = 0.00: 81%|████████▏ | 813/1000 [00:18<00:04, 46.55it/s]
Running chain, α = 0.00: 81%|████████▏ | 813/1000 [00:18<00:04, 46.55it/s]
Running chain, α = 0.00: 81%|████████▏ | 813/1000 [00:18<00:04, 46.55it/s]
Running chain, α = 0.00: 81%|████████▏ | 813/1000 [00:18<00:04, 46.55it/s]
Running chain, α = 0.00: 81%|████████▏ | 813/1000 [00:18<00:04, 46.55it/s]
Running chain, α = 0.00: 82%|████████▏ | 818/1000 [00:18<00:03, 46.53it/s]
Running chain, α = 0.00: 82%|████████▏ | 818/1000 [00:18<00:03, 46.53it/s]
Running chain, α = 0.00: 82%|████████▏ | 818/1000 [00:18<00:03, 46.53it/s]
Running chain, α = 0.00: 82%|████████▏ | 818/1000 [00:18<00:03, 46.53it/s]
Running chain, α = 0.00: 82%|████████▏ | 818/1000 [00:18<00:03, 46.53it/s]
Running chain, α = 0.00: 82%|████████▏ | 818/1000 [00:18<00:03, 46.53it/s]
Running chain, α = 0.00: 82%|████████▏ | 823/1000 [00:18<00:03, 46.18it/s]
Running chain, α = 0.00: 82%|████████▏ | 823/1000 [00:18<00:03, 46.18it/s]
Running chain, α = 0.00: 82%|████████▏ | 823/1000 [00:18<00:03, 46.18it/s]
Running chain, α = 0.00: 82%|████████▏ | 823/1000 [00:18<00:03, 46.18it/s]
Running chain, α = 0.00: 82%|████████▏ | 823/1000 [00:18<00:03, 46.18it/s]
Running chain, α = 0.00: 82%|████████▏ | 823/1000 [00:18<00:03, 46.18it/s]
Running chain, α = 0.00: 83%|████████▎ | 828/1000 [00:18<00:03, 46.27it/s]
Running chain, α = 0.00: 83%|████████▎ | 828/1000 [00:18<00:03, 46.27it/s]
Running chain, α = 0.00: 83%|████████▎ | 828/1000 [00:18<00:03, 46.27it/s]
Running chain, α = 0.00: 83%|████████▎ | 828/1000 [00:18<00:03, 46.27it/s]
Running chain, α = 0.00: 83%|████████▎ | 828/1000 [00:18<00:03, 46.27it/s]
Running chain, α = 0.00: 83%|████████▎ | 828/1000 [00:18<00:03, 46.27it/s]
Running chain, α = 0.00: 83%|████████▎ | 833/1000 [00:18<00:03, 46.29it/s]
Running chain, α = 0.00: 83%|████████▎ | 833/1000 [00:18<00:03, 46.29it/s]
Running chain, α = 0.00: 83%|████████▎ | 833/1000 [00:18<00:03, 46.29it/s]
Running chain, α = 0.00: 83%|████████▎ | 833/1000 [00:18<00:03, 46.29it/s]
Running chain, α = 0.00: 83%|████████▎ | 833/1000 [00:18<00:03, 46.29it/s]
Running chain, α = 0.00: 83%|████████▎ | 833/1000 [00:18<00:03, 46.29it/s]
Running chain, α = 0.00: 84%|████████▍ | 838/1000 [00:18<00:03, 46.27it/s]
Running chain, α = 0.00: 84%|████████▍ | 838/1000 [00:18<00:03, 46.27it/s]
Running chain, α = 0.00: 84%|████████▍ | 838/1000 [00:18<00:03, 46.27it/s]
Running chain, α = 0.00: 84%|████████▍ | 838/1000 [00:18<00:03, 46.27it/s]
Running chain, α = 0.00: 84%|████████▍ | 838/1000 [00:18<00:03, 46.27it/s]
Running chain, α = 0.00: 84%|████████▍ | 838/1000 [00:18<00:03, 46.27it/s]
Running chain, α = 0.00: 84%|████████▍ | 843/1000 [00:18<00:03, 46.43it/s]
Running chain, α = 0.00: 84%|████████▍ | 843/1000 [00:18<00:03, 46.43it/s]
Running chain, α = 0.00: 84%|████████▍ | 843/1000 [00:18<00:03, 46.43it/s]
Running chain, α = 0.00: 84%|████████▍ | 843/1000 [00:18<00:03, 46.43it/s]
Running chain, α = 0.00: 84%|████████▍ | 843/1000 [00:18<00:03, 46.43it/s]
Running chain, α = 0.00: 84%|████████▍ | 843/1000 [00:18<00:03, 46.43it/s]
Running chain, α = 0.00: 85%|████████▍ | 848/1000 [00:18<00:03, 46.41it/s]
Running chain, α = 0.00: 85%|████████▍ | 848/1000 [00:18<00:03, 46.41it/s]
Running chain, α = 0.00: 85%|████████▍ | 848/1000 [00:18<00:03, 46.41it/s]
Running chain, α = 0.00: 85%|████████▍ | 848/1000 [00:19<00:03, 46.41it/s]
Running chain, α = 0.00: 85%|████████▍ | 848/1000 [00:19<00:03, 46.41it/s]
Running chain, α = 0.00: 85%|████████▍ | 848/1000 [00:19<00:03, 46.41it/s]
Running chain, α = 0.00: 85%|████████▌ | 853/1000 [00:19<00:03, 46.50it/s]
Running chain, α = 0.00: 85%|████████▌ | 853/1000 [00:19<00:03, 46.50it/s]
Running chain, α = 0.00: 85%|████████▌ | 853/1000 [00:19<00:03, 46.50it/s]
Running chain, α = 0.00: 85%|████████▌ | 853/1000 [00:19<00:03, 46.50it/s]
Running chain, α = 0.00: 85%|████████▌ | 853/1000 [00:19<00:03, 46.50it/s]
Running chain, α = 0.00: 85%|████████▌ | 853/1000 [00:19<00:03, 46.50it/s]
Running chain, α = 0.00: 86%|████████▌ | 858/1000 [00:19<00:03, 46.51it/s]
Running chain, α = 0.00: 86%|████████▌ | 858/1000 [00:19<00:03, 46.51it/s]
Running chain, α = 0.00: 86%|████████▌ | 858/1000 [00:19<00:03, 46.51it/s]
Running chain, α = 0.00: 86%|████████▌ | 858/1000 [00:19<00:03, 46.51it/s]
Running chain, α = 0.00: 86%|████████▌ | 858/1000 [00:19<00:03, 46.51it/s]
Running chain, α = 0.00: 86%|████████▌ | 858/1000 [00:19<00:03, 46.51it/s]
Running chain, α = 0.00: 86%|████████▋ | 863/1000 [00:19<00:02, 46.05it/s]
Running chain, α = 0.00: 86%|████████▋ | 863/1000 [00:19<00:02, 46.05it/s]
Running chain, α = 0.00: 86%|████████▋ | 863/1000 [00:19<00:02, 46.05it/s]
Running chain, α = 0.00: 86%|████████▋ | 863/1000 [00:19<00:02, 46.05it/s]
Running chain, α = 0.00: 86%|████████▋ | 863/1000 [00:19<00:02, 46.05it/s]
Running chain, α = 0.00: 86%|████████▋ | 863/1000 [00:19<00:02, 46.05it/s]
Running chain, α = 0.00: 87%|████████▋ | 868/1000 [00:19<00:02, 45.99it/s]
Running chain, α = 0.00: 87%|████████▋ | 868/1000 [00:19<00:02, 45.99it/s]
Running chain, α = 0.00: 87%|████████▋ | 868/1000 [00:19<00:02, 45.99it/s]
Running chain, α = 0.00: 87%|████████▋ | 868/1000 [00:19<00:02, 45.99it/s]
Running chain, α = 0.00: 87%|████████▋ | 868/1000 [00:19<00:02, 45.99it/s]
Running chain, α = 0.00: 87%|████████▋ | 868/1000 [00:19<00:02, 45.99it/s]
Running chain, α = 0.00: 87%|████████▋ | 873/1000 [00:19<00:02, 45.75it/s]
Running chain, α = 0.00: 87%|████████▋ | 873/1000 [00:19<00:02, 45.75it/s]
Running chain, α = 0.00: 87%|████████▋ | 873/1000 [00:19<00:02, 45.75it/s]
Running chain, α = 0.00: 87%|████████▋ | 873/1000 [00:19<00:02, 45.75it/s]
Running chain, α = 0.00: 87%|████████▋ | 873/1000 [00:19<00:02, 45.75it/s]
Running chain, α = 0.00: 87%|████████▋ | 873/1000 [00:19<00:02, 45.75it/s]
Running chain, α = 0.00: 88%|████████▊ | 878/1000 [00:19<00:02, 46.00it/s]
Running chain, α = 0.00: 88%|████████▊ | 878/1000 [00:19<00:02, 46.00it/s]
Running chain, α = 0.00: 88%|████████▊ | 878/1000 [00:19<00:02, 46.00it/s]
Running chain, α = 0.00: 88%|████████▊ | 878/1000 [00:19<00:02, 46.00it/s]
Running chain, α = 0.00: 88%|████████▊ | 878/1000 [00:19<00:02, 46.00it/s]
Running chain, α = 0.00: 88%|████████▊ | 878/1000 [00:19<00:02, 46.00it/s]
Running chain, α = 0.00: 88%|████████▊ | 883/1000 [00:19<00:02, 46.24it/s]
Running chain, α = 0.00: 88%|████████▊ | 883/1000 [00:19<00:02, 46.24it/s]
Running chain, α = 0.00: 88%|████████▊ | 883/1000 [00:19<00:02, 46.24it/s]
Running chain, α = 0.00: 88%|████████▊ | 883/1000 [00:19<00:02, 46.24it/s]
Running chain, α = 0.00: 88%|████████▊ | 883/1000 [00:19<00:02, 46.24it/s]
Running chain, α = 0.00: 88%|████████▊ | 883/1000 [00:19<00:02, 46.24it/s]
Running chain, α = 0.00: 89%|████████▉ | 888/1000 [00:19<00:02, 46.33it/s]
Running chain, α = 0.00: 89%|████████▉ | 888/1000 [00:19<00:02, 46.33it/s]
Running chain, α = 0.00: 89%|████████▉ | 888/1000 [00:19<00:02, 46.33it/s]
Running chain, α = 0.00: 89%|████████▉ | 888/1000 [00:19<00:02, 46.33it/s]
Running chain, α = 0.00: 89%|████████▉ | 888/1000 [00:19<00:02, 46.33it/s]
Running chain, α = 0.00: 89%|████████▉ | 888/1000 [00:19<00:02, 46.33it/s]
Running chain, α = 0.00: 89%|████████▉ | 893/1000 [00:19<00:02, 46.23it/s]
Running chain, α = 0.00: 89%|████████▉ | 893/1000 [00:19<00:02, 46.23it/s]
Running chain, α = 0.00: 89%|████████▉ | 893/1000 [00:19<00:02, 46.23it/s]
Running chain, α = 0.00: 89%|████████▉ | 893/1000 [00:19<00:02, 46.23it/s]
Running chain, α = 0.00: 89%|████████▉ | 893/1000 [00:20<00:02, 46.23it/s]
Running chain, α = 0.00: 89%|████████▉ | 893/1000 [00:20<00:02, 46.23it/s]
Running chain, α = 0.00: 90%|████████▉ | 898/1000 [00:20<00:02, 46.32it/s]
Running chain, α = 0.00: 90%|████████▉ | 898/1000 [00:20<00:02, 46.32it/s]
Running chain, α = 0.00: 90%|████████▉ | 898/1000 [00:20<00:02, 46.32it/s]
Running chain, α = 0.00: 90%|████████▉ | 898/1000 [00:20<00:02, 46.32it/s]
Running chain, α = 0.00: 90%|████████▉ | 898/1000 [00:20<00:02, 46.32it/s]
Running chain, α = 0.00: 90%|████████▉ | 898/1000 [00:20<00:02, 46.32it/s]
Running chain, α = 0.00: 90%|█████████ | 903/1000 [00:20<00:02, 46.21it/s]
Running chain, α = 0.00: 90%|█████████ | 903/1000 [00:20<00:02, 46.21it/s]
Running chain, α = 0.00: 90%|█████████ | 903/1000 [00:20<00:02, 46.21it/s]
Running chain, α = 0.00: 90%|█████████ | 903/1000 [00:20<00:02, 46.21it/s]
Running chain, α = 0.00: 90%|█████████ | 903/1000 [00:20<00:02, 46.21it/s]
Running chain, α = 0.00: 90%|█████████ | 903/1000 [00:20<00:02, 46.21it/s]
Running chain, α = 0.00: 91%|█████████ | 908/1000 [00:20<00:01, 46.29it/s]
Running chain, α = 0.00: 91%|█████████ | 908/1000 [00:20<00:01, 46.29it/s]
Running chain, α = 0.00: 91%|█████████ | 908/1000 [00:20<00:01, 46.29it/s]
Running chain, α = 0.00: 91%|█████████ | 908/1000 [00:20<00:01, 46.29it/s]
Running chain, α = 0.00: 91%|█████████ | 908/1000 [00:20<00:01, 46.29it/s]
Running chain, α = 0.00: 91%|█████████ | 908/1000 [00:20<00:01, 46.29it/s]
Running chain, α = 0.00: 91%|█████████▏| 913/1000 [00:20<00:01, 46.33it/s]
Running chain, α = 0.00: 91%|█████████▏| 913/1000 [00:20<00:01, 46.33it/s]
Running chain, α = 0.00: 91%|█████████▏| 913/1000 [00:20<00:01, 46.33it/s]
Running chain, α = 0.00: 91%|█████████▏| 913/1000 [00:20<00:01, 46.33it/s]
Running chain, α = 0.00: 91%|█████████▏| 913/1000 [00:20<00:01, 46.33it/s]
Running chain, α = 0.00: 91%|█████████▏| 913/1000 [00:20<00:01, 46.33it/s]
Running chain, α = 0.00: 92%|█████████▏| 918/1000 [00:20<00:01, 45.81it/s]
Running chain, α = 0.00: 92%|█████████▏| 918/1000 [00:20<00:01, 45.81it/s]
Running chain, α = 0.00: 92%|█████████▏| 918/1000 [00:20<00:01, 45.81it/s]
Running chain, α = 0.00: 92%|█████████▏| 918/1000 [00:20<00:01, 45.81it/s]
Running chain, α = 0.00: 92%|█████████▏| 918/1000 [00:20<00:01, 45.81it/s]
Running chain, α = 0.00: 92%|█████████▏| 918/1000 [00:20<00:01, 45.81it/s]
Running chain, α = 0.00: 92%|█████████▏| 923/1000 [00:20<00:01, 45.94it/s]
Running chain, α = 0.00: 92%|█████████▏| 923/1000 [00:20<00:01, 45.94it/s]
Running chain, α = 0.00: 92%|█████████▏| 923/1000 [00:20<00:01, 45.94it/s]
Running chain, α = 0.00: 92%|█████████▏| 923/1000 [00:20<00:01, 45.94it/s]
Running chain, α = 0.00: 92%|█████████▏| 923/1000 [00:20<00:01, 45.94it/s]
Running chain, α = 0.00: 92%|█████████▏| 923/1000 [00:20<00:01, 45.94it/s]
Running chain, α = 0.00: 93%|█████████▎| 928/1000 [00:20<00:01, 46.14it/s]
Running chain, α = 0.00: 93%|█████████▎| 928/1000 [00:20<00:01, 46.14it/s]
Running chain, α = 0.00: 93%|█████████▎| 928/1000 [00:20<00:01, 46.14it/s]
Running chain, α = 0.00: 93%|█████████▎| 928/1000 [00:20<00:01, 46.14it/s]
Running chain, α = 0.00: 93%|█████████▎| 928/1000 [00:20<00:01, 46.14it/s]
Running chain, α = 0.00: 93%|█████████▎| 928/1000 [00:20<00:01, 46.14it/s]
Running chain, α = 0.00: 93%|█████████▎| 933/1000 [00:20<00:01, 46.32it/s]
Running chain, α = 0.00: 93%|█████████▎| 933/1000 [00:20<00:01, 46.32it/s]
Running chain, α = 0.00: 93%|█████████▎| 933/1000 [00:20<00:01, 46.32it/s]
Running chain, α = 0.00: 93%|█████████▎| 933/1000 [00:20<00:01, 46.32it/s]
Running chain, α = 0.00: 93%|█████████▎| 933/1000 [00:20<00:01, 46.32it/s]
Running chain, α = 0.00: 93%|█████████▎| 933/1000 [00:20<00:01, 46.32it/s]
Running chain, α = 0.00: 94%|█████████▍| 938/1000 [00:20<00:01, 46.34it/s]
Running chain, α = 0.00: 94%|█████████▍| 938/1000 [00:20<00:01, 46.34it/s]
Running chain, α = 0.00: 94%|█████████▍| 938/1000 [00:20<00:01, 46.34it/s]
Running chain, α = 0.00: 94%|█████████▍| 938/1000 [00:20<00:01, 46.34it/s]
Running chain, α = 0.00: 94%|█████████▍| 938/1000 [00:20<00:01, 46.34it/s]
Running chain, α = 0.00: 94%|█████████▍| 938/1000 [00:21<00:01, 46.34it/s]
Running chain, α = 0.00: 94%|█████████▍| 943/1000 [00:21<00:01, 45.54it/s]
Running chain, α = 0.00: 94%|█████████▍| 943/1000 [00:21<00:01, 45.54it/s]
Running chain, α = 0.00: 94%|█████████▍| 943/1000 [00:21<00:01, 45.54it/s]
Running chain, α = 0.00: 94%|█████████▍| 943/1000 [00:21<00:01, 45.54it/s]
Running chain, α = 0.00: 94%|█████████▍| 943/1000 [00:21<00:01, 45.54it/s]
Running chain, α = 0.00: 94%|█████████▍| 943/1000 [00:21<00:01, 45.54it/s]
Running chain, α = 0.00: 95%|█████████▍| 948/1000 [00:21<00:01, 42.27it/s]
Running chain, α = 0.00: 95%|█████████▍| 948/1000 [00:21<00:01, 42.27it/s]
Running chain, α = 0.00: 95%|█████████▍| 948/1000 [00:21<00:01, 42.27it/s]
Running chain, α = 0.00: 95%|█████████▍| 948/1000 [00:21<00:01, 42.27it/s]
Running chain, α = 0.00: 95%|█████████▍| 948/1000 [00:21<00:01, 42.27it/s]
Running chain, α = 0.00: 95%|█████████▍| 948/1000 [00:21<00:01, 42.27it/s]
Running chain, α = 0.00: 95%|█████████▌| 953/1000 [00:21<00:01, 40.29it/s]
Running chain, α = 0.00: 95%|█████████▌| 953/1000 [00:21<00:01, 40.29it/s]
Running chain, α = 0.00: 95%|█████████▌| 953/1000 [00:21<00:01, 40.29it/s]
Running chain, α = 0.00: 95%|█████████▌| 953/1000 [00:21<00:01, 40.29it/s]
Running chain, α = 0.00: 95%|█████████▌| 953/1000 [00:21<00:01, 40.29it/s]
Running chain, α = 0.00: 95%|█████████▌| 953/1000 [00:21<00:01, 40.29it/s]
Running chain, α = 0.00: 96%|█████████▌| 958/1000 [00:21<00:01, 40.23it/s]
Running chain, α = 0.00: 96%|█████████▌| 958/1000 [00:21<00:01, 40.23it/s]
Running chain, α = 0.00: 96%|█████████▌| 958/1000 [00:21<00:01, 40.23it/s]
Running chain, α = 0.00: 96%|█████████▌| 958/1000 [00:21<00:01, 40.23it/s]
Running chain, α = 0.00: 96%|█████████▌| 958/1000 [00:21<00:01, 40.23it/s]
Running chain, α = 0.00: 96%|█████████▌| 958/1000 [00:21<00:01, 40.23it/s]
Running chain, α = 0.00: 96%|█████████▋| 963/1000 [00:21<00:00, 42.08it/s]
Running chain, α = 0.00: 96%|█████████▋| 963/1000 [00:21<00:00, 42.08it/s]
Running chain, α = 0.00: 96%|█████████▋| 963/1000 [00:21<00:00, 42.08it/s]
Running chain, α = 0.00: 96%|█████████▋| 963/1000 [00:21<00:00, 42.08it/s]
Running chain, α = 0.00: 96%|█████████▋| 963/1000 [00:21<00:00, 42.08it/s]
Running chain, α = 0.00: 96%|█████████▋| 963/1000 [00:21<00:00, 42.08it/s]
Running chain, α = 0.00: 97%|█████████▋| 968/1000 [00:21<00:00, 43.42it/s]
Running chain, α = 0.00: 97%|█████████▋| 968/1000 [00:21<00:00, 43.42it/s]
Running chain, α = 0.00: 97%|█████████▋| 968/1000 [00:21<00:00, 43.42it/s]
Running chain, α = 0.00: 97%|█████████▋| 968/1000 [00:21<00:00, 43.42it/s]
Running chain, α = 0.00: 97%|█████████▋| 968/1000 [00:21<00:00, 43.42it/s]
Running chain, α = 0.00: 97%|█████████▋| 968/1000 [00:21<00:00, 43.42it/s]
Running chain, α = 0.00: 97%|█████████▋| 973/1000 [00:21<00:00, 44.42it/s]
Running chain, α = 0.00: 97%|█████████▋| 973/1000 [00:21<00:00, 44.42it/s]
Running chain, α = 0.00: 97%|█████████▋| 973/1000 [00:21<00:00, 44.42it/s]
Running chain, α = 0.00: 97%|█████████▋| 973/1000 [00:21<00:00, 44.42it/s]
Running chain, α = 0.00: 97%|█████████▋| 973/1000 [00:21<00:00, 44.42it/s]
Running chain, α = 0.00: 97%|█████████▋| 973/1000 [00:21<00:00, 44.42it/s]
Running chain, α = 0.00: 98%|█████████▊| 978/1000 [00:21<00:00, 45.05it/s]
Running chain, α = 0.00: 98%|█████████▊| 978/1000 [00:21<00:00, 45.05it/s]
Running chain, α = 0.00: 98%|█████████▊| 978/1000 [00:21<00:00, 45.05it/s]
Running chain, α = 0.00: 98%|█████████▊| 978/1000 [00:21<00:00, 45.05it/s]
Running chain, α = 0.00: 98%|█████████▊| 978/1000 [00:21<00:00, 45.05it/s]
Running chain, α = 0.00: 98%|█████████▊| 978/1000 [00:21<00:00, 45.05it/s]
Running chain, α = 0.00: 98%|█████████▊| 983/1000 [00:21<00:00, 45.62it/s]
Running chain, α = 0.00: 98%|█████████▊| 983/1000 [00:21<00:00, 45.62it/s]
Running chain, α = 0.00: 98%|█████████▊| 983/1000 [00:21<00:00, 45.62it/s]
Running chain, α = 0.00: 98%|█████████▊| 983/1000 [00:22<00:00, 45.62it/s]
Running chain, α = 0.00: 98%|█████████▊| 983/1000 [00:22<00:00, 45.62it/s]
Running chain, α = 0.00: 98%|█████████▊| 983/1000 [00:22<00:00, 45.62it/s]
Running chain, α = 0.00: 99%|█████████▉| 988/1000 [00:22<00:00, 45.99it/s]
Running chain, α = 0.00: 99%|█████████▉| 988/1000 [00:22<00:00, 45.99it/s]
Running chain, α = 0.00: 99%|█████████▉| 988/1000 [00:22<00:00, 45.99it/s]
Running chain, α = 0.00: 99%|█████████▉| 988/1000 [00:22<00:00, 45.99it/s]
Running chain, α = 0.00: 99%|█████████▉| 988/1000 [00:22<00:00, 45.99it/s]
Running chain, α = 0.00: 99%|█████████▉| 988/1000 [00:22<00:00, 45.99it/s]
Running chain, α = 0.00: 99%|█████████▉| 993/1000 [00:22<00:00, 46.25it/s]
Running chain, α = 0.00: 99%|█████████▉| 993/1000 [00:22<00:00, 46.25it/s]
Running chain, α = 0.00: 99%|█████████▉| 993/1000 [00:22<00:00, 46.25it/s]
Running chain, α = 0.00: 99%|█████████▉| 993/1000 [00:22<00:00, 46.25it/s]
Running chain, α = 0.00: 99%|█████████▉| 993/1000 [00:22<00:00, 46.25it/s]
Running chain, α = 0.00: 99%|█████████▉| 993/1000 [00:22<00:00, 46.25it/s]
Running chain, α = 0.00: 100%|█████████▉| 998/1000 [00:22<00:00, 46.35it/s]
Running chain, α = 0.00: 100%|█████████▉| 998/1000 [00:22<00:00, 46.35it/s]
Running chain, α = 0.00: 100%|█████████▉| 998/1000 [00:22<00:00, 46.35it/s]
Running chain, α = 0.00: 100%|██████████| 1000/1000 [00:22<00:00, 44.78it/s]
Out[10]:
{'n_steps': 1001,
'n_dim': 3,
'move_fraction': 0.036,
'hf_call_fraction': 0.001998001998001998,
'n_hf_calls': 2,
'posterior_rmse': 6.881929007119004}
In [11]:
Copied!
samples_single = chain_single.samples
used_hf_single = chain_single.extras.used_hf
fig1, ax1 = plot_chain_2d(
samples_single[:, :2],
used_hf=used_hf_single,
theta_true=theta_true[:2],
title="Single posterior: samples (first two parameters)",
show=True,
)
fig2, ax2 = plot_cumulative_hf_fraction(
used_hf_single,
title="Single posterior: cumulative HF fraction",
show=True,
)
fig_s1, ax_s1 = plot_prediction_at_theta(
model_single.lf_model,
theta=theta_true,
y_obs=y_obs,
t=t,
y_true=hf_forward(theta_true),
title="Surrogate prediction (after MCMC sampling)",
)
samples_single = chain_single.samples
used_hf_single = chain_single.extras.used_hf
fig1, ax1 = plot_chain_2d(
samples_single[:, :2],
used_hf=used_hf_single,
theta_true=theta_true[:2],
title="Single posterior: samples (first two parameters)",
show=True,
)
fig2, ax2 = plot_cumulative_hf_fraction(
used_hf_single,
title="Single posterior: cumulative HF fraction",
show=True,
)
fig_s1, ax_s1 = plot_prediction_at_theta(
model_single.lf_model,
theta=theta_true,
y_obs=y_obs,
t=t,
y_true=hf_forward(theta_true),
title="Surrogate prediction (after MCMC sampling)",
)
Part 2 — DA-MCMC guided active learning with an adaptive subchain (recommended)¶
Now we pass two posteriors [coarse, fine] and sample with sample_adaptive_active_chain.
The adaptive policy updates subchain_length online using LF-HF discrepancy computed at HF calls.
Because tinyDA uses a fixed subsampling_rate per call, we run sampling in chunks and update the
subsampling rate between chunks.
In [12]:
Copied!
theta0 = prior_mean.copy()
proposal = AdaptiveMetropolisShared(
C0=0.1 * prior_cov,
period=100,
share_across_deepcopy=True,
adaptive=True,
sd=1,
)
result_adapt = sample_adaptive_active_chain(
model=model_adapt,
posterior=posterior_adapt,
proposal=proposal,
n_coarse_evals=n_coarse_evals,
initial_parameters=theta0,
chain_key="chain_coarse_0",
config=ChunkedMCMCConfig(chain_key="chain_coarse_0", chunk_size=chunk_size),
store_coarse_chain=True,
)
chain_adapt = result_adapt.chain
chain_adapt.summary(theta_true=theta_true, burn_in=burn_in)
theta0 = prior_mean.copy()
proposal = AdaptiveMetropolisShared(
C0=0.1 * prior_cov,
period=100,
share_across_deepcopy=True,
adaptive=True,
sd=1,
)
result_adapt = sample_adaptive_active_chain(
model=model_adapt,
posterior=posterior_adapt,
proposal=proposal,
n_coarse_evals=n_coarse_evals,
initial_parameters=theta0,
chain_key="chain_coarse_0",
config=ChunkedMCMCConfig(chain_key="chain_coarse_0", chunk_size=chunk_size),
store_coarse_chain=True,
)
chain_adapt = result_adapt.chain
chain_adapt.summary(theta_true=theta_true, burn_in=burn_in)
/opt/hostedtoolcache/Python/3.10.19/x64/lib/python3.10/site-packages/tinyDA/sampler.py:114: UserWarning: subsampling_rate has been deprecated in favour of subchain_length.
Sampling chain 1/1
0%| | 0/12 [00:00<?, ?it/s]
Running chain, α_c = 1.000, α_f = 1.00: 0%| | 0/12 [00:00<?, ?it/s]
Running chain, α_c = 1.000, α_f = 1.00: 8%|▊ | 1/12 [00:00<00:06, 1.59it/s]
Running chain, α_c = 0.273, α_f = 1.00: 8%|▊ | 1/12 [00:00<00:06, 1.59it/s]
Running chain, α_c = 0.273, α_f = 1.00: 17%|█▋ | 2/12 [00:01<00:05, 1.94it/s]
Running chain, α_c = 0.140, α_f = 0.67: 17%|█▋ | 2/12 [00:01<00:05, 1.94it/s]
Running chain, α_c = 0.140, α_f = 0.67: 25%|██▌ | 3/12 [00:01<00:04, 2.16it/s]
Running chain, α_c = 0.094, α_f = 0.50: 25%|██▌ | 3/12 [00:01<00:04, 2.16it/s]
Running chain, α_c = 0.094, α_f = 0.50: 33%|███▎ | 4/12 [00:01<00:03, 2.28it/s]
Running chain, α_c = 0.071, α_f = 0.40: 33%|███▎ | 4/12 [00:01<00:03, 2.28it/s]
Running chain, α_c = 0.071, α_f = 0.40: 42%|████▏ | 5/12 [00:02<00:02, 2.35it/s]
Running chain, α_c = 0.057, α_f = 0.33: 42%|████▏ | 5/12 [00:02<00:02, 2.35it/s]
Running chain, α_c = 0.057, α_f = 0.33: 50%|█████ | 6/12 [00:02<00:02, 2.39it/s]
Running chain, α_c = 0.047, α_f = 0.29: 50%|█████ | 6/12 [00:02<00:02, 2.39it/s]
Running chain, α_c = 0.047, α_f = 0.29: 58%|█████▊ | 7/12 [00:03<00:02, 2.43it/s]
Running chain, α_c = 0.041, α_f = 0.25: 58%|█████▊ | 7/12 [00:03<00:02, 2.43it/s]
Running chain, α_c = 0.041, α_f = 0.25: 67%|██████▋ | 8/12 [00:03<00:01, 2.44it/s]
Running chain, α_c = 0.036, α_f = 0.22: 67%|██████▋ | 8/12 [00:03<00:01, 2.44it/s]
Running chain, α_c = 0.036, α_f = 0.22: 75%|███████▌ | 9/12 [00:03<00:01, 2.45it/s]
Running chain, α_c = 0.032, α_f = 0.20: 75%|███████▌ | 9/12 [00:03<00:01, 2.45it/s]
Running chain, α_c = 0.032, α_f = 0.20: 83%|████████▎ | 10/12 [00:04<00:00, 2.46it/s]
Running chain, α_c = 0.028, α_f = 0.18: 83%|████████▎ | 10/12 [00:04<00:00, 2.46it/s]
Running chain, α_c = 0.028, α_f = 0.18: 92%|█████████▏| 11/12 [00:04<00:00, 2.47it/s]
Running chain, α_c = 0.026, α_f = 0.17: 92%|█████████▏| 11/12 [00:04<00:00, 2.47it/s]
Running chain, α_c = 0.026, α_f = 0.17: 100%|██████████| 12/12 [00:05<00:00, 2.47it/s]
Running chain, α_c = 0.026, α_f = 0.17: 100%|██████████| 12/12 [00:05<00:00, 2.36it/s]
/opt/hostedtoolcache/Python/3.10.19/x64/lib/python3.10/site-packages/tinyDA/sampler.py:114: UserWarning: subsampling_rate has been deprecated in favour of subchain_length.
Sampling chain 1/1
0%| | 0/12 [00:00<?, ?it/s]
Running chain, α_c = 1.000, α_f = 1.00: 0%| | 0/12 [00:00<?, ?it/s]
Running chain, α_c = 1.000, α_f = 1.00: 8%|▊ | 1/12 [00:00<00:05, 2.05it/s]
Running chain, α_c = 0.409, α_f = 0.50: 8%|▊ | 1/12 [00:00<00:05, 2.05it/s]
Running chain, α_c = 0.409, α_f = 0.50: 17%|█▋ | 2/12 [00:01<00:05, 1.97it/s]
Running chain, α_c = 0.372, α_f = 0.33: 17%|█▋ | 2/12 [00:01<00:05, 1.97it/s]
Running chain, α_c = 0.372, α_f = 0.33: 25%|██▌ | 3/12 [00:01<00:04, 2.00it/s]
Running chain, α_c = 0.328, α_f = 0.25: 25%|██▌ | 3/12 [00:01<00:04, 2.00it/s]
Running chain, α_c = 0.328, α_f = 0.25: 33%|███▎ | 4/12 [00:02<00:04, 1.95it/s]
Running chain, α_c = 0.306, α_f = 0.20: 33%|███▎ | 4/12 [00:02<00:04, 1.95it/s]
Running chain, α_c = 0.306, α_f = 0.20: 42%|████▏ | 5/12 [00:02<00:03, 1.95it/s]
Running chain, α_c = 0.302, α_f = 0.17: 42%|████▏ | 5/12 [00:02<00:03, 1.95it/s]
Running chain, α_c = 0.302, α_f = 0.17: 50%|█████ | 6/12 [00:03<00:03, 1.92it/s]
Running chain, α_c = 0.354, α_f = 0.29: 50%|█████ | 6/12 [00:03<00:03, 1.92it/s]
Running chain, α_c = 0.354, α_f = 0.29: 58%|█████▊ | 7/12 [00:03<00:02, 1.86it/s]
Running chain, α_c = 0.331, α_f = 0.25: 58%|█████▊ | 7/12 [00:03<00:02, 1.86it/s]
Running chain, α_c = 0.331, α_f = 0.25: 67%|██████▋ | 8/12 [00:04<00:02, 1.83it/s]
Running chain, α_c = 0.325, α_f = 0.33: 67%|██████▋ | 8/12 [00:04<00:02, 1.83it/s]
Running chain, α_c = 0.325, α_f = 0.33: 75%|███████▌ | 9/12 [00:04<00:01, 1.78it/s]
Running chain, α_c = 0.305, α_f = 0.40: 75%|███████▌ | 9/12 [00:04<00:01, 1.78it/s]
Running chain, α_c = 0.305, α_f = 0.40: 83%|████████▎ | 10/12 [00:05<00:01, 1.77it/s]
Running chain, α_c = 0.289, α_f = 0.36: 83%|████████▎ | 10/12 [00:05<00:01, 1.77it/s]
Running chain, α_c = 0.289, α_f = 0.36: 92%|█████████▏| 11/12 [00:05<00:00, 1.75it/s]
Running chain, α_c = 0.310, α_f = 0.42: 92%|█████████▏| 11/12 [00:05<00:00, 1.75it/s]
Running chain, α_c = 0.310, α_f = 0.42: 100%|██████████| 12/12 [00:06<00:00, 1.72it/s]
Running chain, α_c = 0.310, α_f = 0.42: 100%|██████████| 12/12 [00:06<00:00, 1.82it/s]
/opt/hostedtoolcache/Python/3.10.19/x64/lib/python3.10/site-packages/tinyDA/sampler.py:114: UserWarning: subsampling_rate has been deprecated in favour of subchain_length.
Sampling chain 1/1
0%| | 0/3 [00:00<?, ?it/s]
Running chain, α_c = 1.000, α_f = 1.00: 0%| | 0/3 [00:00<?, ?it/s]
Running chain, α_c = 1.000, α_f = 1.00: 33%|███▎ | 1/3 [00:01<00:03, 1.77s/it]
Running chain, α_c = 0.329, α_f = 1.00: 33%|███▎ | 1/3 [00:01<00:03, 1.77s/it]
Running chain, α_c = 0.329, α_f = 1.00: 67%|██████▋ | 2/3 [00:03<00:01, 1.79s/it]
Running chain, α_c = 0.380, α_f = 1.00: 67%|██████▋ | 2/3 [00:03<00:01, 1.79s/it]
Running chain, α_c = 0.380, α_f = 1.00: 100%|██████████| 3/3 [00:05<00:00, 1.78s/it]
Running chain, α_c = 0.380, α_f = 1.00: 100%|██████████| 3/3 [00:05<00:00, 1.78s/it]
/opt/hostedtoolcache/Python/3.10.19/x64/lib/python3.10/site-packages/tinyDA/sampler.py:114: UserWarning: subsampling_rate has been deprecated in favour of subchain_length.
Sampling chain 1/1
0%| | 0/1 [00:00<?, ?it/s]
Running chain, α_c = 1.000, α_f = 1.00: 0%| | 0/1 [00:00<?, ?it/s]
Running chain, α_c = 1.000, α_f = 1.00: 100%|██████████| 1/1 [00:03<00:00, 3.37s/it]
Running chain, α_c = 1.000, α_f = 1.00: 100%|██████████| 1/1 [00:03<00:00, 3.37s/it]
/opt/hostedtoolcache/Python/3.10.19/x64/lib/python3.10/site-packages/tinyDA/sampler.py:114: UserWarning: subsampling_rate has been deprecated in favour of subchain_length.
Sampling chain 1/1
0%| | 0/1 [00:00<?, ?it/s]
Running chain, α_c = 1.000, α_f = 1.00: 0%| | 0/1 [00:00<?, ?it/s]
Running chain, α_c = 1.000, α_f = 1.00: 100%|██████████| 1/1 [00:02<00:00, 2.58s/it]
Running chain, α_c = 1.000, α_f = 1.00: 100%|██████████| 1/1 [00:02<00:00, 2.58s/it]
Out[12]:
{'n_steps': 1000,
'n_dim': 3,
'move_fraction': 0.22622622622622623,
'hf_call_fraction': 0.023,
'n_hf_calls': 23,
'mean_subchain_length': 118.80597014925372,
'posterior_rmse': 2.488553138656431}
In [13]:
Copied!
samples_adapt = chain_adapt.samples
used_hf_adapt = chain_adapt.extras.used_hf
fig3, ax3 = plot_chain_2d(
samples_adapt[:, :2],
used_hf=used_hf_adapt,
theta_true=theta_true[:2],
title="Adaptive DA-MCMC: samples (first two parameters)",
show=True,
)
fig4, ax4 = plot_cumulative_hf_fraction(
used_hf_adapt,
title="Adaptive DA-MCMC: cumulative HF fraction",
show=True,
)
if chain_adapt.extras.subchain_length is not None:
fig5, ax5 = plot_subchain_length_history(chain_adapt.extras.subchain_length, show=True)
fig_s2, ax_s2 = plot_prediction_at_theta(
model_adapt.lf_model,
theta=theta_true,
y_obs=y_obs,
y_true=hf_forward(theta_true),
t=t,
title="Surrogate prediction (after DA-MCMC sampling)",
)
samples_adapt = chain_adapt.samples
used_hf_adapt = chain_adapt.extras.used_hf
fig3, ax3 = plot_chain_2d(
samples_adapt[:, :2],
used_hf=used_hf_adapt,
theta_true=theta_true[:2],
title="Adaptive DA-MCMC: samples (first two parameters)",
show=True,
)
fig4, ax4 = plot_cumulative_hf_fraction(
used_hf_adapt,
title="Adaptive DA-MCMC: cumulative HF fraction",
show=True,
)
if chain_adapt.extras.subchain_length is not None:
fig5, ax5 = plot_subchain_length_history(chain_adapt.extras.subchain_length, show=True)
fig_s2, ax_s2 = plot_prediction_at_theta(
model_adapt.lf_model,
theta=theta_true,
y_obs=y_obs,
y_true=hf_forward(theta_true),
t=t,
title="Surrogate prediction (after DA-MCMC sampling)",
)
