]>
Commit | Line | Data |
---|---|---|
8b954ab9 | 1 | /* floor function, sparc32 v9 vis2 version. |
f7a9f785 | 2 | Copyright (C) 2013-2016 Free Software Foundation, Inc. |
8b954ab9 DM |
3 | This file is part of the GNU C Library. |
4 | Contributed by David S. Miller <davem@davemloft.net>, 2013. | |
5 | ||
6 | The GNU C Library is free software; you can redistribute it and/or | |
7 | modify it under the terms of the GNU Lesser General Public | |
8 | License as published by the Free Software Foundation; either | |
9 | version 2.1 of the License, or (at your option) any later version. | |
10 | ||
11 | The GNU C Library is distributed in the hope that it will be useful, | |
12 | but WITHOUT ANY WARRANTY; without even the implied warranty of | |
13 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
14 | Lesser General Public License for more details. | |
15 | ||
16 | You should have received a copy of the GNU Lesser General Public | |
17 | License along with the GNU C Library; if not, see | |
18 | <http://www.gnu.org/licenses/>. */ | |
19 | ||
20 | #include <sysdep.h> | |
21 | ||
22 | /* 'siam' (Set Interval Arithmetic Mode) is used to quickly override | |
23 | the rounding mode during this routine. | |
24 | ||
25 | We add then subtract (or subtract than add if the initial | |
26 | value was negative) 2**23 to the value, then subtract it | |
27 | back out. | |
28 | ||
29 | This will clear out the fractional portion of the value and, | |
30 | with suitable 'siam' initiated rouding mode settings, round | |
31 | the final result in the proper direction. */ | |
32 | ||
33 | #define TWO_FIFTYTWO 0x43300000 /* 2**52 */ | |
34 | ||
35 | #define ZERO %f10 /* 0.0 */ | |
36 | #define SIGN_BIT %f12 /* -0.0 */ | |
37 | ||
38 | ENTRY (__floor_vis2) | |
39 | sethi %hi(TWO_FIFTYTWO), %o2 | |
40 | sllx %o0, 32, %o0 | |
41 | or %o0, %o1, %o0 | |
42 | stx %o0, [%sp + 72] | |
43 | sllx %o2, 32, %o2 | |
44 | fzero ZERO | |
45 | ldd [%sp + 72], %f0 | |
46 | fnegd ZERO, SIGN_BIT | |
47 | stx %o2, [%sp + 72] | |
48 | fabsd %f0, %f14 | |
49 | ldd [%sp + 72], %f16 | |
50 | fcmpd %fcc3, %f14, %f16 | |
51 | fmovduge %fcc3, ZERO, %f16 | |
52 | fand %f0, SIGN_BIT, SIGN_BIT | |
53 | for %f16, SIGN_BIT, %f16 | |
54 | siam (1 << 2) | 3 | |
55 | faddd %f0, %f16, %f18 | |
56 | siam (1 << 2) | 0 | |
57 | fsubd %f18, %f16, %f18 | |
58 | siam (0 << 2) | |
59 | retl | |
60 | for %f18, SIGN_BIT, %f0 | |
61 | END (__floor_vis2) |