Function skewness_simd_py 

pub fn skewness_simd_py(
    data: &Bound<'_, PyArray1<f64>>,
    bias: bool,
) -> PyResult<f64>

Compute SIMD-optimized skewness (third standardized moment)