Subversion Repositories shark

Rev

Details | Last modification | View Log | RSS feed

Rev Author Line No. Line
2 pj 1
/* w_atanhf.c -- float version of w_atanh.c.
2
 * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
3
 */
4
 
5
/*
6
 * ====================================================
7
 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
8
 *
9
 * Developed at SunPro, a Sun Microsystems, Inc. business.
10
 * Permission to use, copy, modify, and distribute this
11
 * software is freely granted, provided that this notice
12
 * is preserved.
13
 * ====================================================
14
 */
15
 
16
#ifndef lint
17
static char rcsid[] = "$\Id: w_atanhf.c,v 1.2 1995/05/30 05:50:46 rgrimes Exp $";
18
#endif
19
 
20
/*
21
 * wrapper atanhf(x)
22
 */
23
 
24
#include "math.h"
25
#include "math_private.h"
26
 
27
 
28
#ifdef __STDC__
29
        float atanhf(float x)           /* wrapper atanhf */
30
#else
31
        float atanhf(x)                 /* wrapper atanhf */
32
        float x;
33
#endif
34
{
35
#ifdef _IEEE_LIBM
36
        return __ieee754_atanhf(x);
37
#else
38
        float z,y;
39
        z = __ieee754_atanhf(x);
40
        if(_LIB_VERSION == _IEEE_ || isnanf(x)) return z;
41
        y = fabsf(x);
42
        if(y>=(float)1.0) {
43
            if(y>(float)1.0)
44
                /* atanhf(|x|>1) */
45
                return (float)__kernel_standard((double)x,(double)x,130);
46
            else
47
                /* atanhf(|x|==1) */
48
                return (float)__kernel_standard((double)x,(double)x,131);
49
        } else
50
            return z;
51
#endif
52
}