概述

怎样在MATLAB上做GPU计算呢?主要分为三个步骤:数据的初始化、对GPU数据进行操作、把GPU上的数据回传给CPU

一、数据的初始化

首先要进行数据的初始化。有两种方法可以进行初始化:一是先在CPU上设置好数据然后拷贝到GPU;二是直接在GPU上初始化数据。

1.1.先看看第一种方案:使用gpuArray函数。

N = 6;

M = magic(N);

G = gpuArray(M);

可见把数据从CPU拷贝到GPU上非常简单,只要B = gpuArray (A)就可以了。实际上MATLAB并没有规定一个矩阵定义之后不能改类型,所以以下也是可以的:

M = gpuArray(M);

补充说明: matlab是一种函数,用于产生魔方矩阵,它的每行、列以及对角线的数之和相等。该和的值为1+2+3+.....+n^2的和再除以n,n必须为大于或等于3的整数。例如:

有时候GPU受限于硬件架构,单精度的计算远快于双精度。这时候可以考虑在拷贝的时候顺便转换一下精度:

A = gpuArray(single(B));

1.2.第二种方法是直接在GPU上设置数据:

A = zeros(10, 'gpuArray');

也可以生成一个一维的随机数组:

r = gpuArray.rand(1, 100) % 一行,一百列

运行class函数:

class(r)

得到

ans = gpuArray。可见这是一个在gpu上的数组。

二、对GPU数据进行操作

下面就是对GPU数据进行操作。首先,MATLAB定义了GPU上丰富的库函数,比如快速傅立叶变换:

result = fft(r)

这样result就是另一个GPU上的数组,存储了对r做fft的结果。

三、把GPU上的数据回传给CPU--使用gather函数

最后是如何把GPU上的数据回传给CPU:

B = gather (A);

其中A是GPU上的数据,B是CPU上的数据。B的内容在回传之后等于A。

补充:如何查看自己的电脑能否用MATLAB进行GPU加速计算

并不是所有的电脑都可以用MATLAB进行GPU加速计算。想知道自己的电脑有没有这个能力,运行gpuDevice:

>> gpuDevice

MATLAB目前只支持Nvidia的显卡(即前面显示 CUDADevice with properties)。如果你的显卡是AMD的或者是Intel的,就得考虑另寻它路了。

参考链接:http://blog.sina.com.cn/s/blog_6f062c360102v9ic.html




Matlab+GPU加速学习笔记(一)

在matlab中使用GPU加速,来加速矩阵运算。

首先如前面所说,并不是所有GPU都能在maltab中进行加速的,貌似只有NVDIA的显卡可以吧。

硬件:GeForce GTX 980

软件:Matlab 2015a  (Matlab 2012以后的版本才带有GPU加速的工具箱)

下面开始介绍怎么玩GPU加速

第一步:在matlab命令窗口,运行gpuDevice,查看自己的显卡是否具备GPU加速功能

gpuDevice

正常情况下matlab会输出如下结果代表具备GPU加速功能:

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

CUDADevice with properties:

                    Name: 'GeForce GTX 980'

                   Index: 1

       ComputeCapability: '5.2'

          SupportsDouble: 1

           DriverVersion: 7.5000

          ToolkitVersion: 6.5000

      MaxThreadsPerBlock: 1024

        MaxShmemPerBlock: 49152

      MaxThreadBlockSize: [1024 1024 64]

             MaxGridSize: [2.1475e+09 65535 65535]

               SIMDWidth: 32

             TotalMemory: 4.2950e+09

         AvailableMemory: 3.3203e+09

     MultiprocessorCount: 16

            ClockRateKHz: 1380000

             ComputeMode: 'Default'

    GPUOverlapsTransfers: 1

  KernelExecutionTimeout: 1

        CanMapHostMemory: 1

         DeviceSupported: 1

          DeviceSelected: 1

第二步:CPU和GPU之间的数据交换,大致有如下几个函数和功能:

(1)、将数据从CPU中搬入GPU,函数:gpuArray,用法 M = gpuArray(M);

(2)、在GPU中计算完成后,将数据搬出到CPU存储,函数:gather,用法M=gather(M)。

