def METHOD_NAME(self):
    rmg = RasterModelGrid((4, 5))
    rmg.add_empty("values", at="node")
    node_values = rmg.at_node["values"]
    node_values[:] = np.arange(rmg.number_of_nodes)
    link_values = maps.map_mean_of_link_nodes_to_link(rmg, "values")
    assert_array_equal(
        link_values,
        np.array(
            [
                0.5,
                1.5,
                2.5,
                3.5,
                2.5,
                3.5,
                4.5,
                5.5,
                6.5,
                5.5,
                6.5,
                7.5,
                8.5,
                7.5,
                8.5,
                9.5,
                10.5,
                11.5,
                10.5,
                11.5,
                12.5,
                13.5,
                12.5,
                13.5,
                14.5,
                15.5,
                16.5,
                15.5,
                16.5,
                17.5,
                18.5,
            ]
        ),
    )