fn positive_range(data: &[f64], mask: &[u8]) -> (f64, f64)
Get min/max of positive values within mask