SIMD指令简介

单指令多数据流,即SIMD(Single Instruction, Multiple Data)指一类能够在单个指令周期内同时处理多个数据元素的指令集,利用的是数据级并行来提高运行效率。

不同处理器架构所支持的SIMD指令集各不相同,比如嵌入式处理器ARM Cortex-A系列中支持NEON指令集。而我们常用的ARM Cortex-M系列单片机则使用了ARMv6-M/ARMv7-M架构,指令集较为简化,不支持NEON。其中Cortex-M4和Cortex-M7所使用的指令集也被称为ARMv7E-M,支持部分SIMD指令。对于Cortex-A系列来说,其指令集向下包含,即许多Cortex-M中的SIMD指令,也能够在Cortex-A中使用。

SIMD指令集的最初设计目标就是为了提高多媒体应用的性能(图像相关运算),因为图像数据的像素点都是8位数据,而CPU寄存器通常是32位以上的,因此一条指令只用于计算一个像素点不免显得非常浪费。在32位处理器上,可以用SIMD指令一次性计算4个8位数据,极大地提升了运算效率。

直接使用C语言编写的程序,编译器是不会自动翻译成SIMD指令的。因此为了在嵌入式处理器如单片机上提升图像数据运算速度,需要专门使用SIMD指令进行优化。得益于SIMD指令加速,我们能够在STM32F4上实现简单的光流算法,有兴趣的可以看一下PX4Flow的算法部分。为了能够在优化图像算法,使其能以较快速度运行于嵌入式平台,熟悉掌握SIMD指令是非常有必要的,下面列出一些可以在Cortex-M4和Cortex-M7上使用的SIMD指令(不全面,忽略了部分16位操作数指令),Cortex-A系列处理器基本上都兼容这些指令,当然在那上面还有强大的NEON指令集可以使用。

在STM32的库文件core_cm4_simd.h中可以看到所有SIMD指令的宏定义

部分SIMD指令说明

饱和加法: a+b=c,当计算结果大于c可表示的最大值或者小于c可表示的最小值时,计算结果取值为这个最大值或最小值。

非饱和加法: a+b=c,如果计算结果一出,则直接去掉一出位,剩下的就是结果。

  • __USAD8
    指令说明:无符号值的差的绝对值求和
    指令定义:uint32_t __USAD8 (uint32_t val1, uint32_t val2)
    指令操作:
absdiff1  = val1[7:0]   - val2[7:0]
absdiff2  = val1[15:8]  - val2[15:8]
absdiff3  = val1[23:16] - val2[23:16]
absdiff4  = val1[31:24] - val2[31:24]
res[31:0] = absdiff1 + absdiff2 + absdiff3 + absdiff4
  • __USADA8
    指令说明:无符号值的差的绝对值求和累加
    指令定义:uint32_t __USADA8 (uint32_t val1, uint32_t val2, uint32_t val3)
    指令操作:
absdiff1  = val1[7:0]   - val2[7:0]
absdiff2  = val1[15:8]  - val2[15:8]
absdiff3  = val1[23:16] - val2[23:16]
absdiff4  = val1[31:24] - val2[31:24]
sum       = absdiff1 + absdiff2 + absdiff3 + absdiff4
res[31:0] = sum[31:0] + val3[31:0]

并行加法

  • __SADD8
    指令说明:有符号8位操作数的并行加法
    指令定义:uint32_t __SADD8 (uint32_t val1, uint32_t val2 )
    指令操作:
res[7:0]   = val1[7:0]   + val2[7:0]
res[15:8]  = val1[15:8]  + val2[15:8]
res[23:16] = val1[23:16] + val2[23:16]
res[31:24] = val1[31:24] + val2[31:24]    
  • __QADD8
    指令说明:有符号8位操作数的并行饱和加法
    指令定义:uint32_t __QADD8 (uint32_t val1, uint32_t val2 )
    指令操作:
res[7:0]   = val1[7:0]   + val2[7:0]
res[15:8]  = val1[15:8]  + val2[15:8]
res[23:16] = val1[23:16] + val2[23:16]
res[31:24] = val1[31:24] + val2[31:24]   
  • __SHADD8
    指令说明:有符号8位操作数的并行减半加法
    指令定义:uint32_t __SHADD8 (uint32_t val1, uint32_t val2 )
    指令操作:
