hypotf.s revision ddc0e0b53c661f6e439e3b7072b3ef353eadb4af
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin/*
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin * CDDL HEADER START
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin *
3e14f97f673e8a630f076077de35afdd43dc1587Roger A. Faulkner * The contents of this file are subject to the terms of the
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin * Common Development and Distribution License (the "License").
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin * You may not use this file except in compliance with the License.
7c2fbfb345896881c631598ee3852ce9ce33fb07April Chin *
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin * or http://www.opensolaris.org/os/licensing.
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin * See the License for the specific language governing permissions
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin * and limitations under the License.
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin *
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin * When distributing Covered Code, include this CDDL HEADER in each
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin * If applicable, add the following below this CDDL HEADER, with the
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin * fields enclosed by brackets "[]" replaced with your own identifying
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin * information: Portions Copyright [yyyy] [name of copyright owner]
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin *
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin * CDDL HEADER END
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin */
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin/*
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin * Copyright 2011 Nexenta Systems, Inc. All rights reserved.
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin */
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin/*
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin * Copyright 2005 Sun Microsystems, Inc. All rights reserved.
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin * Use is subject to license terms.
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin */
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin .file "hypotf.s"
7c2fbfb345896881c631598ee3852ce9ce33fb07April Chin
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin#include "libm.h"
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chinLIBM_ANSI_PRAGMA_WEAK(hypotf,function)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin#include "libm_protos.h"
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin ENTRY(hypotf)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin movl 4(%esp),%eax / eax <-- x
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin andl $0x7fffffff,%eax / eax <-- |x|
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin jz .return_abs_y / if x = +/-0, return |y|
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin subl $0x7f800000,%eax / eax <-- |x| - INF
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin jz .return_abs_x / if x = +/-INF, return |x|
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin movl 8(%esp),%eax / eax <-- y
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin andl $0x7fffffff,%eax / eax <-- |y|
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin jz .return_abs_x / if y = +/-0, return |x|
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin subl $0x7f800000,%eax / eax <-- |y| - INF
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin.return_abs_y:
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin flds 8(%esp) / y
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin jz .take_abs / if y = +/-INF, return |y|
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin fmul %st(0),%st / y*y
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin flds 4(%esp) / x,y*y
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin fmul %st(0),%st / x*x,y*y
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin faddp %st,%st(1) / x*x+y*y
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin fsqrt / sqrt(x*x+y*y)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin subl $4,%esp
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin fstps (%esp) / round to single
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin flds (%esp)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin fwait
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin addl $4,%esp
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin ret
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin.return_abs_x:
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin / returns |x|
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin flds 4(%esp)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin.take_abs:
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin fabs
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin ret
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin .align 4
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin SET_SIZE(hypotf)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin