def METHOD_NAME():
    # Recommended in scil_tractogram_math: use precision 0 to manage shifted
    # tractograms. Testing here.
    precision_shifted = 0
    same = sft.streamlines[0]
    shifted_same = same.copy() + 0.5
    different = np.asarray([[1., 0., 0.],
                            [1., 0., 0.],
                            [1., 0., 0.]])
    compared = [same, shifted_same, different]
    # Intersection: same/shifted
    output, indices = perform_tractogram_operation_on_lines(
        intersection_robust, [[same], compared], precision=precision_shifted)
    assert np.array_equal(indices, [0])
    assert len(output) == 1
    # Difference: different
    output, indices = perform_tractogram_operation_on_lines(
        difference_robust, [compared, [same]], precision=precision_shifted)
    logging.warning(indices)
    assert np.array_equal(indices, [2])
    assert len(output) == 1
    # Union: 4 (different, similar/same/shited)
    output, indices = perform_tractogram_operation_on_lines(
        union_robust, [compared, [same]], precision=precision_shifted)
    logging.warning(indices)
    assert len(output) == 2
    assert (indices == [0, 2]).all()