ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/common/portio.c
(Generate patch)

Comparing ray/src/common/portio.c (file contents):
Revision 2.6 by greg, Sat Feb 22 02:07:22 2003 UTC vs.
Revision 2.24 by greg, Fri Feb 19 02:59:32 2021 UTC

# Line 4 | Line 4 | static const char      RCSid[] = "$Id$";
4   /*
5   * Portable i/o for binary files
6   *
7 < * External symbols declared in standard.h
7 > * External symbols declared in rtio.h
8   */
9  
10 < /* ====================================================================
11 < * The Radiance Software License, Version 1.0
12 < *
13 < * Copyright (c) 1990 - 2002 The Regents of the University of California,
14 < * through Lawrence Berkeley National Laboratory.   All rights reserved.
15 < *
16 < * Redistribution and use in source and binary forms, with or without
17 < * modification, are permitted provided that the following conditions
18 < * are met:
19 < *
20 < * 1. Redistributions of source code must retain the above copyright
21 < *         notice, this list of conditions and the following disclaimer.
22 < *
23 < * 2. Redistributions in binary form must reproduce the above copyright
24 < *       notice, this list of conditions and the following disclaimer in
25 < *       the documentation and/or other materials provided with the
26 < *       distribution.
27 < *
28 < * 3. The end-user documentation included with the redistribution,
29 < *           if any, must include the following acknowledgment:
30 < *             "This product includes Radiance software
31 < *                 (http://radsite.lbl.gov/)
32 < *                 developed by the Lawrence Berkeley National Laboratory
33 < *               (http://www.lbl.gov/)."
34 < *       Alternately, this acknowledgment may appear in the software itself,
35 < *       if and wherever such third-party acknowledgments normally appear.
36 < *
37 < * 4. The names "Radiance," "Lawrence Berkeley National Laboratory"
38 < *       and "The Regents of the University of California" must
39 < *       not be used to endorse or promote products derived from this
40 < *       software without prior written permission. For written
41 < *       permission, please contact [email protected].
42 < *
43 < * 5. Products derived from this software may not be called "Radiance",
44 < *       nor may "Radiance" appear in their name, without prior written
45 < *       permission of Lawrence Berkeley National Laboratory.
46 < *
47 < * THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESSED OR IMPLIED
48 < * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
49 < * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
50 < * DISCLAIMED.   IN NO EVENT SHALL Lawrence Berkeley National Laboratory OR
51 < * ITS CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
52 < * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
53 < * LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF
54 < * USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
55 < * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
56 < * OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT
57 < * OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
58 < * SUCH DAMAGE.
59 < * ====================================================================
60 < *
61 < * This software consists of voluntary contributions made by many
62 < * individuals on behalf of Lawrence Berkeley National Laboratory.   For more
63 < * information on Lawrence Berkeley National Laboratory, please see
64 < * <http://www.lbl.gov/>.
65 < */
10 > #include "copyright.h"
11  
12 < #include <stdio.h>
12 > #include "rtio.h"
13  
14 < #ifndef frexp
70 < extern double  frexp();
71 < #endif
72 < #ifndef ldexp
73 < extern double  ldexp();
74 < #endif
14 > #include <math.h>
15  
16  
17   void
18 < putstr(s, fp)                   /* write null-terminated string to fp */
19 < register char  *s;
20 < register FILE  *fp;
18 > putstr(                         /* write null-terminated string to fp */
19 >        char  *s,
20 >        FILE  *fp
21 > )
22   {
23          do
24                  putc(*s, fp);
# Line 86 | Line 27 | register FILE  *fp;
27  
28  
29   void
30 < putint(i, siz, fp)              /* write a siz-byte integer to fp */
31 < long  i;
32 < register int  siz;
33 < register FILE  *fp;
30 > putint(                         /* write a siz-byte integer to fp */
31 >        long  i,
32 >        int  siz,
33 >        FILE  *fp
34 > )
35   {
36 <        while (siz--)
37 <                putc((int)(i>>(siz<<3) & 0xff), fp);
36 >        while (siz > sizeof(long)) {
37 >                putc((i<0)*0xff, fp);
38 >                siz--;
39 >        }
40 >        siz <<= 3;
41 >        while ((siz -= 8) > 0)
42 >                putc((int)(i>>siz & 0xff), fp);
43 >        putc((int)(i & 0xff), fp);
44   }
45  
46  
47   void
48 < putflt(f, fp)                   /* put out floating point number */
49 < double  f;
50 < FILE  *fp;
48 > putflt(                         /* put out floating point number */
49 >        double  f,
50 >        FILE  *fp
51 > )
52   {
53 +        long  m;
54          int  e;
55  
56 <        putint((long)(frexp(f,&e)*0x7fffffff), 4, fp);
57 <        putint((long)e, 1, fp);
56 >        m = frexp(f, &e) * 0x7fffffff;
57 >        if (e > 127) {                  /* overflow */
58 >                m = m > 0 ? (long)0x7fffffff : -(long)0x7fffffff;
59 >                e = 127;
60 >        } else if (e < -128) {          /* underflow */
61 >                m = 0;
62 >                e = 0;
63 >        }
64 >        putint(m, 4, fp);
65 >        putint(e, 1, fp);
66   }
67  
68  
69 + size_t
70 + putbinary(                      /* fwrite() replacement for small objects */
71 +        const void *p,
72 +        size_t elsiz,
73 +        size_t nel,
74 +        FILE *fp)
75 + {
76 +        const char      *s = (const char *)p;
77 +        size_t          nbytes = elsiz*nel;
78 +
79 +        if (nbytes > 128)
80 +                return(fwrite(p, elsiz, nel, fp));
81 +        
82 +        while (nbytes-- > 0)
83 +                if (putc(*s++, fp) == EOF)
84 +                        return((elsiz*nel - nbytes)/elsiz);
85 +
86 +        return(nel);
87 + }
88 +
89 +
90   char *
91 < getstr(s, fp)                   /* get null-terminated string */
92 < char  *s;
93 < register FILE  *fp;
91 > getstr(                         /* get null-terminated string */
92 >        char  *s,
93 >        FILE  *fp
94 > )
95   {
96 <        register char  *cp;
97 <        register int  c;
96 >        char  *cp;
97 >        int  c;
98  
99          cp = s;
100          while ((c = getc(fp)) != EOF)
# Line 126 | Line 106 | register FILE  *fp;
106  
107  
108   long
109 < getint(siz, fp)                 /* get a siz-byte integer */
110 < int  siz;
111 < register FILE  *fp;
109 > getint(                         /* get a siz-byte integer */
110 >        int  siz,
111 >        FILE  *fp
112 > )
113   {
114 <        register int  c;
115 <        register long  r;
114 >        int  c;
115 >        long  r;
116  
117          if ((c = getc(fp)) == EOF)
118                  return(EOF);
119 <        r = 0x80&c ? -1<<8|c : c;               /* sign extend */
119 >        r = c;
120 >        if (c & 0x80)           /* sign extend? */
121 >                r |= ~255L;
122          while (--siz > 0) {
123                  if ((c = getc(fp)) == EOF)
124                          return(EOF);
# Line 147 | Line 130 | register FILE  *fp;
130  
131  
132   double
133 < getflt(fp)                      /* get a floating point number */
134 < FILE  *fp;
133 > getflt(                         /* get a floating point number */
134 >        FILE  *fp
135 > )
136   {
137          long    l;
138          double  d;
139  
140          l = getint(4, fp);
141 +        if (l == EOF && feof(fp))       /* EOF? */
142 +                return((double)EOF);
143 +        if (l == 0) {
144 +                getc(fp);               /* exactly zero -- ignore exponent */
145 +                return(0.0);
146 +        }
147          d = (l + (l > 0 ? .5 : -.5)) * (1./0x7fffffff);
148          return(ldexp(d, (int)getint(1, fp)));
149 + }
150 +
151 +
152 + size_t
153 + getbinary(                      /* fread() replacement for small objects */
154 +        void *p,
155 +        size_t elsiz,
156 +        size_t nel,
157 +        FILE *fp)
158 + {
159 +        char    *s = (char *)p;
160 +        size_t  nbytes = elsiz*nel;
161 +        int     c;
162 +
163 +        if (nbytes > 128)
164 +                return(fread(p, elsiz, nel, fp));
165 +        
166 +        while (nbytes-- > 0) {
167 +                if ((c = getc(fp)) == EOF)
168 +                        return((elsiz*nel - nbytes)/elsiz);
169 +                *s++ = c;
170 +        }
171 +        return(nel);
172   }

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines