Function rstats::vecvecg::partition_indexed

pub fn partition_indexed<T>(
    v: &[T], 
    pivot: T
) -> (Vec<usize, Global>, Vec<usize, Global>) where
    T: PartialOrd<T> + Copy,
    f64: From<T>, 

Partition by pivot gives two sets of indices. Items that are equal to pivot are ignored.