ulps/subr_mpc.c

   1 #include <mpc.h>
   2 #include "mytypes.h"
   3 #include "subr_mpc.h"
   4
   5 #if MPC_VERSION_MAJOR == 0 && MPC_VERSION_MINOR < 9
   6 double complex
   7 mpc_get_dc(mpc_t op, mpfr_rnd_t rnd)
   8 {
   9         mpfr_t mp_real;
  10         mpfr_t mp_imag;
  11         double_complex z;
  12
  13         mpfr_init2(mp_real, 200);
  14         mpfr_init2(mp_imag, 200);
  15
  16         mpc_real(mp_real, op, rnd);
  17         mpc_imag(mp_imag, op, rnd);
  18
  19         z.parts[0] = mpfr_get_d(mp_real, rnd);
  20         z.parts[1] = mpfr_get_d(mp_imag, rnd);
  21
  22         mpfr_clear(mp_real);
  23         mpfr_clear(mp_imag);
  24
  25         return (z.z);
  26 }
  27
  28 long double complex
  29 mpc_get_ldc(mpc_t op, mpfr_rnd_t rnd)
  30 {
  31         mpfr_t mp_real;
  32         mpfr_t mp_imag;
  33         long_double_complex z;
  34
  35         mpfr_init2(mp_real, 200);
  36         mpfr_init2(mp_imag, 200);
  37
  38         mpc_real(mp_real, op, rnd);
  39         mpc_imag(mp_imag, op, rnd);
  40
  41         z.parts[0] = mpfr_get_ld(mp_real, rnd);
  42         z.parts[1] = mpfr_get_ld(mp_imag, rnd);
  43
  44         mpfr_clear(mp_real);
  45         mpfr_clear(mp_imag);
  46
  47         return (z.z);
  48 }
  49 #endif