Print this page
11210 libm should be cstyle(1ONBLD) clean

*** 16,28 **** --- 16,30 ---- * fields enclosed by brackets "[]" replaced with your own identifying * information: Portions Copyright [yyyy] [name of copyright owner] * * CDDL HEADER END */ + /* * Copyright 2011 Nexenta Systems, Inc. All rights reserved. */ + /* * Copyright 2005 Sun Microsystems, Inc. All rights reserved. * Use is subject to license terms. */
*** 34,44 **** #if defined(__i386) && !defined(__amd64) extern int __swapRP(int); #endif fcomplex ! clogf(fcomplex z) { fcomplex ans; float x, y, ax, ay; double dx, dy; int ix, iy, hx, hy; --- 36,47 ---- #if defined(__i386) && !defined(__amd64) extern int __swapRP(int); #endif fcomplex ! clogf(fcomplex z) ! { fcomplex ans; float x, y, ax, ay; double dx, dy; int ix, iy, hx, hy;
*** 49,58 **** --- 52,62 ---- ix = hx & 0x7fffffff; iy = hy & 0x7fffffff; ay = fabsf(y); ax = fabsf(x); F_IM(ans) = atan2f(y, x); + if (ix >= 0x7f800000 || iy >= 0x7f800000) { /* x or y is Inf or NaN */ if (iy == 0x7f800000) F_RE(ans) = ay; else if (ix == 0x7f800000)
*** 63,82 **** --- 67,89 ---- #if defined(__i386) && !defined(__amd64) int rp = __swapRP(fp_extended); #endif dx = (double)ax; dy = (double)ay; + if (ix == 0x3f800000) F_RE(ans) = (float)(0.5 * log1p(dy * dy)); else if (iy == 0x3f800000) F_RE(ans) = (float)(0.5 * log1p(dx * dx)); else if ((ix | iy) == 0) F_RE(ans) = -1.0f / ax; else F_RE(ans) = (float)(0.5 * log(dx * dx + dy * dy)); + #if defined(__i386) && !defined(__amd64) if (rp != fp_extended) (void) __swapRP(rp); #endif } + return (ans); }