下面给大家看看最简单的用法以及处理同样矩阵运算时的时间对比:

clear all
clc
M = rand(2000,2000);            % 生成一个随机矩阵
tic
[A1,B1] = eig(M);               % 求该随机矩阵的特征值和特征向量
t1=toctic
M = gpuArray(M);                % 将数据从CPU中搬到GPU
[A2,B2] = eig(M);               % 求特征值和特征向量
A2 = gather(A2);                % 将数据从GPU中搬到CPU
t2 = toc

输出结果(运行时间对比):

1

2

t1 = 14.1483

t2 =  9.1567

总结:虽然加速不是特别明显,不过还是非常给力的,需要说明的是,一般的GPU对于单精度型的数据运算加速更明显一点,因此我们在编写代码时,最好能够将数据从double型转换为single型,即M=single(M),然后在使用M=gpuArray(M),将数据搬迁到GPU进行计算,下面看看转换为单精度型后的时间对比:

1

2

3

4

5

6

7

8

9

10

11

12

13

clear all

clc

M = rand(2000,2000);            % 生成一个随机矩阵

tic

[A1,B1] = eig(M);                    % 求该随机矩阵的特征值和特征向量

t1=toc

tic

M = single(M);                     % 将数据转换为单精度型

M = gpuArray(M);                % 将数据从CPU中搬到GPU

[A2,B2] = eig(M);                 % 求特征值和特征向量

A2 = gather(A2);                 % 将数据从GPU中搬到CPU

t2 = toc

运行结果:

1

2

t1 =  14.8332

t2 =   5.0963

GPU对eig和svd进行加速,和cpu运行时间对比

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

clear all

A = rand(2000,2000);

tic

A1 = gpuArray(single(A));

[U,S,V] = svd(A1,'econ');

A2 = U*S*V';

A3 = gather(A2);

t = toc

tic

[U1,S1,V1] = svd(A,'econ');

A4 = U1*S1*V1';

t2 = toc

error = norm(A3-A4,'fro')

tic

A1 = gpuArray(single(A));

[vv,dd] = eig(A1);

t3 = toc

结果:

1

2

3

4

5

t =    3.1564

t2 =    7.1644

error =    0.0032

t3 =    5.2244

t4 =   17.6389

  

 

是不是加速更明显了呢?  

注:1、单精度型初始化矩阵函数例:X=rand(10,'single');  %定义在CPU上的一个10x10的随机初始化数组 

2、直接在GPU中初始化矩阵例:GX=rand(10,'gpuArray');   %直接在GPU设备上随机初始化一个10x10的数组 

  3、判断数据DATA是否存在于GPU内存中例:TF=existOnGPU(DATA)

4、计算效率统计例:t=gputimeit(F,N)         %返回执行F操作N次所需的时间,当N=1时可以缺省 




使用CUDA加速MATLAB程序

2018-01-25 17:38:19 qq_32464407 阅读数 6570  收藏

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。

本文链接:https://blog.csdn.net/qq_32464407/article/details/79164305

MATLAB十分擅长矩阵计算,如能借助CUDA工具箱调用NVIDIA GPU加速并行运算,能起到如虎添翼的效果。

最简单的调用方法是使用gpuArray()函数数据放到GPU上运算,再用gather()函数将结果取回(前提是运算支持gpuArray类型的变量,目前尝试过的矩阵乘法、求逆、pdist、pdist2函数均支持,可惜进行二次规划的quadprog并不支持。)。

