scalbn, scalbnf, scalbnl, scalbln, scalblnf, scalblnl

参数

返回值

若未发生错误，则返回 arg 乘以 FLT_RADIX 的 exp 次幂（ arg×FLT_RADIX exp
）。

如果发生因溢出导致的范围错误，将返回 ± HUGE_VAL 、 ±HUGE_VALF 或 ±HUGE_VALL 。

如果发生因下溢导致的范围错误，将返回正确结果（经舍入后）。

错误处理

错误报告方式遵循 math_errhandling 中的规范。

如果实现支持 IEEE 浮点算术 (IEC 60559)，

• 除非发生范围错误，否则永远不会引发 FE_INEXACT （结果是精确的）。
• 除非发生范围错误，否则将忽略 当前舍入模式 。
• 如果 arg 为 ±0，则直接返回该值（不作修改）。
• 如果 arg 为 ±∞，则直接返回该值（不作修改）。
• 如果 exp 为 0，则直接返回 arg （不作修改）。
• 如果 arg 为 NaN，则返回 NaN。

注释

在二进制系统（其中 FLT_RADIX 为 2 ）中， scalbn 等价于 ldexp 。

尽管 scalbn 和 scalbln 被指定为高效执行该操作，但在许多实现中，它们的效率仍低于通过算术运算符进行二的幂次乘除运算。

提供 scalbln 函数的原因是：从最小正浮点数值缩放到最大有限值所需的缩放因子可能超过 32767 ——即标准所保证的 INT_MAX 。特别地，对于80位 long double 类型，该缩放因子为 32828 。

示例

#include <errno.h>
#include <fenv.h>
#include <float.h>
#include <math.h>
#include <stdio.h>
// #pragma STDC FENV_ACCESS ON
int main(void)
{
    printf("scalbn(7, -4) = %f\n", scalbn(7, -4));
    printf("scalbn(1, -1074) = %g (minimum positive subnormal double)\n",
            scalbn(1, -1074));
    printf("scalbn(nextafter(1,0), 1024) = %g (largest finite double)\n",
            scalbn(nextafter(1,0), 1024));
    // 特殊值
    printf("scalbn(-0, 10) = %f\n", scalbn(-0.0, 10));
    printf("scalbn(-Inf, -1) = %f\n", scalbn(-INFINITY, -1));
    // 错误处理
    errno = 0; feclearexcept(FE_ALL_EXCEPT);
    printf("scalbn(1, 1024) = %f\n", scalbn(1, 1024));
    if (errno == ERANGE)
        perror("    errno == ERANGE");
    if (fetestexcept(FE_OVERFLOW))
        puts("    FE_OVERFLOW raised");
}

scalbn(7, -4) = 0.437500
scalbn(1, -1074) = 4.94066e-324 (minimum positive subnormal double)
scalbn(nextafter(1,0), 1024) = 1.79769e+308 (largest finite double)
scalbn(-0, 10) = -0.000000
scalbn(-Inf, -1) = -inf
scalbn(1, 1024) = inf
    errno == ERANGE: Numerical result out of range
    FE_OVERFLOW raised

参考文献

参见