Appease static analyzers by making all code paths which assign
values into p,q sane. Get rid of redundant assignment. Indent
for legibility. NFC.

This doesn't create a functional difference, as all callers
test number >= 0x40000000 anyway.

To see this, note the following:
- consistently, hx is the high bits of x, lx is the low bits,
x is the float.
- & 0x7fffffff zeroes the sign bit, as does fabs.

A case where it isn't easy to see that there's no functional
change is y1, which does:

ix = hx & 0x7fffffff (zero signbit of high bits of x)
y = fabs(x) (this has a zeroed signbit but otherwise same as x)
ix >= 0x40000000
pone(y); qone(y)

qone(x) (also pone) do:
ix = hx & 0x7fffffff

ix in qone and in the calling function are the same number,
and the comparison applies for both, and ix < 0x40000000 isn't
possible.

(Also, no explosions seem to happen when I feed it random numbers)

Add C99 complex support, for double and float.
Most complex function implementations are from the "c9x-complex" library,
originating from the "cephes" math library, see
http://www.netlib.org/cephes/, from Stephen L. Moshier, incorporated and
redistributed with the NetBSD license by permission of the author.

Error behaviour and other boundary conditions (branch cuts)
need to be looked at.

For namespace sanity, I've done the rename/weak alias procedure to
most of the exported functions which are also used internally.
Didn't do so for sin/cos(f) yet because assembler implementations use
them directly, and renaming functions shared between the main libm
and the machine specific "overlay" might raise binary compatibility
issues.

Coverity CID 572: Comment out dead code.

Coverity CID 573: Comment out unreachable code.

__STDC__ is always defined on NetBSD.

More trailing white space.

rcsid facism, WARNSify

Id -> NetBSD

LIBM_SCCS must be defined for the RCS ID to be defined.
This library will soon be used by DJGPP. And since memory is often tight
on DOS machines, there is little benifit to be had from including RCS IDs.

Change int -> int32_t and unsigned int -> u_int32_t to improve portabilty
to machines where int != 32 bits.

Float versions of math functions. From Ian Taylor (ian@cygnus.com), with
minor changes by me.