def METHOD_NAME():
    """Test multiple energy bins MAGIC case"""
    # create a basic meshtally
    coords = [[0, 1, 2], [-1, 3, 4], [10, 12]]
    flux_data = [[1.2, 3.3], [1.6, 1.7], [1.5, 1.4], [2.6, 1.0]]
    flux_error = [[0.11, 0.013], [0.14, 0.19], [0.02, 0.16], [0.04, 0.09]]
    tally = Mesh(structured=True, structured_coords=coords)
    tally.particle = "neutron"
    tally.e_bounds = [0.0, 0.5, 1.0]
    tally.n_flux = NativeMeshTag(2, float)
    tally.n_flux[:] = flux_data
    tally.n_rel_error = NativeMeshTag(2, float)
    tally.n_rel_error[:] = flux_error
    tolerance = 0.15
    null_value = 0.001
    magic(tally, "n_flux", "n_rel_error", tolerance=tolerance, null_value=null_value)
    expected_ww = [
        [0.2307692308, 0.5],
        [0.3076923077, 0.001],
        [0.2884615385, 0.001],
        [0.5, 0.15151515],
    ]
    assert_array_almost_equal(tally.ww_x[:], expected_ww[:])