fn cluster_parts( lmax: u32, k0: f64, radii: &[f64], epsilon: &[Complex], positions: &[[f64; 3]], ) -> Result<ClusterParts>