1 | |
2 | /* |
3 | * Copyright (c) 1998, 2001, Oracle and/or its affiliates. All rights reserved. |
4 | * DO NOT ALTER OR REMOVE COPYRIGHT NOTICES OR THIS FILE HEADER. |
5 | * |
6 | * This code is free software; you can redistribute it and/or modify it |
7 | * under the terms of the GNU General Public License version 2 only, as |
8 | * published by the Free Software Foundation. Oracle designates this |
9 | * particular file as subject to the "Classpath" exception as provided |
10 | * by Oracle in the LICENSE file that accompanied this code. |
11 | * |
12 | * This code is distributed in the hope that it will be useful, but WITHOUT |
13 | * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or |
14 | * FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License |
15 | * version 2 for more details (a copy is included in the LICENSE file that |
16 | * accompanied this code). |
17 | * |
18 | * You should have received a copy of the GNU General Public License version |
19 | * 2 along with this work; if not, write to the Free Software Foundation, |
20 | * Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA. |
21 | * |
22 | * Please contact Oracle, 500 Oracle Parkway, Redwood Shores, CA 94065 USA |
23 | * or visit www.oracle.com if you need additional information or have any |
24 | * questions. |
25 | */ |
26 | /* |
27 | * wrapper atanh(x) |
28 | */ |
29 | |
30 | #include "fdlibm.h" |
31 | |
32 | |
33 | #ifdef __STDC__ |
34 | double atanh(double x) /* wrapper atanh */ |
35 | #else |
36 | double atanh(x) /* wrapper atanh */ |
37 | double x; |
38 | #endif |
39 | { |
40 | #ifdef _IEEE_LIBM |
41 | return __ieee754_atanh(x); |
42 | #else |
43 | double z,y; |
44 | z = __ieee754_atanh(x); |
45 | if(_LIB_VERSION == _IEEE_ || isnan(x)) return z; |
46 | y = fabs(x); |
47 | if(y>=1.0) { |
48 | if(y>1.0) |
49 | return __kernel_standard(x,x,30); /* atanh(|x|>1) */ |
50 | else |
51 | return __kernel_standard(x,x,31); /* atanh(|x|==1) */ |
52 | } else |
53 | return z; |
54 | #endif |
55 | } |
56 | |