不过,我最初在windows版MATLAB下调用gather()函数经常报错,然后便无法再MATLAB中调用GPU,使用reset(gpuDevice(1)))也还是报错。百度一番之后发现原来是windows系统的显卡超时检测和恢复(TDR)功能惹的祸。关闭TDR的方法是在HKLM\System\CurrentControlSet\Control\GraphicsDrivers下创建Dword值TdrLevel,并赋值为0  (https://answers.microsoft.com/zh-hans/windows/forum/windows_7-hardware/win7%E4%B8%AD%E5%A6%82%E4%BD%95%E9%85%8D%E7%BD%AE/69384e71-5075-4afe-a437-372425c0a3bb?auth=1)

用常用的pdist2函数进行测试:

gpuDevice(1)
A=rand(10000,10000);
B=rand(8000,10000);

tic; C = (pdist2((A),(B))); toc
tic; D = gather(pdist2(gpuArray(A),gpuArray(B))); toc;

结果:

ans =

CUDADevice - 属性:

Name: 'GeForce GTX 1050'
                     Index: 1
         ComputeCapability: '6.1'
            SupportsDouble: 1
             DriverVersion: 9.1000
            ToolkitVersion: 8
        MaxThreadsPerBlock: 1024
          MaxShmemPerBlock: 49152
        MaxThreadBlockSize: [1024 1024 64]
               MaxGridSize: [2.1475e+09 65535 65535]
                 SIMDWidth: 32
               TotalMemory: 4.2950e+09
           AvailableMemory: 3.4824e+09
       MultiprocessorCount: 5
              ClockRateKHz: 1493000
               ComputeMode: 'Default'
      GPUOverlapsTransfers: 1
    KernelExecutionTimeout: 0
          CanMapHostMemory: 1
           DeviceSupported: 1
            DeviceSelected: 1

时间已过 161.138259 秒。
时间已过 41.843864 秒。

可见GPU执行pdist2这种高度并行化的运算效率高出不少。




Matlab 并行计算学习初步

2018-08-09 08:19:57 如若明镜 阅读数 18430  收藏 更多

分类专栏: programming

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。

本文链接:https://blog.csdn.net/u011335616/article/details/41929033

Matlab 并行计算学习

1. 简介

高性能计算(High Performance Computing,HPC)是计算机科学的一个分支,研究并行算法和开发相关软件,致力于开发高性能计算机。可见并行计算是高性能计算的不可或缺的重要组成部分。

1.1 并行计算

并行计算(Parallel Computing)是指同时使用多种计算资源解决计算问题的过程,是提高计算机系统计算速度和处理能力的一种有效手段。它的基本思想是用多个处理器来协同求解同一问题,即将被求解的问题分解成若干个部分,各部分均由一个独立的处理机来并行计算。并行计算系统既可以是专门设计的、含有多个处理器的超级计算机,也可以是以某种方式互连的若干台的独立计算机构成的集群。通过并行计算集群完成数据的处理,再将处理的结果返回给用户[1]。

1.2 并行计算平台

平台是并行计算的载体,它决定着你可以用或只能用什么样的技术来实现并行计算。

多核和集群技术的发展,使得并行程序的设计成为提高数值计算效率的主流技术之一。常用的小型计算平台大致分为:由多核和多处理器构建的单计算机平台;由多个计算机组成的集群(Cluster)。前者通过共享内存进行数据交互,后者通过网络进行数据通信。

计算正在从 CPU(中央处理)向 CPU 与 GPU(协同处理)的方向发展。

GPU最早主要应用在图形计算机领域,近年来,它在通用计算机领域得到了迅猛的发展,使用GPU做并行计算已经变得越来越重要和高效。

常用的并行计算技术包括多线程技术、基于共享内存的OpenMP技术,基于集群的MPI技术等。但它们都需要用户处理大量与并行计算算法无关的技术细节,且不提供高效的算法库,与数值计算的关联较为松散。

1.3 Matlab与并行计算

Matlab即是一款数值计算软件,又是一门语言,它已经成为数值计算领域的主流工具。Matlab提供了大量高效的数值计算模块和丰富的数据显示模式,便于用户进行快速算法的研究和科学建模仿真。自Matlab 2009 之后,Matlab推出了并行计算工具箱(Parallel ComputingToolbox,PCT)和并行计算服务(Distributed ComputingServer,DCS),通过PCT和DCS用户可以实现基于多核平台、多处理器平台和集群平台的多种并行计算任务。利用PCT和DCS,用户无需关心多核、多处理器以及集群之间的底层数据通信,而是将主要精力集中于并行算法的设计。此外,PCT还增加了对GPU(Graphics ProcessingUnit)的支持。

2. Matlab并行计算初探

Matlab 并行计算工具箱,使用多核处理器,GPU和计算机集群,解决计算和数据密集型问题。高层次的结构并行for循环,特殊阵列类型和并行数值算法,让你在没有而CUDA或MPI编程基础的条件下并行MATLAB应用程序。您可以使用带有Simulink的工具箱并行运行一个模型的多个仿真。

该工具箱,通过在本地workers(MATLAB计算引擎)上执行应用程序,来使用多核台式机的全部处理能力。不用改变代码,就可以在计算机集群或网格计算服务(使用MATLAB分布式计算服务器™)上运行相同的应用程序。可以交互地或批处理地运行并行应用程序。

包含以下特性:

  • 并行for循环 (parfor),可以在多个处理器上运行任务并行算法。
  • 支持CUDA的NVIDIA GPU。
  • 通过本地workers,充分利用台式机的多核处理器。
  • 计算机集群和网格计算支持(通过MATLAB分布式计算服务)。
  • 以交互方式和批量方式执行并行应用程序。
  • 分布式阵列和SPMD(single-program-multiple-data),用于大型数据集的处理和数据并行算法。

本文仅介绍循环的并行,批处理和GPU并行计算。

2.1 并行池

2.1.1 开启和配置并行池

在使用parfor和spmd之前,首先要配置和开启并行计算池。parpool函数[1]提供并行池的配置和开启功能,其用法如下表 2‑1:

表 2‑1 parpool函数用法

此外,还可以给上述命令指定输出参数,如输入:p = parpool,输出如下图:

图2‑1 输入命令p=parpool时的运行结果

输出参数p为并行池对象,一个池对象有7中属性,上图显示了6种属性。其中,“NumWorkers” 代表worker的个数,也就是说开启了两个并行池;“IdleTimeout”指定了池定时关闭的时长,单位为分钟,图中为30分钟;“SpmdEnabled”指定池是否可以运行SPMD代码,其它属性相见Matlab 并行工具箱手册[2]。

此外,可以通过p.Property 访问和设置属性,如p.IdleTimeout = 10; 可以设置定时关闭池的时长为10分钟。

1.1.1.           获取当前池(gcp)

如果在启动池时没有设置输出参数p,可以使用gcp函数获取当前池的对象,gcp的用法如下表2-2:

表 2‑2 gcp用法

语法和参数

功能

p = gcp

获取当前池对象,如果未开启,则使用默认配置文件创建池。

p = gcp('nocreate')

获取当前池对象,如果未开启,返回空。

1.1.2.           关闭池(delete)

Matlab中,使用delete函数关闭或删除池,用法:delete(poolobj),其中poolobj为要关闭的池对象。

1.2.  循环并行parfor

1.2.1.           Matlab client和Matlab worker

图 2‑2 Matlab client 分配parfor并行任务到多个Matlab

Matlab并行执行parfor循环时,采用client和worker模式。其中client指编写和启动并行代码的Matlab段,worker指并行运行代码的Matlab端。用户可以将Matlab软件理解为一个进程,在同一计算机或网络上的多个计算机上可以运行多个Matlab进程,每个Matlab进程之间通过某种方式进行数据传输,用户首先在client端编写代码,然后再client端运行编写的代码,在运行代码的过程中,如果某个程序需要并行执行,Matlabclient端根据并行代码的关键字及类型将并行代码分配到本机或网络上的其它Matlab进程执行,那些最终执行并行代码的Matlab进程即为worker[3]。右图显示了,Matlab client 对parfor关键字的处理,它将任务分配到多个Matlab workers 上并行执行程序。假设m为Matlab worker的数量,循环次数为n,若n/m为整数的话,则循环将被均匀划分到Matlab worker;否则,循环被非均匀划分,有些Matlab worker的负载会大些。

1.2.2.           并行程序中的循环迭代parfor

应用parfor的一个基本前提是循环可以被等效分解,即分解后的循环的执行顺序不会影响最终结果。parfor的用法与for类似,其语法如下表 2‑3:

表2‑3 parfor的用法

语法和参数

功能

parfor loopvar = initval:endval,statements,end

在并行池中执行循环。其中initvalendval必须为正整数,或者是一个包含连续递增的整数行矢量。parfor i = range; <loop body>; end ,这样单行书写也是可以的,其中的“;”可以是“,”。

parfor (loopvar = initval:endval, M), statements, end

使用M指定执行循环体的Matlab worker的最大数量,M为非负整数,如果M=0,则在client上串行执行。

并行计算间的数据转移大小极限受Java虚拟机(JVM)内存分配的限制。在MJS集群上执行job的client和worker间的单次数据转移量也受此限制。具体大小限制取决于系统结构。如下表所示:

表 2‑4 每次转移数据的最大量与系统结构的关系

系统结构(System Architecture)

每次转移数据的最大数(近似值.)

64-bit

2.0 GB

32-bit

600 MB

1.2.3.           利用parfor并行for循环的步奏

使用parfor并行for循环的基本步骤为:

1)       使用parpool函数配置和开启并行计算池。

2)       将串行循环中的for关键字改为parfor,并注意是否要修改循环体,以满足特定要求,如循环变量的类型要求(参见参考手册)

3)       执行完毕后若不再进行并行计算,使用delete关闭并行池。

1.3.  批处理(batch)

batch函数完成函数或脚本在集群或台式机上的分载执行,也就是实现函数或脚本的并行。这里不再说明语法(详见Matlab 并行工具箱参考手册),而是举两个例子说明它的使用方法(此例子亦来自参考手册)。

1.3.1.           运行批处理任务

创建mywave.m脚本文件,输入如下代码并保存:

for i=1:1024

A(i) = sin(i*2*pi/1024);

end

在命令行窗口输入:job = batch('mywave'),Matlab以下图所示方式并行运行:

图 2‑3 批处理for循环

然后等待作业完成:wait(job),使用load命令将数据从Matlabworker 转移到Matlabclient的workspace,即load(job,’A’)。任务完成后,使用delete(job)永久删除它的数据。

1.3.2.           运行批处理并行循环

打开mywave.m脚本文件,修改for为parfor,保存,代码如下:

parfor i=1:1024

A(i) = sin(i*2*pi/1024);

end

输入:job =batch('mywave','Pool',2),以两个worker来执行parfor并行循环任务,所以此例总共使用3个本地worker,如下图所示:

