英特尔oneAPI - 使用parallel_for完成高斯消去的并行化

原创

已于 2022-05-17 21:34:23 修改 · 1k 阅读

于 2022-05-16 22:28:42 首次发布

本文介绍了如何利用英特尔OneAPI的DPC++实现高斯消去算法的并行化，重点在于第二重循环的并行化处理，通过`parallel_for`简化代码并提升效率。在OneAPI中，通过队列进行主机和设备间的通信，以达到在CPU或GPU上运行并行任务的目的。文章展示了部分关键代码，并分享了在英特尔DevCloud平台上进行开发的经验。

效果

在完成并行程序设计课程作业时，老师向我们介绍了一个平台——英特尔OneAPI。因此简单学习了一下OneAPI中DPC++的使用，并且用它完成了高斯消去并行化的作业，得到了不错的加速比。
在这里插入图片描述

英特尔oneAPI - 使用parallel_for完成高斯消去的并行化

OneAPI介绍

英特尔OneAPI平台是一个用于让程序员只需编写一份代码，就可以在诸如CPU、GPU、FPGA等多种设备、平台上运行的异构计算平台。

在这次高斯消去作业中，我们使用到的是英特尔OneAPI中的DPC++。

DPC++是英特尔对SYCL编程模型的一个实现，它允许我们编写跨平台的并行代码。

高斯消去算法

介绍下高斯消去算法，这个算法是将任意矩阵化为上三角矩阵的一种算法。

高斯消去

如图所示，高斯消去的典型串行算法具体步骤如下：

读入一个N*N的矩阵，设置变量i=0代表当前正在处理的行
通过行变换的方法，用第i行去消除i+1到N行的第i个元素。使得矩阵的i列中，i行以下的元素全为0
i自增，若i>N，则消元结束。

高斯消去算法C++串行实现

跟据这些步骤，我们可以用C++进行串行实现：

int mat[n][n]; // 待处理矩阵
for (int i = 0; i < n; i++) // 当前消去行
{
   
   
        for (int j = i + 1; j < n; j++) // 被消去行
        {
   
      
            float div = mat[j][i] / mat