1N/A/***********************************************************************
1N/A* *
1N/A* This software is part of the ast package *
1N/A* Copyright (c) 1985-2011 AT&T Intellectual Property *
1N/A* and is licensed under the *
1N/A* Common Public License, Version 1.0 *
1N/A* by AT&T Intellectual Property *
1N/A* *
1N/A* A copy of the License is available at *
1N/A* http://www.opensource.org/licenses/cpl1.0.txt *
1N/A* (with md5 checksum 059e8cd6165cb4c31e351f2b69388fd9) *
1N/A* *
1N/A* Information and Software Systems Research *
1N/A* AT&T Research *
1N/A* Florham Park NJ *
1N/A* *
1N/A* Glenn Fowler <gsf@research.att.com> *
1N/A* David Korn <dgk@research.att.com> *
1N/A* Phong Vo <kpv@research.att.com> *
1N/A* *
1N/A***********************************************************************/
1N/A#include "sfhdr.h"
1N/A
1N/A/* Return the length of a double value if coded in a portable format
1N/A**
1N/A** Written by Kiem-Phong Vo
1N/A*/
1N/A
1N/A#if __STD_C
1N/Aint _sfdlen(Sfdouble_t v)
1N/A#else
1N/Aint _sfdlen(v)
1N/ASfdouble_t v;
1N/A#endif
1N/A{
1N/A#define N_ARRAY (16*sizeof(Sfdouble_t))
1N/A reg int n, w;
1N/A Sfdouble_t x;
1N/A int exp;
1N/A
1N/A if(v < 0)
1N/A v = -v;
1N/A
1N/A /* make the magnitude of v < 1 */
1N/A if(v != 0.)
1N/A v = frexpl(v,&exp);
1N/A else exp = 0;
1N/A
1N/A for(w = 1; w <= N_ARRAY; ++w)
1N/A { /* get 2^SF_PRECIS precision at a time */
1N/A n = (int)(x = ldexpl(v,SF_PRECIS));
1N/A v = x-n;
1N/A if(v <= 0.)
1N/A break;
1N/A }
1N/A
1N/A return 1 + sfulen(exp) + w;
1N/A}