Real Dobule

Author: rajithv

ext/symengine/symengine.c

@@ -62,6 +62,10 @@ void Init_symengine() {
     rb_define_method(c_integer, "initialize", cinteger_init, 1);
     rb_define_method(c_integer, "%", cntheory_mod, 1);
 
+    //RealDouble Class
+    c_real_double = rb_define_class_under(m_symengine, "RealDouble", c_basic);
+    rb_define_alloc_func(c_double, cbasic_alloc);
+
     //Rational class
     c_rational = rb_define_class_under(m_symengine, "Rational", c_basic);
     rb_define_alloc_func(c_rational, cbasic_alloc);

ext/symengine/symengine.h

@@ -10,6 +10,7 @@ VALUE m_symengine;
 VALUE c_basic;
 VALUE c_symbol;
 VALUE c_integer;
+VALUE c_real_double;
 VALUE c_rational;
 VALUE c_complex;
 VALUE c_complex_double;

ext/symengine/symengine_utils.c

@@ -6,13 +6,22 @@ void sympify(VALUE operand2, basic_struct *cbasic_operand2) {
     basic_struct *temp;
     VALUE new_operand2, num, den;
     VALUE real, imag;
+    VALUE Rb_Temp_String;
+    char *s;
 
     switch(TYPE(operand2)) {
         case T_FIXNUM:
         case T_BIGNUM:
             GET_SYMINTFROMVAL(operand2, cbasic_operand2);
             break;
 
+        case T_FLOAT:
+            Rb_Temp_String = rb_funcall(operand2, rb_intern("to_s"), 0, NULL);
+            s = StringValueCStr(Rb_Temp_String);
+            real_double_set_d(this, atof(s));
+
+            break;
+
         case T_RATIONAL:
             num = rb_funcall(operand2, rb_intern("numerator"), 0, NULL);
             den = rb_funcall(operand2, rb_intern("denominator"), 0, NULL);
@@ -66,6 +75,8 @@ VALUE Klass_of_Basic(const basic_struct *basic_ptr) {
             return c_symbol;
         case SYMENGINE_INTEGER:
             return c_integer;
+        case SYMENGINE_REAL_DOUBLE:
+            return c_real_double;
         case SYMENGINE_RATIONAL:
             return c_rational;
         case SYMENGINE_COMPLEX: