def test_stat_contour():
dataset = MyDataset()
dataset.models.parameters["x"].frozen = True
fit = Fit(backend="minuit")
fit.optimize([dataset])
result = fit.stat_contour(datasets=[dataset], x="y", y="z")
assert result["success"]
x = result["y"]
assert_allclose(len(x), 10)
assert_allclose(x[0], 299, rtol=1e-5)
assert_allclose(x[-1], 299.133975, rtol=1e-5)
y = result["z"]
assert_allclose(len(y), 10)
assert_allclose(y[0], 0.04, rtol=1e-5)
assert_allclose(y[-1], 0.54, rtol=1e-5)
# Check that original value state wasn't changed
assert_allclose(dataset.models.parameters["y"].value, 300)