图2‑4 批处理并行循环

同样的,依次使用wait(job)load(job,'A'),将数据从worker的workspace转移到client的workspace。最后,使用delet(job)删除数据。

批处理命令还有其它的运行方式,在此不再介绍。

1.4.  MATLAB的GPU计算

Matlab 的GPU计算实现MATLABworkspace和图形处理器(GPU)间的数据传递;在GPU上执行代码。

用户可以使用计算机的GPU处理矩阵运算,在大多数情况下,在GPU上执行要比在CPU上快。

Matlab的GPU计算提供了下表 2‑5所列的14种MATLAB函数和1种C语言函数,表中仅简要介绍函数功能,具体参见Matlab并行工具箱参考手册。

表2‑5 GPU计算函数列表

函数

功能

gpuArray

在GPU上创建阵列。

gather

把分布式阵列或gpuArray转移到Matlab的workspace。

existsOnGPU

确定gpuArray或CUDAKernel在GPU上可用。

gpuDevice

查询或选择GPU设备。

gpuDeviceCount

目前GPU的设备数量。

gputimeit

在GPU上运行函数所需时间。

reset

重置GPU设备并清空它的内存。

wait(GPUDevice)

等待GPU计算完成。

arrayfun

对GPU阵列上的每个元素应用函数。

bsxfun

用于GPU阵列的二元但扩展函数。

pagefun

对GPU阵列的每一页应用函数。

parallel.gpu.CUDAKernel

从PTX和CU代码创建 GPU CUDA kernel 对象。

feval

评估GPU的核。

setConstantMemory

设置GPU的一些常量内存。

mxInitGPU

这是一个C语言函数,在当前选择的设备上初始化Matlab GPU计算库。

Matlab的GPU计算提供了三种Matlab类和一种C类,如下表 2‑6

表2‑6 Matlab GPU计算提供的类

描述

gpuArray

存储在GPU上的阵列

GPUDevice

