d/freebsd-skq
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Implement inline functions to give the complex result x+I*y from float

Browse Source 
or double args x and y.  x+I*y cannot be used directly yet due to compiler
bugs.

Submitted by:	Steve Kargl <sgk@troutmask.apl.washington.edu>
This commit is contained in:
bde 2005-10-29 17:14:11 +00:00
parent bbfb40721e
commit 26610cfe9b
1 changed files with 42 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

42
lib/msun/src/math_private.h
View File

@ -154,6 +154,48 @@ do { \
(d) = sf_u.value; \
} while (0)
#ifdef _COMPLEX_H
/*
* Inline functions that can be used to construct complex values.
*
* The C99 standard intends x+I*y to be used for this, but x+I*y is
* currently unusable in general since gcc introduces many overflow,
* underflow, sign and efficiency bugs by rewriting I*y as
* (0.0+I)*(y+0.0*I) and laboriously computing the full complex product.
* In particular, I*Inf is corrupted to NaN+I*Inf, and I*-0 is corrupted
* to -0.0+I*0.0.
*/
static __inline float complex
cpackf(float x, float y)
{
float complex z;
__real__ z = x;
__imag__ z = y;
return (z);
}
static __inline double complex
cpack(double x, double y)
{
double complex z;
__real__ z = x;
__imag__ z = y;
return (z);
}
static __inline long double complex
cpackl(long double x, long double y)
{
long double complex z;
__real__ z = x;
__imag__ z = y;
return (z);
}
#endif /* _COMPLEX_H */
/*
* ieee style elementary functions
*