1da2e3ebdSchin /***********************************************************************
2da2e3ebdSchin *                                                                      *
3da2e3ebdSchin *               This software is part of the ast package               *
4*b30d1939SAndy Fiddaman *          Copyright (c) 1985-2011 AT&T Intellectual Property          *
5da2e3ebdSchin *                      and is licensed under the                       *
6*b30d1939SAndy Fiddaman *                 Eclipse Public License, Version 1.0                  *
77c2fbfb3SApril Chin *                    by AT&T Intellectual Property                     *
8da2e3ebdSchin *                                                                      *
9da2e3ebdSchin *                A copy of the License is available at                 *
10*b30d1939SAndy Fiddaman *          http://www.eclipse.org/org/documents/epl-v10.html           *
11*b30d1939SAndy Fiddaman *         (with md5 checksum b35adb5213ca9657e911e9befb180842)         *
12da2e3ebdSchin *                                                                      *
13da2e3ebdSchin *              Information and Software Systems Research               *
14da2e3ebdSchin *                            AT&T Research                             *
15da2e3ebdSchin *                           Florham Park NJ                            *
16da2e3ebdSchin *                                                                      *
17da2e3ebdSchin *                 Glenn Fowler <gsf@research.att.com>                  *
18da2e3ebdSchin *                  David Korn <dgk@research.att.com>                   *
19da2e3ebdSchin *                   Phong Vo <kpv@research.att.com>                    *
20da2e3ebdSchin *                                                                      *
21da2e3ebdSchin ***********************************************************************/
22da2e3ebdSchin #include	"sfhdr.h"
23da2e3ebdSchin 
24da2e3ebdSchin /*	Read a portably coded double value
25da2e3ebdSchin **
26da2e3ebdSchin **	Written by Kiem-Phong Vo
27da2e3ebdSchin */
28da2e3ebdSchin 
29da2e3ebdSchin #if __STD_C
sfgetd(Sfio_t * f)30da2e3ebdSchin Sfdouble_t sfgetd(Sfio_t* f)
31da2e3ebdSchin #else
32da2e3ebdSchin Sfdouble_t sfgetd(f)
33da2e3ebdSchin Sfio_t*	f;
34da2e3ebdSchin #endif
35da2e3ebdSchin {
36da2e3ebdSchin 	reg uchar	*s, *ends, c;
37da2e3ebdSchin 	reg int		p, sign, exp;
38da2e3ebdSchin 	Sfdouble_t	v;
39*b30d1939SAndy Fiddaman 	SFMTXDECL(f); /* declare a local stream variable for multithreading */
40da2e3ebdSchin 
417c2fbfb3SApril Chin 	SFMTXENTER(f,-1.);
42da2e3ebdSchin 
43da2e3ebdSchin 	if((sign = sfgetc(f)) < 0 || (exp = (int)sfgetu(f)) < 0)
44da2e3ebdSchin 		SFMTXRETURN(f, -1.);
45da2e3ebdSchin 
46da2e3ebdSchin 	if(f->mode != SF_READ && _sfmode(f,SF_READ,0) < 0)
47da2e3ebdSchin 		SFMTXRETURN(f, -1.);
48da2e3ebdSchin 
49da2e3ebdSchin 	SFLOCK(f,0);
50da2e3ebdSchin 
51da2e3ebdSchin 	v = 0.;
52da2e3ebdSchin 	for(;;)
53da2e3ebdSchin 	{	/* fast read for data */
54da2e3ebdSchin 		if(SFRPEEK(f,s,p) <= 0)
55da2e3ebdSchin 		{	f->flags |= SF_ERROR;
56da2e3ebdSchin 			v = -1.;
57da2e3ebdSchin 			goto done;
58da2e3ebdSchin 		}
59da2e3ebdSchin 
60da2e3ebdSchin 		for(ends = s+p; s < ends; )
61da2e3ebdSchin 		{	c = *s++;
62da2e3ebdSchin 			v += SFUVALUE(c);
63da2e3ebdSchin 			v = ldexpl(v,-SF_PRECIS);
64da2e3ebdSchin 			if(!(c&SF_MORE))
65da2e3ebdSchin 			{	f->next = s;
66da2e3ebdSchin 				goto done;
67da2e3ebdSchin 			}
68da2e3ebdSchin 		}
69da2e3ebdSchin 		f->next = s;
70da2e3ebdSchin 	}
71da2e3ebdSchin 
72da2e3ebdSchin done:
73da2e3ebdSchin 	v = ldexpl(v,(sign&02) ? -exp : exp);
74da2e3ebdSchin 	if(sign&01)
75da2e3ebdSchin 		v = -v;
76da2e3ebdSchin 
77da2e3ebdSchin 	SFOPEN(f,0);
78da2e3ebdSchin 	SFMTXRETURN(f, v);
79da2e3ebdSchin }
80