11dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project/* s_asinhf.c -- float version of s_asinh.c. 21dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com. 31dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project */ 41dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project 51dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project/* 61dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project * ==================================================== 71dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. 81dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project * 91dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project * Developed at SunPro, a Sun Microsystems, Inc. business. 101dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project * Permission to use, copy, modify, and distribute this 111dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project * software is freely granted, provided that this notice 121dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project * is preserved. 131dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project * ==================================================== 141dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project */ 151dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project 16a0ee07829a9ba7e99ef68e8c12551301cc797f0fElliott Hughes#include <sys/cdefs.h> 17a0ee07829a9ba7e99ef68e8c12551301cc797f0fElliott Hughes__FBSDID("$FreeBSD$"); 181dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project 191dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project#include "math.h" 201dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project#include "math_private.h" 211dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project 221dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Projectstatic const float 231dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Projectone = 1.0000000000e+00, /* 0x3F800000 */ 241dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Projectln2 = 6.9314718246e-01, /* 0x3f317218 */ 251dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Projecthuge= 1.0000000000e+30; 261dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project 271dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Projectfloat 281dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Projectasinhf(float x) 291dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project{ 301dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project float t,w; 311dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project int32_t hx,ix; 321dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project GET_FLOAT_WORD(hx,x); 331dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project ix = hx&0x7fffffff; 341dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project if(ix>=0x7f800000) return x+x; /* x is inf or NaN */ 351dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project if(ix< 0x31800000) { /* |x|<2**-28 */ 361dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project if(huge+x>one) return x; /* return x inexact except 0 */ 371dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project } 381dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project if(ix>0x4d800000) { /* |x| > 2**28 */ 391dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project w = __ieee754_logf(fabsf(x))+ln2; 401dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project } else if (ix>0x40000000) { /* 2**28 > |x| > 2.0 */ 411dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project t = fabsf(x); 421dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project w = __ieee754_logf((float)2.0*t+one/(__ieee754_sqrtf(x*x+one)+t)); 431dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project } else { /* 2.0 > |x| > 2**-28 */ 441dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project t = x*x; 451dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project w =log1pf(fabsf(x)+t/(one+__ieee754_sqrtf(one+t))); 461dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project } 471dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project if(hx>0) return w; else return -w; 481dc9e472e19acfe6dc7f41e429236e7eef7ceda1The Android Open Source Project} 49