图形处理单元Graphics processing unit (GPU)

CUDAKernel

执行在GPU上的核

mxGPUArray

它是一个C类,MEX函数可以访问Matlab的gpuArray

关于各类的具体含义和各函数的详细用法请参见MATLAB并行工具箱手册。

下面主要从:GPU设备查询、在GPU上创建数组、在GPU上运行内置函数、在GPU上运行自定义函数四个方面说明Matlab GPU计算支持的技术,此外,用户还可以在GPU上运行CUDA或PTX代码、运行包含CUDA代码的MEX-函数等等,详见MATLAB并行工具箱手册和文献[3]的第10章。

1.4.1.   GPU 设备查询与选择

图 2‑5 GPU设备信息

GPU设备查询与选择

在使用GPU设备之前,可以利用gpuDeviceCount函数查看计算机当前可用GPU设备数量,若存在可用的GPU设备,则输出大于等于1的数。

在Matlab中使用gpuDevice(IDX)来选择第IDX个GPU设备,起始编号为1,如输入gpuDevice(1),输出右图所示结果:

其中,“MaxThreadsPerBlock”表示每个Block的最大线程数;“MaxShmemPerBlock”表示每个Block可用的最大共享内存大小;“TotalMemory”表示GPU拥有的所有内存,图中为2.1475e+09Byte,即2GB;“FreeMemory”代表可用内存;“MultiprocessorCount”表示GPU设备处理器个数;“DeviceSupported”表示本机Matlab支持的GPU设备数。

1.4.2.           在GPU上创建阵列

有两种方式在GPU上创建阵列:从Workspace转移和直接创建。

1)       从Workspace转移到GPU

使用gpuArray函数来把数据从Workspace转移到GPU,如下面的代码所示,将Workspace中的阵列A转移到GPU中,并存为G,原Workspace中的A阵列并不消失。

A = [1 2 3;4 5 6];

G = gpuArray(A);

该函数支持对非稀疏阵列且可以是:'single','double','int8','int16','int32','int64','uint8','uint16','uint32','uint64',or 'logical'类型。

2)       直接创建GPU阵列

gpuArray类提供了一些方法来在GPU上直接创建阵列,这样就不需要从Workspace进行转移。可用的有:

表 2‑7 gpuArray类提供的创建GPU阵列的方法

方法

方法

eye(___,'gpuArray')

rand(___,'gpuArray')

false(___,'gpuArray')

randi(___,'gpuArray')

Inf(___,'gpuArray')

randn(___,'gpuArray')

NaN(___,'gpuArray')

ones(___,'gpuArray')

gpuArray.linspace

true(___,'gpuArray')

gpuArray.logspace

zeros(___,'gpuArray')

gpuArray.colon

如:G = ones(100,100,50,'gpuArray');在GPU上创建一个100-by-100-by-50的阵列。

此外,支持gpuArray阵列的方法还有abs,cos,fft,conv,plot等等共369个方法。若要查看,输入如下命令即可查看:

methods('gpuArray')

如果,要查看某个方法的帮助,输入如下命令查看:

help gpuArray/functionname

其中,functionname为函数名,如:helpgpuArray/dot

3)       从GPU上回收数据

从GPU上回收数据,即把数据从GPU转移到MatlabWorkspace。在Matlab中使用gather函数来取回GPU上的数据并存在Matlab的 Workspace中。可以使用isequal函数来验证转移前后,值是否相同。如下代码将GPU中的阵列G传回到Workspace生成D。

G = gpuArray(ones(100,'uint32'));

D = gather(G);

OK = isequal(D,G)

1.4.3.           在GPU上运行内置函数

如 2)直接创建GPU阵列 所述,Matlab的369个函数都支持对GPU阵列的操作,下面是手册上的一个示例代码:

Ga = rand(1000,'single','gpuArray');

Gfft = fft(Ga);

Gb = (real(Gfft) + Ga) * 6;

G = gather(Gb);

代码主要实现对GPU上的阵列Ga进行FFT变换,并将结果回收到Workspace。可见内置函数的调用与普通串行程序没有大的区别,只是输入阵列的类型为gpuArray。

1.4.4.           在GPU上运行自定义函数

