Commit | Line | Data |
---|---|---|
1da177e4 LT |
1 | /* Machine-dependent software floating-point definitions. |
2 | S/390 kernel version. | |
3 | Copyright (C) 1997,1998,1999 Free Software Foundation, Inc. | |
4 | This file is part of the GNU C Library. | |
5 | Contributed by Richard Henderson (rth@cygnus.com), | |
6 | Jakub Jelinek (jj@ultra.linux.cz), | |
7 | David S. Miller (davem@redhat.com) and | |
8 | Peter Maydell (pmaydell@chiark.greenend.org.uk). | |
9 | ||
10 | The GNU C Library is free software; you can redistribute it and/or | |
11 | modify it under the terms of the GNU Library General Public License as | |
12 | published by the Free Software Foundation; either version 2 of the | |
13 | License, or (at your option) any later version. | |
14 | ||
15 | The GNU C Library is distributed in the hope that it will be useful, | |
16 | but WITHOUT ANY WARRANTY; without even the implied warranty of | |
17 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
18 | Library General Public License for more details. | |
19 | ||
20 | You should have received a copy of the GNU Library General Public | |
21 | License along with the GNU C Library; see the file COPYING.LIB. If | |
22 | not, write to the Free Software Foundation, Inc., | |
23 | 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. */ | |
24 | ||
25 | #ifndef _SFP_MACHINE_H | |
26 | #define _SFP_MACHINE_H | |
27 | ||
1da177e4 LT |
28 | |
29 | #define _FP_W_TYPE_SIZE 32 | |
30 | #define _FP_W_TYPE unsigned long | |
31 | #define _FP_WS_TYPE signed long | |
32 | #define _FP_I_TYPE long | |
33 | ||
34 | #define _FP_MUL_MEAT_S(R,X,Y) \ | |
35 | _FP_MUL_MEAT_1_wide(_FP_WFRACBITS_S,R,X,Y,umul_ppmm) | |
36 | #define _FP_MUL_MEAT_D(R,X,Y) \ | |
37 | _FP_MUL_MEAT_2_wide(_FP_WFRACBITS_D,R,X,Y,umul_ppmm) | |
38 | #define _FP_MUL_MEAT_Q(R,X,Y) \ | |
39 | _FP_MUL_MEAT_4_wide(_FP_WFRACBITS_Q,R,X,Y,umul_ppmm) | |
40 | ||
41 | #define _FP_DIV_MEAT_S(R,X,Y) _FP_DIV_MEAT_1_udiv(S,R,X,Y) | |
42 | #define _FP_DIV_MEAT_D(R,X,Y) _FP_DIV_MEAT_2_udiv(D,R,X,Y) | |
43 | #define _FP_DIV_MEAT_Q(R,X,Y) _FP_DIV_MEAT_4_udiv(Q,R,X,Y) | |
44 | ||
45 | #define _FP_NANFRAC_S ((_FP_QNANBIT_S << 1) - 1) | |
46 | #define _FP_NANFRAC_D ((_FP_QNANBIT_D << 1) - 1), -1 | |
47 | #define _FP_NANFRAC_Q ((_FP_QNANBIT_Q << 1) - 1), -1, -1, -1 | |
48 | #define _FP_NANSIGN_S 0 | |
49 | #define _FP_NANSIGN_D 0 | |
50 | #define _FP_NANSIGN_Q 0 | |
51 | ||
52 | #define _FP_KEEPNANFRACP 1 | |
53 | ||
54 | /* | |
55 | * If one NaN is signaling and the other is not, | |
56 | * we choose that one, otherwise we choose X. | |
57 | */ | |
58 | #define _FP_CHOOSENAN(fs, wc, R, X, Y, OP) \ | |
59 | do { \ | |
60 | if ((_FP_FRAC_HIGH_RAW_##fs(X) & _FP_QNANBIT_##fs) \ | |
61 | && !(_FP_FRAC_HIGH_RAW_##fs(Y) & _FP_QNANBIT_##fs)) \ | |
62 | { \ | |
63 | R##_s = Y##_s; \ | |
64 | _FP_FRAC_COPY_##wc(R,Y); \ | |
65 | } \ | |
66 | else \ | |
67 | { \ | |
68 | R##_s = X##_s; \ | |
69 | _FP_FRAC_COPY_##wc(R,X); \ | |
70 | } \ | |
71 | R##_c = FP_CLS_NAN; \ | |
72 | } while (0) | |
73 | ||
74 | /* Some assembly to speed things up. */ | |
75 | #define __FP_FRAC_ADD_3(r2,r1,r0,x2,x1,x0,y2,y1,y0) ({ \ | |
76 | unsigned int __r2 = (x2) + (y2); \ | |
77 | unsigned int __r1 = (x1); \ | |
78 | unsigned int __r0 = (x0); \ | |
79 | __asm__ (" alr %2,%3\n" \ | |
80 | " brc 12,0f\n" \ | |
81 | " lhi 0,1\n" \ | |
82 | " alr %1,0\n" \ | |
83 | " brc 12,0f\n" \ | |
84 | " alr %0,0\n" \ | |
85 | "0:" \ | |
86 | : "+&d" (__r2), "+&d" (__r1), "+&d" (__r0) \ | |
87 | : "d" (y0), "i" (1) : "cc", "0" ); \ | |
88 | __asm__ (" alr %1,%2\n" \ | |
89 | " brc 12,0f\n" \ | |
90 | " ahi %0,1\n" \ | |
91 | "0:" \ | |
92 | : "+&d" (__r2), "+&d" (__r1) \ | |
93 | : "d" (y1) : "cc" ); \ | |
94 | (r2) = __r2; \ | |
95 | (r1) = __r1; \ | |
96 | (r0) = __r0; \ | |
97 | }) | |
98 | ||
99 | #define __FP_FRAC_SUB_3(r2,r1,r0,x2,x1,x0,y2,y1,y0) ({ \ | |
100 | unsigned int __r2 = (x2) - (y2); \ | |
101 | unsigned int __r1 = (x1); \ | |
102 | unsigned int __r0 = (x0); \ | |
103 | __asm__ (" slr %2,%3\n" \ | |
104 | " brc 3,0f\n" \ | |
105 | " lhi 0,1\n" \ | |
106 | " slr %1,0\n" \ | |
107 | " brc 3,0f\n" \ | |
108 | " slr %0,0\n" \ | |
109 | "0:" \ | |
110 | : "+&d" (__r2), "+&d" (__r1), "+&d" (__r0) \ | |
111 | : "d" (y0) : "cc", "0" ); \ | |
112 | __asm__ (" slr %1,%2\n" \ | |
113 | " brc 3,0f\n" \ | |
114 | " ahi %0,-1\n" \ | |
115 | "0:" \ | |
116 | : "+&d" (__r2), "+&d" (__r1) \ | |
117 | : "d" (y1) : "cc" ); \ | |
118 | (r2) = __r2; \ | |
119 | (r1) = __r1; \ | |
120 | (r0) = __r0; \ | |
121 | }) | |
122 | ||
123 | #define __FP_FRAC_DEC_3(x2,x1,x0,y2,y1,y0) __FP_FRAC_SUB_3(x2,x1,x0,x2,x1,x0,y2,y1,y0) | |
124 | ||
125 | /* Obtain the current rounding mode. */ | |
126 | #define FP_ROUNDMODE mode | |
127 | ||
128 | /* Exception flags. */ | |
129 | #define FP_EX_INVALID 0x800000 | |
130 | #define FP_EX_DIVZERO 0x400000 | |
131 | #define FP_EX_OVERFLOW 0x200000 | |
132 | #define FP_EX_UNDERFLOW 0x100000 | |
133 | #define FP_EX_INEXACT 0x080000 | |
134 | ||
135 | /* We write the results always */ | |
136 | #define FP_INHIBIT_RESULTS 0 | |
137 | ||
138 | #endif |