### Summary The history matching method `runs` only runs one wave; it would be useful to extend the functionality to: - Have multiple waves as part of the workflow - Allow an early stopping option with wave < n_waves if the number of emulated samples falls within the region of interest (e.g. 90%). ### Basic Example For an HM object `hmw`, we could: history_matching_results = hmw.run_waves(n_waves=10, n_simulations=100, n_test_samples=10000) history_matching_results will be a list of parameters and implausibility scores for each wave. ### Drawbacks Needs to be well tested. ### Unresolved questions _No response_ ### Implementation PR This is being work at #632 ### Reference Issues _No response_