Matlab中使用arrayfunbsxfun来执行自定义的函数,即:

result = arrayfun(@myFunction,arg1,arg2);

下面举例来说明,创建myfun.m文件,打开输入如下代码:

function [ absga1,absga2,maxga ] = myfun( ga1,ga2 )

absga1 = abs(ga1);

absga2 = abs(ga2);

maxga = max(absga1,absga2);

该代码实现对输入阵列ga1,ga2的每个元素取绝对值,然后取出两者中对应最大的。为了运行并查看结果,创建一个脚本,输入如下代码:

ga1 = [-1 5 -3]

ga2 = gpuArray([4 -6 -2])

[absga1,absga2,maxga] = arrayfun(@myfun,ga1,ga2)

absda1 = gather(absga1)

absda2 = gather(absga2)

maxda = gather(maxga)

通过上述程序,可以看出,在自定义函数myfun中,只要有一个输入参数的类型为gpuArray,则整个函数myfun都将运行在GPU中,对于非gpuArray类型的阵列,Matlab自动进行转换。

除了可以调用Matlab自身支持的函数外,可以自己编写支持GPU的Matlab函数。自定义的函数中可以调用的函数和操作如下表 2 8:

表2‑8 自定义函数中支持gpuArray的函数和操作

abs

and

acos

acosh

acot

acoth

acsc

acsch

asec

asech

asin

asinh

atan

atan2

atanh

beta

betaln

bitand

bitcmp

bitget

bitor

bitset

bitshift

bitxor

ceil

complex

conj

cos

cosh

cot

coth

csc

csch

double

eps

eq

erf

erfc

erfcinv

erfcx

erfinv

exp

expm1

false

fix

floor

gamma

gammaln

ge

gt

hypot

imag

Inf

int8

int16

int32

int64

intmax

intmin

isfinite

isinf

isnan

ldivide

le

log

log2

log10

log1p

logical

lt

max

min

minus

mod

NaN

ne

not

or

pi

plus

pow2

power

rand

randi

randn

rdivide

real

reallog

realmax

realmin

realpow

realsqrt

rem

round

sec

sech

sign

sin

single

sinh

sqrt

tan

tanh

times

true

uint8

uint16

uint32

uint64

xor

+

-

.*

./

.\

.^

==

~=

<

<=

>

>=

&

|

~

&&

||

以下的标量扩展版:

*

/

\

^

分支指令:

break

continue

else

elseif

for

if

return

while

2.       总结

本文简要介绍了Matlab并行工具箱及其用法,包括多处理器的并行和GPU的并行,具体包括:并行计算池的开启,查询与关闭;for循环的并行;批处理;GPU计算。由于时间关系并未介绍分布式阵列和SPMD。

尝试并行化自己改进的Hough变换算法和FCM,可是并行后的效果很差(这和使用的电脑的配置有关),本文中并未给出,将来会进一步研究。

之所以选择高性能计算这门课,是因为在研究学习中,发现程序的运行速度是一个很重要和亟待解决的问题。已经尽量采用矢量化编程实现,可是它的本质是牺牲空间换时间,仍不能满足需求。高性能计算让我找到了突破点,将来会继续开展相关的工作。

[1] http://baike.baidu.com/view/1666.htm?fr=aladdin

[2] 并行工具箱在线手册:http://cn.mathworks.com/help/distcomp/index.html

[3] 实战 Matlab之并行程序设计[J]. 2012.

