Function gaussian_smooth_3d_masked 

fn gaussian_smooth_3d_masked(
    data: &[f64],
    mask: &[u8],
    nx: usize,
    ny: usize,
    nz: usize,
    sigmas: &[f64; 3],
) -> Vec<f64>

3D Gaussian smoothing with anisotropic sigma