1/*********************************************************************** 2* * 3* This software is part of the ast package * 4* Copyright (c) 1985-2011 AT&T Intellectual Property * 5* and is licensed under the * 6* Common Public License, Version 1.0 * 7* by AT&T Intellectual Property * 8* * 9* A copy of the License is available at * 10* http://www.opensource.org/licenses/cpl1.0.txt * 11* (with md5 checksum 059e8cd6165cb4c31e351f2b69388fd9) * 12* * 13* Information and Software Systems Research * 14* AT&T Research * 15* Florham Park NJ * 16* * 17* Glenn Fowler <gsf@research.att.com> * 18* David Korn <dgk@research.att.com> * 19* Phong Vo <kpv@research.att.com> * 20* * 21***********************************************************************/ 22#include "sfhdr.h" 23 24/* Return the length of a double value if coded in a portable format 25** 26** Written by Kiem-Phong Vo 27*/ 28 29#if __STD_C 30int _sfdlen(Sfdouble_t v) 31#else 32int _sfdlen(v) 33Sfdouble_t v; 34#endif 35{ 36#define N_ARRAY (16*sizeof(Sfdouble_t)) 37 reg int n, w; 38 Sfdouble_t x; 39 int exp; 40 41 if(v < 0) 42 v = -v; 43 44 /* make the magnitude of v < 1 */ 45 if(v != 0.) 46 v = frexpl(v,&exp); 47 else exp = 0; 48 49 for(w = 1; w <= N_ARRAY; ++w) 50 { /* get 2^SF_PRECIS precision at a time */ 51 n = (int)(x = ldexpl(v,SF_PRECIS)); 52 v = x-n; 53 if(v <= 0.) 54 break; 55 } 56 57 return 1 + sfulen(exp) + w; 58} 59