isnan implementation for Power5.
[kopensolaris-gnu/glibc.git] / sysdeps / powerpc / powerpc32 / power5 / fpu / s_isnan.S
1 /* isnan().  PowerPC32 version.
2    Copyright (C) 2008 Free Software Foundation, Inc.
3    This file is part of the GNU C Library.
4
5    The GNU C Library is free software; you can redistribute it and/or
6    modify it under the terms of the GNU Lesser General Public
7    License as published by the Free Software Foundation; either
8    version 2.1 of the License, or (at your option) any later version.
9
10    The GNU C Library is distributed in the hope that it will be useful,
11    but WITHOUT ANY WARRANTY; without even the implied warranty of
12    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
13    Lesser General Public License for more details.
14
15    You should have received a copy of the GNU Lesser General Public
16    License along with the GNU C Library; if not, write to the Free
17    Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
18    02111-1307 USA.  */
19
20 #include <sysdep.h>
21 #include <math_ldbl_opt.h>
22
23 /* int __isnan(x)  */
24         .machine power5
25 EALIGN (__isnan, 4, 0)
26         stwu    r1,-32(r1)
27         cfi_adjust_cfa_offset (32)
28         ori     r1,r1,0
29         stfd    fp1,24(r1)      /* copy FPR to GPR */
30         ori     r1,r1,0
31         lwz     r4,24(r1)
32         lwz     r5,28(r1)
33         lis     r0,0x7ff0       /* const long r0 0x7ff00000 00000000 */
34         clrlwi  r4,r4,1         /* x = fabs(x) */
35         cmpw    cr7,r4,r0       /* if (fabs(x) =< inf) */
36         cmpwi   cr6,r5,0
37         li      r3,0            /* then return 0 */
38         addi    r1,r1,32
39         cfi_adjust_cfa_offset (-32)
40         bltlr+  cr7
41         bgt-    cr7,L(NaN)
42         beqlr+  cr6
43 L(NaN):
44         li      r3,1            /* else return 1 */
45         blr
46         END (__isnan)
47
48 hidden_def (__isnan)
49 weak_alias (__isnan, isnan)
50
51 #ifdef NO_LONG_DOUBLE
52 strong_alias (__isnan, __isnanl)
53 weak_alias (__isnan, isnanl)
54 #endif
55
56 #ifndef IS_IN_libm
57 # if LONG_DOUBLE_COMPAT(libc, GLIBC_2_0)
58 compat_symbol (libc, __isnan, __isnanl, GLIBC_2_0);
59 compat_symbol (libc, isnan, isnanl, GLIBC_2_0);
60 # endif
61 #endif
62