MATLAB上的GPU加速计算相关推荐

  1. MATLAB上的GPU加速计算——学习笔记

    MATLAB目前只支持Nvidia的显卡.如果你的显卡是AMD的或者是Intel的,就得考虑另寻它路了. MATLAB可谓工程计算中的神器,一方面它自带丰富的函数库,另一方面它所有的数据都是内建的矩阵 ...

  2. 如何在MATLAB上用GPU加速计算

    (1)首先,GPU设备确认,Matlab目前只支持Nvidia的显卡.想知道自己的电脑有没有这个能力,在Matlab中运行 gpuDevice. 只要数据格式是gpuArray格式的,那么计算过程会自 ...

  3. matlab 中使用 GPU 加速运算

    为了提高大规模数据处理的能力,matlab 的 GPU 并行计算,本质上是在 cuda 的基础上开发的 wrapper,也就是说 matlab 目前只支持 NVIDIA 的显卡. 1. GPU 硬件支 ...

  4. 如何在MATLAB上使用GUP加速跑代码

    如何在MATLAB上使用GUP加速跑代码 CPU和GPU的主要区别 查看CUDA版本并下载安装 怎么检查CUDA是否安装成功 确认MATLAB与cuda对应版本 在MATLAB查看GPU版本 测试gp ...

  5. Pytorch之GPU加速计算问题以及model=model.to(device)

    文章目录 model=model.to(device) GPU加速计算问题 技巧一 model=model.to(device) model=model.to(device) 这个有什么用,其实就是把 ...

  6. tensorflow的GPU加速计算

    参考    tensorflow的GPU加速计算 - 云+社区 - 腾讯云 一.概述 tensorflow程序可以通过tf.device函数来指定运行每一个操作的设备,这个设备可以是本地的CPU或者G ...

  7. float gpu 加速_(总结篇)使用 MATLAB GPU 加速计算|MATLAB 并行计算与分布式服务器|MATLAB技术论坛...

    本帖最后由 蓝云风翼 于 2013-12-18 17:28 编辑 注: 利用gpu加速有一下工具 1.JACKET 可从帖子中寻找 2.MATLAB a.并行计算工具箱 gpuArray,查看支持gp ...

  8. 光线追踪渲染实战(三):OpenGL 光线追踪,用 GPU 加速计算!

    项目代码仓库: GitHub:https://github.com/AKGWSB/EzRT gitee:https://gitee.com/AKGWSB/EzRT 目录 前言 0. 前置知识 1. 布 ...

  9. MATLAB使用GPU加速计算

    先上结论 1.对于特征值运算eig()函数来说,GPU的加速效果是很明显的 2.如果要求精度不高,采用单精度计算,加速效果更加明显 首先查看自己的电脑是不是支持GPU计算 在matlab的终端中输入: ...

最新文章

  1. pytorch筛选统计
  2. 专家手把手教你写出高水平个人简历
  3. java如何获得键值_如何在java中取map中的键值 的两种方法
  4. python 投票脚本_使用Python脚本拉取2014 CSDN博客之星投票情况
  5. spring 事务-使用@Transactional 注解(事务隔离级别)
  6. Java关键字final使用详解
  7. C语言---二进制和文本文件的备份
  8. 【POJ - 1269 】Intersecting Lines (计算几何,直线间的位置关系)
  9. 机器学习知识总结系列- 特征工程(1-1)
  10. CrawlSpider 详解
  11. EntityFramework 实体拆分与表拆分
  12. 【英语学习】【Level 07】U01 Making friends L2 A new friend from the past
  13. iOS-UICollectionView快速构造/拖拽重排/轮播实现
  14. js中java式的类成员
  15. 算法笔记_面试题_14. strStr 长字符串中查找短字符串
  16. Ubuntu/Linux备份/还原系统
  17. python add sheet_Python系列(6)——使用openpyxl将三种类型的数据添加进新sheet及已有的sheet中且不覆盖原数据...
  18. 如何模拟自动鼠标点击或记录鼠标和回放
  19. 20130723 上海OOW第二日
  20. memory management 蓝屏win10,要怎么解决这个蓝屏问题

热门文章

  1. 基本软件项目管理考试题目
  2. 第十五章 文正学院成绩录入系统
  3. Opencv学习笔记完整版
  4. JPG、GIF、PNG和BMP格式的图片各有什么优点和缺点
  5. 云效x钉钉:让研发工作更简单
  6. LINK : fatal error LNK1181: cannot open input file “libcurl.lib“
  7. AcWing-C/C++语法基础【合集1】
  8. php本地安装mysql报错443_Mysql性能剖析工具--Anemometer部署
  9. STM32单片机RGB红蓝调光植物补光系统红光蓝光PWM调色调节亮度
  10. sap服务器迁移性能问题,专家详解SAP数据迁移的六个方法