pub fn kernel_matrix(xs: &[Vec<f32>]) -> QuantumResult<Vec<Vec<f32>>>
Compute the full kernel matrix K[i,j] = k(xs[i], xs[j]).