DSP代码优化

最新推荐文章于 2022-10-06 11:11:36 发布
最新推荐文章于 2022-10-06 11:11:36 发布
阅读量1.1k 收藏 6
点赞数 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”

墨墨无文
关注
  • 2
    点赞
  • 踩
  • 6
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
TMS320C66x汇编语言的DSP代码优化技术
07-26
随着主频的不断增加及多核处理器的出现，DSP运算量也不断地提高，但是从系统集成度的提高和成本角度考虑，仍然需要对C语言代码进行优化以提高效率。介绍了TMS320C66x汇编和流水线技术，阐述了在CCS编译器的实时运行环境下，编程需要注意的有关问题，以及如何利用流水线技术提高DSP的处理效率。
DSP C6000代码优化
LWN7321的博客
04-22 687
忙完项目的事情，对DSP的学习又近了一步。 在此介绍一下C6000代码优化的问题： 1.手动汇编优化 软件流水 1.1.更短的执行时间 1.2.更小的代码空间 两部分需要权衡 线性汇编比汇编更加简单 编译器把高级语言编译成汇编，汇编器是吧汇编编译成.obj的二进制代码，连接器把所有文件连接到一起生成可执行文件 1.31.先指令 后-O（n）选项，还不行就线性汇编，很少到手动汇编 2. 项目经验 首先得评估 代码运行时间，这个在DSP DEBUG -TOOL -clock，利用该工具可以很好的评估到
dsp中c代码优化
04-02
dsp中c代码优化方法，具有一定参考价值
DSP代码优化方法
07-14
DSP代码优化方法
论文研究-基于TMS320C6000系列DSP代码优化方法研究 .pdf
08-15
基于TMS320C6000系列DSP代码优化方法研究，邱春武，余文，随着DSP任务系统复杂度的增大，为了充分利用TMS320C6000系列DSP丰富的内部资源，需要对程序进行优化。程序的优化一般分三个阶段进行，�
restrict #pragma MUST_ITERATE( , , ) _nassert
05-08 8166
restrict的作用就是限制一个指针对一块内存的访问，进一步说就是如果一块内存区域通过一个受限制指针访问，那么它就不能通过另一个受限指针访问。 引入restrict的目的是确保同一块内存上没有其它引用，让编译器更好地优化指令，生成更有效的汇编代码。 把#pragma MUST_ITERATE( , , )放在循环体之前，告知开发板循环次数，改善软件流水 MUST_ITERATE
C6000系列DSP优化总结
是哥哥呀的专栏
03-09 1176
1.对指针使用restrict来限定memory访问,可以减少pipe 2.improve the balance of resources提高资源平衡,A-side和B-side,因为C6000有两个.D单元和两个.T地址路径，在任何给定的循环（A侧和B侧）上都可用，所以编译器必须在一侧（A侧）对至少两个操作进行分区,使用#pragma MUST_ITERATE(min, max, multiple);来强制迭代 注意事项: #pragma MUST_ITERATE(5) =#pragma MU..
DSP优化
audio_algorithm的博客
05-20 1577
转自：算法的嵌入式移植（二）DSP优化 – Deep Studio 前面一个篇章：算法的嵌入式移植（一）C代码优化部分 侧重语言层面的优化，而更多针对DSP特性进行的优化将在本篇章进行介绍 ，内容大多是摘抄书本和网络资源，因为是整理以往的笔记，来源已经不可靠了，不过不影响其中内容的经典。 一． 前序 1.图像算法在嵌入式移植时(主要针对DSP芯片)优化的原则和步骤？ 原则：算法效果达到预期之前最好不要做过多的优化 步骤：windows下的算法级优化—>C语言的优化—>DSP下的C编译器
C6000 DSP优化技术入门
小裘的博客
10-06 1559
DSP优化技术：restrict 指针；循环展开；地址对齐；程序级优化；线性汇编。
数据打包技术
寻找强烈无比的存在感！！！
04-02 1595
转载自http://dsalli0927.blog.163.com/blog/static/888076072008715535584/ C6000访问存储器是很费时的，要提高C6000的数据处理率，应该使1个Load/Store指令能够访问多个数据。当程序需要对一连传的短型数据进行操作时，可使用字（整型）一次访问2个短型数据；然后用C6000的相应指令，如同时进行2个16位的加法指令，
编译指示_#pragma在DSP中的编程技巧
dbl240906551的博客
02-03 1218
转载自：https://blog.csdn.net/ranyeliang/article/details/102454462 编译指示(Pragma Directives)可能是所有的预处理指令中最复杂的了，它的作用是设定编译器的状态或者是指示编译器完成一些特定的动作。#pragma指令对编译器给出了如何处理特定的函数、对象和代码段的方法,在保持与C/C++语言完全兼容的情况下,给出主机(比如C28x)或操作系统(比如DSP/BIOS)专有的特征。这些编译指示的使用较为复杂，但是我们还必须要了解它们，因为它
DSP编程技巧之：不得不看的编译指示_#pragma
u013749583的博客
10-19 5980
编译指示(Pragma Directives)可能是所有的预处理指令中最复杂的了，它的作用是设定编译器的状态或者是指示编译器完成一些特定的动作。#pragma指令对编译器给出了如何处理特定的函数、对象和代码段的方法,在保持与C/C++语言完全兼容的情况下,给出主机(比如C28x)或操作系统(比如DSP/BIOS)专有的特征。这些编译指示的使用较为复杂，但是我们还必须要了解它们，因为它们是程序中必不
#pragma STRUCT_ALIGN And #pragma pack明辨
dinkysummydinky的专栏
03-13 993
Posted by Bernie Thompson TI  on Feb 03 2010 09:46 AM Verified Answer Verified by George Mock Mastermind41700 points For what you are trying to do it sounds like you
DSP移植和优化之“预编译伪指令: DATA_ALIGN和STRUCT_ALIGN
寻找强烈无比的存在感！！！
05-31 5532
(1)#pragma DATA_ALIGN 语法： #pragma DATA_ALIGN( symbol，constant） 作用:排列symbol到constant 指点的列边界，constant为2的n次方，n=0,1,2.. 说白了就是首地址的字节对齐，首地址字对齐的最大好处就是可以告诉编译器充分利用行或块拷贝。 假设一个CACHE 的line 是128 字节【
DSP编程优化总结
nuaahz的博客
04-03 3945
一、编程优化 1、双重循环、多重循环优化 （1）多重循环拆成单层循环，减少循环层数； 例如，双重循环内一个cycle只使用了一个乘法器，拆成单层循环后，一个cycle可使用2个乘法器，充分利用DSP乘法器资源，同时运算速度也会加快； （2）循环次数少的放在外层循环，循环次数多的放在内存循环； （3）二维数组的二重循环：二维数组的行循环放在外层循环，列循环放在内层循环； ...
DSP程序优化方法
yuchen的博客
07-25 4622
对程序进行优化前后，程序在DM6437平台上的运行效率能有明显的提升，常用的优化方法一般从三个方面展开。 1 对程序编译工具进行合理配置 单击project然后选择build options打开配置界面，如下图所示。 这些编译选项分为以下几类： Basic：所包含的具体选项有：-g，用于选择是否生成全面的调试诊断信息；-ms(n)，在进行减小代码尺寸与提升程序执行效率之间权衡；-o(n...
DSP程序结构优化技巧
怒放的生命的博客
01-05 2776
       DSP的内部结构具有特殊性，如果程序的编写能够在一定程度上与DSP的结构相适应，利用C编译器的优化功能就可以在程序级的基础上大大提高算法的性能[1～5]。具体如下：        (1) 尽量使用int 类型中间变量        在图像处理程序中，数据通常为8bit的变量，而C6XDSP的内部寄存器和数据通道为32bit，在图像处理的滤波、卷积等过程中，需要较多的中间变量，如果...
PI控制器DSP代码
最新发布
10-25
下面是一个简单的PI控制器DSP代码示例： ``` float kp = 0.5; // 比例参数 float ki = 0.1; // 积分参数 float error = 0; // 误差 float integral = 0; // 积分项 float output = 0; // 控制量 while(1) { float...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值