1da2e3ebdSchin /***********************************************************************
2da2e3ebdSchin * *
3da2e3ebdSchin * This software is part of the ast package *
4*3e14f97fSRoger A. Faulkner * Copyright (c) 1985-2010 AT&T Intellectual Property *
5da2e3ebdSchin * and is licensed under the *
6da2e3ebdSchin * Common Public License, Version 1.0 *
77c2fbfb3SApril Chin * by AT&T Intellectual Property *
8da2e3ebdSchin * *
9da2e3ebdSchin * A copy of the License is available at *
10da2e3ebdSchin * http://www.opensource.org/licenses/cpl1.0.txt *
11da2e3ebdSchin * (with md5 checksum 059e8cd6165cb4c31e351f2b69388fd9) *
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;
397c2fbfb3SApril Chin SFMTXDECL(f);
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