Function intersection_indices 

pub fn intersection_indices<P: TrxScalar>(
    a: &TrxFile<P>,
    b: &TrxFile<P>,
) -> Vec<usize>

Compute the intersection: streamlines present in both a and b. Returns indices into a.