def test_nichita_problem_5(model, expected_phase_split):
P = convert_bar_to_Pa(76)
T = 270
result = stability_test(model, P, T, model.mixture.z)
assert result.phase_split == expected_phase_split
def test_nichita_problem_2(model, pressure, expected_phase_split):
P = convert_bar_to_Pa(pressure)
T = 277.6
result = stability_test(model, P, T, model.mixture.z)
assert result.phase_split == expected_phase_split
def test_stochastic_consistency(sample_model):
num_of_realizations = 10
P = sample_model.P
T = sample_model.T
z = sample_model.z
for run in range(num_of_realizations):
result = stability_test(sample_model, P, T, z)
assert result.x.sum() == pytest.approx(1.0, rel=1e-5)
def test_nichita_problem_5(model, expected_phase_split):
P = convert_bar_to_Pa(76)
T = 270
result = stability_test(model,
P,
T,
model.mixture.z,
solver_args=PygmoSelfAdaptiveDESettings(20,
500,
seed=seed))
assert result.phase_split == expected_phase_split