1112158Sdas/****************************************************************
2112158Sdas
3112158SdasThe author of this software is David M. Gay.
4112158Sdas
5112158SdasCopyright (C) 1998, 2000 by Lucent Technologies
6112158SdasAll Rights Reserved
7112158Sdas
8112158SdasPermission to use, copy, modify, and distribute this software and
9112158Sdasits documentation for any purpose and without fee is hereby
10112158Sdasgranted, provided that the above copyright notice appear in all
11112158Sdascopies and that both that the copyright notice and this
12112158Sdaspermission notice and warranty disclaimer appear in supporting
13112158Sdasdocumentation, and that the name of Lucent or any of its entities
14112158Sdasnot be used in advertising or publicity pertaining to
15112158Sdasdistribution of the software without specific, written prior
16112158Sdaspermission.
17112158Sdas
18112158SdasLUCENT DISCLAIMS ALL WARRANTIES WITH REGARD TO THIS SOFTWARE,
19112158SdasINCLUDING ALL IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS.
20112158SdasIN NO EVENT SHALL LUCENT OR ANY OF ITS ENTITIES BE LIABLE FOR ANY
21112158SdasSPECIAL, INDIRECT OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
22112158SdasWHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER
23112158SdasIN AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION,
24112158SdasARISING OUT OF OR IN CONNECTION WITH THE USE OR PERFORMANCE OF
25112158SdasTHIS SOFTWARE.
26112158Sdas
27112158Sdas****************************************************************/
28112158Sdas
29165743Sdas/* Please send bug reports to David M. Gay (dmg at acm dot org,
30165743Sdas * with " at " changed at "@" and " dot " changed to ".").	*/
31112158Sdas
32112158Sdas#include "gdtoaimp.h"
33112158Sdas
34112158Sdas void
35112158Sdas#ifdef KR_headers
36112158SdasULtof(L, bits, exp, k) ULong *L; ULong *bits; Long exp; int k;
37112158Sdas#else
38112158SdasULtof(ULong *L, ULong *bits, Long exp, int k)
39112158Sdas#endif
40112158Sdas{
41112158Sdas	switch(k & STRTOG_Retmask) {
42112158Sdas	  case STRTOG_NoNumber:
43112158Sdas	  case STRTOG_Zero:
44112158Sdas		*L = 0;
45112158Sdas		break;
46112158Sdas
47112158Sdas	  case STRTOG_Normal:
48112158Sdas	  case STRTOG_NaNbits:
49219557Sdas		L[0] = (bits[0] & 0x7fffff) | ((exp + 0x7f + 23) << 23);
50112158Sdas		break;
51112158Sdas
52112158Sdas	  case STRTOG_Denormal:
53112158Sdas		L[0] = bits[0];
54112158Sdas		break;
55112158Sdas
56112158Sdas	  case STRTOG_Infinite:
57112158Sdas		L[0] = 0x7f800000;
58112158Sdas		break;
59112158Sdas
60112158Sdas	  case STRTOG_NaN:
61165743Sdas		L[0] = f_QNAN;
62112158Sdas	  }
63112158Sdas	if (k & STRTOG_Neg)
64112158Sdas		L[0] |= 0x80000000L;
65112158Sdas	}
66112158Sdas
67112158Sdas int
68112158Sdas#ifdef KR_headers
69112158Sdasstrtorf(s, sp, rounding, f) CONST char *s; char **sp; int rounding; float *f;
70112158Sdas#else
71112158Sdasstrtorf(CONST char *s, char **sp, int rounding, float *f)
72112158Sdas#endif
73112158Sdas{
74112158Sdas	static FPI fpi0 = { 24, 1-127-24+1,  254-127-24+1, 1, SI };
75112158Sdas	FPI *fpi, fpi1;
76112158Sdas	ULong bits[1];
77112158Sdas	Long exp;
78112158Sdas	int k;
79112158Sdas
80112158Sdas	fpi = &fpi0;
81112158Sdas	if (rounding != FPI_Round_near) {
82112158Sdas		fpi1 = fpi0;
83112158Sdas		fpi1.rounding = rounding;
84112158Sdas		fpi = &fpi1;
85112158Sdas		}
86112158Sdas	k = strtodg(s, sp, fpi, &exp, bits);
87112158Sdas	ULtof((ULong*)f, bits, exp, k);
88112158Sdas	return k;
89112158Sdas	}
90