res[7:0]   = val1[7:0]   + val2[7:0]  >> 1
res[15:8]  = val1[15:8]  + val2[15:8] >> 1
res[23:16] = val1[23:16] + val2[23:16] >> 1
res[31:24] = val1[31:24] + val2[31:24] >> 1  
  • __UADD8
    指令说明:无符号8位操作数的并行加法
    指令定义:uint32_t __UADD8 (uint32_t val1, uint32_t val2 )
    指令操作:
res[7:0]   = val1[7:0]   + val2[7:0]
res[15:8]  = val1[15:8]  + val2[15:8]
res[23:16] = val1[23:16] + val2[23:16]
res[31:24] = val1[31:24] + val2[31:24]    
  • __UQADD8
    指令说明:无符号8位操作数的并行饱和加法
    指令定义:uint32_t __UQADD8 (uint32_t val1, uint32_t val2 )
    指令操作:
res[7:0]   = val1[7:0]   + val2[7:0]
res[15:8]  = val1[15:8]  + val2[15:8]
res[23:16] = val1[23:16] + val2[23:16]
res[31:24] = val1[31:24] + val2[31:24]    
  • __UHADD8
    指令说明:无符号8位操作数的并行减半加法
    指令定义:uint32_t __UHADD8 (uint32_t val1, uint32_t val2 )
    指令操作:
res[7:0]   = val1[7:0]   + val2[7:0]  >> 1
res[15:8]  = val1[15:8]  + val2[15:8] >> 1
res[23:16] = val1[23:16] + val2[23:16] >> 1
res[31:24] = val1[31:24] + val2[31:24] >> 1  

并行减法

  • __SSUB8
    指令说明:有符号8位操作数的并行减法
    指令定义:uint32_t __SSUB8 (uint32_t val1, uint32_t val2 )
    指令操作:
res[7:0]   = val1[7:0]   - val2[7:0]
res[15:8]  = val1[15:8]  - val2[15:8]
res[23:16] = val1[23:16] - val2[23:16]
res[31:24] = val1[31:24] - val2[31:24]
  • __QSUB8
    指令说明:有符号8位操作数的并行饱和减法
    指令定义:uint32_t __QSUB8 (uint32_t val1, uint32_t val2 )
    指令操作:
res[7:0]   = val1[7:0]   - val2[7:0]
res[15:8]  = val1[15:8]  - val2[15:8]
res[23:16] = val1[23:16] - val2[23:16]
res[31:24] = val1[31:24] - val2[31:24]
  • __SHSUB8
    指令说明:有符号8位操作数的并行减半减法
    指令定义:uint32_t __SHSUB8 (uint32_t val1, uint32_t val2 )
    指令操作:
res[7:0]   = val1[7:0]   - val2[7:0]   >> 1
res[15:8]  = val1[15:8]  - val2[15:8]  >> 1
res[23:16] = val1[23:16] - val2[23:16] >> 1
res[31:24] = val1[31:24] - val2[31:24] >> 1
  • __USUB8
    指令说明:无符号8位操作数的并行减法
    指令定义:uint32_t __USUB8 (uint32_t val1, uint32_t val2 )
    指令操作:
res[7:0]   = val1[7:0]   - val2[7:0]
res[15:8]  = val1[15:8]  - val2[15:8]
res[23:16] = val1[23:16] - val2[23:16]
res[31:24] = val1[31:24] - val2[31:24]
  • __UQSUB8
    指令说明:无符号8位操作数的并行饱和减法
    指令定义:uint32_t __UQSUB8 (uint32_t val1, uint32_t val2 )
    指令操作:
res[7:0]   = val1[7:0]   - val2[7:0]
res[15:8]  = val1[15:8]  - val2[15:8]
res[23:16] = val1[23:16] - val2[23:16]
res[31:24] = val1[31:24] - val2[31:24]
  • __UHSUB8
    指令说明:无符号8位操作数的并行减半减法
    指令定义:uint32_t __UHSUB8 (uint32_t val1, uint32_t val2 )
    指令操作:
res[7:0]   = val1[7:0]   - val2[7:0]   >> 1
res[15:8]  = val1[15:8]  - val2[15:8]  >> 1
res[23:16] = val1[23:16] - val2[23:16] >> 1
res[31:24] = val1[31:24] - val2[31:24] >> 1

指令集详细参考

