Vdsp(bf561)中的浮点运算(10):fract16类型表示

快乐虾

http://blog.csdn.net/lights_joy/

lights@hb165.com  

本文适用于

ADSP-BF561

Visual DSP++ 5.0 (update 6)

Vdsp dual processor simulate

  

欢迎转载,但请保留作者信息

 

 

 

fract16VDSP的另一种浮点类型,它使用1.15的格式:

 

fract16 data type

 

最高位为符号位,其余为数据位,没有指数和尾数的概念。

最大值:+0.999969482421875

最小值:-1

分辨率:2-15 = 0.000030517578125

blackfin里面使用fract类型进行计算可以获得最好的性能,但是要求在使用之前对需要进行计算的数据进行归一化,而且必须保证在计算过程中不会产生溢出。

在对数据归一化的时候,通常是需要进行浮点计算的,因此要产生很多额外的开销,如果使用fract进行的计算量小于归一化的计算量,那就没有必要使用它。比如vdsp文档里面提及的那个例子:

/* return a*b calculated using fract implementation */

float mul_fp(float a, float b) {

   int sign_a, sign_b, sign_res;

   float scaled_a, scaled_b, fract_div_res, result;

   int exp_a, exp_b, exp_res;

   fract32 fract_a, fract_b, fract_res;

   fract32 fract_exp_a, fract_exp_b, fract_exp_res;

   fract16 fract_reshi, fract_reslo;

  

   /* if either input is 0, return 0. */

   if (a == 0.0 || b == 0.0)

      return 0.0;

 

   /* get sign and take absolute of inputs */

   if (*(unsigned int *)&a & 0x80000000) {

      sign_a=-1;

      a = fabs(a);

   } else

      sign_a=1;

 

   if (*(unsigned int *)&b & 0x80000000) {

      sign_b=-1;

      b = fabs(b);

   } else

      sign_b=1;

 

   /* compute sign of result */

   sign_res = sign_a * sign_b;

 

   /* scale inputs */

   scaled_a = frexpf(a, &exp_a);

   scaled_b = frexpf(b, &exp_b);

 

   /* convert scaled inputs to fract */

   fract_a = float_to_fr32(scaled_a);

   fract_b = float_to_fr32(scaled_b);

 

   /* compose extended precision ETSI inputs */

   fract_exp_a = L_Comp(extract_h(fract_a),

                        extract_l(fract_a));

   fract_exp_b = L_Comp(extract_h(fract_b),

                        extract_l(fract_b));

 

   /* do fractional multiplication in extended precision */

   fract_res = Mpy_32(extract_h(fract_exp_a),

                      extract_l(fract_exp_a),

                      extract_h(fract_exp_b),

                      extract_l(fract_exp_b));

 

   /* multiply exponents by adding */

   exp_res = exp_a + exp_b;

 

   /* convert mul result back to float */

   fract_div_res = fr32_to_float(fract_res);

 

   /* compose the floating-point result */

   result = ldexpf(fract_div_res, exp_res);

 

   /* negate result if necessary */

   result = result * sign_res;

   /* return result */

   return result;

}  /* mul_fp */

这个例子用于完成两个浮点数的乘法,它的乘法计算使用fract类型来完成,但是它需要将float类型转换为fract类型,计算完成后再把fract类型转换为float类型,整个过程需要743cycle,远远大于直接使用float进行计算的92cycle,想必没有人会愿意使用这样的方法吧。

使用这种类型的数据,由于它不使用指数位,自然也没有naninf这样的概念,使用时更加需要小心!

vdsp的编译器里面,fract16不是作为一个内置类型出现的,而是使用了一个typedef来完成:

typedef  short    fract16;

由此造成的问题是在vdsp下无法直接查看fract16类型变量的值,非得通过寄存器窗口手动将数据类型设置为fract16,很是有点别扭!

注意fract16的表示,它的最高位表示  -20,也就是说,如果这一位为1,那么整个数将在-1的基础上加上后面的那一串值,这与float的最高位单纯的表示符号位不同。因此在float类型中,只要简单改变最高位的状态即可改变浮点数的符号,但是对于fract这却是非常错误的做法。

比如对于0.2,以float表示时的十六进制值为0x3e4c cccd,而-0.2的值以十六进制表示则为0xbe4c cccd,即它们只差一个最高位。

0.2如果表示为fract16的值为0x1999,但是-0.2的值却是0xe666,它们的值相加刚好是0xffff

使用negate_fr1x16可以改变一个fract16数值的符号位:

/* Returns the 16-bit result of the negation of the input parameter (-_x).  If

 * the input is 0x8000, saturation occurs and 0x7fff is returned. */

#pragma inline

#pragma always_inline

static fract16 negate_fr1x16(fract16 _x)

  { return (fract16)__builtin_negate_fr2x16((int)_x); }

展开__builtin_negate_fr2x16可以发现它其它就是这样的操作:

       R1 = - R0 (V);

       R1 = R1.L (X);

直接使用整数的取反。

 

 

1       参考资料

Vdsp(bf561)中的浮点运算(9):long doublefloat的比较(2009-8-14)

Vdsp(bf561)中的浮点运算(8):float除法运算(2009-8-14)

Vdsp(bf561)中的浮点运算(7):float乘法运算(2009-8-13)

Vdsp(bf561)中的浮点运算(6):float加减运算(2009-8-13)

Vdsp(bf561)中的浮点运算(5):float类型表示总结(2009-8-12)

Vdsp(bf561)中的浮点运算(4):FLT_MAX(2009-8-12)

Vdsp(bf561)中的浮点运算(3):FLT_MIN(2008-12-19)

Vdsp(bf561)中的浮点运算(2):float的疑问(2008-12-18)

Vdsp(bf561)中的浮点运算(1):文档的说法(2008-12-16)

 

 

 

 

 

 

 

 

 

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

嵌云阁主

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值