bitonic双调排序c代码和verilog实现

2024-04-04 18:18

本文主要是介绍bitonic双调排序c代码和verilog实现,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

这个算法除了两两比较和换位没有其他复杂操作,很适合在fpga上实现。

在fpga上如果需要排序的点比较多,实际应用就不能把这些点放到reg,只能放bram,而用bram,每个周期只能读写bram中的一项,现有网上的例子几乎所有都是资源不限制,纯仿真用不能实用的代码。

测试工程里点的个数是2048点,bram一项存两个点,每个点位宽是50bit,使用16个点的reg作为缓存,如果为加快速度,可以增加1项bram存的点个数,增加缓存点数,以资源为代价加快速度。

    实际结果完成2048个点排序需要大约3万周期,在时钟跑150M情况下,大约需要200us。

可以在我的github页面下载

C程序,VS2013的工程:https://github.com/tishi43/bitonic_my

包括bitonic的两种实现函数,bitonic()和bitonic2(),以及bitonic_fpga,模拟fpga读入16个数到缓存。

verilog的modelsim仿真工程:https://github.com/tishi43/bitonic_verilog

这个链接跳转有问题,直接把上面这个https地址拷贝到浏览器的地址栏去访问

资源使用情况,用到4238个LUT和1860个reg,使用performance high的综合策略。

+-------------------------+------+-------+-----------+-------+

|        Site Type        | Used | Fixed | Available | Util% |

+-------------------------+------+-------+-----------+-------+

| Slice LUTs*             | 4238 |     0 |    171900 |  2.47 |

|   LUT as Logic          | 4238 |     0 |    171900 |  2.47 |

|   LUT as Memory         |    0 |     0 |     70400 |  0.00 |

| Slice Registers         | 1860 |     0 |    343800 |  0.54 |

|   Register as Flip Flop | 1860 |     0 |    343800 |  0.54 |

|   Register as Latch     |    0 |     0 |    343800 |  0.00 |

| F7 Muxes                |    0 |     0 |    109300 |  0.00 |

| F8 Muxes                |    0 |     0 |     54650 |  0.00 |

+-------------------------+------+-------+-----------+-------+

最后贴两段C代码,两种bitonic的实现,看注释就知道原理了


//        5, 7,  15, 4,   0, 3,   11, 9,  12, 8,  1, 14, 13, 2, 6, 10//step 1: [5 7]  [4 15]   [0 3]   [9 11]  [8 12]  [1 14] [2 13] [6 10]//step 1之后相邻两个一组,两两排序,两个里面都正序排列//step 2,   4个一组,//round 1,每4个1组里,第一个与最后一个比较,第二个与倒数第二个比较,如此,5和15比,7和4比,//round 2, 再两个1组两两比较//step 2之后相邻每四个一组,组内正序排列//          ______//         |      |//round 1 [5 4 7 15]  [0 3 9 11]   [8 1 12 14]    [2 6 13 10]//           |_|//round 2 [4 5][7 15] [0 3][9 11]  [1 8] [12 14]  [2 6] [10 13]//step 3,round 1,8个1组比较,round 2, 4个1组比较,round 3,两个1组比较//round 1 第1个和倒数第一个比较,第二个和倒数第二个比较,round 2开始是间隔1个比较//step 3之后,相邻每8个一组,组内正序排列,//           _______________//          |               |//round 1  [4 5 3 0 15 7 9 11]         [1 8 6 2  14 12 10 13]//            |__________|//           ___//          |   |//round 2  [3 0 4 5]    [9 7 15 11]    [1 2 6 8]   [10 12 14 13]//            |___|//round 3  [0 3] [4 5]  [7 9] [11 15]  [1 2] [6 8] [10 12] [13 14]//step 4, round 1, 16个1组比较,round 2,8个1组比较,round 3,4个1组比较,round 4,2个1组比较//round 1 第1个和倒数第一个比较,第二个和倒数第二个比较,round 2开始是间隔1个比较//            ____________________________________________________________//           |                                                            |//round 1:  [ 0  3  4  5   7  6  2   1    15    11  9  8   10  12   13   14 ]//               |___________________________________________________|//              _____________//             |             |//round 2    [ 0  3 2 1      7  6 4 5 ]  [10 11 9 8  15 12 13 14]//                |_____________|//round 3    [0  1 2 3]    [4  5 7  6]  [9 8  10 11]   [13 12 15 14]//round 4    [0 1] [2 3]   [4 5] [6 7]  [8 9] [10 11]  [12 13] [14 15]void bitonic(unsigned int *data, int N)
{int i;//遍历stepint j;//遍历roundint k; //遍历groupint l; //遍历group里每对数据int ii;int steps = log2(N);int groups;int rounds;int pairs; //一组有几对数据int M; //一组有几个数据,M=2*pairsfor (i = 1; i <= steps; i++){rounds = i;for (j = 0; j < rounds; j++){//step1,rounds=1,groups=1,//step2,rounds=2,j=0,groups=4,j=1,groups=8//step3,rounds=3,j=0,groups=2,j=1,groups=4,j=2,groups=8groups = N/(1<<(rounds-j));M = N/groups;pairs = M/2;for (k = 0; k < groups; k++){for (l = 0; l < pairs; l++){if (j == 0){if (data[k*M + l] >= data[(k + 1)*M - l-1]) //只有round 1 是组内第一点和最后一点,第二点和倒数第二点这样比较{int temp = data[k*M + l];data[k*M + l] = data[(k + 1)*M - l-1];data[(k + 1)*M - l-1] = temp;}}else{if (data[k*M + l] >= data[k*M +M/2+l]) //round 2之后都是间隔M/2-1点之间的比较{int temp = data[k*M + l];data[k*M + l] = data[k*M + M / 2 + l];data[k*M + M / 2 + l] = temp;}}}}printf("step %2d round %2d: ",i,j+1);for (ii = 0; ii < N; ii++){printf("%4d ",data[ii]);}printf("\n");}}
}


