DSP代码优化

最新推荐文章于 2024-02-07 17:33:58 发布
最新推荐文章于 2024-02-07 17:33:58 发布
阅读量1.2k 收藏 7
点赞数 2
分类专栏： DSP开发与技巧 文章标签： DSP优化

1 注意使用ccs自带的优化工具
1.1 选择恰当的编译器选项

  • 必须要用的选项 –O[2|3]
  • 可以使用-mt（要确保写的数据和读的数据在内存空间上没有重合）
  • -mh<num>  Specify speculative load byte count threshold
  • 如果源代码里含有永远不会执行的代码，使用选项-mo Place each function in a separate subsection
  • 如果考虑可执行程序的大小，加上-ms[0-3]。
  • 不要加上-g –gp –ss –ml3 –mu
  • -s[-k -al]、–mw 、(-on2 -o3)、-consultant 可以在产生分析信息的同时不影响生成代码的性能 

1.2 确保循环中的次数变量（一般for(i; i<N; i++), 此处的i）是有符号整数。

1.3 给编译器提供尽可能多的信息，如限定符restrict， 编译指示MUST_ITERATE、_nasserts()等。

1.4 内存中数据结构的对齐，可使用DATA_ALIGN、 STRUCT_ALIGN等

1.5 参考Compiler Consultant和.nfo文件的建议信息

2 利用优化器的意见
当编译选项中有-s时，在生成的*.asm文件中会有优化器的意见。如"C:/CCStudio_v3.3/C6000/cgtools/bin/cl6x" -g -k -s -on2 -o3 -mt -mw -mv6400+ --mem_model:data=near --consultant -@"Debug.lkf" "lesson_c.c"。其中，
void lesson_c(short *xptr, short *yptr, short *zptr, short *w_sum, int N)
{
    int i, w_vec1, w_vec2;
    short w1,w2;
    w1 = zptr[0];
    w2 = zptr[1];
    for (i = 0; i < N; i++) {
        w_vec1 =  xptr[i] * w1;
        w_vec2 =  yptr[i] * w2;
        w_sum[i] = (w_vec1 + w_vec2) >> 15;
    }
}

生成的优化器意见为：
;** --------------------------------------------------------------------------*
;** 27    -----------------------    w1 = *zptr;
;** 28    -----------------------    w2 = zptr[1];
;** 29    -----------------------    if ( N <= 0 ) goto g4;
;** --------------------------------------------------------------------------*
;**      -----------------------    U$17 = xptr;
;**      -----------------------    U$20 = yptr;
;**      -----------------------    U$26 = w_sum;
;** 31    -----------------------    L$1 = N;
;**      -----------------------    #pragma MUST_ITERATE(1, 1099511627775, 1)
;**      -----------------------    #pragma LOOP_FLAGS(4096u)
;**    -----------------------g3:
;** 31    -----------------------    *U$26++ = _mpy(*U$17++, w1)+_mpy(*U$20++, w2)>>15;
;** 29    -----------------------    if ( --L$1 ) goto g3;
;**    -----------------------g4:
;**      -----------------------    return;
从中可以看出，编译器自动加入了N是否为0的判断。如果改为：
void lesson_c(short *xptr, short *yptr, short *zptr, short *w_sum, int N)
{
    int i, w_vec1, w_vec2;
    short w1,w2;
    w1 = zptr[0];
    w2 = zptr[1];
    #pragma MUST_ITERATE(1)  //告诉编译器至少循环一次(编译器就不会再判断N是否为0)
    for (i = 0; i < N; i++)
    {
        w_vec1 =  xptr[i] * w1;
        w_vec2 =  yptr[i] * w2;
        w_sum[i] = (w_vec1 + w_vec2) >> 15;
    }
}
如下意见中已经没有了对N是否为0的判断：
;** --------------------------------------------------------------------------*
;** 27    -----------------------    w1 = *zptr;
;** 28    -----------------------    w2 = zptr[1];
;**      -----------------------    U$15 = xptr;
;**      -----------------------    U$18 = yptr;
;**      -----------------------    U$24 = w_sum;
;** 32    -----------------------    L$1 = N;
;**      -----------------------    #pragma MUST_ITERATE(1, 4294967295, 1)
;**      -----------------------    #pragma LOOP_FLAGS(4096u)
;**    -----------------------g2:
;** 32    -----------------------    *U$24++ = _mpy(*U$15++, w1)+_mpy(*U$18++, w2)>>15;
;** 30    -----------------------    if ( --L$1 ) goto g2;
;**      -----------------------    return;

