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 /* Write out a floating point value in a portable format
25da2e3ebdSchin **
26da2e3ebdSchin ** Written by Kiem-Phong Vo.
27da2e3ebdSchin */
28da2e3ebdSchin
29da2e3ebdSchin #if __STD_C
_sfputd(Sfio_t * f,Sfdouble_t v)30da2e3ebdSchin int _sfputd(Sfio_t* f, Sfdouble_t v)
31da2e3ebdSchin #else
32da2e3ebdSchin int _sfputd(f,v)
33da2e3ebdSchin Sfio_t* f;
34da2e3ebdSchin Sfdouble_t v;
35da2e3ebdSchin #endif
36da2e3ebdSchin {
37da2e3ebdSchin #define N_ARRAY (16*sizeof(Sfdouble_t))
38da2e3ebdSchin reg ssize_t n, w;
39da2e3ebdSchin reg uchar *s, *ends;
40da2e3ebdSchin int exp;
41da2e3ebdSchin uchar c[N_ARRAY];
42da2e3ebdSchin Sfdouble_t x;
437c2fbfb3SApril Chin SFMTXDECL(f);
44da2e3ebdSchin
457c2fbfb3SApril Chin SFMTXENTER(f,-1);
46da2e3ebdSchin
47da2e3ebdSchin if(f->mode != SF_WRITE && _sfmode(f,SF_WRITE,0) < 0)
48da2e3ebdSchin SFMTXRETURN(f, -1);
49da2e3ebdSchin SFLOCK(f,0);
50da2e3ebdSchin
51da2e3ebdSchin /* get the sign of v */
52da2e3ebdSchin if(v < 0.)
53da2e3ebdSchin { v = -v;
54da2e3ebdSchin n = 1;
55da2e3ebdSchin }
56da2e3ebdSchin else n = 0;
57da2e3ebdSchin
58da2e3ebdSchin /* make the magnitude of v < 1 */
59da2e3ebdSchin if(v != 0.)
60da2e3ebdSchin v = frexpl(v,&exp);
61da2e3ebdSchin else exp = 0;
62da2e3ebdSchin
63da2e3ebdSchin /* code the sign of v and exp */
64da2e3ebdSchin if((w = exp) < 0)
65da2e3ebdSchin { n |= 02;
66da2e3ebdSchin w = -w;
67da2e3ebdSchin }
68da2e3ebdSchin
69da2e3ebdSchin /* write out the signs and the exp */
70da2e3ebdSchin SFOPEN(f,0);
71da2e3ebdSchin if(sfputc(f,n) < 0 || (w = sfputu(f,w)) < 0)
72da2e3ebdSchin SFMTXRETURN(f, -1);
73da2e3ebdSchin SFLOCK(f,0);
74da2e3ebdSchin w += 1;
75da2e3ebdSchin
76da2e3ebdSchin s = (ends = &c[0])+sizeof(c);
77da2e3ebdSchin while(s > ends)
78da2e3ebdSchin { /* get 2^SF_PRECIS precision at a time */
79da2e3ebdSchin n = (int)(x = ldexpl(v,SF_PRECIS));
80da2e3ebdSchin *--s = n|SF_MORE;
81da2e3ebdSchin v = x-n;
82da2e3ebdSchin if(v <= 0.)
83da2e3ebdSchin break;
84da2e3ebdSchin }
85da2e3ebdSchin
86da2e3ebdSchin /* last byte is not SF_MORE */
87da2e3ebdSchin ends = &c[0] + sizeof(c) -1;
88da2e3ebdSchin *ends &= ~SF_MORE;
89da2e3ebdSchin
90da2e3ebdSchin /* write out coded bytes */
91da2e3ebdSchin n = ends - s + 1;
92da2e3ebdSchin w = SFWRITE(f,(Void_t*)s,n) == n ? w+n : -1;
93da2e3ebdSchin
94da2e3ebdSchin SFOPEN(f,0);
95da2e3ebdSchin SFMTXRETURN(f,w);
96da2e3ebdSchin }
97