<b>Vdsp(bf561)中的浮点运算（13）:fract16乘法运算</b>[VC/C++编程]

vdsp供应了三种差别的方法举行fract16的乘法运算.

1.1 mult_fr1x16函数

这个函数定义为:

#pragma inline
#pragma always_inline
static fract16　 mult_fr1x16(fract16　 __a, fract16　 __b) {
fract16　 __rval = __builtin_mult_fr1x16(__a, __b);
return __rval;
}

从这里可以看出我们实际可以利用__builtin_mult_fr1x16这一函数调用.

写一个很简单的程序:

typedef fract16 ft;

ft calc(ft x, ft y)
{
ft r;
r = __builtin_mult_fr1x16(x, y);
return r;
}

__builtin_mult_fr1x16展开后的汇编代码为:

// line 29
R0.L = R0.L * R1.L (T);
R0 = R0.L (X);
W[FP + 16] = R0;

因而完成这样一个运算将只需求一个cycle的时间.

这里的乘法运算利用了(T)尾缀,文档里这样注释:

Signed fraction with truncation. Truncate Accumulator 9.31 format value at bit 16. (Perform no rounding.) Saturate the result to 1.15 precision in destination register half. Result is between minimum -1 and maximum 1-2-15 (or, expressed in hex, between minimum 0x8000 and maximum 0x7FFF).

这种计算方法直接将累加器里的数举行截断而不举行任何舍入的处理.

1.2 multr_fr1x16

这个函数定义为:

/* Performs a 16-bit fractional multiplication of the two input
** parameters. The result is rounded to 16 bits. Whether the
** rounding is biased or unbiased depends what the RND_MOD bit
** in the ASTAT register is set to.
*/

#pragma inline
#pragma always_inline
static fract16　 multr_fr1x16(fract16　 __a, fract16　 __b) {
fract16　 __rval = __builtin_multr_fr1x16(__a, __b);
return __rval;
}

　　以上是“<b>Vdsp(bf561)中的浮点运算（13）:fract16乘法运算</b>[VC/C++编程]”的内容，如果你对以上该文章内容感兴趣，你可以看看七道奇为您推荐以下文章：