diff options
-rw-r--r-- | include/complex.h | 11 | ||||
-rw-r--r-- | src/complex/cimag.c | 3 | ||||
-rw-r--r-- | src/complex/cimagf.c | 3 | ||||
-rw-r--r-- | src/complex/cimagl.c | 3 | ||||
-rw-r--r-- | src/complex/creal.c | 4 | ||||
-rw-r--r-- | src/complex/crealf.c | 4 | ||||
-rw-r--r-- | src/complex/creall.c | 4 | ||||
-rw-r--r-- | src/internal/libm.h | 8 |
8 files changed, 20 insertions, 20 deletions
diff --git a/include/complex.h b/include/complex.h index 8ee70575..90496bd5 100644 --- a/include/complex.h +++ b/include/complex.h @@ -97,6 +97,17 @@ double creal(double complex); float crealf(float complex); long double creall(long double complex); +#define __CREALIMAG(x, t, i) \ + ((union { _Complex t __z; t __xy[2]; }){(_Complex t)(x)}.__xy[i]) + +#define creal(x) __CREALIMAG(x, double, 0) +#define crealf(x) __CREALIMAG(x, float, 0) +#define creall(x) __CREALIMAG(x, long double, 0) + +#define cimag(x) __CREALIMAG(x, double, 1) +#define cimagf(x) __CREALIMAG(x, float, 1) +#define cimagl(x) __CREALIMAG(x, long double, 1) + #ifdef __cplusplus } #endif diff --git a/src/complex/cimag.c b/src/complex/cimag.c index 5e1ad46b..00955641 100644 --- a/src/complex/cimag.c +++ b/src/complex/cimag.c @@ -2,6 +2,5 @@ double (cimag)(double complex z) { - union dcomplex u = {z}; - return u.a[1]; + return cimag(z); } diff --git a/src/complex/cimagf.c b/src/complex/cimagf.c index 99fffc58..f7bcd76e 100644 --- a/src/complex/cimagf.c +++ b/src/complex/cimagf.c @@ -2,6 +2,5 @@ float (cimagf)(float complex z) { - union fcomplex u = {z}; - return u.a[1]; + return cimagf(z); } diff --git a/src/complex/cimagl.c b/src/complex/cimagl.c index d9a0780c..9ec24eee 100644 --- a/src/complex/cimagl.c +++ b/src/complex/cimagl.c @@ -2,6 +2,5 @@ long double (cimagl)(long double complex z) { - union lcomplex u = {z}; - return u.a[1]; + return cimagl(z); } diff --git a/src/complex/creal.c b/src/complex/creal.c index 2bb91812..f6703040 100644 --- a/src/complex/creal.c +++ b/src/complex/creal.c @@ -1,6 +1,6 @@ #include <complex.h> -double creal(double complex z) +double (creal)(double complex z) { - return z; + return creal(z); } diff --git a/src/complex/crealf.c b/src/complex/crealf.c index 078232f0..5dc3ff1d 100644 --- a/src/complex/crealf.c +++ b/src/complex/crealf.c @@ -1,6 +1,6 @@ #include <complex.h> -float crealf(float complex z) +float (crealf)(float complex z) { - return z; + return crealf(z); } diff --git a/src/complex/creall.c b/src/complex/creall.c index 56e64097..fd9dc347 100644 --- a/src/complex/creall.c +++ b/src/complex/creall.c @@ -1,6 +1,6 @@ #include <complex.h> -long double creall(long double complex z) +long double (creall)(long double complex z) { - return z; + return creall(z); } diff --git a/src/internal/libm.h b/src/internal/libm.h index 67c42b98..8c5474a8 100644 --- a/src/internal/libm.h +++ b/src/internal/libm.h @@ -173,14 +173,6 @@ union lcomplex { long double a[2]; }; -// FIXME: move to complex.h ? -#define creal(z) ((double)(z)) -#define crealf(z) ((float)(z)) -#define creall(z) ((long double)(z)) -#define cimag(z) ((union dcomplex){(z)}.a[1]) -#define cimagf(z) ((union fcomplex){(z)}.a[1]) -#define cimagl(z) ((union lcomplex){(z)}.a[1]) - /* x + y*I is not supported properly by gcc */ #define cpack(x,y) ((union dcomplex){.a={(x),(y)}}.z) #define cpackf(x,y) ((union fcomplex){.a={(x),(y)}}.z) |