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.23 by greg, Fri Jan 15 18:31:38 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 >        siz <<= 3;
37 >        while ((siz -= 8) > 0)
38 >                putc((int)(i>>siz & 0xff), fp);
39 >        putc((int)(i & 0xff), fp);
40   }
41  
42  
43   void
44 < putflt(f, fp)                   /* put out floating point number */
45 < double  f;
46 < FILE  *fp;
44 > putflt(                         /* put out floating point number */
45 >        double  f,
46 >        FILE  *fp
47 > )
48   {
49 +        long  m;
50          int  e;
51  
52 <        putint((long)(frexp(f,&e)*0x7fffffff), 4, fp);
53 <        putint((long)e, 1, fp);
52 >        m = frexp(f, &e) * 0x7fffffff;
53 >        if (e > 127) {                  /* overflow */
54 >                m = m > 0 ? (long)0x7fffffff : -(long)0x7fffffff;
55 >                e = 127;
56 >        } else if (e < -128) {          /* underflow */
57 >                m = 0;
58 >                e = 0;
59 >        }
60 >        putint(m, 4, fp);
61 >        putint(e, 1, fp);
62   }
63  
64  
65 + size_t
66 + putbinary(                      /* fwrite() replacement for small objects */
67 +        const void *p,
68 +        size_t elsiz,
69 +        size_t nel,
70 +        FILE *fp)
71 + {
72 +        const char      *s = (const char *)p;
73 +        size_t          nbytes = elsiz*nel;
74 +
75 +        if (nbytes > 128)
76 +                return(fwrite(p, elsiz, nel, fp));
77 +        
78 +        while (nbytes-- > 0)
79 +                if (putc(*s++, fp) == EOF)
80 +                        return((elsiz*nel - nbytes)/elsiz);
81 +
82 +        return(nel);
83 + }
84 +
85 +
86   char *
87 < getstr(s, fp)                   /* get null-terminated string */
88 < char  *s;
89 < register FILE  *fp;
87 > getstr(                         /* get null-terminated string */
88 >        char  *s,
89 >        FILE  *fp
90 > )
91   {
92 <        register char  *cp;
93 <        register int  c;
92 >        char  *cp;
93 >        int  c;
94  
95          cp = s;
96          while ((c = getc(fp)) != EOF)
# Line 126 | Line 102 | register FILE  *fp;
102  
103  
104   long
105 < getint(siz, fp)                 /* get a siz-byte integer */
106 < int  siz;
107 < register FILE  *fp;
105 > getint(                         /* get a siz-byte integer */
106 >        int  siz,
107 >        FILE  *fp
108 > )
109   {
110 <        register int  c;
111 <        register long  r;
110 >        int  c;
111 >        long  r;
112  
113          if ((c = getc(fp)) == EOF)
114                  return(EOF);
# Line 147 | Line 124 | register FILE  *fp;
124  
125  
126   double
127 < getflt(fp)                      /* get a floating point number */
128 < FILE  *fp;
127 > getflt(                         /* get a floating point number */
128 >        FILE  *fp
129 > )
130   {
131          long    l;
132          double  d;
133  
134          l = getint(4, fp);
135 +        if (l == EOF && feof(fp))       /* EOF? */
136 +                return((double)EOF);
137 +        if (l == 0) {
138 +                getc(fp);               /* exactly zero -- ignore exponent */
139 +                return(0.0);
140 +        }
141          d = (l + (l > 0 ? .5 : -.5)) * (1./0x7fffffff);
142          return(ldexp(d, (int)getint(1, fp)));
143 + }
144 +
145 +
146 + size_t
147 + getbinary(                      /* fread() replacement for small objects */
148 +        void *p,
149 +        size_t elsiz,
150 +        size_t nel,
151 +        FILE *fp)
152 + {
153 +        char    *s = (char *)p;
154 +        size_t  nbytes = elsiz*nel;
155 +        int     c;
156 +
157 +        if (nbytes > 128)
158 +                return(fread(p, elsiz, nel, fp));
159 +        
160 +        while (nbytes-- > 0) {
161 +                if ((c = getc(fp)) == EOF)
162 +                        return((elsiz*nel - nbytes)/elsiz);
163 +                *s++ = c;
164 +        }
165 +        return(nel);
166   }

Diff Legend

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