Handle +inf/+inf.
[kopensolaris-gnu/glibc.git] / sysdeps / i386 / fpu / s_fdim.S
1 /* Compute positive difference.
2    Copyright (C) 1997, 1999, 2004 Free Software Foundation, Inc.
3    This file is part of the GNU C Library.
4    Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
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, write to the Free
18    Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
19    02111-1307 USA.  */
20
21 #include <sysdep.h>
22
23         .text
24 ENTRY(__fdim)
25         fldl    4(%esp)         // x
26         fldl    12(%esp)        // x : y
27
28         fucom   %st(1)
29         fnstsw
30         sahf
31         jp      1f
32
33         jc      3f
34
35         fstp    %st(1)
36         fldz
37         jmp     2f
38
39 3:      fsubrp  %st, %st(1)
40         ret
41
42 1:      fxam
43         fnstsw
44         andb    $0x45, %ah
45         cmpb    $0x01, %ah
46         je      2f
47
48         fxch
49 2:      fstp    %st(1)
50         ret
51 END(__fdim)
52 weak_alias (__fdim, fdim)