1 1.1 jtc /* w_hypotf.c -- float version of w_hypot.c. 2 1.1 jtc * Conversion to float by Ian Lance Taylor, Cygnus Support, ian (at) cygnus.com. 3 1.1 jtc */ 4 1.1 jtc 5 1.1 jtc /* 6 1.1 jtc * ==================================================== 7 1.1 jtc * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. 8 1.1 jtc * 9 1.1 jtc * Developed at SunPro, a Sun Microsystems, Inc. business. 10 1.1 jtc * Permission to use, copy, modify, and distribute this 11 1.5 simonb * software is freely granted, provided that this notice 12 1.1 jtc * is preserved. 13 1.1 jtc * ==================================================== 14 1.1 jtc */ 15 1.1 jtc 16 1.4 lukem #include <sys/cdefs.h> 17 1.2 jtc #if defined(LIBM_SCCS) && !defined(lint) 18 1.7 drochner __RCSID("$NetBSD: w_hypotf.c,v 1.7 2007/08/10 21:20:36 drochner Exp $"); 19 1.1 jtc #endif 20 1.1 jtc 21 1.1 jtc /* 22 1.1 jtc * wrapper hypotf(x,y) 23 1.1 jtc */ 24 1.1 jtc 25 1.7 drochner #include "namespace.h" 26 1.1 jtc #include "math.h" 27 1.1 jtc #include "math_private.h" 28 1.1 jtc 29 1.7 drochner #ifdef __weak_alias 30 1.7 drochner __weak_alias(hypotf, _hypotf) 31 1.7 drochner #endif 32 1.1 jtc 33 1.6 wiz float 34 1.6 wiz hypotf(float x, float y) /* wrapper hypotf */ 35 1.1 jtc { 36 1.1 jtc #ifdef _IEEE_LIBM 37 1.1 jtc return __ieee754_hypotf(x,y); 38 1.1 jtc #else 39 1.1 jtc float z; 40 1.1 jtc z = __ieee754_hypotf(x,y); 41 1.1 jtc if(_LIB_VERSION == _IEEE_) return z; 42 1.1 jtc if((!finitef(z))&&finitef(x)&&finitef(y)) 43 1.1 jtc /* hypot overflow */ 44 1.1 jtc return (float)__kernel_standard((double)x,(double)y,104); 45 1.1 jtc else 46 1.1 jtc return z; 47 1.1 jtc #endif 48 1.1 jtc } 49