1#ifndef _MATH_PRIVATE_H_
2#error "Never use <math_ldbl.h> directly; include <math_private.h> instead."
3#endif
4
5/* A union which permits us to convert between a long double and
6   three 32 bit ints.  */
7
8#if __FLOAT_WORD_ORDER == BIG_ENDIAN
9
10typedef union
11{
12  long double value;
13  struct
14  {
15    int sign_exponent:16;
16    unsigned int empty:16;
17    u_int32_t msw;
18    u_int32_t lsw;
19  } parts;
20} ieee_long_double_shape_type;
21
22#endif
23
24#if __FLOAT_WORD_ORDER == LITTLE_ENDIAN
25
26typedef union
27{
28  long double value;
29  struct
30  {
31    u_int32_t lsw;
32    u_int32_t msw;
33    int sign_exponent:16;
34    unsigned int empty:16;
35  } parts;
36} ieee_long_double_shape_type;
37
38#endif
39
40/* Get three 32 bit ints from a double.  */
41
42#define GET_LDOUBLE_WORDS(exp,ix0,ix1,d)			\
43do {								\
44  ieee_long_double_shape_type ew_u;				\
45  ew_u.value = (d);						\
46  (exp) = ew_u.parts.sign_exponent;				\
47  (ix0) = ew_u.parts.msw;					\
48  (ix1) = ew_u.parts.lsw;					\
49} while (0)
50
51/* Set a double from two 32 bit ints.  */
52
53#define SET_LDOUBLE_WORDS(d,exp,ix0,ix1)			\
54do {								\
55  ieee_long_double_shape_type iw_u;				\
56  iw_u.parts.sign_exponent = (exp);				\
57  iw_u.parts.msw = (ix0);					\
58  iw_u.parts.lsw = (ix1);					\
59  (d) = iw_u.value;						\
60} while (0)
61
62/* Get the more significant 32 bits of a long double mantissa.  */
63
64#define GET_LDOUBLE_MSW(v,d)					\
65do {								\
66  ieee_long_double_shape_type sh_u;				\
67  sh_u.value = (d);						\
68  (v) = sh_u.parts.msw;						\
69} while (0)
70
71/* Set the more significant 32 bits of a long double mantissa from an int.  */
72
73#define SET_LDOUBLE_MSW(d,v)					\
74do {								\
75  ieee_long_double_shape_type sh_u;				\
76  sh_u.value = (d);						\
77  sh_u.parts.msw = (v);						\
78  (d) = sh_u.value;						\
79} while (0)
80
81/* Get int from the exponent of a long double.  */
82
83#define GET_LDOUBLE_EXP(exp,d)					\
84do {								\
85  ieee_long_double_shape_type ge_u;				\
86  ge_u.value = (d);						\
87  (exp) = ge_u.parts.sign_exponent;				\
88} while (0)
89
90/* Set exponent of a long double from an int.  */
91
92#define SET_LDOUBLE_EXP(d,exp)					\
93do {								\
94  ieee_long_double_shape_type se_u;				\
95  se_u.value = (d);						\
96  se_u.parts.sign_exponent = (exp);				\
97  (d) = se_u.value;						\
98} while (0)
99