| 280 | |
| 281 | @pytest.mark.parametrize("T, m", test_data) |
| 282 | def test_mstump_include_discords(T, m): |
| 283 | for width in range(T.shape[0]): |
| 284 | for i in range(T.shape[0] - width): |
| 285 | include = np.asarray(range(i, i + width + 1)) |
| 286 | |
| 287 | excl_zone = int(np.ceil(m / 4)) |
| 288 | |
| 289 | ref_P, ref_I = naive.mstump(T, m, excl_zone, include, discords=True) |
| 290 | comp_P, comp_I = mstump(T, m, include, discords=True) |
| 291 | |
| 292 | npt.assert_almost_equal(ref_P, comp_P) |
| 293 | npt.assert_almost_equal(ref_I, comp_I) |
| 294 | |
| 295 | |
| 296 | @pytest.mark.parametrize("T, m", test_data) |