def METHOD_NAME(self):
    net = gen.cubic([4, 4, 1])
    net['edge.label'] = np.ones(24, dtype=bool)
    assert net['node.coords'].shape[0] == 16
    assert net['edge.conns'].shape[0] == 24
    assert np.all(net['edge.conns'][0] == [0, 1])
    net = ops.trim_nodes(net, 0)
    assert net['node.coords'].shape[0] == 15
    assert net['edge.conns'].shape[0] == 22
    assert net['edge.label'].shape[0] == 22
    assert np.all(net['edge.conns'][0] == [0, 1])
    net = ops.trim_nodes(net, [5, 10])
    assert net['node.coords'].shape[0] == 13
    assert net['edge.conns'].shape[0] == 15
    assert net['edge.label'].shape[0] == 15
    # ax = plot_edges(net['edge.conns'], net['node.coords'])
    # ax = plot_nodes(net['node.coords'], ax=ax)