def METHOD_NAME(self):
    g = cubic(shape=[10, 10, 10])
    hits = tools.find_coincident_nodes(g)
    assert hits == []
    g['node.coords'][1, :] = g['node.coords'][0, :]
    g['node.coords'][2, :] = g['node.coords'][0, :]
    g['node.coords'][4, :] = g['node.coords'][0, :]
    hits = tools.find_coincident_nodes(g)
    assert len(hits) == 1
    assert np.all(hits[0] == [0, 1, 2, 4])
    g['node.coords'][10, :] = g['node.coords'][20, :]
    hits = tools.find_coincident_nodes(g)
    assert len(hits) == 2
    # TODO: The following line seems to fail on the CI, as if the entries
    # in hits are in a different order. There is no reason I can see why
    # this would happen though. Need to figure this out.
    # assert np.all(hits[1] == [10, 20])