写点什么

用 MXnet 实战深度学习之一:安装 GPU 版 mxnet 并跑一个 MNIST 手写数字识别

2015 年 11 月 24 日

【编者的话】本文是一系列深度学习的简单实战教程的第一篇,用 MXnet 做实现平台的实例代码,简单讲解深度学习常用的一些技术方向和实战样例。来自 Phunter 的博客,Phunter 的职业是数据科学家,DMLC ( http://dmlc.ml/ ) 工具包贡献者之一。

用 MXnet 实战深度学习之一:安装 GPU 版 mxnet 并跑一个 MNIST 手写数字识别

我想写一系列深度学习的简单实战教程,用 MXnet 做实现平台的实例代码简单讲解深度学习常用的一些技术方向和实战样例。这一系列的主要内容偏向于讲解实际的例子,从样例和代码里中学习解决实际问题。我会默认读者有一定神经网络和深度学习的基础知识,读者在这里不会看到大段推导和理论阐述。基础理论知识十分重要,如果读者对理论知识有兴趣,可以参看已有的深度学习教程补充和巩固理论基础,这里 http://deeplearning.net/reading-list/tutorials/ 有一些不错的理论教程,相关的理论知识在此不赘述。

MXnet: 轻量化分布式可移植深度学习计算平台

MXnet 是一群聪明勇敢勤劳的年轻计算机科学家实现的开源深度学习计算平台,它是 DMLC 分布式机器学习通用工具包 http://dmlc.ml/ 的重要部分(如果你知道 xgboost https://github.com/dmlc/xgboost 这个并行 GBT 的实现,应该对 DMLC 也不会陌生)。MXnet 的优点是,轻量化、可移植性高、也可轻松分布式并行,并且高效利用显存,更可以灵活的运行在移动设备上。它的代码和使用方法也简洁明了,适合学习实战。这么有意思的深度学习工具平台,大家快去点这个 github 连接给它加个星加个 fork 吧,传送门: https://github.com/dmlc/mxnet

安装 MXnet

MXnet 支持 Linux,Windows 和 Mac 平台。本文使用的主要平台是 ubuntu 14.04 LTS。提醒注意,这一些系列教程使用 CUDA 平台做 GPU 运算,而在本文写作的时候 CUDA 暂时还不支持最新的 ubuntu 15.10 版本的环境和编译器(主要是 gcc 5.2 的编译器),所以强烈建议坚守 14.04 LTS 版本或者是最多到 15.04 版。

安装环境可以是带 nVidia 显卡的实体机器或者是带 GPU 的云服务器。如果选择实体机,请不要通过虚拟机安装,比如原生 Windows 下面跑个虚拟的 Linux,因为多数虚拟机软件不支持直接调用本机显卡。如果选择云服务器,请一定选择 GPU instance 比如 AWS 的 g2.2xlarge 或 g2.8xlarge,或者是 terminal.com 的 GPU instance。注意:terminal.com 号称运行时可以改虚拟机的类型,但是纯 CPU 的虚拟机在运行时不能无缝切换到 GPU,建议一开始就选择 GPU instance。

以下安装步骤参考于官方文档: http://mxnt.ml/en/latest/build.html#building-on-linux ,本文根据 CUDA 的安装和实际操作略有修改。

基本依赖的安装

MXnet 的另一个优点就是它只需要很少的第三方包,它基本只需要 gcc 的编译器,BLAS 以及可选安装 OpenCV。这里如果还没有安装 git 可以顺道安装一下。

复制代码
sudo apt-get update
sudo apt-get install -y build-essential git libblas-dev libopencv-dev

下载 MXnet

复制代码
git clone --recursive https://github.com/dmlc/mxnet

这里提醒注意一定不要忘记–recursive 参数,因为 mxnet 依赖于 DMLC 通用工具包 http://dmlc.ml/ ,–recursive 参数可以自动加载 mshadow 等依赖。这里暂时不要着急编译,我们还要装一下 CUDA。

安装 CUDA

这里提到的 CUDA 安装方法也适用于除 MXnet 之外的其他深度学习软件包。我们通过 nVidia 官方链接下载安装 CUDA 驱动和工具包,请前往 https://developer.nvidia.com/cuda-downloads 选择对应的安装方式。国内读者建议网络安装方式 deb(network),这样 ubuntu 会选择就近的国内的源安装,速度可能比较快。

如果用 ubuntu 14.04,不用去官网,直接运行以下这些命令也可以调用官网下载(安装包较大需要耐心等待):

复制代码
wget http://developer.download.nvidia.com/compute/cuda/repos/ubuntu1404/x86_64/cuda-repo-ubuntu1404_7.5-18_amd64.deb
sudo dpkg -i cuda-repo-ubuntu1404_7.5-18_amd64.deb
sudo apt-get update
sudo apt-get install cuda

如果一切安装成功,可以用 nvidia-smi 命令查看你的显卡使用情况,一般空闲的显卡状态是这个样子的:

显卡型号取决于个人经济能力,不过 mxnet 的显存利用率高,一般一个 4G 的显卡就足够处理多数别的工具包要很多显存的问题。

可选安装:MXnet 也支持 cuDNN,它是 nVidia 推出的深度学习加速工具包,能高效实现一些卷积等深度学习常用操作,在内存使用和计算速度上面能有所提高。大家可以到这里 https://developer.nvidia.com/cudnn 申请开发者项目,如果批准通过可以下载安装 cuDNN 工具包,具体请参照 nVidia 官方教程。

编译支持 GPU 的 MXnet

MXnet 需要打开一个编译和链接选项来支持 CUDA。在前一步 git clone 得到的 mxnet/ 目录里找到 mxnet/make/ 子目录,把该目录下的 config.mk 复制到 mxnet/ 目录,用文本编辑器打开,找到并修改以下三行:

复制代码
USE_CUDA = 1
USE_CUDA_PATH = /usr/local/cuda
USE_BLAS = atlas

其中第二行是 CUDA 的安装目录。如果选择默认安装方式,它会在 /usr/local/cuda 或者是类似 /usr/local/cuda-7.5 这样的原始安装目录,如果是自定义目录的安装,请自行修改本条。这里的 BLAS 使用的是 ubuntu 的 atlas 实现,也请一并修改。

修改之后,在 mxnet/ 目录下编译(-j4 是可选参数表示用 4 线程编译):

复制代码
make -j4

注意:如果没有 CUDA 支持的显卡(比如 Intel 的 Iris 显卡或者 AMD 的 R 系列显卡)或者没有显卡,安装和编译 GPU 版本的 mxnet 会出错。解决方法是,把 USE_CUDA = 1 改回 USE_CUDA = 0,并确保 USE_OPENMP = 1,mxnet 会自动编译 CPU 版本并使用 OpenMP 进行多核 CPU 计算。根据问题的不同,GPU 版本对比 CPU 版一般会有 20-30 倍左右的加速。

安装 Python 支持

MXnet 支持 python 调用。简单来说就这么安装:

复制代码
cd python; python setup.py install

建议使用 python 2.7 版本,需要预先安装 setuptools 和 numpy。如果你的系统安装 Numpy 有些困难,可以考虑安装 Anaconda 或者 Miniconda 之类的 python 发行版:

复制代码
wget https://repo.continuum.io/miniconda/Miniconda-latest-Linux-x86_64.sh
bash Miniconda-latest-Linux-x86_64.sh
(确认回答若干安装问题后)
conda install numpy

运行 MNIST 手写数字识别

当 MXnet 一切安装好之后,可以试试看一下最简单的例子,MNIST 手写数字识别。MNIST 数据集包含 6 万个手写数字的训练数据集以及 1 万个测试数据集,每个图片是 28x28 的灰度图。在 mxnet/example/mnist 里可以找到 MXnet 自带 MNIST 的识别样例,我们可以先运行一下试试:

复制代码
cd mxnet/example/mnist
python mlp.py

mlp.py 会自动下载 MNIST 数据集,在第一次运行的时候耐心等待一下。

注意:mlp.py 默认使用 CPU,训练过程可以跑起来但是很慢。我们已经安装了 GPU,只需要修改一行代码,把 FeedForward 调用的 CPU 部分改成 GPU 即可让 MXnet 运行在 GPU 上:

复制代码
model = mx.model.FeedForward(
ctx = mx.cpu(), symbol = mlp, num_epoch = 20,
learning_rate = 0.1, momentum = 0.9, wd = 0.00001)

变成:

复制代码
model = mx.model.FeedForward(
ctx = mx.gpu(), symbol = mlp, num_epoch = 20,
learning_rate = 0.1, momentum = 0.9, wd = 0.00001)

再运行一下,是不是快多了呢?MXnet 的优点就是接口简洁。运行的时候,nvidia-smi 命令查看显卡状态差不多是这个样子的:

可以看到 python 进程在使用 GPU,因为这是个比较小的问题同时 MXnet 的显存优化较好,GPU 使用率 30% 到 40% 之间,显存占用 67MB。

可能出现的问题

运行 GPU 例子的时候可能会遇到这样的问题:

复制代码
ImportError: libcudart.so.7.0: cannot open shared object file: No such file 

这是因为没有把 CUDA 的动态链接库加入 PATH 里,解决方法是,可以在./bashrc 里面加入:

复制代码
export LD_LIBRARY_PATH=/usr/local/cuda-7.5/targets/x86_64-linux/lib/:$LD_LIBRARY_PATH

或者是在编译 MXnet 的时候,在 config.mk 里的

复制代码
ADD_LDFLAGS = -I/usr/local/cuda-7.5/targets/x86_64-linux/lib/
ADD_CFLAGS =-I/usr/local/cuda-7.5/targets/x86_64-linux/lib/

MNIST 代码简单讲解:设计一个最简单的多层神经网络

mlp.py 实现的是一个多层感知器网络(multilayer perceptron (MLP) )或者叫多层神经网络。在 MXnet 里,实现一个 MLP 首先需要定义一下这个 MLP 的结构,比如在代码里一个三层网络的 MLP 就是这样的:

复制代码
data = mx.symbol.Variable('data')
fc1 = mx.symbol.FullyConnected(data = data, name='fc1', num_hidden=128)
act1 = mx.symbol.Activation(data = fc1, name='relu1', act_type="relu")
fc2 = mx.symbol.FullyConnected(data = act1, name = 'fc2', num_hidden = 64)
act2 = mx.symbol.Activation(data = fc2, name='relu2', act_type="relu")
fc3 = mx.symbol.FullyConnected(data = act2, name='fc3', num_hidden=10)
mlp = mx.symbol.Softmax(data = fc3, name = 'mlp')

简单解释一下这几行代码:MNIST 的数据集每组数据是 28x28 的灰度图像,差不多如下图:

每组数据就可以表示成一个长度为 28x28=784 的一维数组,数组的每个元素是这个像素的灰度值。MLP 的每一层需要定义这一层节点的样式,比如 fc1 就是接受输入的第一层,它定义为一个全链接层 mx.symbol.FullyConnected,通过 data 接受输入,这一层包含了 128 个节点(num_hidden)。每一层也需要定义激活函数 Activation,比如第一层到第二层之间的激活函数就是 relu(代表 rectified linear unit 或者叫 Rectifier)ReLu 是深度神经网络里最常见的一个激活函数,主要因为计算函数相对容易和梯度下降不会发散,并且由于 MNIST 的问题比较稀疏更适合 ReLU。限于这里篇幅主要是为了介绍实现一个网络,关于 ReLU 的相关背景知识请参考 wikipedia 和其他相关教程。第二层网络 fc2 和第一层相似,它接受 fc1 的数据作为输入,输出给第三层。第三层网络 fc3 和前两层类似,不一样的是它是个结果输出层,产生的是输入图片对应于 0-9 总共 10 个数字里每个数字的概率,所以它的 num_hidden=10。

设计好了网络结构之后,MXnet 需要声明输入 feature 的格式,因为每个图片都是 28x28 大小,按照每个像素的灰度值展开成一列向量就是 784 维,我们可以告诉 mxnet 数据的输入尺寸是 784,mnist_iterator 是一个 python generator 一次提供 100 组数据给我们刚刚设计的 MLP,参见同目录的 data.py:

复制代码
train, val = mnist_iterator(batch_size=100, input_shape = (784,))

接下来就让 MXnet 建立并运行这个一个模型,就是这样简单,如果你会 scikit-learn 会感到很亲切,对不对(记得刚刚修改的指定 GPU 运行的那一行么?):

复制代码
model = mx.model.FeedForward(
ctx = mx.gpu(), symbol = mlp, num_epoch = 20,
learning_rate = 0.1, momentum = 0.9, wd = 0.00001)
model.fit(X=train, eval_data=val)

到这里,大家就基本会实现一个多层感知器 MLP,恭喜你们这是掌握深度学习的第一步。MXnet 的方式比 Caffe 等其他工具要写个配置文件简单的多了。工业界和学术界的多数深度学习的实际问题都是围绕着设计多层感知器展开,在结构设计激活函数设计等方面有很多有意思的问题。

有读者会问,MLP 是不是非要像 MNIST 手写数字识别这么设计。不是的,这个三层网络只是一个最简单的 MLP 的例子,这里每一层并不一定需要这样。设计一个更好更高效的多层神经网络和艺术一样没有止境。比如在 MNIST 同一个目录下的 lenet.py 就是用 Yann Lecun 设计的卷积网络实现数字识别,每层网络需要做的是 ConvolutionActivation 和 Pooling(如果想知道这三个具体是什么,请参看他的深度学习教程,以后的文章里面可能也会提到。

当做课后作业,读者可以自己试试调一下mlp.py 里不同的节点数和激活函数看看对数字识别率有什么提升,也可以增加num_epoch 调整learning_rate 等参数,在转发、评论或留言写下你们的设计方法和识别准确度(并没有奖励,嗯)。Kaggle 针对MNIST 数据集有一个教学比赛,读者可以用MXnet 训练一个自己的MNIST 模型,把结果提交上去比一比,记得说你是用MXnet 做的哟,传送门: https://www.kaggle.com/c/digit-recognizer

后记

这篇文章是这一系列的第一篇,我本意是想写个 MXnet 的 GPU 安装方法,后来想想加个例子讲解一下各种模型顺便当做另外一种深度学习入门教程吧。后续的一些文章会挑选 mxnet 自带的例子,介绍一些常见的有意思的深度学习模型,比如 RNN,LSTM,以及它们在 MXnet 里的实现,比如写个自动作词机模仿汪峰老师作词之类的。MXnet 这么有意思的深度学习工具平台,大家快去这个 github 连接给它加个星加个 fork 吧,传送门: https://github.com/dmlc/mxnet

编后语

《他山之石》是 InfoQ 中文站新推出的一个专栏,精选来自国内外技术社区和个人博客上的技术文章,让更多的读者朋友受益,本栏目转载的内容都经过原作者授权。文章推荐可以发送邮件到 editors@cn.infoq.com。


感谢杜小芳对本文的审校。

给InfoQ 中文站投稿或者参与内容翻译工作,请邮件至 editors@cn.infoq.com 。也欢迎大家通过新浪微博( @InfoQ @丁晓昀),微信(微信号: InfoQChina )关注我们,并与我们的编辑和其他读者朋友交流(欢迎加入 InfoQ 读者交流群(已满),InfoQ 读者交流群(#2))。

2015 年 11 月 24 日 17:1122231

评论

发布
暂无评论
发现更多内容

游戏夜读 | 游戏关卡设计师

game1night

架构演化

满山李子

大型互联网架构与集群技术

cxy

大型互联网应用系统使用技术方案和手段

如何进行高效学习

淡蓝色

深度思考 方法论 感悟 随笔杂谈

奈学:数据湖有哪些缺点?

奈学教育

数据湖

奈学:数据湖有哪些缺点?

古月木易

数据湖

如何学 Java,我说点不太一样的学习方式

四猿外

学习 程序员 个人成长 程序员成长

深入理解Kubernetes的Service:回归本源的场景需求

韩超

Kubernetes 微服务 服务

嗨,兄弟,别担心,这年头谁还没有一点焦虑!

攀鱼飞岩

管理 程序员人生 成长 个人感想 程序员素养

奈学:数据湖和数据仓库的区别有哪些?

古月木易

数据仓库 数据湖

读闲书自由和财务自由

池建强

读书 财务自由

聊一聊程序员如何增加收入

张小方

程序员 互联网 面试 副业赚钱 薪资

架构师训练营第 4 周——学习总结

在野

极客大学架构师训练营

轻松上手promise原理(2):then的简单实现

前端小帅

计算机操作系统基础(五)---Linux的进程管理

书旅

php 线程 多线程 操作系统 进程

实现简单的"纤程"

Near

谈谈架构和微服务<一>

Gabriel

架构 微服务 微服务架构 领域驱动设计 软件设计

SpringBatch系列之Remote-Chunking

稻草鸟人

大数据 Spring Boot SpringBatch 批量任务

ARTS 打卡 Week 05

teoking

创新管理体系标准ISO56002介绍

涛哥

数字化转型 创新

架构师训练营第四周课后作业

竹森先生

极客时间 极客大学架构师训练营

极客大学架构师训练营第四周学习总结

竹森先生

极客大学 极客大学架构师训练营

如何成为一名合格的 C/C++ 开发者?

张小方

c++ Linux 编程语言 架构设计 后端开发

我写了一本操作系统词典送给你

cxuan

操作系统 计算机

奈学:数据湖和数据仓库的区别有哪些?

奈学教育

数据仓库 数据湖

架构师训练营第 4 周作业

在野

极客大学架构师训练营

典型的大型互联网应用系统

Z冰红茶

第四章总结

MySQL 实战 45 讲笔记(2)-查询优化

程序员老王

MySQL

动态规划算法重点在于找上一个的公式,Google Code Review,John 易筋 ARTS 打卡 Week 06

John(易筋)

ARTS 打卡计划

用MXnet实战深度学习之一:安装GPU版mxnet并跑一个MNIST手写数字识别-InfoQ