Vdsp(bf561)中的浮点运算（16）：fract2x16类型

快乐虾

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

[e-mail]

本文适用于

ADSP-BF561

Visual DSP++ 5.0 (update 6)

Vdsp dual processor simulate

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

由于BF561内部带有两个16位的MAC，因此它将可以在一个周期内进行两个fract16类型的运算。

为适应这种特性，vdsp引入了一个称之为fract2x16的类型。它其实是定义为一个int类型的整数，但是其实际意义却是要用高低16位分别来表示两个fract16类型。

typedef int _raw32;

typedef _raw32 raw2x16;

typedef raw2x16 fract2x16;

要查看fract2x16类型的值还是只能通过data register窗口，手动将类型改为fract16，这样在寄存器的高16位和低16位就能分别看见这两个值了。

使用compose_fr2x16函数可以构造一个fr2x16数据：

The notation used to represent two fract16 values packed into a fract2x16 is {a,b}, where “a” is the fract16 packed into the high half, and “b” is the fract16 packed into the low half.

fract2x16 compose_fr2x16(fract16 f1, fract16 f2)

Takes two fract16 values, and returns a fract2x16 value.

直接看看它在头文件的定义：

/* Takes two fract16 values, and returns a fract2x16 value.

* Input: two fract16 values

* Returns: {_x,_y} */

#pragma inline

#pragma always_inline

static fract2x16 compose_fr2x16(fract16 _x, fract16 _y) {

return compose_2x16(_x,_y);

}

/* Composes a packed integer from two short inputs.

*/

#pragma inline

#pragma always_inline

static int compose_2x16(short __a, short __b) {

int __rval = __builtin_compose_2x16(__a, __b);

return __rval;

}

所以咱直接使用__builtin_compose_2x16得了。

result = __builtin_compose_2x16(0x7fff, 0x8000);

这行语句生成汇编就成了：

R0 = -32768 /* 2147450880 */;

R0.H = 32767 /* 2147450880 */;

[FP + 16] = R0;

从中可以很清楚看出组合的过程。

要从fr2x16得到高半部分可以用这个调用：

/* Takes a fract2x16 and returns the 'high half' fract16.

* Input: _x{a,b}

* Returns: a.

*/

#pragma inline

#pragma always_inline

static fract16 high_of_fr2x16(fract2x16 _x) {

return high_of_2x16(_x);

}

#pragma inline

#pragma always_inline

static _raw16 high_of_2x16(raw2x16 _x) {

return __builtin_extract_hi(_x);

}

取低半部分可以用下面的函数：

/* Takes a fract2x16 and returns the 'low half' fract16

* Input: _x{a,b}

* Returns: b

*/

#pragma inline

#pragma always_inline

static fract16 low_of_fr2x16(fract2x16 _x) {

return low_of_2x16(_x);

}

#pragma inline

#pragma always_inline

static _raw16 low_of_2x16(raw2x16 _x) {

return __builtin_extract_lo(_x);

}

1 参考资料

Vdsp(bf561)中的浮点运算（15）：vdsp库的一个BUG(2009-8-19)

Vdsp(bf561)中的浮点运算（14）：fract16除法(2009-8-18)

Vdsp(bf561)中的浮点运算（13）：fract16乘法运算(2009-8-18)

Vdsp(bf561)中的浮点运算（12）：fract16加减运算(2009-8-17)

Vdsp(bf561)中的浮点运算（11）：fract16与float的转换(2009-8-17)

Vdsp(bf561)中的浮点运算（10）：fract16类型表示(2009-8-17)

Vdsp(bf561)中的浮点运算（9）：long double和float的比较(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)