11.1. 数组 API 支持(实验性)#
一些主要依赖 NumPy(而不是使用 Cython)来实现其 fit
、 predict
或 transform
方法的算法逻辑的 scikit-learn 估计器可以配置为接受任何兼容 Array API 的输入数据结构,并自动将操作分派到基础命名空间,而不是依赖于 NumPy。
在这个阶段,此支持被视为 实验性的,必须显式启用,如下所述。
Note
目前,只有 cupy.array_api
、 array-api-strict
、 cupy
和 PyTorch
已知可以与 scikit-learn 的估计器一起工作。
11.1.1. 示例用法#
以下是一个代码片段示例,演示如何使用 CuPy 在 GPU 上运行 LinearDiscriminantAnalysis
>>> from sklearn.datasets import make_classification
>>> from sklearn import config_context
>>> from sklearn.discriminant_analysis import LinearDiscriminantAnalysis
>>> import cupy
>>> X_np, y_np = make_classification(random_state=0)
>>> X_cu = cupy.asarray(X_np)
>>> y_cu = cupy.asarray(y_np)
>>> X_cu.device
<CUDA Device 0>
>>> with config_context(array_api_dispatch=True):
... lda = LinearDiscriminantAnalysis()
... X_trans = lda.fit_transform(X_cu, y_cu)
>>> X_trans.device
<CUDA Device 0>
模型训练后,数组形式的拟合属性也将来自与训练数据相同的 Array API 命名空间。例如,如果使用 CuPy 的
数组API命名空间用于训练,那么拟合的属性将在GPU上。我们提供了一个实验性的 _estimator_with_converted_arrays
工具,用于将估计器的属性从数组API转移到ndarray:
>>> from sklearn.utils._array_api import _estimator_with_converted_arrays
>>> cupy_to_ndarray = lambda array : array.get()
>>> lda_np = _estimator_with_converted_arrays(lda, cupy_to_ndarray)
>>> X_trans = lda_np.transform(X_np)
>>> type(X_trans)
<class 'numpy.ndarray'>
11.1.1.1. PyTorch支持#
通过设置 array_api_dispatch=True
并直接传递张量,支持PyTorch张量:
>>> import torch
>>> X_torch = torch.asarray(X_np, device="cuda", dtype=torch.float32)
>>> y_torch = torch.asarray(y_np, device="cuda", dtype=torch.float32)
>>> with config_context(array_api_dispatch=True):
... lda = LinearDiscriminantAnalysis()
... X_trans = lda.fit_transform(X_torch, y_torch)
>>> type(X_trans)
<class 'torch.Tensor'>
>>> X_trans.device.type
'cuda'
11.1.2. 支持 数组API
兼容输入#
scikit-learn中的估计器和其他工具支持数组API兼容输入。
11.1.2.1. 估计器#
decomposition.PCA
(使用svd_solver="full"
、svd_solver="randomized"
和power_iteration_normalizer="QR"
)linear_model.Ridge
(使用solver="svd"
)discriminant_analysis.LinearDiscriminantAnalysis
(使用solver="svd"
)preprocessing.KernelCenterer
preprocessing.MaxAbsScaler
preprocessing.MinMaxScaler
preprocessing.Normalizer
11.1.2.2. 元估计器#
接受数组API输入的元估计器,前提是基础估计器也支持:
model_selection.GridSearchCV
model_selection.RandomizedSearchCV
model_selection.HalvingGridSearchCV
model_selection.HalvingRandomSearchCV
11.1.2.3. 指标#
sklearn.metrics.cluster.entropy
11.1.2.4. 工具#
model_selection.train_test_split
覆盖范围预计会随时间增长。请关注GitHub上的专门 meta-issue 以跟踪进展。
11.1.2.5. 返回值类型和拟合属性#
当使用兼容Array API的输入调用函数或方法时,约定是返回与输入数据相同数组容器类型和设备的数组值。
同样,当使用兼容Array API的输入拟合估计器时,拟合属性将是与输入相同的库中的数组,并存储在同一设备上。随后的 predict
和 transform
方法期望输入来自与传递给 fit
方法的数据相同的数组库和设备。
但请注意,返回标量值的评分函数返回Python标量(通常是 float
实例)而不是数组标量值。
11.1.3. 常见估计器检查#
将 array_api_support
标签添加到估计器的标签集中,以表明它支持Array API。这将启用专门的检查作为常见测试的一部分,以验证在使用普通NumPy和Array API输入时估计器结果是否相同。
要运行这些检查,你需要在你的测试环境中安装
array_api_compat 。要运行完整的检查集,你需要安装 PyTorch 和 CuPy 并且拥有一个GPU。无法执行或缺少依赖项的检查将自动跳过。因此,使用
-v
标志运行测试以查看哪些检查被跳过是很重要的:
pip install array-api-compat # 以及根据需要安装其他库
pytest -k "array_api" -v
11.1.3.1. 关于 MPS 设备支持的说明#
在 macOS 上,PyTorch 可以使用 Metal Performance Shaders (MPS) 来访问硬件加速器(例如 M1 或 M2 芯片的内部 GPU 组件)。然而,在撰写本文时,PyTorch 的 MPS 设备支持还不完整。有关更多详细信息,请参阅以下 github 问题:
要在 PyTorch 中启用 MPS 支持,请在运行测试之前设置环境变量 PYTORCH_ENABLE_MPS_FALLBACK=1
:
PYTORCH_ENABLE_MPS_FALLBACK=1 pytest -k "array_api" -v
在撰写本文时,所有 scikit-learn 测试应该都能通过,然而,计算速度不一定比使用 CPU 设备更好。