DSP

优化饱和加法运算

2019-07-13 16:46发布

通常计算机上的加法运算是循环的,例如 0x0001 + 0x7fff = 0x8000, 其中0x7fffshort类型的最大正值,0x8000是最小负值。假设我们是把两个声音相加,由于结果的符号反转,将出现不连续,如果我们把上面的数当做Q15定点小数来看的话,本来是正1(0x7fff),由于加了个很小的数(0x0001),结果却变为了负1(0x8000)。 因此DSP芯片都支持饱和加法,当结果超出范围的时候,就取范围的上下限为结果:0x0001 + 0x7fff = 0x7fff。例如在5510芯片中,为了使用饱和加法只需要把状态位:SATASATD设为1即可(SATASATD的具体含义请参考手册)。这种位操作只能在汇编下进行,因此为了方便C语言编程,编译器提供了伪函数_sadd_ssub。使用这些伪函数的时候,编译器会自动添加设置和清除饱和位的语句,例如: a=_sadd(a,b)可能会编译成 BSET ST3_SATA ;设置饱和位
ADD T0, T1
BCLR ST3_SATA ;
清除饱和位 当几个_sadd连续使用的时候,编译器只设置和清除一次饱和位。但是当几个_sadd调用中间还有一般的加法的话,就会设置和清除好几次饱和位了。 例如 c= _sadd(c,b); // 设置并清除一次
c= a+c;
c= _sadd(c,b); //
再设置并清除一次 如果在一个循环中有多次这样的交叉调用的话,饱和位将被设置和清除很多次,这样显然是浪费DSP资源的。因此尽量把饱和加减的调用集中,为了减少设置清除饱和位的次数,一些不需要饱和运算的加减也可以使用饱和运算。如果可以把整个函数都改写为饱和加减运算的话,那么干脆在函数最前面添加asm(" BSET ST3_SATA");手工设置饱和位,这样程序中可以直接是用加减号运算。