Print this page
11210 libm should be cstyle(1ONBLD) clean
*** 20,37 ****
*/
/*
* Copyright 2011 Nexenta Systems, Inc. All rights reserved.
*/
/*
* Copyright 2006 Sun Microsystems, Inc. All rights reserved.
* Use is subject to license terms.
*/
#pragma weak __atanh = atanh
! /* INDENT OFF */
/*
* atanh(x)
* Code originated from 4.3bsd.
* Modified by K.C. Ng for SUN 4.0 libm.
* Method :
--- 20,38 ----
*/
/*
* Copyright 2011 Nexenta Systems, Inc. All rights reserved.
*/
+
/*
* Copyright 2006 Sun Microsystems, Inc. All rights reserved.
* Use is subject to license terms.
*/
#pragma weak __atanh = atanh
!
/*
* atanh(x)
* Code originated from 4.3bsd.
* Modified by K.C. Ng for SUN 4.0 libm.
* Method :
*** 46,70 ****
* Special cases:
* atanh(x) is NaN if |x| > 1 with signal;
* atanh(NaN) is that NaN with no signal;
* atanh(+-1) is +-INF with signal.
*/
- /* INDENT ON */
#include "libm.h"
#include "libm_protos.h"
#include <math.h>
double
! atanh(double x) {
double t;
if (isnan(x))
return (x * x); /* switched from x + x for Cheetah */
t = fabs(x);
if (t > 1.0)
return (_SVID_libm_err(x, x, 30)); /* sNaN */
if (t == 1.0)
return (_SVID_libm_err(x, x, 31)); /* x/0; */
t = t / (1.0 - t);
return (copysign(0.5, x) * log1p(t + t));
}
--- 47,75 ----
* Special cases:
* atanh(x) is NaN if |x| > 1 with signal;
* atanh(NaN) is that NaN with no signal;
* atanh(+-1) is +-INF with signal.
*/
#include "libm.h"
#include "libm_protos.h"
#include <math.h>
double
! atanh(double x)
! {
double t;
if (isnan(x))
return (x * x); /* switched from x + x for Cheetah */
+
t = fabs(x);
+
if (t > 1.0)
return (_SVID_libm_err(x, x, 30)); /* sNaN */
+
if (t == 1.0)
return (_SVID_libm_err(x, x, 31)); /* x/0; */
+
t = t / (1.0 - t);
return (copysign(0.5, x) * log1p(t + t));
}