上一篇文章中我们熟悉了浮点数在内存中的存储方式,这篇文章中我们来分析一下数学库中提供的函数。数学库中包含了很多函数,比如取整函数、三角函数、指数对数函数、贝塞尔函数、伽马函数等等,很多函数是对具体一个算法或数学问题的描述。我们从中挑选了一个简单的函数,目的是分析清楚数学库中浮点数的处理过程,同时又不要陷入数学问题当中。我们挑选的函数是truncf(),这个函数声明如下:
float truncf(float x);
这是一个取整函数,向绝对值小的方向取整,比如truncf(3.2) = 3.0,truncf(3.7) = 3.0。
我们分析的代码取自glibc,在分析truncf()之前,我们先看几个数据结构和宏:
typedef unsigned int u_int32_t;
typedef union
{
float value;
u_int32_t word;
} ieee_float_shape_type;
# define GET_FLOAT_WORD(i,d) \
do { \
ieee_float_shape_type gf_u; \
gf_u.value = (d);