//另一种实现,和第一种区别是组内没有第一个与最后一个比较,第二个与倒数第二个比较,取数规律都是一样的,这种适合fpga实现//但是比较有正序和逆序//        5, 7,  15, 4,   0, 3,   11, 9,  12, 8,  1, 14, 13, 2, 6, 10//         正      逆       正       逆    正      逆       正    逆//step 1: [5 7]  [15 4]   [0 3]   [11 9]  [8 12]  [14 1] [2 13] [10 6]//step 2,   4个数正,4个数逆,交替//           正          逆           正              逆//round 1 [5 4 15 7]  [11 9 0 3]   [8 1 14 12]    [10 13 2 6]//         正   正      逆    逆    正     正       逆    逆//round 2 [4 5][7 15] [11 9][3 0]  [1 8] [12 14]  [13 10] [6 2]//step 3,8个数正,8个数逆//                  正                       逆//round 1  [4 5 3 0    11 9 7  15]    [13 10 12 14    1 8 6 2]//            正          正                逆          逆//round 2  [3 0 4 5]   [7 9 11 15]    [13 14 12 10]   [6 8 1 2]//          正   正      正   正        逆     逆       逆  逆//round 3  [0 3][4 5]  [7 9][11 15]  [14 13][12 10]  [8 6][2 1]//step 4, 16个数正,16个数逆,也就是全部正//round 1:  [ 0  3  4  5   7  6  2 1    14  13  12  10  8  9  11 15 ]//round 2    [ 0  3 2 1    7  6 4 5 ]  [8 9 11 10    14 13 12 15]//round 3    [0  1 2 3]    [4  5 7  6]  [8 9  11 10]  [12 13 14 15]//round 4    [0 1] [2 3]   [4 5] [6 7]  [8 9] [10 11]  [12 13] [14 15]void bitonic2(unsigned int *data, int N)
{int i;//遍历stepint j;//遍历roundint k; //遍历groupint l; //遍历group里每对数据int ii;int steps = log2(N);int groups;int rounds;int pairs;         //一组有几对数据int M;             //一组有几个数据,M=2*pairsint ascend;        //1=升序 0=降序for (i = 1; i <= steps; i++){rounds = i;//step 1,最大组2个数据,//step 2,最大组4个数据//...for (j = 0; j < rounds; j++){//step 1,rounds=1,groups=1,//step 2,rounds=2,j=0,groups=4,j=1,groups=8//step 3,rounds=3,j=0,groups=2,j=1,groups=4,j=2,groups=8groups = N / (1 << (rounds - j));M = N / groups;pairs = M / 2;for (k = 0; k < groups; k++){//round 1, k第0bit决定升降//round 2, k第1bit决定升降//...ascend = ((k >> j) & 0x1) ==0;for (l = 0; l < pairs; l++){int swap = ascend ? data[k*M + l] >= data[k*M + M / 2 + l] :data[k*M + l] <= data[k*M + M / 2 + l];if (swap){int temp = data[k*M + l];data[k*M + l] = data[k*M + M / 2 + l];data[k*M + M / 2 + l] = temp;}}}printf("step %2d round %2d: ", i, j + 1);for (ii = 0; ii < N; ii++){printf("%4d ", data[ii]);}printf("\n");}}
}

