1cc36ccd1SDavid Schultz /****************************************************************
2cc36ccd1SDavid Schultz
3cc36ccd1SDavid Schultz The author of this software is David M. Gay.
4cc36ccd1SDavid Schultz
5cc36ccd1SDavid Schultz Copyright (C) 1998, 2000 by Lucent Technologies
6cc36ccd1SDavid Schultz All Rights Reserved
7cc36ccd1SDavid Schultz
8cc36ccd1SDavid Schultz Permission to use, copy, modify, and distribute this software and
9cc36ccd1SDavid Schultz its documentation for any purpose and without fee is hereby
10cc36ccd1SDavid Schultz granted, provided that the above copyright notice appear in all
11cc36ccd1SDavid Schultz copies and that both that the copyright notice and this
12cc36ccd1SDavid Schultz permission notice and warranty disclaimer appear in supporting
13cc36ccd1SDavid Schultz documentation, and that the name of Lucent or any of its entities
14cc36ccd1SDavid Schultz not be used in advertising or publicity pertaining to
15cc36ccd1SDavid Schultz distribution of the software without specific, written prior
16cc36ccd1SDavid Schultz permission.
17cc36ccd1SDavid Schultz
18cc36ccd1SDavid Schultz LUCENT DISCLAIMS ALL WARRANTIES WITH REGARD TO THIS SOFTWARE,
19cc36ccd1SDavid Schultz INCLUDING ALL IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS.
20cc36ccd1SDavid Schultz IN NO EVENT SHALL LUCENT OR ANY OF ITS ENTITIES BE LIABLE FOR ANY
21cc36ccd1SDavid Schultz SPECIAL, INDIRECT OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
22cc36ccd1SDavid Schultz WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER
23cc36ccd1SDavid Schultz IN AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION,
24cc36ccd1SDavid Schultz ARISING OUT OF OR IN CONNECTION WITH THE USE OR PERFORMANCE OF
25cc36ccd1SDavid Schultz THIS SOFTWARE.
26cc36ccd1SDavid Schultz
27cc36ccd1SDavid Schultz ****************************************************************/
28cc36ccd1SDavid Schultz
29c88250a5SDavid Schultz /* Please send bug reports to David M. Gay (dmg at acm dot org,
30c88250a5SDavid Schultz * with " at " changed at "@" and " dot " changed to "."). */
31cc36ccd1SDavid Schultz
32cc36ccd1SDavid Schultz #include "gdtoaimp.h"
33cc36ccd1SDavid Schultz
34cc36ccd1SDavid Schultz int
35cc36ccd1SDavid Schultz #ifdef KR_headers
strtopf(s,sp,f)36cc36ccd1SDavid Schultz strtopf(s, sp, f) CONST char *s; char **sp; float *f;
37cc36ccd1SDavid Schultz #else
38cc36ccd1SDavid Schultz strtopf(CONST char *s, char **sp, float *f)
39cc36ccd1SDavid Schultz #endif
40cc36ccd1SDavid Schultz {
414848dd08SDavid Schultz static FPI fpi0 = { 24, 1-127-24+1, 254-127-24+1, 1, SI };
42cc36ccd1SDavid Schultz ULong bits[1], *L;
43cc36ccd1SDavid Schultz Long exp;
44cc36ccd1SDavid Schultz int k;
454848dd08SDavid Schultz #ifdef Honor_FLT_ROUNDS
464848dd08SDavid Schultz #include "gdtoa_fltrnds.h"
474848dd08SDavid Schultz #else
484848dd08SDavid Schultz #define fpi &fpi0
494848dd08SDavid Schultz #endif
50cc36ccd1SDavid Schultz
514848dd08SDavid Schultz k = strtodg(s, sp, fpi, &exp, bits);
52cc36ccd1SDavid Schultz L = (ULong*)f;
53cc36ccd1SDavid Schultz switch(k & STRTOG_Retmask) {
54cc36ccd1SDavid Schultz case STRTOG_NoNumber:
55cc36ccd1SDavid Schultz case STRTOG_Zero:
56cc36ccd1SDavid Schultz L[0] = 0;
57cc36ccd1SDavid Schultz break;
58cc36ccd1SDavid Schultz
59cc36ccd1SDavid Schultz case STRTOG_Normal:
60cc36ccd1SDavid Schultz case STRTOG_NaNbits:
61*50dad48bSDavid Schultz L[0] = (bits[0] & 0x7fffff) | ((exp + 0x7f + 23) << 23);
62cc36ccd1SDavid Schultz break;
63cc36ccd1SDavid Schultz
64cc36ccd1SDavid Schultz case STRTOG_Denormal:
65cc36ccd1SDavid Schultz L[0] = bits[0];
66cc36ccd1SDavid Schultz break;
67cc36ccd1SDavid Schultz
68cc36ccd1SDavid Schultz case STRTOG_Infinite:
69cc36ccd1SDavid Schultz L[0] = 0x7f800000;
70cc36ccd1SDavid Schultz break;
71cc36ccd1SDavid Schultz
72cc36ccd1SDavid Schultz case STRTOG_NaN:
73c88250a5SDavid Schultz L[0] = f_QNAN;
74cc36ccd1SDavid Schultz }
75cc36ccd1SDavid Schultz if (k & STRTOG_Neg)
76cc36ccd1SDavid Schultz L[0] |= 0x80000000L;
77cc36ccd1SDavid Schultz return k;
78cc36ccd1SDavid Schultz }
79