Function remove_masked_mean 

fn remove_masked_mean(values: &mut [f64], mask: &[u8])

Subtract the mean over mask from values inside mask (outside is left untouched).