pub fn argsort_radix_u64(data: &[u64], descending: bool) -> Vec<usize>
LSD Radix argsort for u64 - O(n·k) where k=8 (bytes)