14887Schin /***********************************************************************
24887Schin * *
34887Schin * This software is part of the ast package *
4*12068SRoger.Faulkner@Oracle.COM * Copyright (c) 1985-2010 AT&T Intellectual Property *
54887Schin * and is licensed under the *
64887Schin * Common Public License, Version 1.0 *
78462SApril.Chin@Sun.COM * by AT&T Intellectual Property *
84887Schin * *
94887Schin * A copy of the License is available at *
104887Schin * http://www.opensource.org/licenses/cpl1.0.txt *
114887Schin * (with md5 checksum 059e8cd6165cb4c31e351f2b69388fd9) *
124887Schin * *
134887Schin * Information and Software Systems Research *
144887Schin * AT&T Research *
154887Schin * Florham Park NJ *
164887Schin * *
174887Schin * Glenn Fowler <gsf@research.att.com> *
184887Schin * David Korn <dgk@research.att.com> *
194887Schin * Phong Vo <kpv@research.att.com> *
204887Schin * *
214887Schin ***********************************************************************/
224887Schin #include "sfhdr.h"
234887Schin
244887Schin /* Read a portably coded double value
254887Schin **
264887Schin ** Written by Kiem-Phong Vo
274887Schin */
284887Schin
294887Schin #if __STD_C
sfgetd(Sfio_t * f)304887Schin Sfdouble_t sfgetd(Sfio_t* f)
314887Schin #else
324887Schin Sfdouble_t sfgetd(f)
334887Schin Sfio_t* f;
344887Schin #endif
354887Schin {
364887Schin reg uchar *s, *ends, c;
374887Schin reg int p, sign, exp;
384887Schin Sfdouble_t v;
398462SApril.Chin@Sun.COM SFMTXDECL(f);
404887Schin
418462SApril.Chin@Sun.COM SFMTXENTER(f,-1.);
424887Schin
434887Schin if((sign = sfgetc(f)) < 0 || (exp = (int)sfgetu(f)) < 0)
444887Schin SFMTXRETURN(f, -1.);
454887Schin
464887Schin if(f->mode != SF_READ && _sfmode(f,SF_READ,0) < 0)
474887Schin SFMTXRETURN(f, -1.);
484887Schin
494887Schin SFLOCK(f,0);
504887Schin
514887Schin v = 0.;
524887Schin for(;;)
534887Schin { /* fast read for data */
544887Schin if(SFRPEEK(f,s,p) <= 0)
554887Schin { f->flags |= SF_ERROR;
564887Schin v = -1.;
574887Schin goto done;
584887Schin }
594887Schin
604887Schin for(ends = s+p; s < ends; )
614887Schin { c = *s++;
624887Schin v += SFUVALUE(c);
634887Schin v = ldexpl(v,-SF_PRECIS);
644887Schin if(!(c&SF_MORE))
654887Schin { f->next = s;
664887Schin goto done;
674887Schin }
684887Schin }
694887Schin f->next = s;
704887Schin }
714887Schin
724887Schin done:
734887Schin v = ldexpl(v,(sign&02) ? -exp : exp);
744887Schin if(sign&01)
754887Schin v = -v;
764887Schin
774887Schin SFOPEN(f,0);
784887Schin SFMTXRETURN(f, v);
794887Schin }
80