3 利用软件流水信息优化循环
void lesson_c(short *xptr, short *yptr, short *zptr, short *w_sum, int N)
{
    int i, w_vec1, w_vec2;
    short w1,w2;
    w1 = zptr[0];
    w2 = zptr[1];
    for (i = 0; i < N; i++)
    {
        w_vec1 =  xptr[i] * w1;
        w_vec2 =  yptr[i] * w2;
        w_sum[i] = (w_vec1 + w_vec2) >> 15;
    }
}
上述代码的编译软件流水信息为：
;*----------------------------------------------------------------------------*
;*   SOFTWARE PIPELINE INFORMATION
;*
;*      Loop source line                 : 30                        ;循环开始的行数
;*      Loop opening brace source line   : 31
;*      Loop closing brace source line   : 35
;*      Known Minimum Trip Count         : 1                ;循环最小次数
;*      Known Max Trip Count Factor      : 1                ;循环的因子，循环次数是循环因子的倍数，如果知道循环因子的话，便于编译器自动铺开（unroll）代码
;*      Loop Carried Dependency Bound(^) : 0           ;内存读写瓶颈，如果有的话，后面的汇编代码注释里相应语句含有^标志
;*      Unpartitioned Resource Bound     : 2               ;资源瓶颈
;*      Partitioned Resource Bound(*)    : 2
;*      Resource Partition:
;*                                   A-side   B-side
;*      .L units                         0        0     
;*      .S units                        1        0     
;*      .D units                        2*       1     
;*      .M units                       1        1     
;*      .X cross paths             1        0     
;*      .T address paths         2*       1     
;*      Long read paths          0        0     
;*      Long write paths          0        0     
;*      Logical  ops (.LS)         0        0     (.L or .S unit)
;*      Addition ops (.LSD)      1        0     (.L or .S or .D unit)
;*      Bound(.L .S .LS)           1        0     
;*      Bound(.L .S .D .LS .LSD)     2*       1                 ;资源使用不平衡（没有完全利用可用的计算能力）
;*
;*      Searching for software pipeline schedule at ...
;*      ii = 2  Schedule found with 6 iterations in parallel
;*
;*      Register Usage Table:
;*          +-----------------------------------------------------------------+
;*          |AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA|BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB|
;*          |00000000001111111111222222222233|00000000001111111111222222222233|
;*          |01234567890123456789012345678901|01234567890123456789012345678901|
;*          |--------------------------------+--------------------------------|
;*       0: |   ******                       |    * **                        |
;*       1: |   *  ***                       |    ****                        |
;*          +-----------------------------------------------------------------+
;*
;*      Done
;*
;*      Loop will be splooped
;*      Collapsed epilog stages     : 0
;*      Collapsed prolog stages     : 0
;*      Minimum required memory pad : 0 bytes
;*
;*      Minimum safe trip count     : 1
;*----------------------------------------------------------------------------*
;*        SINGLE SCHEDULED ITERATION                    ;需要加-mw选项
;*
;*        $C$C23:
;*   0              LDH     .D2T2   *B6++,B5          ; |32|              ;一次装载16bit，浪费了剩余16bit带宽（假设数据带宽为32bit）
;*   1              LDH     .D1T1   *A6++,A5          ; |32| 
;*   2              NOP             3
;*   5              MPY     .M2     B5,B7,B4          ; |32| 
;*   6              MPY     .M1     A5,A8,A4          ; |32| 
;*   7              NOP             1
;*   8              ADD     .L1X    B4,A4,A3          ; |32| 
;*   9              SHR     .S1     A3,15,A3          ; |32| 
;*  10              STH     .D1T1   A3,*A7++          ; |32| 
;*         ||           SPBR            $C$C23
;*  11              NOP             1
;*  12              ; BRANCHCC OCCURS {$C$C23}        ; |30|          ;一次循环需要12时钟周期（没有展开）
;*----------------------------------------------------------------------------*

