#include <stdio.h>
#include <ruby.h>
#include <mpfr.h>
#include <mpfi.h>
#include <mpfi_io.h>
#include "ruby_mpfi.h"

VALUE r_mpfi_complex;

typedef struct __MPFIComplex{
  MPFI *re;
  MPFI *im;
} MPFIComplex;

/* Although x->re and x->im can have different precision from each other for MPFIComplex x, */
/* we implement mpfi_complex_get_prec assuming that they have same presition. */
#define mpfi_complex_get_prec(x) mpfi_get_prec(x->re)

void mpfi_complex_init(MPFIComplex *x);
void mpfi_complex_clear(MPFIComplex *x);
void mpfi_complex_set_zeros(MPFIComplex *x);
void mpfi_complex_set_real_part(MPFIComplex *x, MPFI *a);
void mpfi_complex_set_imaginary_part(MPFIComplex *x, MPFI *a);
void mpfi_complex_set(MPFIComplex *new, MPFIComplex *x);
void mpfi_complex_conjugate(MPFIComplex *new, MPFIComplex *x);
void mpfi_complex_add(MPFIComplex *new, MPFIComplex *x, MPFIComplex *y);
void mpfi_complex_sub(MPFIComplex *new, MPFIComplex *x, MPFIComplex *y);
void mpfi_complex_mul(MPFIComplex *new, MPFIComplex *x, MPFIComplex *y);
void mpfi_complex_div(MPFIComplex *new, MPFIComplex *x, MPFIComplex *y);

/* void mpfi_complex_mul_real(MPFIComplex *new, MPFIComplex *x, MPFI *y); */
/* void mpfi_complex_mul_pure_imaginary(MPFIComplex *new, MPFIComplex *x, MPFI *y); */
/* void mpfi_complex_div_real(MPFIComplex *new, MPFIComplex *x, MPFI *y); */
/* void mpfi_complex_div_pure_imaginary(MPFIComplex *new, MPFIComplex *x, MPFI *y); */
/* void mpfi_complex_abs(MPFI *new, MPFIComplex *x); */