[资料分享] c6000系列的C代码优化(三)

Jacktang   2017-11-9 21:33 楼主
三、16位变为 32 位操作,使用 intrinsi函数,用 const等。   
        1、源代码:  
  • Word32 L_mpy_ll(Word32 L_var1, Word32 L_var2)
  • {
  •    double  aReg;
  •     Word32 lvar;
  •     /* (unsigned)low1 * (unsigned)low1 */
  •         aReg = (double)(0xffff & L_var1) * (double)(0xffff & L_var2) * 2.0;
  •       /* >> 16 */
  •         aReg = (aReg / 65536);
  •     aReg = floor(aReg);
  •     /* (unsigned)low1 * (signed)high2 */
  •         aReg += (double)(0xffff & L_var1) * ((double)L_shr(L_var2,16)) * 2.0 ;
  •     /* (unsigned)low2 * (signed)high1 */
  •         aReg += (double)(0xffff & L_var2) * ((double)L_shr(L_var1,16)) * 2.0 ;
  •       /* >> 16 */
  •         aReg = (aReg / 65536);
  •     aReg = floor(aReg);
  •       /* (signed)high1 * (signed)high2 */
  •         aReg += (double)(L_shr(L_var1,16)) * (double)(L_shr(L_var2,16)) * 2. 0;
  •     /* saturate result.. */
  •     lvar = L_saturate(aReg);
  •     return(lvar);
  • }

[color=rgb(51, 102, 153) !important]复制代码

       2、改编后的代码:  
  • static inline Word32 L_mpy_ll(Word32 L_var1, Word32 L_var2)
  • {
  •     Word32 aReg_hh;
  •     Word40 aReg,aReg_ll,aReg_lh,aReg_hl;
  •     aReg_ll = (Word40)_mpyu(L_var1, L_var2)>>16;
  •     aReg_lh = (Word40)_mpyluhs(L_var1, L_var2);
  •     aReg_hl = (Word40)_mpyhslu(L_var1, L_var2);
  •     aReg_hh = _smpyh(L_var1, L_var2);
  •     aReg    = _lsadd(aReg_ll, _lsadd(aReg_lh, aReg_hl));
  •     aReg    = _lsadd(aReg>>15, aReg_hh);
  •     return(_sat(aReg));
  • }

[color=rgb(51, 102, 153) !important]复制代码

       3、优化方法说明:  
        C6000 编译器提供的 intrinsic  可快速优化 C代码,intrinsic 用前下划线表示同调用函数一样可以调用它, 即直接内联为 C6000 的函数。 例如, 在上例的源代码中没有使用 intrinsics,每一行 C代码需多个指令周期,在改编后的代码中,每一行代码仅需一个指令周期。例如,“aReg_ll = (Word40)_mpyu(L_var1, L_var2)>>16”中“_mpyu”就是一个 intrinsics 函数,它表示两个无符号数的高 16 位相乘,结果返回。C6000 支持的所有 intrinsics指令及其功能参见《TMS320C6000 系列 DSP 的原理与应用》一书的第 265、266 页,该书还提供了另外的例子。这些内联函数定义在 CCS所在的C6000\CGTOOLS\Include 目录下的 C6X.h 文件中。下面这个例子是 C6000的 “Programmer's Guide” 上提取的使用 intrinsics优化 C代码的例子。
        源代码:  
  • int dotprod(const short *a, const short *b, unsigned int N)
  • {
  •         int i, sum = 0;
  •         for (i = 0; i < N; i++)
  •         sum += a * b;
  •     return sum;
  • }

[color=rgb(51, 102, 153) !important]复制代码

       改编后代码:  
  • int dotprod(const int *a, const int *b, unsigned int N)
  • {
  •         int i, sum1 = 0, sum2 = 0;
  •         for (i = 0; i < (N >> 1); i++)
  •     {
  •         sum1 += _mpy (a, b);
  •         sum2 += _mpyh(a, b);
  •     }
  •      return sum1 + sum2;
  • }

[color=rgb(51, 102, 153) !important]复制代码


回复评论

暂无评论,赶紧抢沙发吧
电子工程世界版权所有 京B2-20211791 京ICP备10001474号-1 京公网安备 11010802033920号
    写回复