pub fn scrambled_radical_inverse(base_index: u16, a: u64, perm: &[u16]) -> Float
Compute the radical inverse, but put each pixel through the permutation table for the given base.