修改为下面的代码时：
#define WORD_ALIGNED(x)  (_nassert(((int)(x) & 0x3) == 0))
#define DWORD_ALIGNED(x) (_nassert(((int)(x) & 0x7) == 0))
void lesson3_c(short * restrict xptr, short * restrict yptr, short *zptr, short *w_sum, int N)
{
    int i, w_vec1, w_vec2;
    short w1, w2;
    WORD_ALIGNED(xptr); //保证内存装载的带宽
    WORD_ALIGNED(yptr);                                   
    w1 = zptr[0];
    w2 = zptr[1];
    #pragma MUST_ITERATE(48, , 2);   //已知最小循环次数=48，循环因子factor=2, 可以铺开代码
    for (i = 0; i < N; i++)
    {
        w_vec1 =  xptr[i] * w1;
        w_vec2 =  yptr[i] * w2;
        w_sum[i] = (w_vec1 + w_vec2) >> 15;
    }
}
上述代码的编译软件流水信息为：
;*----------------------------------------------------------------------------*
;*   SOFTWARE PIPELINE INFORMATION
;*
;*      Loop source line                 : 59
;*      Loop opening brace source line   : 60
;*      Loop closing brace source line   : 64
;*      Loop Unroll Multiple             : 4x                       ;循环铺开的次数
;*      Known Minimum Trip Count         : 12                    
;*      Known Max Trip Count Factor      : 1
;*      Loop Carried Dependency Bound(^) : 0
;*      Unpartitioned Resource Bound     : 4
;*      Partitioned Resource Bound(*)    : 4
;*      Resource Partition:
;*                                A-side   B-side
;*      .L units                     0        0     
;*      .S units                     2        2     
;*      .D units                     4*       4*                     ;铺开循环保证了资源使用的平衡
;*      .M units                     2        2     
;*      .X cross paths               2        2     
;*      .T address paths             4*       4*    
;*      Long read paths              0        0     
;*      Long write paths             0        0     
;*      Logical  ops (.LS)           0        0     (.L or .S unit)
;*      Addition ops (.LSD)          2        2     (.L or .S or .D unit)
;*      Bound(.L .S .LS)             1        1     
;*      Bound(.L .S .D .LS .LSD)     3        3     
;*
;*      Searching for software pipeline schedule at ...
;*         ii = 4  Schedule found with 4 iterations in parallel
;*
;*      Register Usage Table:
;*          +-----------------------------------------------------------------+
;*          |AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA|BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB|
;*          |00000000001111111111222222222233|00000000001111111111222222222233|
;*          |01234567890123456789012345678901|01234567890123456789012345678901|
;*          |--------------------------------+--------------------------------|
;*       0: |    *****        ***            |    ******      ****            |
;*       1: |   *******         *            |    ******      ** *            |
;*       2: |   *******      ****            |    * ****      **              |
;*       3: |   *******      ****            |    ******      **              |
;*          +-----------------------------------------------------------------+
;*
;*      Done
;*
;*      Loop will be splooped
;*      Collapsed epilog stages     : 0
;*      Collapsed prolog stages     : 0
;*      Minimum required memory pad : 0 bytes
;*
;*      Minimum safe trip count     : 1 (after unrolling)
;*----------------------------------------------------------------------------*
;*        SINGLE SCHEDULED ITERATION
;*
;*        $C$C24:
;*   0              LDW     .D1T1   *A7++(8),A9       ; |61|             ;每次装载一个字
;*     ||           LDW     .D2T2   *B6++(8),B17      ; |61| 
;*   1              LDW     .D1T1   *A4++(8),A3       ; |61| 
;*   2              NOP             1
;*   3              LDW     .D2T2   *B16++(8),B4      ; |61| 
;*   4              NOP             2
;*   6              MPY2    .M1     A9,A8,A17:A16     ; |61| 
;*   7              MPY2    .M1     A3,A8,A19:A18     ; |61| 
;*     ||           MPY2    .M2     B17,B7,B5:B4      ; |61| 
;*   8              MPY2    .M2     B4,B7,B19:B18     ; |61| 
;*   9              NOP             2
;*  11              ADD     .L2X    B4,A16,B17        ; |61| 
;*  12              ADD     .L2X    B18,A18,B5        ; |61| 
;*     ||           SHR     .S2     B17,15,B4         ; |61| 
;*     ||           ADD     .L1X    B5,A17,A3         ; |61| 
;*  13              SHR     .S2     B5,15,B4          ; |61| 
;*     ||           ADD     .L1X    B19,A19,A19       ; |61| 
;*     ||           STH     .D2T2   B4,*B9++(8)       ; |61| 
;*     ||           SHR     .S1     A3,15,A18         ; |61| 
;*  14              STH     .D2T2   B4,*B8++(8)       ; |61| 
;*     ||           SHR     .S1     A19,15,A9         ; |61| 
;*     ||           STH     .D1T1   A18,*A5++(8)      ; |61| 
;*  15              STH     .D1T1   A9,*A6++(8)       ; |61| 
;*     ||           SPBR            $C$C24
;*  16              ; BRANCHCC OCCURS {$C$C24}        ; |59|       ;4次循环使用16个时钟周期（循环展开，一次循环展开为4次）
;*----------------------------------------------------------------------------*

