Print this page
11210 libm should be cstyle(1ONBLD) clean
*** 20,29 ****
--- 20,30 ----
*/
/*
* Copyright 2011 Nexenta Systems, Inc. All rights reserved.
*/
+
/*
* Copyright 2006 Sun Microsystems, Inc. All rights reserved.
* Use is subject to license terms.
*/
*** 31,61 ****
#include "libm.h" /* fabsl/isnanl/isinfl/signbitl */
#include "complex_wrapper.h"
#include "longdouble.h"
! /* INDENT OFF */
/*
* ldcomplex cacoshl(ldcomplex z);
* cacosh z = +-i cacos z .
* In order to make conj(cacosh(z))=cacosh(conj(z)),
* we define
* cacosh z = sign(Im(z))*i cacos z .
*
*/
- /* INDENT ON */
ldcomplex
! cacoshl(ldcomplex z) {
ldcomplex w, ans;
long double x, y;
w = cacosl(z);
x = LD_RE(z);
y = LD_IM(z);
if (isnanl(y)) {
LD_IM(ans) = y + y;
if (isinfl(x))
LD_RE(ans) = fabsl(x);
else
LD_RE(ans) = y;
} else if (signbitl(y) == 0) {
--- 32,64 ----
#include "libm.h" /* fabsl/isnanl/isinfl/signbitl */
#include "complex_wrapper.h"
#include "longdouble.h"
!
/*
* ldcomplex cacoshl(ldcomplex z);
* cacosh z = +-i cacos z .
* In order to make conj(cacosh(z))=cacosh(conj(z)),
* we define
* cacosh z = sign(Im(z))*i cacos z .
*
*/
ldcomplex
! cacoshl(ldcomplex z)
! {
ldcomplex w, ans;
long double x, y;
w = cacosl(z);
x = LD_RE(z);
y = LD_IM(z);
+
if (isnanl(y)) {
LD_IM(ans) = y + y;
+
if (isinfl(x))
LD_RE(ans) = fabsl(x);
else
LD_RE(ans) = y;
} else if (signbitl(y) == 0) {
*** 63,69 ****
--- 66,73 ----
LD_IM(ans) = LD_RE(w);
} else {
LD_RE(ans) = LD_IM(w);
LD_IM(ans) = -LD_RE(w);
}
+
return (ans);
}