pytorch 算子调用kernel示例(MINIST)

pytorch 算子调用kernel示例(MINIST)

当进行 MNIST 分类任务时,PyTorch 中的每一个算子会根据设备类型(CPU 或 CUDA)自动选择合适的内核(kernel)进行计算。本文以GPU为例,介绍算子调用kernel的过程。

1. 模型定义与前向传播

python 复制代码
import torch
import torch.nn as nn
import torch.optim as optim
import torch.nn.functional as F

class SimpleCNN(nn.Module):
    def __init__(self):
        super(SimpleCNN, self).__init__()
        self.conv1 = nn.Conv2d(1, 32, kernel_size=3, stride=1, padding=1)  # 卷积层1
        self.conv2 = nn.Conv2d(32, 64, kernel_size=3, stride=1, padding=1)  # 卷积层2
        self.fc1 = nn.Linear(64 * 7 * 7, 128)  # 全连接层1
        self.fc2 = nn.Linear(128, 10)  # 全连接层2
    
    def forward(self, x):
        x = F.relu(self.conv1(x))  # 激活函数 ReLU
        x = F.max_pool2d(x, 2)  # 最大池化
        x = F.relu(self.conv2(x))  # 激活函数 ReLU
        x = F.max_pool2d(x, 2)  # 最大池化
        x = x.view(-1, 64 * 7 * 7)  # 张量展平
        x = F.relu(self.fc1(x))  # 全连接层激活
        x = self.fc2(x)  # 输出层
        return F.log_softmax(x, dim=1)  # 计算 softmax 概率

2. 数据加载与模型放置到 GPU

python 复制代码
from torchvision import datasets, transforms
from torch.utils.data import DataLoader

# 加载MNIST数据集
transform = transforms.Compose([transforms.ToTensor()])
train_dataset = datasets.MNIST('./data', train=True, download=True, transform=transform)
train_loader = DataLoader(train_dataset, batch_size=64, shuffle=True)

# 将模型和数据移动到 GPU
device = torch.device('cuda' if torch.cuda.is_available() else 'cpu')
model = SimpleCNN().to(device)  # 模型加载到 GPU
optimizer = optim.SGD(model.parameters(), lr=0.01, momentum=0.9)

# 训练过程
for epoch in range(1, 2):  # 运行1个 epoch
    model.train()
    for data, target in train_loader:
        data, target = data.to(device), target.to(device)  # 数据加载到 GPU
        optimizer.zero_grad()
        output = model(data)
        loss = F.nll_loss(output, target)
        loss.backward()
        optimizer.step()

3. 算子调用与 GPU 内核函数

当模型和数据都在 GPU 上时,PyTorch 的每一个算子会通过其调度机制(dispatch mechanism)调用相应的 CUDA 内核来加速计算。

(1) Conv2d (卷积层) 的 GPU 调用

在卷积操作中,PyTorch 会调用 Conv2d 算子。此时,设备被指定为 CUDA,调度系统会调用 GPU 专用的内核。

  • 算子调用:

    python 复制代码
    x = F.relu(self.conv1(x))
  • 内核调用

    • 对于 CUDA:调用 conv2d_cuda_kernel
  • CUDA 实现的原型代码conv2d_cuda_kernel 使用高度优化的 cuDNN 库来执行卷积操作,具体原型可能如下:

c++ 复制代码
Tensor conv2d_cuda(
    const Tensor& input,
    const Tensor& weight,
    const Tensor& bias,
    IntArrayRef stride,
    IntArrayRef padding,
    IntArrayRef dilation,
    int64_t groups) {
    
    // 使用 cuDNN 卷积函数来加速计算
    return at::cudnn_convolution(input, weight, bias, padding, stride, dilation, groups);
}

at::cudnn_convolution 是一个封装好的接口,用于调用 cuDNN 库中的卷积操作。cuDNN 提供了一系列高度优化的卷积算法,能够根据输入数据的大小和 GPU 的架构,自动选择最优的计算方式。

(2) ReLU (激活函数) 的 GPU 调用

ReLU 激活函数会根据设备类型,调用 CUDA 内核来进行计算。

  • 算子调用

    python 复制代码
    x = F.relu(self.conv1(x))
  • 内核调用

    • 对于 CUDA:调用 relu_cuda_kernel
  • CUDA 实现的原型代码

    cpp 复制代码
    Tensor relu_cuda(const Tensor& self) {
      // 使用 CUDA 进行并行计算 ReLU
      return at::clamp_min_cuda(self, 0);
    }
(3) MaxPool2d (池化层) 的 GPU 调用
  • 算子调用

    python 复制代码
    x = F.max_pool2d(x, 2)
  • 内核调用

    • 对于 CUDA:调用 max_pool2d_cuda_kernel
  • CUDA 实现的原型代码

    cpp 复制代码
    Tensor max_pool2d_cuda(
        const Tensor& self,
        IntArrayRef kernel_size,
        IntArrayRef stride,
        IntArrayRef padding,
        bool ceil_mode) {
      // 使用 CUDA 并行池化运算
      ...
    }
(4) Linear (全连接层) 的 GPU 调用
  • 算子调用

    python 复制代码
    x = F.relu(self.fc1(x))
  • 内核调用

    • 对于 CUDA:调用 linear_cuda_kernel
  • CUDA 实现的原型代码

    cpp 复制代码
    Tensor linear_cuda(
        const Tensor& input,
        const Tensor& weight,
        const Tensor& bias) {
      // 使用 CUDA 进行矩阵乘法和偏置加法
      return at::addmm_cuda(bias, input, weight.t());
    }
(5) Softmax (输出层) 的 GPU 调用
  • 算子调用

    python 复制代码
    return F.log_softmax(x, dim=1)
  • 内核调用

    • 对于 CUDA:调用 log_softmax_cuda_kernel
  • CUDA 实现的原型代码

    cpp 复制代码
    Tensor log_softmax_cuda(const Tensor& self, int64_t dim, bool half_to_float) {
      // 使用 CUDA 并行计算 softmax
      ...
    }

Reference:

  1. https://github.com/pytorch/pytorch
  2. https://docs.nvidia.com/deeplearning/cudnn/api/index.html
  3. https://github.com/pytorch/pytorch
相关推荐
bst@微胖子26 分钟前
Python高级语法之selenium
开发语言·python·selenium
Luis Li 的猫猫1 小时前
深度学习中的知识蒸馏
人工智能·经验分享·深度学习·学习·算法
查理零世1 小时前
【蓝桥杯集训·每日一题2025】 AcWing 6118. 蛋糕游戏 python
python·算法·蓝桥杯
带娃的IT创业者2 小时前
机器学习实战(6):支持向量机(SVM)——强大的非线性分类器
算法·机器学习·支持向量机
魔尔助理顾问3 小时前
一个简洁高效的Flask用户管理示例
后端·python·flask
java1234_小锋3 小时前
一周学会Flask3 Python Web开发-request请求对象与url传参
开发语言·python·flask·flask3
木觞清3 小时前
PyTorch与TensorFlow的对比:哪个框架更适合你的项目?
人工智能·pytorch·tensorflow
诚信爱国敬业友善6 小时前
常见排序方法的总结归类
开发语言·python·算法
wyg_0311137 小时前
用deepseek学大模型04-模型可视化与数据可视化
人工智能·机器学习·信息可视化
架构默片7 小时前
【JAVA工程师从0开始学AI】,第五步:Python类的“七十二变“——当Java的铠甲遇见Python的液态金属
java·开发语言·python