src/common/portio.c

#ifndef lint
static const char       RCSid[] = "$Id: portio.c,v 2.14 2006/12/23 17:27:45 greg Exp $";
#endif
/*
 * Portable i/o for binary files
 *
 * External symbols declared in rtio.h
 */

#include "copyright.h"

#include "rtio.h"

#include <math.h>

#ifdef getc_unlocked            /* avoid horrendous overhead of flockfile */
#undef getc
#undef putc
#define getc    getc_unlocked
#define putc    putc_unlocked
#endif


void
putstr(                         /* write null-terminated string to fp */
        char  *s,
        FILE  *fp
)
{
        do
                putc(*s, fp);
        while (*s++);
}


void
putint(                         /* write a siz-byte integer to fp */
        long  i,
        int  siz,
        FILE  *fp
)
{
        while (siz--)
                putc((int)(i>>(siz<<3) & 0xff), fp);
}


void
putflt(                         /* put out floating point number */
        double  f,
        FILE  *fp
)
{
        long  m;
        int  e;

        m = frexp(f, &e) * 0x7fffffff;
        if (e > 127) {                  /* overflow */
                m = m > 0 ? (long)0x7fffffff : -(long)0x7fffffff;
                e = 127;
        } else if (e < -128) {          /* underflow */
                m = 0;
                e = 0;
        }
        putint(m, 4, fp);
        putint((long)e, 1, fp);
}


char *
getstr(                         /* get null-terminated string */
        char  *s,
        FILE  *fp
)
{
        char  *cp;
        int  c;

        cp = s;
        while ((c = getc(fp)) != EOF)
                if ((*cp++ = c) == '\0')
                        return(s);
                
        return(NULL);
}


long
getint(                         /* get a siz-byte integer */
        int  siz,
        FILE  *fp
)
{
        int  c;
        long  r;

        if ((c = getc(fp)) == EOF)
                return(EOF);
        r = 0x80&c ? -1<<8|c : c;               /* sign extend */
        while (--siz > 0) {
                if ((c = getc(fp)) == EOF)
                        return(EOF);
                r <<= 8;
                r |= c;
        }
        return(r);
}


double
getflt(                         /* get a floating point number */
        FILE  *fp
)
{
        long    l;
        double  d;

        l = getint(4, fp);
        if (l == EOF && feof(fp))       /* EOF? */
                return((double)EOF);
        if (l == 0) {
                getc(fp);               /* exactly zero -- ignore exponent */
                return(0.0);
        }
        d = (l + (l > 0 ? .5 : -.5)) * (1./0x7fffffff);
        return(ldexp(d, (int)getint(1, fp)));
}
Revision:	2.15
Committed:	Fri Oct 19 01:40:33 2012 UTC (11 years, 6 months ago) by greg
Content type:	text/plain
Branch:	MAIN
CVS Tags:	rad4R2P2, rad5R0, rad4R2, rad4R2P1
Changes since 2.14:	+31 -23 lines
Log Message:	Added check for EOF in getflt() and completed ANSIficiation
#	User	Rev	Content
1	greg	2.1	#ifndef lint
2	greg	2.15	static const char RCSid[] = "$Id: portio.c,v 2.14 2006/12/23 17:27:45 greg Exp $";
3	greg	2.1	#endif
4			/*
5			* Portable i/o for binary files
6	greg	2.6	*
7	greg	2.14	* External symbols declared in rtio.h
8	greg	2.6	*/
9
10	greg	2.7	#include "copyright.h"
11	greg	2.1
12	schorsch	2.8	#include "rtio.h"
13	greg	2.1
14	greg	2.9	#include <math.h>
15	greg	2.10
16			#ifdef getc_unlocked /* avoid horrendous overhead of flockfile */
17	greg	2.11	#undef getc
18			#undef putc
19	greg	2.10	#define getc getc_unlocked
20			#define putc putc_unlocked
21			#endif
22	greg	2.1
23	greg	2.3
24	greg	2.6	void
25	greg	2.15	putstr( /* write null-terminated string to fp */
26			char *s,
27			FILE *fp
28			)
29	greg	2.1	{
30			do
31			putc(*s, fp);
32			while (*s++);
33			}
34
35
36	greg	2.6	void
37	greg	2.15	putint( /* write a siz-byte integer to fp */
38			long i,
39			int siz,
40			FILE *fp
41			)
42	greg	2.1	{
43			while (siz--)
44	greg	2.2	putc((int)(i>>(siz<<3) & 0xff), fp);
45	greg	2.1	}
46
47
48	greg	2.6	void
49	greg	2.15	putflt( /* put out floating point number */
50			double f,
51			FILE *fp
52			)
53	greg	2.1	{
54	greg	2.13	long m;
55	greg	2.1	int e;
56
57	greg	2.13	m = frexp(f, &e) * 0x7fffffff;
58			if (e > 127) { /* overflow */
59			m = m > 0 ? (long)0x7fffffff : -(long)0x7fffffff;
60			e = 127;
61			} else if (e < -128) { /* underflow */
62			m = 0;
63			e = 0;
64			}
65			putint(m, 4, fp);
66	greg	2.1	putint((long)e, 1, fp);
67			}
68
69
70			char *
71	greg	2.15	getstr( /* get null-terminated string */
72			char *s,
73			FILE *fp
74			)
75	greg	2.1	{
76	greg	2.15	char *cp;
77			int c;
78	greg	2.1
79			cp = s;
80			while ((c = getc(fp)) != EOF)
81			if ((*cp++ = c) == '\0')
82			return(s);
83
84			return(NULL);
85			}
86
87
88			long
89	greg	2.15	getint( /* get a siz-byte integer */
90			int siz,
91			FILE *fp
92			)
93	greg	2.1	{
94	greg	2.15	int c;
95			long r;
96	greg	2.1
97			if ((c = getc(fp)) == EOF)
98			return(EOF);
99			r = 0x80&c ? -1<<8\|c : c; /* sign extend */
100			while (--siz > 0) {
101			if ((c = getc(fp)) == EOF)
102			return(EOF);
103			r <<= 8;
104			r \|= c;
105			}
106			return(r);
107			}
108
109
110			double
111	greg	2.15	getflt( /* get a floating point number */
112			FILE *fp
113			)
114	greg	2.1	{
115	greg	2.5	long l;
116	greg	2.2	double d;
117	greg	2.1
118	greg	2.5	l = getint(4, fp);
119	greg	2.15	if (l == EOF && feof(fp)) /* EOF? */
120			return((double)EOF);
121	greg	2.12	if (l == 0) {
122			getc(fp); /* exactly zero -- ignore exponent */
123			return(0.0);
124			}
125	greg	2.5	d = (l + (l > 0 ? .5 : -.5)) * (1./0x7fffffff);
126	greg	2.1	return(ldexp(d, (int)getint(1, fp)));
127			}