4 consultant advice 和 *.nfo文件
当编译时加上--consultant和–on2 –o3，可以查看相应的consultant advice 和*.nfo文件。打开profile，运行程序，这是可以查看viewer的consultant内容。*.nfo文件编译时便可以生成，内容没有consultant全。下面是lesson3_c.nfo的内容

TMS320C6x C/C++ Optimizer               v6.0.8
Build Number 1GKUL-JA0KH827-RSAQQ-TAV-ZAZG_W_Q_Y

        ======File-level Analysis Summary======

extern void _lesson3_c() is called from 0 sites in this file.
    It appears to be inlineable (size = 58 units)
    It calls these functions:
    <NONE>

        ======= End file-level Analysis =======

extern void _lesson3_c() is called from 0 sites in this file.
    It appears to be inlineable (size = 58 units)
    It calls these functions:
    <NONE>

ADVICE: In function lesson3_c()
    in the 'for' loop with loop variable 'i' at lines 39-44
    for the statement w_sum[i] = _mpy(xptr[i], w1)+_mpy(yptr[i], w2)>>15; at line 41

    The address of w_sum[i] for the first iteration of the loop is &w_sum[0].
    This pointer is aligned to a 16 bit boundary.

    Consider adding an assertion just before the loop:

        _nassert( ((int)w_sum % 4) == 0 );  /* 32-bit aligned */
       or    _nassert( ((int)w_sum % 8) == 0 );  /* 64-bit aligned */

    to specify that multiple elements of w_sum[i]
    may be accessed in parallel.

ADVICE: In function lesson3_c()
    in the 'for' loop with loop variable 'i' at lines 39-44
    for the statement w_sum[i] = _mpy(xptr[i], w1)+_mpy(yptr[i], w2)>>15; at line 41

    The address of yptr[i] for the first iteration of the loop is &yptr[0].
    This pointer is aligned to a 32 bit boundary.

    Consider adding an assertion just before the loop:

        _nassert( ((int)yptr % 8) == 0 );  /* 64-bit aligned */

    to specify that multiple elements of yptr[i]
    may be accessed in parallel.

ADVICE: In function lesson3_c()
    in the 'for' loop with loop variable 'i' at lines 39-44
    for the statement w_sum[i] = _mpy(xptr[i], w1)+_mpy(yptr[i], w2)>>15; at line 41

    The address of xptr[i] for the first iteration of the loop is &xptr[0].
    This pointer is aligned to a 32 bit boundary.

    Consider adding an assertion just before the loop:

        _nassert( ((int)xptr % 8) == 0 );  /* 64-bit aligned */

    to specify that multiple elements of xptr[i]
    may be accessed in parallel.
<<NULL MIX DOMAIN>>

== END OF INFO OUTPUT==

转载于“https://blog.csdn.net/tercel_zhang/article/details/52944286”

