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

*** 20,45 **** */ /* * Copyright 2011 Nexenta Systems, Inc. All rights reserved. */ /* * Copyright 2006 Sun Microsystems, Inc. All rights reserved. * Use is subject to license terms. */ #pragma weak __csqrtf = csqrtf #include "libm.h" /* sqrt/fabsf/sqrtf */ #include "complex_wrapper.h" - /* INDENT OFF */ static const float zero = 0.0F; ! /* INDENT ON */ fcomplex ! csqrtf(fcomplex z) { fcomplex ans; double dt, dx, dy; float x, y, t, ax, ay, w; int ix, iy, hx, hy; --- 20,46 ---- */ /* * Copyright 2011 Nexenta Systems, Inc. All rights reserved. */ + /* * Copyright 2006 Sun Microsystems, Inc. All rights reserved. * Use is subject to license terms. */ #pragma weak __csqrtf = csqrtf #include "libm.h" /* sqrt/fabsf/sqrtf */ #include "complex_wrapper.h" static const float zero = 0.0F; ! fcomplex ! csqrtf(fcomplex z) ! { fcomplex ans; double dt, dx, dy; float x, y, t, ax, ay, w; int ix, iy, hx, hy;
*** 49,93 **** hy = THE_WORD(y); ix = hx & 0x7fffffff; iy = hy & 0x7fffffff; ay = fabsf(y); ax = fabsf(x); if (ix >= 0x7f800000 || iy >= 0x7f800000) { /* x or y is Inf or NaN */ ! if (iy == 0x7f800000) F_IM(ans) = F_RE(ans) = ay; ! else if (ix == 0x7f800000) { if (hx > 0) { F_RE(ans) = ax; F_IM(ans) = ay * zero; } else { F_RE(ans) = ay * zero; F_IM(ans) = ax; } ! } else F_IM(ans) = F_RE(ans) = ax + ay; } else if (iy == 0) { if (hx >= 0) { F_RE(ans) = sqrtf(ax); F_IM(ans) = zero; } else { F_IM(ans) = sqrtf(ax); F_RE(ans) = zero; } } else { ! dx = (double) ax; ! dy = (double) ay; dt = sqrt(0.5 * (sqrt(dx * dx + dy * dy) + dx)); ! t = (float) dt; ! w = (float) (dy / (dt + dt)); if (hx >= 0) { F_RE(ans) = t; F_IM(ans) = w; } else { F_IM(ans) = t; F_RE(ans) = w; } } if (hy < 0) F_IM(ans) = -F_IM(ans); return (ans); } --- 50,99 ---- hy = THE_WORD(y); ix = hx & 0x7fffffff; iy = hy & 0x7fffffff; ay = fabsf(y); ax = fabsf(x); + if (ix >= 0x7f800000 || iy >= 0x7f800000) { /* x or y is Inf or NaN */ ! if (iy == 0x7f800000) { F_IM(ans) = F_RE(ans) = ay; ! } else if (ix == 0x7f800000) { if (hx > 0) { F_RE(ans) = ax; F_IM(ans) = ay * zero; } else { F_RE(ans) = ay * zero; F_IM(ans) = ax; } ! } else { F_IM(ans) = F_RE(ans) = ax + ay; + } } else if (iy == 0) { if (hx >= 0) { F_RE(ans) = sqrtf(ax); F_IM(ans) = zero; } else { F_IM(ans) = sqrtf(ax); F_RE(ans) = zero; } } else { ! dx = (double)ax; ! dy = (double)ay; dt = sqrt(0.5 * (sqrt(dx * dx + dy * dy) + dx)); ! t = (float)dt; ! w = (float)(dy / (dt + dt)); ! if (hx >= 0) { F_RE(ans) = t; F_IM(ans) = w; } else { F_IM(ans) = t; F_RE(ans) = w; } } + if (hy < 0) F_IM(ans) = -F_IM(ans); + return (ans); }