def METHOD_NAME():
    same = sft.streamlines[0]
    different = np.asarray([[1., 0., 0.],
                            [1., 0., 0.],
                            [1., 0., 0.]])
    similar = same + 0.0001
    compared = [same, different, similar]
    # Intersection: should find 2 similar.
    output, indices = perform_tractogram_operation_on_lines(
        intersection, [[same], compared])
    assert len(output) == 1
    # Intersection less precise. Should find 3 similar.
    # (but can't be tested now; returns the rounded unique streamline)
    output, indices = perform_tractogram_operation_on_lines(
        intersection, [[same], compared], precision=1)
    assert len(output) == 1
    # Difference: A - B: should return 0
    output, indices = perform_tractogram_operation_on_lines(
        difference, [[same], compared])
    assert len(output) == 0
    # Difference: B - A: should return 2
    output, indices = perform_tractogram_operation_on_lines(
        difference, [compared, [same]])
    assert len(output) == 2
    # Difference: B - A less precise: should return 1
    output, indices = perform_tractogram_operation_on_lines(
        difference, [compared, [same]], precision=1)
    assert len(output) == 1
    # Union: should combine the two same,
    output, indices = perform_tractogram_operation_on_lines(
        union, [[same], compared])
    assert len(output) == 3
    # Union less precise: should combine the similar too
    output, indices = perform_tractogram_operation_on_lines(
        union, [[same], compared], precision=1)
    assert len(output) == 2