Function difference_indices 

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

Compute the difference: streamlines in a but not in b. Returns indices into a.