这篇关于bitonic双调排序c代码和verilog实现的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/876507

相关文章

Python实现终端清屏的几种方式详解

《Python实现终端清屏的几种方式详解》在使用Python进行终端交互式编程时,我们经常需要清空当前终端屏幕的内容,本文为大家整理了几种常见的实现方法,有需要的小伙伴可以参考下... 目录方法一:使用 `os` 模块调用系统命令方法二:使用 `subprocess` 模块执行命令方法三:打印多个换行符模拟

SpringBoot+EasyPOI轻松实现Excel和Word导出PDF

《SpringBoot+EasyPOI轻松实现Excel和Word导出PDF》在企业级开发中,将Excel和Word文档导出为PDF是常见需求,本文将结合​​EasyPOI和​​Aspose系列工具实... 目录一、环境准备与依赖配置1.1 方案选型1.2 依赖配置(商业库方案)二、Excel 导出 PDF

Python实现MQTT通信的示例代码

《Python实现MQTT通信的示例代码》本文主要介绍了Python实现MQTT通信的示例代码,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一... 目录1. 安装paho-mqtt库‌2. 搭建MQTT代理服务器(Broker)‌‌3. pytho

使用zip4j实现Java中的ZIP文件加密压缩的操作方法

《使用zip4j实现Java中的ZIP文件加密压缩的操作方法》本文介绍如何通过Maven集成zip4j1.3.2库创建带密码保护的ZIP文件,涵盖依赖配置、代码示例及加密原理,确保数据安全性,感兴趣的... 目录1. zip4j库介绍和版本1.1 zip4j库概述1.2 zip4j的版本演变1.3 zip4

MySQL进行数据库审计的详细步骤和示例代码

《MySQL进行数据库审计的详细步骤和示例代码》数据库审计通过触发器、内置功能及第三方工具记录和监控数据库活动,确保安全、完整与合规,Java代码实现自动化日志记录,整合分析系统提升监控效率,本文给大... 目录一、数据库审计的基本概念二、使用触发器进行数据库审计1. 创建审计表2. 创建触发器三、Java

python生成随机唯一id的几种实现方法

《python生成随机唯一id的几种实现方法》在Python中生成随机唯一ID有多种方法,根据不同的需求场景可以选择最适合的方案,文中通过示例代码介绍的非常详细,需要的朋友们下面随着小编来一起学习学习... 目录方法 1:使用 UUID 模块(推荐)方法 2:使用 Secrets 模块(安全敏感场景)方法

Spring StateMachine实现状态机使用示例详解

《SpringStateMachine实现状态机使用示例详解》本文介绍SpringStateMachine实现状态机的步骤,包括依赖导入、枚举定义、状态转移规则配置、上下文管理及服务调用示例,重点解... 目录什么是状态机使用示例什么是状态机状态机是计算机科学中的​​核心建模工具​​,用于描述对象在其生命

Spring Boot 结合 WxJava 实现文章上传微信公众号草稿箱与群发

《SpringBoot结合WxJava实现文章上传微信公众号草稿箱与群发》本文将详细介绍如何使用SpringBoot框架结合WxJava开发工具包,实现文章上传到微信公众号草稿箱以及群发功能,... 目录一、项目环境准备1.1 开发环境1.2 微信公众号准备二、Spring Boot 项目搭建2.1 创建

IntelliJ IDEA2025创建SpringBoot项目的实现步骤

《IntelliJIDEA2025创建SpringBoot项目的实现步骤》本文主要介绍了IntelliJIDEA2025创建SpringBoot项目的实现步骤,文中通过示例代码介绍的非常详细,对大家... 目录一、创建 Spring Boot 项目1. 新建项目2. 基础配置3. 选择依赖4. 生成项目5.

Linux下删除乱码文件和目录的实现方式

《Linux下删除乱码文件和目录的实现方式》:本文主要介绍Linux下删除乱码文件和目录的实现方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录linux下删除乱码文件和目录方法1方法2总结Linux下删除乱码文件和目录方法1使用ls -i命令找到文件或目录