number 0.9.7 → 0.9.9

data/ext/{compare.c → number/compare.c}

@@ -657,7 +657,7 @@ int complex_contain_zero_p (Complex* complex)
 
 int complex_intersect_p (Complex* a, Complex* b)
 {
-    return (complex_re_intersect_p(a, b) && complex_im_intersect_p(a, b)) || complex_contain_ll_p(a, b) || complex_contain_lu_p(a, b) || complex_contain_ul_p(a, b) || complex_contain_uu_p(a, b) || (complex_re_enclose_p(a, b) && complex_im_enclose_p(b, a) || (complex_im_enclose_p(a, b) && complex_re_enclose_p(b, a)));
+    return ( complex_re_intersect_p(a, b) && complex_im_intersect_p(a, b)) || complex_contain_ll_p(a, b) || complex_contain_lu_p(a, b) || complex_contain_ul_p(a, b) || complex_contain_uu_p(a, b) || ((complex_re_enclose_p(a, b) && complex_im_enclose_p(b, a)) || (complex_im_enclose_p(a, b) && complex_re_enclose_p(b, a)));
 }
 
 int complex_disjoint_p (Complex* a, Complex* b)

data/ext/{complex.c → number/complex.c}

@@ -111,7 +111,6 @@ void mpc_set_r_r (mpc_t num, Real* re, Real* im)
  */
 void mpfr_set_r (mpfr_t num, Real* real)
 {
-    char sign;
     char buf[10000];
     char real_buf[10000];
     
@@ -155,7 +154,7 @@ void complex_free (Complex* complex)
  */
 long complex_strlen (Complex* complex)
 {
-    return interval_strlen(complex->re) + interval_strlen(complex->im) + 10;
+    return (interval_strlen(complex->re) + interval_strlen(complex->im) + 10);
 }
 
 /*
@@ -189,7 +188,7 @@ void complex_print (char* buf, Complex* complex)
 void complex_print_int (char* buf, Complex* complex)
 {
     Real* rounded;
-    rounded = (Real*)real_round_origin(complex->re->n);
+    rounded = real_round_origin(complex->re->n);
     
     if (!REAL(rounded))
     {
@@ -238,7 +237,7 @@ void complex_print_decimal (char* buf, Complex* complex)
 Complex* complex_new (Interval* re, Interval* im)
 {
     Complex* complex;
-    complex = (Complex*)malloc(sizeof(Complex));
+    complex = malloc(sizeof(Complex));
     
     if (complex == NULL)
     {
@@ -246,8 +245,8 @@ Complex* complex_new (Interval* re, Interval* im)
     }
     
     memset(complex, 0, sizeof(Complex));
-    complex->re = (Interval*)re;
-    complex->im = (Interval*)im;
+    complex->re = re;
+    complex->im = im;
     
     return complex;
 }
@@ -259,7 +258,13 @@ Complex* complex_new (Interval* re, Interval* im)
  */
 Complex* complex_from_real (Real* real)
 {
-    return complex_new((Interval*)interval_from_real(real), (Interval*)interval_zero());
+    Interval* re;
+    Interval* im;
+    
+    re = interval_from_real(real);
+    im = interval_zero();
+    
+    return complex_new(re, im);
 }
 
 /*
@@ -271,8 +276,8 @@ Complex* complex_from_int_str (char* buf)
     Real* real;
     Complex* result;
     
-    real = (Real*)real_from_str(buf, strlen(buf) - (*buf == '-'));
-    result = (Complex*)complex_from_real(real);
+    real = real_from_str(buf, strlen(buf) - (*buf == '-'));
+    result = complex_from_real(real);
     
     real_free(real);
     
@@ -290,8 +295,8 @@ Complex* complex_from_float_str (char* buf)
     
     mpfr_set_str(tmp_complex_from_float_str, buf, 10, MPFR_RNDN);
     
-    real = (Real*)real_from_mpfr(tmp_complex_from_float_str, ROUND_NEAREST);
-    result = (Complex*)complex_from_real(real);
+    real = real_from_mpfr(tmp_complex_from_float_str, ROUND_NEAREST);
+    result = complex_from_real(real);
     
     real_free(real);
     
@@ -303,7 +308,7 @@ Complex* complex_from_float_str (char* buf)
  */
 Complex* complex_nan ()
 {
-    return complex_new((Interval*)interval_nan(), (Interval*)interval_zero());
+    return complex_new(interval_nan(), interval_zero());
 }
 
 /*
@@ -311,7 +316,7 @@ Complex* complex_nan ()
  */
 Complex* complex_pos_inf ()
 {
-    return complex_new((Interval*)interval_pos_inf(), (Interval*)interval_zero());
+    return complex_new(interval_pos_inf(), interval_zero());
 }
 
 /*
@@ -319,7 +324,7 @@ Complex* complex_pos_inf ()
  */
 Complex* complex_neg_inf ()
 {
-    return complex_new((Interval*)interval_neg_inf(), (Interval*)interval_zero());
+    return complex_new(interval_neg_inf(), interval_zero());
 }
 
 /*
@@ -327,7 +332,7 @@ Complex* complex_neg_inf ()
  */
 Complex* complex_zero ()
 {
-    return complex_new((Interval*)interval_zero(), (Interval*)interval_zero());
+    return complex_new(interval_zero(), interval_zero());
 }
 
 /*
@@ -335,7 +340,7 @@ Complex* complex_zero ()
  */
 Complex* complex_one ()
 {
-    return complex_new((Interval*)interval_one(), (Interval*)interval_zero());
+    return complex_new(interval_one(), interval_zero());
 }
 
 /*
@@ -344,7 +349,7 @@ Complex* complex_one ()
  */
 Complex* complex_pi ()
 {
-    return complex_new((Interval*)interval_pi(), (Interval*)interval_zero());
+    return complex_new(interval_pi(), interval_zero());
 }
 
 /*
@@ -353,7 +358,7 @@ Complex* complex_pi ()
  */
 Complex* complex_e ()
 {
-    return complex_new((Interval*)interval_e(), (Interval*)interval_zero());
+    return complex_new(interval_e(), interval_zero());
 }
 
 /*
@@ -361,7 +366,7 @@ Complex* complex_e ()
  */
 Complex* complex_i ()
 {
-    return complex_new((Interval*)interval_zero(), (Interval*)interval_one());
+    return complex_new(interval_zero(), interval_one());
 }
 
 /*
@@ -370,7 +375,7 @@ Complex* complex_i ()
  */
 Complex* complex_dup (Complex* complex)
 {
-    return complex_new((Interval*)interval_dup(complex->re), (Interval*)interval_dup(complex->im));
+    return complex_new(interval_dup(complex->re), interval_dup(complex->im));
 }
 
 /*
@@ -381,8 +386,8 @@ Complex* complex_add (Complex* a, Complex* b)
     Interval* re;
     Interval* im;
     
-    re = (Interval*)interval_add(a->re, b->re);
-    im = (Interval*)interval_add(a->im, b->im);
+    re = interval_add(a->re, b->re);
+    im = interval_add(a->im, b->im);
     
     return complex_new(re, im);
 }
@@ -395,8 +400,8 @@ Complex* complex_subtract (Complex* a, Complex* b)
     Interval* re;
     Interval* im;
     
-    re = (Interval*)interval_subtract(a->re, b->re);
-    im = (Interval*)interval_subtract(a->im, b->im);
+    re = interval_subtract(a->re, b->re);
+    im = interval_subtract(a->im, b->im);
     
     return complex_new(re, im);
 }
@@ -413,20 +418,20 @@ Complex* complex_multiply (Complex* a, Complex* b)
     Interval* re;
     Interval* im;
     
-    re_re = (Interval*)interval_multiply(a->re, b->re);
-    re_im = (Interval*)interval_multiply(a->re, b->im);
-    im_re = (Interval*)interval_multiply(a->im, b->re);
-    im_im = (Interval*)interval_multiply(a->im, b->im);
+    re_re = interval_multiply(a->re, b->re);
+    re_im = interval_multiply(a->re, b->im);
+    im_re = interval_multiply(a->im, b->re);
+    im_im = interval_multiply(a->im, b->im);
     
-    re = (Interval*)interval_subtract(re_re, im_im);
-    im = (Interval*)interval_add(re_im, im_re);
+    re = interval_subtract(re_re, im_im);
+    im = interval_add(re_im, im_re);
     
     interval_free(re_re);
     interval_free(re_im);
     interval_free(im_re);
     interval_free(im_im);
     
-    return (Complex*)complex_new(re, im);
+    return complex_new(re, im);
 }
 
 /*
@@ -446,9 +451,9 @@ Complex* complex_divide (Complex* a, Complex* b)
     Interval* re;
     Interval* im;
     
-    denominator_re = (Interval*)interval_multiply(b->re, b->re);
-    denominator_im = (Interval*)interval_multiply(b->im, b->im);
-    denominator = (Interval*)interval_add(denominator_re, denominator_im);
+    denominator_re = interval_multiply(b->re, b->re);
+    denominator_im = interval_multiply(b->im, b->im);
+    denominator = interval_add(denominator_re, denominator_im);
     
     if (real_zero_p(denominator->l) || real_zero_p(denominator->u) || interval_span_zero_p(denominator))
     {
@@ -466,21 +471,21 @@ Complex* complex_divide (Complex* a, Complex* b)
         mpc_real(tmp_divide_re, tmp_divide_quotient, MPFR_RNDN);
         mpc_imag(tmp_divide_im, tmp_divide_quotient, MPFR_RNDN);
         
-        re = (Interval*)interval_new((Real*)real_neg_inf(), true, (Real*)real_from_mpfr(tmp_divide_re, ROUND_NEAREST), true, (Real*)real_pos_inf());
-        im = (Interval*)interval_new((Real*)real_neg_inf(), true, (Real*)real_from_mpfr(tmp_divide_im, ROUND_NEAREST), true, (Real*)real_pos_inf());
+        re = interval_new(real_neg_inf(), true, real_from_mpfr(tmp_divide_re, ROUND_NEAREST), true, real_pos_inf());
+        im = interval_new(real_neg_inf(), true, real_from_mpfr(tmp_divide_im, ROUND_NEAREST), true, real_pos_inf());
     }
     else
     {
-        re_re = (Interval*)interval_multiply(a->re, b->re);
-        re_im = (Interval*)interval_multiply(a->re, b->im);
-        im_re = (Interval*)interval_multiply(a->im, b->re);
-        im_im = (Interval*)interval_multiply(a->im, b->im);
+        re_re = interval_multiply(a->re, b->re);
+        re_im = interval_multiply(a->re, b->im);
+        im_re = interval_multiply(a->im, b->re);
+        im_im = interval_multiply(a->im, b->im);
         
-        numerator_re = (Interval*)interval_add(re_re, im_im);
-        numerator_im = (Interval*)interval_subtract(im_re, re_im);
+        numerator_re = interval_add(re_re, im_im);
+        numerator_im = interval_subtract(im_re, re_im);
         
-        re = (Interval*)interval_divide(numerator_re, denominator);
-        im = (Interval*)interval_divide(numerator_im, denominator);
+        re = interval_divide(numerator_re, denominator);
+        im = interval_divide(numerator_im, denominator);
         
         interval_free(re_re);
         interval_free(re_im);
@@ -494,7 +499,7 @@ Complex* complex_divide (Complex* a, Complex* b)
     interval_free(denominator_im);
     interval_free(denominator);
     
-    return (Complex*)complex_new(re, im);
+    return complex_new(re, im);
 }
 
 /*
@@ -505,10 +510,10 @@ Complex* complex_modulo (Complex* a, Complex* b)
     Interval* re;
     Interval* im;
     
-    re = (Interval*)interval_zero();
-    im = (Interval*)interval_zero();
+    re = interval_zero();
+    im = interval_zero();
     
-    return (Complex*)complex_new(re, im);
+    return complex_new(re, im);
 }
 
 /*
@@ -553,7 +558,7 @@ Complex* complex_abs (Complex* complex)
     ul_mpc_ternary = mpc_abs(tmp_abs_ul_re, tmp_abs_ul, MPC_RNDUD);
     uu_mpc_ternary = mpc_abs(tmp_abs_uu_re, tmp_abs_uu, MPC_RNDUU);
     
-    bounds = (Bounds*)bounds_new();
+    bounds = bounds_new();
     bounds_update_re(bounds, tmp_abs_ll_re, complex->re->L, MPC_INEX_RE(ll_mpc_ternary) == 0);
     bounds_update_re(bounds, tmp_abs_lu_re, complex->re->L, MPC_INEX_RE(lu_mpc_ternary) == 0);
     bounds_update_re(bounds, tmp_abs_ul_re, complex->re->R, MPC_INEX_RE(ul_mpc_ternary) == 0);
@@ -561,8 +566,8 @@ Complex* complex_abs (Complex* complex)
     
     if (interval_span_zero_p(complex->im))
     {
-        l_abs = (Real*)real_abs(complex->re->l);
-        u_abs = (Real*)real_abs(complex->re->u);
+        l_abs = real_abs(complex->re->l);
+        u_abs = real_abs(complex->re->u);
         
         mpfr_set_r(tmp_abs_re_abs, real_lt_p(l_abs, u_abs) ? l_abs : u_abs);
         bounds_update_re(bounds, tmp_abs_re_abs, true, true);
@@ -573,8 +578,8 @@ Complex* complex_abs (Complex* complex)
     
     if (interval_span_zero_p(complex->re))
     {
-        l_abs = (Real*)real_abs(complex->im->l);
-        u_abs = (Real*)real_abs(complex->im->u);
+        l_abs = real_abs(complex->im->l);
+        u_abs = real_abs(complex->im->u);
         
         mpfr_set_r(tmp_abs_im_abs, real_lt_p(l_abs, u_abs) ? l_abs : u_abs);
         bounds_update_re(bounds, tmp_abs_im_abs, true, true);
@@ -589,11 +594,11 @@ Complex* complex_abs (Complex* complex)
         bounds_update_re(bounds, tmp_abs_zero, true, true);
     }
     
-    re = (Interval*)interval_new(real_from_mpfr(bounds->re_min, ROUND_DOWN), bounds->re_min_closed, real_from_mpfr(tmp_abs_n_re, ROUND_NEAREST), bounds->re_max_closed, real_from_mpfr(bounds->re_max, bounds->re_max_exact ? ROUND_UP : ROUND_UP1));
+    re = interval_new(real_from_mpfr(bounds->re_min, ROUND_DOWN), bounds->re_min_closed, real_from_mpfr(tmp_abs_n_re, ROUND_NEAREST), bounds->re_max_closed, real_from_mpfr(bounds->re_max, bounds->re_max_exact ? ROUND_UP : ROUND_UP1));
     
     bounds_free(bounds);
     
-    return (Complex*)complex_new(re, (Interval*)interval_zero());
+    return complex_new(re, interval_zero());
 }
 
 /*
@@ -607,10 +612,10 @@ Complex* complex_area (Complex* complex)
     Real* area;
     Complex* result;
     
-    re_width = (Real*)real_subtract(complex->re->u, complex->re->l);
-    im_width = (Real*)real_subtract(complex->im->u, complex->im->l);
-    area = (Real*)real_multiply(re_width, im_width);
-    result = (Complex*)complex_from_real(area);
+    re_width = real_subtract(complex->re->u, complex->re->l);
+    im_width = real_subtract(complex->im->u, complex->im->l);
+    area = real_multiply(re_width, im_width);
+    result = complex_from_real(area);
     
     real_free(re_width);
     real_free(im_width);
@@ -632,20 +637,19 @@ Complex* complex_diagonal (Complex* complex)
     Real* im_width2;
     Real* diag2;
     Interval* re;
-    Complex* result;
     
-    re_width = (Real*)real_subtract(complex->re->u, complex->re->l);
-    im_width = (Real*)real_subtract(complex->im->u, complex->im->l);
+    re_width = real_subtract(complex->re->u, complex->re->l);
+    im_width = real_subtract(complex->im->u, complex->im->l);
     
-    re_width2 = (Real*)real_multiply(re_width, re_width);
-    im_width2 = (Real*)real_multiply(im_width, im_width);
+    re_width2 = real_multiply(re_width, re_width);
+    im_width2 = real_multiply(im_width, im_width);
     
-    diag2 = (Real*)real_add(re_width2, im_width2);
+    diag2 = real_add(re_width2, im_width2);
     
     mpfr_set_r(tmp_diagonal, diag2);
     mpfr_ternary = mpfr_sqrt(tmp_diagonal, tmp_diagonal, MPFR_RNDN);
     
-    re = (Interval*)interval_new((Real*)real_from_mpfr(tmp_diagonal, ROUND_DOWN), true, (Real*)real_from_mpfr(tmp_diagonal, ROUND_NEAREST), true, (Real*)real_from_mpfr(tmp_diagonal, (mpfr_ternary == 0) ? ROUND_UP : ROUND_UP1));
+    re = interval_new(real_from_mpfr(tmp_diagonal, ROUND_DOWN), true, real_from_mpfr(tmp_diagonal, ROUND_NEAREST), true, real_from_mpfr(tmp_diagonal, (mpfr_ternary == 0) ? ROUND_UP : ROUND_UP1));
     
     real_free(re_width);
     real_free(im_width);
@@ -653,7 +657,7 @@ Complex* complex_diagonal (Complex* complex)
     real_free(im_width2);
     real_free(diag2);
     
-    return (Complex*)complex_new(re, (Interval*)interval_zero());
+    return complex_new(re, interval_zero());
 }
 
 /*
@@ -662,7 +666,7 @@ Complex* complex_diagonal (Complex* complex)
  */
 Complex* complex_conjugate (Complex* complex)
 {
-    return (Complex*)complex_new((Interval*)interval_dup(complex->re), (Interval*)interval_negate(complex->im));
+    return complex_new(interval_dup(complex->re), interval_negate(complex->im));
 }
 
 /*
@@ -670,7 +674,7 @@ Complex* complex_conjugate (Complex* complex)
  */
 Complex* complex_negate (Complex* complex)
 {
-    return (Complex*)complex_new((Interval*)interval_negate(complex->re), (Interval*)interval_negate(complex->im));
+    return complex_new(interval_negate(complex->re), interval_negate(complex->im));
 }
 
 /*
@@ -679,7 +683,7 @@ Complex* complex_negate (Complex* complex)
  */
 Complex* complex_reflect (Complex* complex)
 {
-    return (Complex*)complex_new((Interval*)interval_dup(complex->im), (Interval*)interval_dup(complex->re));
+    return complex_new(interval_dup(complex->im), interval_dup(complex->re));
 }
 
 /*
@@ -720,7 +724,7 @@ Complex* complex_sqrt (Complex* complex)
     ul_mpc_ternary = mpc_sqrt(tmp_sqrt_ul, tmp_sqrt_ul, MPC_RNDUD);
     uu_mpc_ternary = mpc_sqrt(tmp_sqrt_uu, tmp_sqrt_uu, MPC_RNDUU);
     
-    bounds = (Bounds*)bounds_new();
+    bounds = bounds_new();
     bounds_update(bounds, tmp_sqrt_ll, complex->re->L, complex->im->L, ll_mpc_ternary);
     bounds_update(bounds, tmp_sqrt_lu, complex->re->L, complex->im->R, lu_mpc_ternary);
     bounds_update(bounds, tmp_sqrt_ul, complex->re->R, complex->im->L, ul_mpc_ternary);
@@ -734,10 +738,10 @@ Complex* complex_sqrt (Complex* complex)
         bounds_update(bounds, tmp_sqrt_re_sqrt, complex->re->L, true, re_sqrt_mpc_ternary);
     }
     
-    re = (Interval*)interval_new((Real*)real_from_mpfr(bounds->re_min, ROUND_DOWN), bounds->re_min_closed, (Real*)real_from_mpfr(tmp_sqrt_n_re, ROUND_NEAREST), bounds->re_max_closed, (Real*)real_from_mpfr(bounds->re_max, bounds->re_max_exact ? ROUND_UP : ROUND_UP1));
-    im = (Interval*)interval_new((Real*)real_from_mpfr(bounds->im_min, ROUND_DOWN), bounds->im_min_closed, (Real*)real_from_mpfr(tmp_sqrt_n_im, ROUND_NEAREST), bounds->im_max_closed, (Real*)real_from_mpfr(bounds->im_max, bounds->im_max_exact ? ROUND_UP : ROUND_UP1));
+    re = interval_new(real_from_mpfr(bounds->re_min, ROUND_DOWN), bounds->re_min_closed, real_from_mpfr(tmp_sqrt_n_re, ROUND_NEAREST), bounds->re_max_closed, real_from_mpfr(bounds->re_max, bounds->re_max_exact ? ROUND_UP : ROUND_UP1));
+    im = interval_new(real_from_mpfr(bounds->im_min, ROUND_DOWN), bounds->im_min_closed, real_from_mpfr(tmp_sqrt_n_im, ROUND_NEAREST), bounds->im_max_closed, real_from_mpfr(bounds->im_max, bounds->im_max_exact ? ROUND_UP : ROUND_UP1));
     
     bounds_free(bounds);
     
-    return (Complex*)complex_new(re, im);
+    return complex_new(re, im);
 }

data/ext/number/decompose.c

@@ -0,0 +1,231 @@
+#include "number.h"
+
+Complex* complex_im (Complex* complex)
+{
+    return complex_new(interval_dup(complex->im), interval_zero());
+}
+
+Complex* complex_im_l (Complex* complex)
+{
+    return complex_from_real(complex->im->l);
+}
+
+Complex* complex_im_n (Complex* complex)
+{
+    return complex_from_real(complex->im->n);
+}
+
+Complex* complex_im_u (Complex* complex)
+{
+    return complex_from_real(complex->im->u);
+}
+
+Complex* complex_im_magnitude (Complex* complex)
+{
+    Real* l_abs;
+    Real* u_abs;
+    Complex* result;
+    
+    l_abs = real_abs(complex->im->l);
+    u_abs = real_abs(complex->im->u);
+    
+    result = complex_from_real(real_gt_p(l_abs, u_abs) ? l_abs : u_abs);
+    
+    real_free(l_abs);
+    real_free(u_abs);
+    
+    return result;
+}
+
+Complex* complex_im_mignitude (Complex* complex)
+{
+    Real* l_abs;
+    Real* u_abs;
+    Real* zero;
+    Complex* result;
+    
+    l_abs = real_abs(complex->im->l);
+    u_abs = real_abs(complex->im->u);
+    zero = real_zero();
+    
+    result = complex_from_real(interval_span_zero_p(complex->im) ? zero : (real_lt_p(l_abs, u_abs) ? l_abs : u_abs));
+    
+    real_free(l_abs);
+    real_free(u_abs);
+    real_free(zero);
+    
+    return result;
+}
+
+Complex* complex_im_width (Complex* complex)
+{
+    Real* width;
+    Complex* result;
+    
+    width = real_subtract(complex->im->u, complex->im->l);
+    result = complex_from_real(width);
+    
+    real_free(width);
+    
+    return result;
+}
+
+Complex* complex_re (Complex* complex)
+{
+    return complex_new(interval_dup(complex->re), interval_zero());
+}
+
+Complex* complex_re_l (Complex* complex)
+{
+    return complex_from_real(complex->re->l);
+}
+
+Complex* complex_re_n (Complex* complex)
+{
+    return complex_from_real(complex->re->n);
+}
+
+Complex* complex_re_u (Complex* complex)
+{
+    return complex_from_real(complex->re->u);
+}
+
+Complex* complex_re_magnitude (Complex* complex)
+{
+    Real* l_abs;
+    Real* u_abs;
+    Complex* result;
+    
+    l_abs = real_abs(complex->re->l);
+    u_abs = real_abs(complex->re->u);
+    
+    result = complex_from_real(real_gt_p(l_abs, u_abs) ? l_abs : u_abs);
+    
+    real_free(l_abs);
+    real_free(u_abs);
+    
+    return result;
+}
+
+Complex* complex_re_mignitude (Complex* complex)
+{
+    Real* l_abs;
+    Real* u_abs;
+    Real* zero;
+    Complex* result;
+    
+    l_abs = real_abs(complex->re->l);
+    u_abs = real_abs(complex->re->u);
+    zero = real_zero();
+    
+    result = complex_from_real(interval_span_zero_p(complex->re) ? zero : (real_lt_p(l_abs, u_abs) ? l_abs : u_abs));
+    
+    real_free(l_abs);
+    real_free(u_abs);
+    real_free(zero);
+    
+    return result;
+}
+
+Complex* complex_re_width (Complex* complex)
+{
+    Real* width;
+    Complex* result;
+    
+    width = real_subtract(complex->re->u, complex->re->l);
+    result = complex_from_real(width);
+    
+    real_free(width);
+    
+    return result;
+}
+
+Complex* complex_ll (Complex* complex)
+{
+    return complex_new(interval_from_real(complex->re->l), interval_from_real(complex->im->l));
+}
+
+Complex* complex_lu (Complex* complex)
+{
+    return complex_new(interval_from_real(complex->re->l), interval_from_real(complex->im->u));
+}
+
+Complex* complex_n (Complex* complex)
+{
+    return complex_new(interval_from_real(complex->re->n), interval_from_real(complex->im->n));
+}
+
+Complex* complex_ul (Complex* complex)
+{
+    return complex_new(interval_from_real(complex->re->u), interval_from_real(complex->im->l));
+}
+
+Complex* complex_uu (Complex* complex)
+{
+    return complex_new(interval_from_real(complex->re->u), interval_from_real(complex->im->u));
+}
+
+Complex* complex_magnitude (Complex* complex)
+{
+    Real* re_l_abs;
+    Real* re_u_abs;
+    Real* im_l_abs;
+    Real* im_u_abs;
+    Complex* result;
+    
+    re_l_abs = real_abs(complex->re->l);
+    re_u_abs = real_abs(complex->re->u);
+    im_l_abs = real_abs(complex->im->l);
+    im_u_abs = real_abs(complex->im->u);
+    
+    result = complex_new(interval_from_real(real_gt_p(re_l_abs, re_u_abs) ? re_l_abs : re_u_abs), interval_from_real(real_gt_p(im_l_abs, im_u_abs) ? im_l_abs : im_u_abs));
+    
+    real_free(re_l_abs);
+    real_free(re_u_abs);
+    real_free(im_l_abs);
+    real_free(im_u_abs);
+    
+    return result;
+}
+
+Complex* complex_mignitude (Complex* complex)
+{
+    Real* re_l_abs;
+    Real* re_u_abs;
+    Real* im_l_abs;
+    Real* im_u_abs;
+    Real* zero;
+    Complex* result;
+    
+    re_l_abs = real_abs(complex->re->l);
+    re_u_abs = real_abs(complex->re->u);
+    im_l_abs = real_abs(complex->im->l);
+    im_u_abs = real_abs(complex->im->u);
+    
+    zero = real_zero();
+    result = complex_new(interval_from_real(interval_span_zero_p(complex->re) ? zero : (real_lt_p(re_l_abs, re_u_abs) ? re_l_abs : re_u_abs)), interval_from_real(interval_span_zero_p(complex->im) ? zero : (real_lt_p(im_l_abs, im_u_abs) ? im_l_abs : im_u_abs)));
+    
+    real_free(re_l_abs);
+    real_free(re_u_abs);
+    real_free(im_l_abs);
+    real_free(im_u_abs);
+    real_free(zero);
+    
+    return result;
+}
+
+Complex* complex_width (Complex* complex)
+{
+    Real* re_width;
+    Real* im_width;
+    Complex* result;
+    
+    re_width = real_subtract(complex->re->u, complex->re->l);
+    im_width = real_subtract(complex->im->u, complex->im->l);
+    result = complex_new(interval_from_real(re_width), interval_from_real(im_width));
+    
+    real_free(re_width);
+    real_free(im_width);
+    
+    return result;
+}