Imported Upstream version 93u+upstream

1 files changed, 96 insertions, 0 deletions

diff --git a/src/lib/libast/sfio/_sfputd.c b/src/lib/libast/sfio/_sfputd.c
new file mode 100644
index 0000000..733caa6
--- /dev/null
+++ b/src/lib/libast/sfio/_sfputd.c
@@ -0,0 +1,96 @@
+/***********************************************************************
+*                                                                      *
+*               This software is part of the ast package               *
+*          Copyright (c) 1985-2011 AT&T Intellectual Property          *
+*                      and is licensed under the                       *
+*                 Eclipse Public License, Version 1.0                  *
+*                    by AT&T Intellectual Property                     *
+*                                                                      *
+*                A copy of the License is available at                 *
+*          http://www.eclipse.org/org/documents/epl-v10.html           *
+*         (with md5 checksum b35adb5213ca9657e911e9befb180842)         *
+*                                                                      *
+*              Information and Software Systems Research               *
+*                            AT&T Research                             *
+*                           Florham Park NJ                            *
+*                                                                      *
+*                 Glenn Fowler <gsf@research.att.com>                  *
+*                  David Korn <dgk@research.att.com>                   *
+*                   Phong Vo <kpv@research.att.com>                    *
+*                                                                      *
+***********************************************************************/
+#include	"sfhdr.h"
+
+/*	Write out a floating point value in a portable format
+**
+**	Written by Kiem-Phong Vo.
+*/
+
+#if __STD_C
+int _sfputd(Sfio_t* f, Sfdouble_t v)
+#else
+int _sfputd(f,v)
+Sfio_t*		f;
+Sfdouble_t	v;
+#endif
+{
+#define N_ARRAY		(16*sizeof(Sfdouble_t))
+	reg ssize_t	n, w;
+	reg uchar	*s, *ends;
+	int		exp;
+	uchar		c[N_ARRAY];
+	Sfdouble_t	x;
+	SFMTXDECL(f);
+
+	SFMTXENTER(f,-1);
+
+	if(f->mode != SF_WRITE && _sfmode(f,SF_WRITE,0) < 0)
+		SFMTXRETURN(f, -1);
+	SFLOCK(f,0);
+
+	/* get the sign of v */
+	if(v < 0.)
+	{	v = -v;
+		n = 1;
+	}
+	else	n = 0;
+
+	/* make the magnitude of v < 1 */
+	if(v != 0.)
+		v = frexpl(v,&exp);
+	else	exp = 0;
+
+	/* code the sign of v and exp */
+	if((w = exp) < 0)
+	{	n |= 02;
+		w = -w;
+	}
+
+	/* write out the signs and the exp */
+	SFOPEN(f,0);
+	if(sfputc(f,n) < 0 || (w = sfputu(f,w)) < 0)
+		SFMTXRETURN(f, -1);
+	SFLOCK(f,0);
+	w += 1;
+
+	s = (ends = &c[0])+sizeof(c);
+	while(s > ends)
+	{	/* get 2^SF_PRECIS precision at a time */
+		n = (int)(x = ldexpl(v,SF_PRECIS));
+		*--s = n|SF_MORE;
+		v = x-n;
+		if(v <= 0.)
+			break;
+	}
+
+	/* last byte is not SF_MORE */
+	ends = &c[0] + sizeof(c) -1;
+	*ends &= ~SF_MORE;
+
+	/* write out coded bytes */
+	n = ends - s + 1;
+	w = SFWRITE(f,(Void_t*)s,n) == n ? w+n : -1;
+
+	SFOPEN(f,0);
+	SFMTXRETURN(f,w);
+}