number 0.9.0

data/ext/number.c

@@ -0,0 +1,1651 @@
+#include "number.h"
+
+long digs;
+VALUE rb_cNumber;
+
+VALUE rb_Integer_to_number (VALUE integer)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_from_int_str(RSTRING_PTR(rb_funcall(integer, rb_intern("to_s"), 0)));
+    
+    return result;
+}
+
+VALUE rb_Float_to_number (VALUE flo)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_from_float_str(RSTRING_PTR(rb_funcall(rb_str_new2("%.15e"), rb_intern("%"), 1, flo)));
+    
+    return result;
+}
+
+VALUE rb_Complex_to_number (VALUE complex)
+{
+    Complex* re;
+    Complex* im;
+    
+    re = (Complex*)DATA_PTR(rb_convert_type(rb_funcall(complex, rb_intern("real"), 0), T_DATA, "Number", "to_number"));
+    im = (Complex*)DATA_PTR(rb_convert_type(rb_funcall(complex, rb_intern("imag"), 0), T_DATA, "Number", "to_number"));
+    
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_new((Interval*)interval_dup(re->re), (Interval*)interval_dup(im->re));
+    
+    return result;
+}
+
+VALUE rb_f_Number (int argc, VALUE* argv, VALUE obj)
+{
+    int L;
+    int R;
+    Complex* tmp1;
+    Complex* tmp2;
+    Complex* tmp3;
+    VALUE arg1;
+    VALUE arg2;
+    VALUE arg3;
+    VALUE arg4;
+    VALUE arg5;
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    
+    switch (rb_scan_args(argc, argv, "14", &arg1, &arg2, &arg3, &arg4, &arg5))
+    {
+        case 1:
+            DATA_PTR(result) = (Complex*)complex_dup(DATA_PTR(rb_convert_type(arg1, T_DATA, "Number", "to_number")));
+            break;
+        
+        case 2:
+            tmp1 = (Complex*)DATA_PTR(rb_convert_type(arg1, T_DATA, "Number", "to_number"));
+            tmp2 = (Complex*)DATA_PTR(rb_convert_type(arg2, T_DATA, "Number", "to_number"));
+            
+            DATA_PTR(result) = (Complex*)complex_new((Interval*)interval_dup(tmp1->re), (Interval*)interval_dup(tmp2->re));
+            break;
+        
+        default:
+            L = (argc == 3) ? true : RTEST(arg4);
+            R = (argc == 5) ? RTEST(arg5) : L;
+            
+            tmp1 = (Complex*)DATA_PTR(rb_convert_type(arg1, T_DATA, "Number", "to_number"));
+            tmp2 = (Complex*)DATA_PTR(rb_convert_type(arg2, T_DATA, "Number", "to_number"));
+            tmp3 = (Complex*)DATA_PTR(rb_convert_type(arg3, T_DATA, "Number", "to_number"));
+            
+            DATA_PTR(result) = (Complex*)complex_new((Interval*)interval_new((Real*)real_dup(tmp1->re->n), L, (Real*)real_dup(tmp2->re->n), R, (Real*)real_dup(tmp3->re->n)), (Interval*)interval_zero());
+            break;
+    }
+    
+    return result;
+}
+
+VALUE rb_Number_s_allocate (VALUE klass)
+{
+    void complex_free (Complex*);
+    return Data_Wrap_Struct(klass, 0, complex_free, 0);
+}
+
+VALUE rb_Number_s_digits (VALUE klass)
+{
+    return LONG2FIX(digs);
+}
+
+VALUE rb_Number_s_set_digits (VALUE klass, VALUE n)
+{
+    digs = FIX2LONG(n);
+    return n;
+}
+
+VALUE rb_Number_s_nan (VALUE klass)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_nan();
+    
+    return result;
+}
+
+VALUE rb_Number_s_pos_inf (VALUE klass)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_pos_inf();
+    
+    return result;
+}
+
+VALUE rb_Number_s_neg_inf (VALUE klass)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_neg_inf();
+    
+    return result;
+}
+
+VALUE rb_Number_s_zero (VALUE klass)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_zero();
+    
+    return result;
+}
+
+VALUE rb_Number_s_one (VALUE klass)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_one();
+    
+    return result;
+}
+
+VALUE rb_Number_s_pi (VALUE klass)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_pi();
+    
+    return result;
+}
+
+VALUE rb_Number_s_e (VALUE klass)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_e();
+    
+    return result;
+}
+
+VALUE rb_Number_s_i (VALUE klass)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_i();
+    
+    return result;
+}
+
+VALUE rb_Number_singleton_method_added (VALUE number, VALUE name)
+{
+    ID method_id = rb_to_id(name);
+    
+    rb_remove_method_id(rb_singleton_class(number), method_id);
+    rb_raise(rb_eTypeError, "can't define singleton method \"%s\" for %s", rb_id2name(method_id), rb_obj_classname(number));
+    
+    return Qnil; /* not reached */
+}
+
+VALUE rb_Number_initialize_copy (VALUE number, VALUE copy)
+{
+    rb_raise(rb_eTypeError, "can't copy %s", rb_obj_classname(number));
+    return Qnil; /* not reached */
+}
+
+VALUE rb_Number_to_s (VALUE number)
+{
+    char buf[1000];
+    VALUE str;
+    
+    complex_print(buf, DATA_PTR(number));
+    
+    return rb_str_new_cstr(buf);
+}
+
+VALUE rb_Number_dummy (VALUE number)
+{
+    return number;
+}
+
+VALUE rb_Number_coerce (VALUE a, VALUE b)
+{
+    return rb_assoc_new(rb_convert_type(b, T_DATA, "Number", "to_number"), a);
+}
+
+VALUE rb_Number_add (VALUE a, VALUE b)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_add(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number")));
+    
+    return result;
+}
+
+VALUE rb_Number_subtract (VALUE a, VALUE b)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_subtract(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number")));
+    
+    return result;
+}
+
+VALUE rb_Number_multiply (VALUE a, VALUE b)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_multiply(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number")));
+    
+    return result;
+}
+
+VALUE rb_Number_divide (VALUE a, VALUE b)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_divide(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number")));
+    
+    return result;
+}
+
+VALUE rb_Number_modulo (VALUE a, VALUE b)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_modulo(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number")));
+    
+    return result;
+}
+
+VALUE rb_Number_abs (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_abs(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_area (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_area(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_complex_p (VALUE number)
+{
+    return complex_complex_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_conjugate (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_conjugate(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_diagonal (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_diagonal(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_imag (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_im(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_negate (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_negate(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_real (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_re(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_real_p (VALUE number)
+{
+    return complex_real_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_reflect (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_reflect(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_round (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_round(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_sqrt (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_sqrt(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_round_c_c (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_round_c_c(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_round_c_f (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_round_c_f(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_round_f_c (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_round_f_c(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_round_f_f (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_round_f_f(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_round_i_i (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_round_i_i(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_round_i_o (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_round_i_o(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_round_o_i (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_round_o_i(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_round_o_o (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_round_o_o(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_complex_ll (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_ll(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_complex_lu (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_lu(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_complex_n (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_n(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_complex_ul (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_ul(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_complex_uu (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_uu(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_complex_magnitude (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_magnitude(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_complex_mignitude (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_mignitude(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_complex_width (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_width(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_imag_l (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_im_l(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_imag_n (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_im_n(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_imag_u (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_im_u(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_imag_magnitude (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_im_magnitude(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_imag_mignitude (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_im_mignitude(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_imag_width (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_im_width(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_real_l (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_re_l(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_real_n (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_re_n(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_real_u (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_re_u(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_real_magnitude (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_re_magnitude(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_real_mignitude (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_re_mignitude(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_real_width (VALUE number)
+{
+    VALUE result = rb_obj_alloc(rb_cNumber);
+    DATA_PTR(result) = (Complex*)complex_re_width(DATA_PTR(number));
+    
+    return result;
+}
+
+VALUE rb_Number_real_intersect_p (VALUE a, VALUE b)
+{
+    return complex_re_intersect_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_disjoint_p (VALUE a, VALUE b)
+{
+    return complex_re_disjoint_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_bounded_p (VALUE number)
+{
+    return complex_re_unbounded_p(DATA_PTR(number)) ? Qfalse : Qtrue;
+}
+
+VALUE rb_Number_real_left_bounded_p (VALUE number)
+{
+    return complex_re_left_unbounded_p(DATA_PTR(number)) ? Qfalse : Qtrue;
+}
+
+VALUE rb_Number_real_right_bounded_p (VALUE number)
+{
+    return complex_re_right_unbounded_p(DATA_PTR(number)) ? Qfalse : Qtrue;
+}
+
+VALUE rb_Number_real_unbounded_p (VALUE number)
+{
+    return complex_re_unbounded_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_left_unbounded_p (VALUE number)
+{
+    return complex_re_left_unbounded_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_right_unbounded_p (VALUE number)
+{
+    return complex_re_right_unbounded_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_contain_p (VALUE a, VALUE b)
+{
+    return complex_re_contain_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_contain_l_p (VALUE a, VALUE b)
+{
+    return complex_re_contain_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_contain_n_p (VALUE a, VALUE b)
+{
+    return complex_re_contain_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_contain_u_p (VALUE a, VALUE b)
+{
+    return complex_re_contain_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_contain_zero_p (VALUE number)
+{
+    return complex_re_contain_zero_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_eq_p (VALUE a, VALUE b)
+{
+    return complex_re_eq_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_bounds_eq_p (VALUE a, VALUE b)
+{
+    return complex_re_bounds_eq_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_normal_p (VALUE number)
+{
+    return complex_re_normal_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_zero_p (VALUE number)
+{
+    return complex_re_zero_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_nan_p (VALUE number)
+{
+    return complex_re_nan_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_l_eq_l_p (VALUE a, VALUE b)
+{
+    return complex_re_l_eq_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_l_eq_n_p (VALUE a, VALUE b)
+{
+    return complex_re_l_eq_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_l_eq_u_p (VALUE a, VALUE b)
+{
+    return complex_re_l_eq_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_n_eq_l_p (VALUE a, VALUE b)
+{
+    return complex_re_n_eq_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_n_eq_n_p (VALUE a, VALUE b)
+{
+    return complex_re_n_eq_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_n_eq_u_p (VALUE a, VALUE b)
+{
+    return complex_re_n_eq_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_u_eq_l_p (VALUE a, VALUE b)
+{
+    return complex_re_u_eq_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_u_eq_n_p (VALUE a, VALUE b)
+{
+    return complex_re_u_eq_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_u_eq_u_p (VALUE a, VALUE b)
+{
+    return complex_re_u_eq_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_l_gt_l_p (VALUE a, VALUE b)
+{
+    return complex_re_l_gt_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_l_gt_n_p (VALUE a, VALUE b)
+{
+    return complex_re_l_gt_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_l_gt_u_p (VALUE a, VALUE b)
+{
+    return complex_re_l_gt_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_n_gt_l_p (VALUE a, VALUE b)
+{
+    return complex_re_n_gt_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_n_gt_n_p (VALUE a, VALUE b)
+{
+    return complex_re_n_gt_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_n_gt_u_p (VALUE a, VALUE b)
+{
+    return complex_re_n_gt_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_u_gt_l_p (VALUE a, VALUE b)
+{
+    return complex_re_u_gt_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_u_gt_n_p (VALUE a, VALUE b)
+{
+    return complex_re_u_gt_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_u_gt_u_p (VALUE a, VALUE b)
+{
+    return complex_re_u_gt_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_l_ge_l_p (VALUE a, VALUE b)
+{
+    return complex_re_l_ge_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_l_ge_n_p (VALUE a, VALUE b)
+{
+    return complex_re_l_ge_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_l_ge_u_p (VALUE a, VALUE b)
+{
+    return complex_re_l_ge_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_n_ge_l_p (VALUE a, VALUE b)
+{
+    return complex_re_n_ge_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_n_ge_n_p (VALUE a, VALUE b)
+{
+    return complex_re_n_ge_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_n_ge_u_p (VALUE a, VALUE b)
+{
+    return complex_re_n_ge_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_u_ge_l_p (VALUE a, VALUE b)
+{
+    return complex_re_u_ge_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_u_ge_n_p (VALUE a, VALUE b)
+{
+    return complex_re_u_ge_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_u_ge_u_p (VALUE a, VALUE b)
+{
+    return complex_re_u_ge_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_l_lt_l_p (VALUE a, VALUE b)
+{
+    return complex_re_l_lt_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_l_lt_n_p (VALUE a, VALUE b)
+{
+    return complex_re_l_lt_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_l_lt_u_p (VALUE a, VALUE b)
+{
+    return complex_re_l_lt_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_n_lt_l_p (VALUE a, VALUE b)
+{
+    return complex_re_n_lt_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_n_lt_n_p (VALUE a, VALUE b)
+{
+    return complex_re_n_lt_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_n_lt_u_p (VALUE a, VALUE b)
+{
+    return complex_re_n_lt_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_u_lt_l_p (VALUE a, VALUE b)
+{
+    return complex_re_u_lt_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_u_lt_n_p (VALUE a, VALUE b)
+{
+    return complex_re_u_lt_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_u_lt_u_p (VALUE a, VALUE b)
+{
+    return complex_re_u_lt_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_l_le_l_p (VALUE a, VALUE b)
+{
+    return complex_re_l_le_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_l_le_n_p (VALUE a, VALUE b)
+{
+    return complex_re_l_le_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_l_le_u_p (VALUE a, VALUE b)
+{
+    return complex_re_l_le_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_n_le_l_p (VALUE a, VALUE b)
+{
+    return complex_re_n_le_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_n_le_n_p (VALUE a, VALUE b)
+{
+    return complex_re_n_le_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_n_le_u_p (VALUE a, VALUE b)
+{
+    return complex_re_n_le_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_u_le_l_p (VALUE a, VALUE b)
+{
+    return complex_re_u_le_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_u_le_n_p (VALUE a, VALUE b)
+{
+    return complex_re_u_le_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_real_u_le_u_p (VALUE a, VALUE b)
+{
+    return complex_re_u_le_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_intersect_p (VALUE a, VALUE b)
+{
+    return complex_im_intersect_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_disjoint_p (VALUE a, VALUE b)
+{
+    return complex_im_disjoint_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_bounded_p (VALUE number)
+{
+    return complex_im_unbounded_p(DATA_PTR(number)) ? Qfalse : Qtrue;
+}
+
+VALUE rb_Number_imag_left_bounded_p (VALUE number)
+{
+    return complex_im_left_unbounded_p(DATA_PTR(number)) ? Qfalse : Qtrue;
+}
+
+VALUE rb_Number_imag_right_bounded_p (VALUE number)
+{
+    return complex_im_right_unbounded_p(DATA_PTR(number)) ? Qfalse : Qtrue;
+}
+
+VALUE rb_Number_imag_unbounded_p (VALUE number)
+{
+    return complex_im_unbounded_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_left_unbounded_p (VALUE number)
+{
+    return complex_im_left_unbounded_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_right_unbounded_p (VALUE number)
+{
+    return complex_im_right_unbounded_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_contain_p (VALUE a, VALUE b)
+{
+    return complex_im_contain_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_contain_l_p (VALUE a, VALUE b)
+{
+    return complex_im_contain_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_contain_n_p (VALUE a, VALUE b)
+{
+    return complex_im_contain_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_contain_u_p (VALUE a, VALUE b)
+{
+    return complex_im_contain_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_contain_zero_p (VALUE number)
+{
+    return complex_im_contain_zero_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_eq_p (VALUE a, VALUE b)
+{
+    return complex_im_eq_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_bounds_eq_p (VALUE a, VALUE b)
+{
+    return complex_im_bounds_eq_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_nan_p (VALUE number)
+{
+    return complex_im_nan_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_normal_p (VALUE number)
+{
+    return complex_im_normal_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_zero_p (VALUE number)
+{
+    return complex_im_zero_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_l_eq_l_p (VALUE a, VALUE b)
+{
+    return complex_im_l_eq_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_l_eq_n_p (VALUE a, VALUE b)
+{
+    return complex_im_l_eq_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_l_eq_u_p (VALUE a, VALUE b)
+{
+    return complex_im_l_eq_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_n_eq_l_p (VALUE a, VALUE b)
+{
+    return complex_im_n_eq_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_n_eq_n_p (VALUE a, VALUE b)
+{
+    return complex_im_n_eq_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_n_eq_u_p (VALUE a, VALUE b)
+{
+    return complex_im_n_eq_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_u_eq_l_p (VALUE a, VALUE b)
+{
+    return complex_im_u_eq_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_u_eq_n_p (VALUE a, VALUE b)
+{
+    return complex_im_u_eq_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_u_eq_u_p (VALUE a, VALUE b)
+{
+    return complex_im_u_eq_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_l_gt_l_p (VALUE a, VALUE b)
+{
+    return complex_im_l_gt_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_l_gt_n_p (VALUE a, VALUE b)
+{
+    return complex_im_l_gt_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_l_gt_u_p (VALUE a, VALUE b)
+{
+    return complex_im_l_gt_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_n_gt_l_p (VALUE a, VALUE b)
+{
+    return complex_im_n_gt_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_n_gt_n_p (VALUE a, VALUE b)
+{
+    return complex_im_n_gt_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_n_gt_u_p (VALUE a, VALUE b)
+{
+    return complex_im_n_gt_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_u_gt_l_p (VALUE a, VALUE b)
+{
+    return complex_im_u_gt_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_u_gt_n_p (VALUE a, VALUE b)
+{
+    return complex_im_u_gt_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_u_gt_u_p (VALUE a, VALUE b)
+{
+    return complex_im_u_gt_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_l_ge_l_p (VALUE a, VALUE b)
+{
+    return complex_im_l_ge_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_l_ge_n_p (VALUE a, VALUE b)
+{
+    return complex_im_l_ge_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_l_ge_u_p (VALUE a, VALUE b)
+{
+    return complex_im_l_ge_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_n_ge_l_p (VALUE a, VALUE b)
+{
+    return complex_im_n_ge_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_n_ge_n_p (VALUE a, VALUE b)
+{
+    return complex_im_n_ge_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_n_ge_u_p (VALUE a, VALUE b)
+{
+    return complex_im_n_ge_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_u_ge_l_p (VALUE a, VALUE b)
+{
+    return complex_im_u_ge_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_u_ge_n_p (VALUE a, VALUE b)
+{
+    return complex_im_u_ge_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_u_ge_u_p (VALUE a, VALUE b)
+{
+    return complex_im_u_ge_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_l_lt_l_p (VALUE a, VALUE b)
+{
+    return complex_im_l_lt_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_l_lt_n_p (VALUE a, VALUE b)
+{
+    return complex_im_l_lt_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_l_lt_u_p (VALUE a, VALUE b)
+{
+    return complex_im_l_lt_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_n_lt_l_p (VALUE a, VALUE b)
+{
+    return complex_im_n_lt_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_n_lt_n_p (VALUE a, VALUE b)
+{
+    return complex_im_n_lt_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_n_lt_u_p (VALUE a, VALUE b)
+{
+    return complex_im_n_lt_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_u_lt_l_p (VALUE a, VALUE b)
+{
+    return complex_im_u_lt_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_u_lt_n_p (VALUE a, VALUE b)
+{
+    return complex_im_u_lt_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_u_lt_u_p (VALUE a, VALUE b)
+{
+    return complex_im_u_lt_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_l_le_l_p (VALUE a, VALUE b)
+{
+    return complex_im_l_le_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_l_le_n_p (VALUE a, VALUE b)
+{
+    return complex_im_l_le_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_l_le_u_p (VALUE a, VALUE b)
+{
+    return complex_im_l_le_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_n_le_l_p (VALUE a, VALUE b)
+{
+    return complex_im_n_le_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_n_le_n_p (VALUE a, VALUE b)
+{
+    return complex_im_n_le_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_n_le_u_p (VALUE a, VALUE b)
+{
+    return complex_im_n_le_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_u_le_l_p (VALUE a, VALUE b)
+{
+    return complex_im_u_le_l_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_u_le_n_p (VALUE a, VALUE b)
+{
+    return complex_im_u_le_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_imag_u_le_u_p (VALUE a, VALUE b)
+{
+    return complex_im_u_le_u_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_complex_intersect_p (VALUE a, VALUE b)
+{
+    return complex_intersect_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_complex_disjoint_p (VALUE a, VALUE b)
+{
+    return complex_disjoint_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_complex_bounded_p (VALUE number)
+{
+    return complex_unbounded_p(DATA_PTR(number)) ? Qfalse : Qtrue;
+}
+
+VALUE rb_Number_complex_unbounded_p (VALUE number)
+{
+    return complex_unbounded_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_complex_contain_p (VALUE a, VALUE b)
+{
+    return complex_contain_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_complex_contain_ll_p (VALUE a, VALUE b)
+{
+    return complex_contain_ll_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_complex_contain_lu_p (VALUE a, VALUE b)
+{
+    return complex_contain_lu_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_complex_contain_n_p (VALUE a, VALUE b)
+{
+    return complex_contain_n_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_complex_contain_ul_p (VALUE a, VALUE b)
+{
+    return complex_contain_ul_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_complex_contain_uu_p (VALUE a, VALUE b)
+{
+    return complex_contain_uu_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_complex_contain_zero_p (VALUE number)
+{
+    return complex_contain_zero_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_complex_nan_p (VALUE number)
+{
+    return complex_nan_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_complex_normal_p (VALUE number)
+{
+    return complex_normal_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_complex_zero_p (VALUE number)
+{
+    return complex_zero_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_complex_real_p (VALUE number)
+{
+    return complex_real_p(DATA_PTR(number)) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_complex_eq_p (VALUE a, VALUE b)
+{
+    return complex_eq_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+VALUE rb_Number_complex_bounds_eq_p (VALUE a, VALUE b)
+{
+    return complex_bounds_eq_p(DATA_PTR(a), DATA_PTR(rb_convert_type(b, T_DATA, "Number", "to_number"))) ? Qtrue : Qfalse;
+}
+
+void Init_number ()
+{
+    init_tmp_nums();
+    
+    digs = 17;
+    
+    rb_define_method(rb_cInteger, "to_number", rb_Integer_to_number, 0);
+    rb_define_method(rb_cFloat, "to_number", rb_Float_to_number, 0);
+    rb_define_method(rb_cComplex, "to_number", rb_Complex_to_number, 0);
+    
+    rb_cNumber = rb_define_class("Number", rb_cObject);
+    rb_undef_method(CLASS_OF(rb_cNumber), "new");
+    rb_define_alloc_func(rb_cNumber, rb_Number_s_allocate);
+    rb_define_global_function("Number", rb_f_Number, -1);
+    
+    rb_define_singleton_method(rb_cNumber, "digits", rb_Number_s_digits, 0);
+    rb_define_singleton_method(rb_cNumber, "digits=", rb_Number_s_set_digits, 1);
+    
+    rb_define_singleton_method(rb_cNumber, "nan", rb_Number_s_nan, 0);
+    rb_define_singleton_method(rb_cNumber, "pos_inf", rb_Number_s_pos_inf, 0);
+    rb_define_singleton_method(rb_cNumber, "neg_inf", rb_Number_s_neg_inf, 0);
+    rb_define_singleton_method(rb_cNumber, "zero", rb_Number_s_zero, 0);
+    rb_define_singleton_method(rb_cNumber, "one", rb_Number_s_one, 0);
+    rb_define_singleton_method(rb_cNumber, "pi", rb_Number_s_pi, 0);
+    rb_define_singleton_method(rb_cNumber, "e", rb_Number_s_e, 0);
+    rb_define_singleton_method(rb_cNumber, "i", rb_Number_s_i, 0);
+    
+    rb_define_method(rb_cNumber, "singleton_method_added", rb_Number_singleton_method_added, 1);
+    rb_define_method(rb_cNumber, "initialize_copy", rb_Number_initialize_copy, 1);
+    
+    rb_define_method(rb_cNumber, "to_s", rb_Number_to_s, 0);
+    rb_define_method(rb_cNumber, "to_number", rb_Number_dummy, 0);
+    rb_define_method(rb_cNumber, "coerce", rb_Number_coerce, 1);
+    
+    rb_define_method(rb_cNumber, "+@", rb_Number_dummy, 0);
+    rb_define_method(rb_cNumber, "-@", rb_Number_negate, 0);
+    
+    rb_define_method(rb_cNumber, "+", rb_Number_add, 1);
+    rb_define_method(rb_cNumber, "-", rb_Number_subtract, 1);
+    rb_define_method(rb_cNumber, "*", rb_Number_multiply, 1);
+    rb_define_method(rb_cNumber, "/", rb_Number_divide, 1);
+    rb_define_method(rb_cNumber, "%", rb_Number_modulo, 1);
+    
+    rb_define_method(rb_cNumber, "<", rb_Number_real_n_lt_n_p, 1);
+    rb_define_method(rb_cNumber, "<=", rb_Number_real_n_le_n_p, 1);
+    rb_define_method(rb_cNumber, ">", rb_Number_real_n_gt_n_p, 1);
+    rb_define_method(rb_cNumber, ">=", rb_Number_real_n_ge_n_p, 1);
+    rb_define_method(rb_cNumber, "==", rb_Number_real_n_eq_n_p, 1);
+    rb_define_method(rb_cNumber, "===", rb_Number_complex_eq_p, 1);
+    
+    rb_define_method(rb_cNumber, "abs", rb_Number_abs, 0);
+    rb_define_method(rb_cNumber, "area", rb_Number_area, 0);
+    rb_define_method(rb_cNumber, "complex?", rb_Number_complex_p, 0);
+    rb_define_method(rb_cNumber, "conjugate", rb_Number_conjugate, 0);
+    rb_define_method(rb_cNumber, "diagonal", rb_Number_diagonal, 0);
+    rb_define_method(rb_cNumber, "imag", rb_Number_imag, 0);
+    rb_define_method(rb_cNumber, "negate", rb_Number_negate, 0);
+    rb_define_method(rb_cNumber, "reflect", rb_Number_reflect, 0);
+    rb_define_method(rb_cNumber, "round", rb_Number_round, 0);
+    rb_define_method(rb_cNumber, "real", rb_Number_real, 0);
+    rb_define_method(rb_cNumber, "real?", rb_Number_real_p, 0);
+    rb_define_method(rb_cNumber, "sqrt", rb_Number_sqrt, 0);
+    
+    rb_define_method(rb_cNumber, "ceil", rb_Number_round_c_c, 0);
+    rb_define_method(rb_cNumber, "floor", rb_Number_round_f_f, 0);
+    rb_define_method(rb_cNumber, "round_c_c", rb_Number_round_c_c, 0);
+    rb_define_method(rb_cNumber, "round_c_f", rb_Number_round_c_f, 0);
+    rb_define_method(rb_cNumber, "round_f_c", rb_Number_round_f_c, 0);
+    rb_define_method(rb_cNumber, "round_f_f", rb_Number_round_f_f, 0);
+    
+    rb_define_method(rb_cNumber, "round_i_i", rb_Number_round_i_i, 0);
+    rb_define_method(rb_cNumber, "round_i_o", rb_Number_round_i_o, 0);
+    rb_define_method(rb_cNumber, "round_o_i", rb_Number_round_o_i, 0);
+    rb_define_method(rb_cNumber, "round_o_o", rb_Number_round_o_o, 0);
+    
+    rb_define_method(rb_cNumber, "infimum", rb_Number_real_l, 0);
+    rb_define_method(rb_cNumber, "l", rb_Number_real_l, 0);
+    rb_define_method(rb_cNumber, "real_l", rb_Number_real_l, 0);
+    rb_define_method(rb_cNumber, "imag_l", rb_Number_imag_l, 0);
+    rb_define_method(rb_cNumber, "complex_ll", rb_Number_complex_ll, 0);
+    rb_define_method(rb_cNumber, "complex_lu", rb_Number_complex_lu, 0);
+    
+    rb_define_method(rb_cNumber, "magnitude", rb_Number_real_magnitude, 0);
+    rb_define_method(rb_cNumber, "real_magnitude", rb_Number_real_magnitude, 0);
+    rb_define_method(rb_cNumber, "imag_magnitude", rb_Number_imag_magnitude, 0);
+    rb_define_method(rb_cNumber, "complex_magnitude", rb_Number_complex_magnitude, 0);
+    
+    rb_define_method(rb_cNumber, "mignitude", rb_Number_real_mignitude, 0);
+    rb_define_method(rb_cNumber, "real_mignitude", rb_Number_real_mignitude, 0);
+    rb_define_method(rb_cNumber, "imag_mignitude", rb_Number_imag_mignitude, 0);
+    rb_define_method(rb_cNumber, "complex_mignitude", rb_Number_complex_mignitude, 0);
+    
+    rb_define_method(rb_cNumber, "nominal", rb_Number_real_n, 0);
+    rb_define_method(rb_cNumber, "n", rb_Number_real_n, 0);
+    rb_define_method(rb_cNumber, "real_n", rb_Number_real_n, 0);
+    rb_define_method(rb_cNumber, "imag_n", rb_Number_imag_n, 0);
+    rb_define_method(rb_cNumber, "complex_n", rb_Number_complex_n, 0);
+    
+    rb_define_method(rb_cNumber, "supremum", rb_Number_real_u, 0);
+    rb_define_method(rb_cNumber, "u", rb_Number_real_u, 0);
+    rb_define_method(rb_cNumber, "real_u", rb_Number_real_u, 0);
+    rb_define_method(rb_cNumber, "imag_u", rb_Number_imag_u, 0);
+    rb_define_method(rb_cNumber, "complex_ul", rb_Number_complex_uu, 0);
+    rb_define_method(rb_cNumber, "complex_uu", rb_Number_complex_ul, 0);
+    
+    rb_define_method(rb_cNumber, "width", rb_Number_real_width, 0);
+    rb_define_method(rb_cNumber, "real_width", rb_Number_real_width, 0);
+    rb_define_method(rb_cNumber, "imag_width", rb_Number_imag_width, 0);
+    rb_define_method(rb_cNumber, "complex_width", rb_Number_complex_width, 0);
+    
+    rb_define_method(rb_cNumber, "intersect?", rb_Number_real_intersect_p, 1);
+    rb_define_method(rb_cNumber, "real_intersect?", rb_Number_real_intersect_p, 1);
+    rb_define_method(rb_cNumber, "imag_intersect?", rb_Number_imag_intersect_p, 1);
+    rb_define_method(rb_cNumber, "complex_intersect?", rb_Number_complex_intersect_p, 1);
+    
+    rb_define_method(rb_cNumber, "disjoint?", rb_Number_real_disjoint_p, 1);
+    rb_define_method(rb_cNumber, "real_disjoint?", rb_Number_real_disjoint_p, 1);
+    rb_define_method(rb_cNumber, "imag_disjoint?", rb_Number_imag_disjoint_p, 1);
+    rb_define_method(rb_cNumber, "complex_disjoint?", rb_Number_complex_disjoint_p, 1);
+    
+    rb_define_method(rb_cNumber, "contain?", rb_Number_real_contain_p, 1);
+    rb_define_method(rb_cNumber, "real_contain?", rb_Number_real_contain_p, 1);
+    rb_define_method(rb_cNumber, "imag_contain?", rb_Number_imag_contain_p, 1);
+    rb_define_method(rb_cNumber, "complex_contain?", rb_Number_complex_contain_p, 1);
+    
+    rb_define_method(rb_cNumber, "contain_l?", rb_Number_real_contain_l_p, 1);
+    rb_define_method(rb_cNumber, "real_contain_l?", rb_Number_real_contain_l_p, 1);
+    rb_define_method(rb_cNumber, "imag_contain_l?", rb_Number_imag_contain_l_p, 1);
+    rb_define_method(rb_cNumber, "complex_contain_ll?", rb_Number_complex_contain_ll_p, 1);
+    rb_define_method(rb_cNumber, "complex_contain_lu?", rb_Number_complex_contain_lu_p, 1);
+    
+    rb_define_method(rb_cNumber, "contain_n?", rb_Number_real_contain_n_p, 1);
+    rb_define_method(rb_cNumber, "real_contain_n?", rb_Number_real_contain_n_p, 1);
+    rb_define_method(rb_cNumber, "imag_contain_n?", rb_Number_imag_contain_n_p, 1);
+    rb_define_method(rb_cNumber, "complex_contain_n?", rb_Number_complex_contain_n_p, 1);
+    
+    rb_define_method(rb_cNumber, "contain_u?", rb_Number_real_contain_u_p, 1);
+    rb_define_method(rb_cNumber, "real_contain_u?", rb_Number_real_contain_u_p, 1);
+    rb_define_method(rb_cNumber, "imag_contain_u?", rb_Number_imag_contain_u_p, 1);
+    rb_define_method(rb_cNumber, "complex_contain_ul?", rb_Number_complex_contain_ul_p, 1);
+    rb_define_method(rb_cNumber, "complex_contain_uu?", rb_Number_complex_contain_uu_p, 1);
+    
+    rb_define_method(rb_cNumber, "contain_zero?", rb_Number_real_contain_zero_p, 0);
+    rb_define_method(rb_cNumber, "real_contain_zero?", rb_Number_real_contain_zero_p, 0);
+    rb_define_method(rb_cNumber, "imag_contain_zero?", rb_Number_imag_contain_zero_p, 0);
+    rb_define_method(rb_cNumber, "complex_contain_zero?", rb_Number_complex_contain_zero_p, 0);
+    
+    rb_define_method(rb_cNumber, "eq?", rb_Number_real_eq_p, 1);
+    rb_define_method(rb_cNumber, "real_eq?", rb_Number_real_eq_p, 1);
+    rb_define_method(rb_cNumber, "imag_eq?", rb_Number_imag_eq_p, 1);
+    rb_define_method(rb_cNumber, "complex_eq?", rb_Number_complex_eq_p, 1);
+    
+    rb_define_method(rb_cNumber, "bounds_eq?", rb_Number_real_bounds_eq_p, 1);
+    rb_define_method(rb_cNumber, "real_bounds_eq?", rb_Number_real_bounds_eq_p, 1);
+    rb_define_method(rb_cNumber, "imag_bounds_eq?", rb_Number_imag_bounds_eq_p, 1);
+    rb_define_method(rb_cNumber, "complex_bounds_eq?", rb_Number_complex_bounds_eq_p, 1);
+    
+    rb_define_method(rb_cNumber, "bounded?", rb_Number_real_bounded_p, 0);
+    rb_define_method(rb_cNumber, "real_bounded?", rb_Number_real_bounded_p, 0);
+    rb_define_method(rb_cNumber, "imag_bounded?", rb_Number_imag_bounded_p, 0);
+    rb_define_method(rb_cNumber, "complex_bounded?", rb_Number_complex_bounded_p, 0);
+    
+    rb_define_method(rb_cNumber, "unbounded?", rb_Number_real_unbounded_p, 0);
+    rb_define_method(rb_cNumber, "real_unbounded?", rb_Number_real_unbounded_p, 0);
+    rb_define_method(rb_cNumber, "imag_unbounded?", rb_Number_imag_unbounded_p, 0);
+    rb_define_method(rb_cNumber, "complex_unbounded?", rb_Number_complex_unbounded_p, 0);
+    
+    rb_define_method(rb_cNumber, "left_bounded?", rb_Number_real_left_bounded_p, 0);
+    rb_define_method(rb_cNumber, "real_left_bounded?", rb_Number_real_left_bounded_p, 0);
+    rb_define_method(rb_cNumber, "imag_left_bounded?", rb_Number_imag_left_bounded_p, 0);
+    
+    rb_define_method(rb_cNumber, "left_unbounded?", rb_Number_real_left_unbounded_p, 0);
+    rb_define_method(rb_cNumber, "real_left_unbounded?", rb_Number_real_left_unbounded_p, 0);
+    rb_define_method(rb_cNumber, "imag_left_unbounded?", rb_Number_imag_left_unbounded_p, 0);
+    
+    rb_define_method(rb_cNumber, "right_bounded?", rb_Number_real_right_bounded_p, 0);
+    rb_define_method(rb_cNumber, "real_right_bounded?", rb_Number_real_right_bounded_p, 0);
+    rb_define_method(rb_cNumber, "imag_right_bounded?", rb_Number_imag_right_bounded_p, 0);
+    
+    rb_define_method(rb_cNumber, "right_unbounded?", rb_Number_real_right_unbounded_p, 0);
+    rb_define_method(rb_cNumber, "real_right_unbounded?", rb_Number_real_right_unbounded_p, 0);
+    rb_define_method(rb_cNumber, "imag_right_unbounded?", rb_Number_imag_right_unbounded_p, 0);
+    
+    rb_define_method(rb_cNumber, "normal?", rb_Number_real_normal_p, 0);
+    rb_define_method(rb_cNumber, "real_normal?", rb_Number_real_normal_p, 0);
+    rb_define_method(rb_cNumber, "imag_normal?", rb_Number_imag_normal_p, 0);
+    rb_define_method(rb_cNumber, "complex_normal?", rb_Number_complex_normal_p, 0);
+    
+    rb_define_method(rb_cNumber, "nan?", rb_Number_real_nan_p, 0);
+    rb_define_method(rb_cNumber, "real_nan?", rb_Number_real_nan_p, 0);
+    rb_define_method(rb_cNumber, "imag_nan?", rb_Number_imag_nan_p, 0);
+    rb_define_method(rb_cNumber, "complex_nan?", rb_Number_complex_nan_p, 0);
+    
+    //rb_define_method(rb_cNumber, "finite?", rb_Number_real_finite_p, 0);
+    //rb_define_method(rb_cNumber, "real_finite?", rb_Number_real_finite_p, 0);
+    //rb_define_method(rb_cNumber, "imag_finite?", rb_Number_imag_finite_p, 0);
+    //rb_define_method(rb_cNumber, "complex_finite?", rb_Number_complex_finite_p, 0);
+    //
+    //rb_define_method(rb_cNumber, "infinite?", rb_Number_real_infinite_p, 0);
+    //rb_define_method(rb_cNumber, "real_infinite?", rb_Number_real_infinite_p, 0);
+    //rb_define_method(rb_cNumber, "imag_infinite?", rb_Number_imag_infinite_p, 0);
+    //rb_define_method(rb_cNumber, "complex_infinite?", rb_Number_complex_infinite_p, 0);
+    
+    rb_define_method(rb_cNumber, "zero?", rb_Number_real_zero_p, 0);
+    rb_define_method(rb_cNumber, "real_zero?", rb_Number_real_zero_p, 0);
+    rb_define_method(rb_cNumber, "imag_zero?", rb_Number_imag_zero_p, 0);
+    rb_define_method(rb_cNumber, "complex_zero?", rb_Number_complex_zero_p, 0);
+    
+    rb_define_method(rb_cNumber, "l_eq_l?", rb_Number_real_l_eq_l_p, 1);
+    rb_define_method(rb_cNumber, "real_l_eq_l?", rb_Number_real_l_eq_l_p, 1);
+    rb_define_method(rb_cNumber, "imag_l_eq_l?", rb_Number_imag_l_eq_l_p, 1);
+    
+    rb_define_method(rb_cNumber, "l_eq_n?", rb_Number_real_l_eq_n_p, 1);
+    rb_define_method(rb_cNumber, "real_l_eq_n?", rb_Number_real_l_eq_n_p, 1);
+    rb_define_method(rb_cNumber, "imag_l_eq_n?", rb_Number_imag_l_eq_n_p, 1);
+    
+    rb_define_method(rb_cNumber, "l_eq_u?", rb_Number_real_l_eq_u_p, 1);
+    rb_define_method(rb_cNumber, "real_l_eq_u?", rb_Number_real_l_eq_u_p, 1);
+    rb_define_method(rb_cNumber, "imag_l_eq_u?", rb_Number_imag_l_eq_u_p, 1);
+    
+    rb_define_method(rb_cNumber, "n_eq_l?", rb_Number_real_n_eq_l_p, 1);
+    rb_define_method(rb_cNumber, "real_n_eq_l?", rb_Number_real_n_eq_l_p, 1);
+    rb_define_method(rb_cNumber, "imag_n_eq_l?", rb_Number_imag_n_eq_l_p, 1);
+    
+    rb_define_method(rb_cNumber, "n_eq_n?", rb_Number_real_n_eq_n_p, 1);
+    rb_define_method(rb_cNumber, "real_n_eq_n?", rb_Number_real_n_eq_n_p, 1);
+    rb_define_method(rb_cNumber, "imag_n_eq_n?", rb_Number_imag_n_eq_n_p, 1);
+    
+    rb_define_method(rb_cNumber, "n_eq_u?", rb_Number_real_n_eq_u_p, 1);
+    rb_define_method(rb_cNumber, "real_n_eq_u?", rb_Number_real_n_eq_u_p, 1);
+    rb_define_method(rb_cNumber, "imag_n_eq_u?", rb_Number_imag_n_eq_u_p, 1);
+    
+    rb_define_method(rb_cNumber, "u_eq_l?", rb_Number_real_u_eq_l_p, 1);
+    rb_define_method(rb_cNumber, "real_u_eq_l?", rb_Number_real_u_eq_l_p, 1);
+    rb_define_method(rb_cNumber, "imag_u_eq_l?", rb_Number_imag_u_eq_l_p, 1);
+    
+    rb_define_method(rb_cNumber, "u_eq_n?", rb_Number_real_u_eq_n_p, 1);
+    rb_define_method(rb_cNumber, "real_u_eq_n?", rb_Number_real_u_eq_n_p, 1);
+    rb_define_method(rb_cNumber, "imag_u_eq_n?", rb_Number_imag_u_eq_n_p, 1);
+    
+    rb_define_method(rb_cNumber, "u_eq_u?", rb_Number_real_u_eq_u_p, 1);
+    rb_define_method(rb_cNumber, "real_u_eq_u?", rb_Number_real_u_eq_u_p, 1);
+    rb_define_method(rb_cNumber, "imag_u_eq_u?", rb_Number_imag_u_eq_u_p, 1);
+    
+    rb_define_method(rb_cNumber, "l_gt_l?", rb_Number_real_l_gt_l_p, 1);
+    rb_define_method(rb_cNumber, "real_l_gt_l?", rb_Number_real_l_gt_l_p, 1);
+    rb_define_method(rb_cNumber, "imag_l_gt_l?", rb_Number_imag_l_gt_l_p, 1);
+    
+    rb_define_method(rb_cNumber, "l_gt_n?", rb_Number_real_l_gt_n_p, 1);
+    rb_define_method(rb_cNumber, "real_l_gt_n?", rb_Number_real_l_gt_n_p, 1);
+    rb_define_method(rb_cNumber, "imag_l_gt_n?", rb_Number_imag_l_gt_n_p, 1);
+    
+    rb_define_method(rb_cNumber, "l_gt_u?", rb_Number_real_l_gt_u_p, 1);
+    rb_define_method(rb_cNumber, "real_l_gt_u?", rb_Number_real_l_gt_u_p, 1);
+    rb_define_method(rb_cNumber, "imag_l_gt_u?", rb_Number_imag_l_gt_u_p, 1);
+    
+    rb_define_method(rb_cNumber, "n_gt_l?", rb_Number_real_n_gt_l_p, 1);
+    rb_define_method(rb_cNumber, "real_n_gt_l?", rb_Number_real_n_gt_l_p, 1);
+    rb_define_method(rb_cNumber, "imag_n_gt_l?", rb_Number_imag_n_gt_l_p, 1);
+    
+    rb_define_method(rb_cNumber, "n_gt_n?", rb_Number_real_n_gt_n_p, 1);
+    rb_define_method(rb_cNumber, "real_n_gt_n?", rb_Number_real_n_gt_n_p, 1);
+    rb_define_method(rb_cNumber, "imag_n_gt_n?", rb_Number_imag_n_gt_n_p, 1);
+    
+    rb_define_method(rb_cNumber, "n_gt_u?", rb_Number_real_n_gt_u_p, 1);
+    rb_define_method(rb_cNumber, "real_n_gt_u?", rb_Number_real_n_gt_u_p, 1);
+    rb_define_method(rb_cNumber, "imag_n_gt_u?", rb_Number_imag_n_gt_u_p, 1);
+    
+    rb_define_method(rb_cNumber, "u_gt_l?", rb_Number_real_u_gt_l_p, 1);
+    rb_define_method(rb_cNumber, "real_u_gt_l?", rb_Number_real_u_gt_l_p, 1);
+    rb_define_method(rb_cNumber, "imag_u_gt_l?", rb_Number_imag_u_gt_l_p, 1);
+    
+    rb_define_method(rb_cNumber, "u_gt_n?", rb_Number_real_u_gt_n_p, 1);
+    rb_define_method(rb_cNumber, "real_u_gt_n?", rb_Number_real_u_gt_n_p, 1);
+    rb_define_method(rb_cNumber, "imag_u_gt_n?", rb_Number_imag_u_gt_n_p, 1);
+    
+    rb_define_method(rb_cNumber, "u_gt_u?", rb_Number_real_u_gt_u_p, 1);
+    rb_define_method(rb_cNumber, "real_u_gt_u?", rb_Number_real_u_gt_u_p, 1);
+    rb_define_method(rb_cNumber, "imag_u_gt_u?", rb_Number_imag_u_gt_u_p, 1);
+    
+    rb_define_method(rb_cNumber, "l_ge_l?", rb_Number_real_l_ge_l_p, 1);
+    rb_define_method(rb_cNumber, "real_l_ge_l?", rb_Number_real_l_ge_l_p, 1);
+    rb_define_method(rb_cNumber, "imag_l_ge_l?", rb_Number_imag_l_ge_l_p, 1);
+    
+    rb_define_method(rb_cNumber, "l_ge_n?", rb_Number_real_l_ge_n_p, 1);
+    rb_define_method(rb_cNumber, "real_l_ge_n?", rb_Number_real_l_ge_n_p, 1);
+    rb_define_method(rb_cNumber, "imag_l_ge_n?", rb_Number_imag_l_ge_n_p, 1);
+    
+    rb_define_method(rb_cNumber, "l_ge_u?", rb_Number_real_l_ge_u_p, 1);
+    rb_define_method(rb_cNumber, "real_l_ge_u?", rb_Number_real_l_ge_u_p, 1);
+    rb_define_method(rb_cNumber, "imag_l_ge_u?", rb_Number_imag_l_ge_u_p, 1);
+    
+    rb_define_method(rb_cNumber, "n_ge_l?", rb_Number_real_n_ge_l_p, 1);
+    rb_define_method(rb_cNumber, "real_n_ge_l?", rb_Number_real_n_ge_l_p, 1);
+    rb_define_method(rb_cNumber, "imag_n_ge_l?", rb_Number_imag_n_ge_l_p, 1);
+    
+    rb_define_method(rb_cNumber, "n_ge_n?", rb_Number_real_n_ge_n_p, 1);
+    rb_define_method(rb_cNumber, "real_n_ge_n?", rb_Number_real_n_ge_n_p, 1);
+    rb_define_method(rb_cNumber, "imag_n_ge_n?", rb_Number_imag_n_ge_n_p, 1);
+    
+    rb_define_method(rb_cNumber, "n_ge_u?", rb_Number_real_n_ge_u_p, 1);
+    rb_define_method(rb_cNumber, "real_n_ge_u?", rb_Number_real_n_ge_u_p, 1);
+    rb_define_method(rb_cNumber, "imag_n_ge_u?", rb_Number_imag_n_ge_u_p, 1);
+    
+    rb_define_method(rb_cNumber, "u_ge_l?", rb_Number_real_u_ge_l_p, 1);
+    rb_define_method(rb_cNumber, "real_u_ge_l?", rb_Number_real_u_ge_l_p, 1);
+    rb_define_method(rb_cNumber, "imag_u_ge_l?", rb_Number_imag_u_ge_l_p, 1);
+    
+    rb_define_method(rb_cNumber, "u_ge_n?", rb_Number_real_u_ge_n_p, 1);
+    rb_define_method(rb_cNumber, "real_u_ge_n?", rb_Number_real_u_ge_n_p, 1);
+    rb_define_method(rb_cNumber, "imag_u_ge_n?", rb_Number_imag_u_ge_n_p, 1);
+    
+    rb_define_method(rb_cNumber, "u_ge_u?", rb_Number_real_u_ge_u_p, 1);
+    rb_define_method(rb_cNumber, "real_u_ge_u?", rb_Number_real_u_ge_u_p, 1);
+    rb_define_method(rb_cNumber, "imag_u_ge_u?", rb_Number_imag_u_ge_u_p, 1);
+    
+    rb_define_method(rb_cNumber, "l_lt_l?", rb_Number_real_l_lt_l_p, 1);
+    rb_define_method(rb_cNumber, "real_l_lt_l?", rb_Number_real_l_lt_l_p, 1);
+    rb_define_method(rb_cNumber, "imag_l_lt_l?", rb_Number_imag_l_lt_l_p, 1);
+    
+    rb_define_method(rb_cNumber, "l_lt_n?", rb_Number_real_l_lt_n_p, 1);
+    rb_define_method(rb_cNumber, "real_l_lt_n?", rb_Number_real_l_lt_n_p, 1);
+    rb_define_method(rb_cNumber, "imag_l_lt_n?", rb_Number_imag_l_lt_n_p, 1);
+    
+    rb_define_method(rb_cNumber, "l_lt_u?", rb_Number_real_l_lt_u_p, 1);
+    rb_define_method(rb_cNumber, "real_l_lt_u?", rb_Number_real_l_lt_u_p, 1);
+    rb_define_method(rb_cNumber, "imag_l_lt_u?", rb_Number_imag_l_lt_u_p, 1);
+    
+    rb_define_method(rb_cNumber, "n_lt_l?", rb_Number_real_n_lt_l_p, 1);
+    rb_define_method(rb_cNumber, "real_n_lt_l?", rb_Number_real_n_lt_l_p, 1);
+    rb_define_method(rb_cNumber, "imag_n_lt_l?", rb_Number_imag_n_lt_l_p, 1);
+    
+    rb_define_method(rb_cNumber, "n_lt_n?", rb_Number_real_n_lt_n_p, 1);
+    rb_define_method(rb_cNumber, "real_n_lt_n?", rb_Number_real_n_lt_n_p, 1);
+    rb_define_method(rb_cNumber, "imag_n_lt_n?", rb_Number_imag_n_lt_n_p, 1);
+    
+    rb_define_method(rb_cNumber, "n_lt_u?", rb_Number_real_n_lt_u_p, 1);
+    rb_define_method(rb_cNumber, "real_n_lt_u?", rb_Number_real_n_lt_u_p, 1);
+    rb_define_method(rb_cNumber, "imag_n_lt_u?", rb_Number_imag_n_lt_u_p, 1);
+    
+    rb_define_method(rb_cNumber, "u_lt_l?", rb_Number_real_u_lt_l_p, 1);
+    rb_define_method(rb_cNumber, "real_u_lt_l?", rb_Number_real_u_lt_l_p, 1);
+    rb_define_method(rb_cNumber, "imag_u_lt_l?", rb_Number_imag_u_lt_l_p, 1);
+    
+    rb_define_method(rb_cNumber, "u_lt_n?", rb_Number_real_u_lt_n_p, 1);
+    rb_define_method(rb_cNumber, "real_u_lt_n?", rb_Number_real_u_lt_n_p, 1);
+    rb_define_method(rb_cNumber, "imag_u_lt_n?", rb_Number_imag_u_lt_n_p, 1);
+    
+    rb_define_method(rb_cNumber, "u_lt_u?", rb_Number_real_u_lt_u_p, 1);
+    rb_define_method(rb_cNumber, "real_u_lt_u?", rb_Number_real_u_lt_u_p, 1);
+    rb_define_method(rb_cNumber, "imag_u_lt_u?", rb_Number_imag_u_lt_u_p, 1);
+    
+    rb_define_method(rb_cNumber, "l_le_l?", rb_Number_real_l_le_l_p, 1);
+    rb_define_method(rb_cNumber, "real_l_le_l?", rb_Number_real_l_le_l_p, 1);
+    rb_define_method(rb_cNumber, "imag_l_le_l?", rb_Number_imag_l_le_l_p, 1);
+    
+    rb_define_method(rb_cNumber, "l_le_n?", rb_Number_real_l_le_n_p, 1);
+    rb_define_method(rb_cNumber, "real_l_le_n?", rb_Number_real_l_le_n_p, 1);
+    rb_define_method(rb_cNumber, "imag_l_le_n?", rb_Number_imag_l_le_n_p, 1);
+    
+    rb_define_method(rb_cNumber, "l_le_u?", rb_Number_real_l_le_u_p, 1);
+    rb_define_method(rb_cNumber, "real_l_le_u?", rb_Number_real_l_le_u_p, 1);
+    rb_define_method(rb_cNumber, "imag_l_le_u?", rb_Number_imag_l_le_u_p, 1);
+    
+    rb_define_method(rb_cNumber, "n_le_l?", rb_Number_real_n_le_l_p, 1);
+    rb_define_method(rb_cNumber, "real_n_le_l?", rb_Number_real_n_le_l_p, 1);
+    rb_define_method(rb_cNumber, "imag_n_le_l?", rb_Number_imag_n_le_l_p, 1);
+    
+    rb_define_method(rb_cNumber, "n_le_n?", rb_Number_real_n_le_n_p, 1);
+    rb_define_method(rb_cNumber, "real_n_le_n?", rb_Number_real_n_le_n_p, 1);
+    rb_define_method(rb_cNumber, "imag_n_le_n?", rb_Number_imag_n_le_n_p, 1);
+    
+    rb_define_method(rb_cNumber, "n_le_u?", rb_Number_real_n_le_u_p, 1);
+    rb_define_method(rb_cNumber, "real_n_le_u?", rb_Number_real_n_le_u_p, 1);
+    rb_define_method(rb_cNumber, "imag_n_le_u?", rb_Number_imag_n_le_u_p, 1);
+    
+    rb_define_method(rb_cNumber, "u_le_l?", rb_Number_real_u_le_l_p, 1);
+    rb_define_method(rb_cNumber, "real_u_le_l?", rb_Number_real_u_le_l_p, 1);
+    rb_define_method(rb_cNumber, "imag_u_le_l?", rb_Number_imag_u_le_l_p, 1);
+    
+    rb_define_method(rb_cNumber, "u_le_n?", rb_Number_real_u_le_n_p, 1);
+    rb_define_method(rb_cNumber, "real_u_le_n?", rb_Number_real_u_le_n_p, 1);
+    rb_define_method(rb_cNumber, "imag_u_le_n?", rb_Number_imag_u_le_n_p, 1);
+    
+    rb_define_method(rb_cNumber, "u_le_u?", rb_Number_real_u_le_u_p, 1);
+    rb_define_method(rb_cNumber, "real_u_le_u?", rb_Number_real_u_le_u_p, 1);
+    rb_define_method(rb_cNumber, "imag_u_le_u?", rb_Number_imag_u_le_u_p, 1);
+}