(T, m)
| 182 | |
| 183 | @pytest.mark.parametrize("T, m", test_data) |
| 184 | def test_maamp(T, m): |
| 185 | excl_zone = int(np.ceil(m / 4)) |
| 186 | |
| 187 | for p in [1.0, 2.0, 3.0]: |
| 188 | ref_P, ref_I = naive.maamp(T, m, excl_zone, p=p) |
| 189 | comp_P, comp_I = maamp(T, m, p=p) |
| 190 | |
| 191 | npt.assert_almost_equal(ref_P, comp_P) |
| 192 | npt.assert_almost_equal(ref_I, comp_I) |
| 193 | |
| 194 | |
| 195 | @pytest.mark.parametrize("T, m", test_data) |