xref: /freebsd-src/contrib/gdtoa/strtorf.c (revision 50dad48bb740a8e56d185d9e8c165e0758f46e25)
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  void
35cc36ccd1SDavid Schultz #ifdef KR_headers
ULtof(L,bits,exp,k)36cc36ccd1SDavid Schultz ULtof(L, bits, exp, k) ULong *L; ULong *bits; Long exp; int k;
37cc36ccd1SDavid Schultz #else
38cc36ccd1SDavid Schultz ULtof(ULong *L, ULong *bits, Long exp, int k)
39cc36ccd1SDavid Schultz #endif
40cc36ccd1SDavid Schultz {
41cc36ccd1SDavid Schultz 	switch(k & STRTOG_Retmask) {
42cc36ccd1SDavid Schultz 	  case STRTOG_NoNumber:
43cc36ccd1SDavid Schultz 	  case STRTOG_Zero:
44cc36ccd1SDavid Schultz 		*L = 0;
45cc36ccd1SDavid Schultz 		break;
46cc36ccd1SDavid Schultz 
47cc36ccd1SDavid Schultz 	  case STRTOG_Normal:
48cc36ccd1SDavid Schultz 	  case STRTOG_NaNbits:
49*50dad48bSDavid Schultz 		L[0] = (bits[0] & 0x7fffff) | ((exp + 0x7f + 23) << 23);
50cc36ccd1SDavid Schultz 		break;
51cc36ccd1SDavid Schultz 
52cc36ccd1SDavid Schultz 	  case STRTOG_Denormal:
53cc36ccd1SDavid Schultz 		L[0] = bits[0];
54cc36ccd1SDavid Schultz 		break;
55cc36ccd1SDavid Schultz 
56cc36ccd1SDavid Schultz 	  case STRTOG_Infinite:
57cc36ccd1SDavid Schultz 		L[0] = 0x7f800000;
58cc36ccd1SDavid Schultz 		break;
59cc36ccd1SDavid Schultz 
60cc36ccd1SDavid Schultz 	  case STRTOG_NaN:
61c88250a5SDavid Schultz 		L[0] = f_QNAN;
62cc36ccd1SDavid Schultz 	  }
63cc36ccd1SDavid Schultz 	if (k & STRTOG_Neg)
64cc36ccd1SDavid Schultz 		L[0] |= 0x80000000L;
65cc36ccd1SDavid Schultz 	}
66cc36ccd1SDavid Schultz 
67cc36ccd1SDavid Schultz  int
68cc36ccd1SDavid Schultz #ifdef KR_headers
strtorf(s,sp,rounding,f)69cc36ccd1SDavid Schultz strtorf(s, sp, rounding, f) CONST char *s; char **sp; int rounding; float *f;
70cc36ccd1SDavid Schultz #else
71cc36ccd1SDavid Schultz strtorf(CONST char *s, char **sp, int rounding, float *f)
72cc36ccd1SDavid Schultz #endif
73cc36ccd1SDavid Schultz {
74cc36ccd1SDavid Schultz 	static FPI fpi0 = { 24, 1-127-24+1,  254-127-24+1, 1, SI };
75cc36ccd1SDavid Schultz 	FPI *fpi, fpi1;
76cc36ccd1SDavid Schultz 	ULong bits[1];
77cc36ccd1SDavid Schultz 	Long exp;
78cc36ccd1SDavid Schultz 	int k;
79cc36ccd1SDavid Schultz 
80cc36ccd1SDavid Schultz 	fpi = &fpi0;
81cc36ccd1SDavid Schultz 	if (rounding != FPI_Round_near) {
82cc36ccd1SDavid Schultz 		fpi1 = fpi0;
83cc36ccd1SDavid Schultz 		fpi1.rounding = rounding;
84cc36ccd1SDavid Schultz 		fpi = &fpi1;
85cc36ccd1SDavid Schultz 		}
86cc36ccd1SDavid Schultz 	k = strtodg(s, sp, fpi, &exp, bits);
87cc36ccd1SDavid Schultz 	ULtof((ULong*)f, bits, exp, k);
88cc36ccd1SDavid Schultz 	return k;
89cc36ccd1SDavid Schultz 	}
90