pub fn simd_binarize<F>(data: &Array2<F>, threshold: F) -> Result<Array2<F>>where F: Float + NumCast + SimdUnifiedOps,
SIMD-accelerated binarization with validation