3
/* Return the length of a double value if coded in a portable format
5
** Written by Kiem-Phong Vo
9
int _sfdlen(Sfdouble_t v)
15
#define N_ARRAY (16*sizeof(Sfdouble_t))
23
/* make the magnitude of v < 1 */
28
for(w = 1; w <= N_ARRAY; ++w)
29
{ /* get 2^SF_PRECIS precision at a time */
30
n = (int)(x = ldexp(v,SF_PRECIS));
36
return 1 + sfulen(exp) + w;