Vdsp(bf561)中的浮点运算（10）：fract16范例暗示
副标题#e#
fract16是VDSP的另一种浮点范例，它利用1.15的名目：
最高位为标记位，其余为数据位，没有指数和尾数的观念。
最大值：+0.999969482421875
最小值：-1
判别率：2-15 = 0.000030517578125
在blackfin内里利用fract范例举办计较可以得到最好的机能，可是要求在利用之前对需要举办计较的数据举办归一化，并且必需担保在计较进程中不会发生溢出。
在对数据归一化的时候，凡是是需要举办浮点计较的，因此要发生许多特另外开销，假如利用fract举办的计较劲小于归一化的计较劲，那就没有须要利用它。好比vdsp文档内里提及的谁人例子：
#p#副标题#e#
/* 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范例，整个进程需要743个cycle，远远大于直接利用float举办计较的92个cycle，想必没有人会愿意利用这样的要领吧。
利用这种范例的数据，由于它不利用指数位，自然也没有nan和inf这样的观念，利用时越发需要小心！
在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);
直接利用整数的取反。
