Aren/NetBSD
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Spare ourselves a fabs call. We already check the sign later.

Browse Source 
w = r + y*I is the same as w = r because this is the y == 0 case.

no functional change.
This commit is contained in:
maya 2016-12-31 15:33:03 +00:00
parent 0e2548254a
commit 03a73f9282
1 changed files with 9 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
lib/libm/complex/csqrt.c
View File

@ -1,4 +1,4 @@
/* $NetBSD: csqrt.c,v 1.1 2007/08/20 16:01:37 drochner Exp $ */
/* $NetBSD: csqrt.c,v 1.2 2016/12/31 15:33:03 maya Exp $ */
/*-
* Copyright (c) 2007 The NetBSD Foundation, Inc.
@ -45,23 +45,24 @@ csqrt(double complex z)
if (x == 0.0) {
w = 0.0 + y * I;
} else {
r = fabs(x);
r = sqrt(r);
if (x < 0.0) {
r = sqrt(-x);
w = 0.0 + r * I;
} else {
w = r + y * I;
r = sqrt(x);
w = r;
}
}
return w;
}
if (x == 0.0) {
r = fabs(y);
r = sqrt(0.5 * r);
if (y > 0)
if (y > 0) {
r = sqrt(0.5 * y);
w = r + r * I;
else
} else {
r = sqrt(-0.5 * y);
w = r - r * I;
}
return w;
}
/* Rescale to avoid internal overflow or underflow. */