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 2005 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 "hypotf.s" 30*5d9d9091SRichard Lowe 31*5d9d9091SRichard Lowe#include "libm.h" 32*5d9d9091SRichard LoweLIBM_ANSI_PRAGMA_WEAK(hypotf,function) 33*5d9d9091SRichard Lowe#include "libm_protos.h" 34*5d9d9091SRichard Lowe 35*5d9d9091SRichard Lowe ENTRY(hypotf) 36*5d9d9091SRichard Lowe movl 4(%esp),%eax / eax <-- x 37*5d9d9091SRichard Lowe andl $0x7fffffff,%eax / eax <-- |x| 38*5d9d9091SRichard Lowe jz .return_abs_y / if x = +/-0, return |y| 39*5d9d9091SRichard Lowe subl $0x7f800000,%eax / eax <-- |x| - INF 40*5d9d9091SRichard Lowe jz .return_abs_x / if x = +/-INF, return |x| 41*5d9d9091SRichard Lowe movl 8(%esp),%eax / eax <-- y 42*5d9d9091SRichard Lowe andl $0x7fffffff,%eax / eax <-- |y| 43*5d9d9091SRichard Lowe jz .return_abs_x / if y = +/-0, return |x| 44*5d9d9091SRichard Lowe subl $0x7f800000,%eax / eax <-- |y| - INF 45*5d9d9091SRichard Lowe.return_abs_y: 46*5d9d9091SRichard Lowe flds 8(%esp) / y 47*5d9d9091SRichard Lowe jz .take_abs / if y = +/-INF, return |y| 48*5d9d9091SRichard Lowe fmul %st(0),%st / y*y 49*5d9d9091SRichard Lowe flds 4(%esp) / x,y*y 50*5d9d9091SRichard Lowe fmul %st(0),%st / x*x,y*y 51*5d9d9091SRichard Lowe faddp %st,%st(1) / x*x+y*y 52*5d9d9091SRichard Lowe fsqrt / sqrt(x*x+y*y) 53*5d9d9091SRichard Lowe subl $4,%esp 54*5d9d9091SRichard Lowe fstps (%esp) / round to single 55*5d9d9091SRichard Lowe flds (%esp) 56*5d9d9091SRichard Lowe fwait 57*5d9d9091SRichard Lowe addl $4,%esp 58*5d9d9091SRichard Lowe ret 59*5d9d9091SRichard Lowe 60*5d9d9091SRichard Lowe.return_abs_x: 61*5d9d9091SRichard Lowe / returns |x| 62*5d9d9091SRichard Lowe flds 4(%esp) 63*5d9d9091SRichard Lowe.take_abs: 64*5d9d9091SRichard Lowe fabs 65*5d9d9091SRichard Lowe ret 66*5d9d9091SRichard Lowe .align 4 67*5d9d9091SRichard Lowe SET_SIZE(hypotf) 68