pub fn find_interval_range( range: (usize, usize), v: &[Float], pred: &dyn Fn(&[Float], usize) -> bool, ) -> usize