Cortex-M4和Cortex-M7中的SIMD指令相关推荐

  1. 使用SIMD指令实现数据库运算(论文阅读笔记)

    摘要 现代cpu的指令允许在多个数据元素上并行执行基本操作.这些指令称为SIMD指令,因为它们将一条指令应用于多个数据元素.SIMD技术最初被内置到商业处理器中,以加速多媒体应用程序的性能.SIMD指 ...

  2. ST-LINK下载程序出现flash download failed - ‘cortex m4‘

    ST-LINK下载程序出现flash download failed - 'cortex m4' 可能原因:项目测试时使用ST-link下载程序调试,之前都好好的可以下载程序,之后可能做掉电存储等一些 ...

  3. ARM Cortex M4使用浮点运算单元(FPU)

    1.ARM Cortex M4   ARM Cortex-M4处理器是由ARM专门开发的最新嵌入式处理器,在M3的基础上强化了运算能力,新加了浮点.DSP.并行计算等.Cortex-M4处理器的最大亮 ...

  4. Cortex M4 中断向量表

    前言 重新学习Cortex M4的中断向量表内容.中断向量表是Cortex M4自带的功能,厂家可以根据需要自行定义,所以需要学习该内容,直接看Cortex M4内核相关即可.我这查看的是<Co ...

  5. Helix MP3 解码器在cortex m4应用

    Helix MP3 解码器在cortex m4应用 Helix MP3 解码器可浮点和定点实现.该算法可运行在任意 32 位定点处理器上,并完全使用 C 语言编码,可选择用优化的汇编指令替换某些代码段 ...

  6. 国内LORA首颗SOC芯片ASR6601ARM Cortex M4

    ASR6601是一种通用的LPWAN无线通信芯片组,集成射频收发器.调制解调器和32位RISC MCU.MCU使用ARM Cortex M4,频率为48 MHz工作频率.射频150MHZ-960MHZ ...

  7. stm32mp1 Cortex M4开发篇12:扩展板震动马达控制

    写在前面: 本文章为<ARM Cortex-M4裸机开发篇>系列中的一篇,,全系列总计14篇.笔者使用的开发平台为华清远见FS-MP1A开发板(STM32MP157开发板),Cortex- ...

  8. GCC中SIMD指令的应用方法

    X86架构上的多媒体应用开发,如果能够使用SIMD指令进行优化, 性能将大大提高.目前,IA-32的SIMD指令包括MMX,SSE,SSE2等几级. 在GCC的开发环境中,有几种使用SIMD指令的方式 ...

  9. 数字信号处理中的 SIMD

    系列文章目录 数字信号处理中的 SIMD Neon intrinsics 简明教程 文章目录 系列文章目录 0. 1. 前言 2. SIMD 是什么 3. SIMD 伪代码示例 4. SIMD 是如何 ...

最新文章

  1. 平方变换载波同步 matlab,matlab源码-costas载波同步环.docx
  2. 有道概率题:一个有趣的抛硬币问题
  3. 模块化和组件化的定义以及两者的区别
  4. WebRTC / Jitsi / 架构
  5. 配置Windows Server 2008群集
  6. 4-2cglib 动态代理
  7. 单元测试 : Googel test测试框架
  8. 30.Linux/Unix 系统编程手册(上) -- 线程:线程同步
  9. 破解 找回 lockdir 加密的文件
  10. 拉普拉斯平滑(Laplacian smoothing)
  11. JAVA自学知识点评定标准--自尚学堂马士兵
  12. msxml6_x86.msi和msxml6_ia64.msi和msxml6_x64.msi的选择
  13. Android studio Gradle:Executing tasks - 一直卡在Gradle:Build Running
  14. 【AI数学原理】函数求导(精髓篇)
  15. TS 版 Promise 详解
  16. python驱动photoshop_用Python用Photoshop打开PDF
  17. 线段树 (更新区间查询点)秋实大哥与小朋友
  18. 谷歌浏览器如何设置搜索结果的语言
  19. asp net夜话之八 数据绑定控件
  20. 清理Windows系统的C盘空间

热门文章

  1. excel 外部链接 乱码_在Excel文件中查找外部链接
  2. 初入职常见问题:弱网络环境构造、压力测试、内存泄漏、性能持续优化工具
  3. Robust Consistent Video Depth Estimation_具有鲁棒一致性的视频深度估计
  4. Arthas Spring Boot Starter工程启动报错
  5. 洛谷P1425小鱼的游泳时间【c语言】
  6. 网络编程培训之一 编程实现IP/TCP/UDP报文
  7. 博途PLC 1200/1500PLC 工艺对象PID PID_Compact详细解读
  8. mdx词典包_欧路词典PC端 词库安装 渲染
  9. 计算机 竞赛 甘肃,甘肃省第三届创新杯计算机应用能力大赛获奖名单
  10. 李笑来和 stormzhang,其实是一类人