Vdsp(bf561)中的浮点运算（16）：fract2x16类型
快乐虾
http://blog.csdn.net/lights_joy/
本文适用于
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)