墨墨无文
关注
专栏目录
TI DSP TMS320C66x学习笔记之DSP优化经验（八-1）
Mars_WH的专栏
05-03 1万+
为了实现循环的pipeline。我们需要把循环体进行拆分，拆分时要考虑以下几点： （1）、拆分成几个循环比较合适？在各个循环能pipeline的前提下，拆开的循环个数越少越好。这就要求尽可能让各个循环的运算量接近。 （2）考虑在什么地方把程序拆开比较合适？循环体里的数据流往往并不是单一的，在拆开的断点处势必要用中间变量保存上次的循环运算结果，供以后的循环用。适当的拆开循环体，使所需的中间变量越少越好。 （3）循环体中的函数调用必须定义成内嵌形式，含有函数调用的循环系统是无法使之pipeline的
DSP程序优化方法
左超
09-14 3629
1、选择合适的算法和数据结构     选择一种合适的数据结构很重要，如果在一堆随机存放的数中使用了大量的插入和删除指令，那使用链表要快得多。数组与指针语句具有十分密切的关系，一般来说，指针比较灵活简洁，而数组则比较直观，容易理解。对于大部分的编译器，使用指针比使用数组生成的代码更短，执行效率更高。     在许多种情况下，可以用指针运算代替数组索引，这样做常常能产生又快又短的代码。与数组索引相比，
dsp中c代码优化
04-02
dsp中c代码优化方法，具有一定参考价值
DSP代码优化方法
07-14
DSP代码优化方法
dsp 优化心得
06-15 2852
1.多重if else 嵌套需要用函数指针来优化！2.ldbu  .D   *A_s++[A_lx],A_sreg1 //是先取A_s的值，然后在加偏移A_lx四、       1、源代码：void fir_fxd1(short input[], short coefs[], short out[]){    int i, j;    for (i = 
DSP的代码优化求余运算，用移位实现乘除法运算，延时函数
07-31
位操作只需一个指令周期即可完成，而大部分的 C 编译器的“ % ”运算均是调用子程序来完成，代码长、执行速度慢。
TMS320C66x汇编语言的DSP代码优化技术
07-26
随着主频的不断增加及多核处理器的出现，DSP运算量也不断地提高，但是从系统集成度的提高和成本角度考虑，仍然需要对C语言代码进行优化以提高效率。介绍了TMS320C66x汇编和流水线技术，阐述了在CCS编译器的实时运行...
TI DSP代码优化指南
01-16
### TI DSP代码优化指南 #### 一、概述 在TI（德州仪器）CCS（Code Composer Studio）开发环境中，对DSP（数字信号处理器）程序进行优化是提高性能的关键步骤之一。本文将详细介绍如何利用CCS中的Compiler ...
DSP程序代码优化方法，初学者不要错过，值得开发人员参考
12-03
DSP开发者通常面临着C语言代码的优化问题，此文章详细介绍了DSP代码的优化措施。
DSP_TMS320F28377D_算法加速方法4_代码本身的优化
heqiunong的专栏
09-04 1517
前面3篇的优化思路是从硬件本身和函数库这些方向去加速， 本文则仅从去考虑加速的方法，下面提供7种代码本身的优化思路，并验证其是否有效。
DSP C66 代码优化基本介绍
xrh003的博客
05-08 2663
C66 代码优化基本介绍C6000系列处理器的闪光之处就是它可以通过循环提高运行速度。这在以循环为中心的数字信号处理、图像处理和其他数学程序方面有着非常明显的优势。“软件流水”的技术对提高循环代码的性能做出的贡献最大。软件流水只有在使用-o2 或 -o3 编译选项时，才会被启用。如下图示，如果不使用软件流水，循环就会在循环体 i 完成后再开始循环体 i+1。软件流水技术允许循环体出现重叠。因此，只...
6678DSP多核处理和优化
最新发布
langzi113的专栏
02-07 1401
本文旨在纪录6678多核处理的调试过程。
DSP程序结构优化技巧
怒放的生命的博客
01-05 2889
       DSP的内部结构具有特殊性，如果程序的编写能够在一定程度上与DSP的结构相适应，利用C编译器的优化功能就可以在程序级的基础上大大提高算法的性能[1～5]。具体如下：        (1) 尽量使用int 类型中间变量        在图像处理程序中，数据通常为8bit的变量，而C6XDSP的内部寄存器和数据通道为32bit，在图像处理的滤波、卷积等过程中，需要较多的中间变量，如果...
DSP优化
audio_algorithm的博客
05-20 1836
转自：算法的嵌入式移植（二）DSP优化 – Deep Studio 前面一个篇章：算法的嵌入式移植（一）C代码优化部分 侧重语言层面的优化，而更多针对DSP特性进行的优化将在本篇章进行介绍 ，内容大多是摘抄书本和网络资源，因为是整理以往的笔记，来源已经不可靠了，不过不影响其中内容的经典。 一． 前序 1.图像算法在嵌入式移植时(主要针对DSP芯片)优化的原则和步骤？ 原则：算法效果达到预期之前最好不要做过多的优化 步骤：windows下的算法级优化—>C语言的优化—>DSP下的C编译器
ccs代码优化
Focustc
08-25 8586
进行代码优化，先要找出程序的瓶颈，即占用CPU时间较多的代码，然后对其进行有针对性的优化。使用CCS提供的代码剖析工具Profile可以统计显示出程序中各个重要段和函数的运行时间，找出运算量较大的程序段，优化这些程序段，对于提高算法的性能有巨大影响。 　　⑴ 联合使用-pm和-03编译选项，对代码进行项目级的优化：CCS提供了强大功能的编译选项，从-O0到-O3共四级优化。-O3编译选项使能软件
TMS320C6000 DSP代码优化与性能提升技术
“TMS320C6000代码优化技术主要探讨了在TMS320C6000 DSP平台上如何进行高效的代码优化，以提高指令并行执行度和代码执行效率，充分利用系统资源。文章提到了C6x代码的三种形式——C代码、线性汇编代码和汇编代码，并...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值