complex.h

Go to the documentation of this file.

// mightymandel -- GPU-based Mandelbrot Set explorer
// Copyright (C) 2012,2013,2014 Claude Heiland-Allen
// License GPL3+ http://www.gnu.org/licenses/gpl.html

#ifndef COMPLEX_H
#define COMPLEX_H 1

#include <mpfr.h>

typedef mpfr_t R;

typedef struct {
  R x; 
  R y; 
} C;

#define c_set_prec(o, p) do{\
  mpfr_set_prec(o.x, p);\
  mpfr_set_prec(o.y, p);\
}while(0)

#define c_set(o, l) do{\
  mpfr_set(o.x, l.x, MPFR_RNDN);\
  mpfr_set(o.y, l.y, MPFR_RNDN);\
}while(0)

#define c_mag2(o, l, t1, t2) do{\
  mpfr_sqr(t1, l.x, MPFR_RNDN);\
  mpfr_sqr(t2, l.y, MPFR_RNDN);\
  mpfr_add(o, t1, t2, MPFR_RNDN);\
}while(0)

#define c_add(o, l, r) do{\
  mpfr_add(o.x, l.x, r.x, MPFR_RNDN);\
  mpfr_add(o.y, l.y, r.y, MPFR_RNDN);\
}while(0)

#define c_sqr(o, l, t1, t2, t3) do{\
  mpfr_sqr(t1, l.x, MPFR_RNDN);\
  mpfr_sqr(t2, l.y, MPFR_RNDN);\
  mpfr_mul(t3, l.x, l.y, MPFR_RNDN);\
  mpfr_sub(o.x, t1, t2, MPFR_RNDN);\
  mpfr_mul_2ui(o.y, t3, 1, MPFR_RNDN);\
}while(0)

#define c_mul(o, l, r, t1, t2, t3, t4) do{\
  mpfr_mul(t1, l.x, r.x, MPFR_RNDN);\
  mpfr_mul(t2, l.y, r.y, MPFR_RNDN);\
  mpfr_mul(t3, l.x, r.y, MPFR_RNDN);\
  mpfr_mul(t4, l.y, r.x, MPFR_RNDN);\
  mpfr_sub(o.x, t1, t2, MPFR_RNDN);\
  mpfr_add(o.y, t3, t4, MPFR_RNDN);\
}while(0)

#define c_mul_2ui(o, l, r) do{\
  mpfr_mul_2ui(o.x, l.x, r, MPFR_RNDN);\
  mpfr_mul_2ui(o.y, l.y, r, MPFR_RNDN);\
}while(0)

#endif