4 * The contents of this file are subject to the terms of the
5 * Common Development and Distribution License (the "License").
6 * You may not use this file except in compliance with the License.
8 * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
9 * or http://www.opensolaris.org/os/licensing.
10 * See the License for the specific language governing permissions
11 * and limitations under the License.
13 * When distributing Covered Code, include this CDDL HEADER in each
14 * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
15 * If applicable, add the following below this CDDL HEADER, with the
16 * fields enclosed by brackets "[]" replaced with your own identifying
17 * information: Portions Copyright [yyyy] [name of copyright owner]
22 * Copyright 2011 Nexenta Systems, Inc. All rights reserved.
25 * Copyright 2005 Sun Microsystems, Inc. All rights reserved.
26 * Use is subject to license terms.
32 LIBM_ANSI_PRAGMA_WEAK
(hypotf
,function
)
33 #include "libm_protos.h"
36 movl
4(%esp
),%eax
/ eax
<-- x
37 andl $
0x7fffffff,%eax
/ eax
<-- |x|
38 jz
.return_abs_y / if x = +/-0, return |y|
39 subl $
0x7f800000,%eax
/ eax
<-- |x|
- INF
40 jz
.return_abs_x / if x = +/-INF, return |x|
41 movl
8(%esp
),%eax
/ eax
<-- y
42 andl $
0x7fffffff,%eax
/ eax
<-- |y|
43 jz
.return_abs_x / if y = +/-0, return |x|
44 subl $
0x7f800000,%eax
/ eax
<-- |y|
- INF
47 jz
.take_abs / if y = +/-INF, return |y|
50 fmul %st(0),%st / x
*x
,y
*y
51 faddp
%st,%st(1) / x
*x+y
*y
54 fstps
(%esp
) / round to single