xref: /illumos-gate/usr/src/lib/libm/amd64/src/asinl.S (revision 5d9d9091f564c198a760790b0bfa72c44e17912b)
1*5d9d9091SRichard Lowe/*
2*5d9d9091SRichard Lowe * CDDL HEADER START
3*5d9d9091SRichard Lowe *
4*5d9d9091SRichard Lowe * The contents of this file are subject to the terms of the
5*5d9d9091SRichard Lowe * Common Development and Distribution License (the "License").
6*5d9d9091SRichard Lowe * You may not use this file except in compliance with the License.
7*5d9d9091SRichard Lowe *
8*5d9d9091SRichard Lowe * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
9*5d9d9091SRichard Lowe * or http://www.opensolaris.org/os/licensing.
10*5d9d9091SRichard Lowe * See the License for the specific language governing permissions
11*5d9d9091SRichard Lowe * and limitations under the License.
12*5d9d9091SRichard Lowe *
13*5d9d9091SRichard Lowe * When distributing Covered Code, include this CDDL HEADER in each
14*5d9d9091SRichard Lowe * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
15*5d9d9091SRichard Lowe * If applicable, add the following below this CDDL HEADER, with the
16*5d9d9091SRichard Lowe * fields enclosed by brackets "[]" replaced with your own identifying
17*5d9d9091SRichard Lowe * information: Portions Copyright [yyyy] [name of copyright owner]
18*5d9d9091SRichard Lowe *
19*5d9d9091SRichard Lowe * CDDL HEADER END
20*5d9d9091SRichard Lowe */
21*5d9d9091SRichard Lowe/*
22*5d9d9091SRichard Lowe * Copyright 2011 Nexenta Systems, Inc.  All rights reserved.
23*5d9d9091SRichard Lowe */
24*5d9d9091SRichard Lowe/*
25*5d9d9091SRichard Lowe * Copyright 2006 Sun Microsystems, Inc.  All rights reserved.
26*5d9d9091SRichard Lowe * Use is subject to license terms.
27*5d9d9091SRichard Lowe */
28*5d9d9091SRichard Lowe
29*5d9d9091SRichard Lowe	.file	"asinl.s"
30*5d9d9091SRichard Lowe
31*5d9d9091SRichard Lowe#include "libm.h"
32*5d9d9091SRichard LoweLIBM_ANSI_PRAGMA_WEAK(asinl,function)
33*5d9d9091SRichard Lowe
34*5d9d9091SRichard Lowe	ENTRY(asinl)
35*5d9d9091SRichard Lowe	fldt	8(%rsp)			/ push x
36*5d9d9091SRichard Lowe	fld1				/ push 1
37*5d9d9091SRichard Lowe	fld	%st(1)			/ x , 1 , x
38*5d9d9091SRichard Lowe	fabs				/ |x| , 1 , x
39*5d9d9091SRichard Lowe	fucomip	%st(1),%st
40*5d9d9091SRichard Lowe	ja	9f
41*5d9d9091SRichard Lowe	fadd	%st(1),%st		/ 1+x,x
42*5d9d9091SRichard Lowe	fld1				/ 1,1+x,x
43*5d9d9091SRichard Lowe	fsub	%st(2),%st		/ 1-x,1+x,x
44*5d9d9091SRichard Lowe	fmulp	%st,%st(1)		/ (1-x)*(1+x),x
45*5d9d9091SRichard Lowe	fsqrt				/ sqrt((1-x)*(1+x)),x
46*5d9d9091SRichard Lowe	fpatan				/ atan(x/sqrt((1-x)*(1+x)))
47*5d9d9091SRichard Lowe	ret
48*5d9d9091SRichard Lowe9:
49*5d9d9091SRichard Lowe	/ |x| > 1
50*5d9d9091SRichard Lowe	fstp	%st(0)			/ x
51*5d9d9091SRichard Lowe	fsub	%st,%st(0)		/ +/-0 or NaN+invalid
52*5d9d9091SRichard Lowe	fdiv	%st,%st(0)		/ NaN+invalid or NaN
53*5d9d9091SRichard Lowe	ret
54*5d9d9091SRichard Lowe	.align	16
55*5d9d9091SRichard Lowe	SET_SIZE(asinl)
56