ruby-mpfr 0.0.7 → 0.0.8

data/History.txt

@@ -1,3 +1,13 @@
+=== 0.0.8
+* 2 major enhancements:
+  * Support MPFR v3.
+    * Add some methods.
+    * Note that ruby-mpfr does not support MPFR v2.
+  * Support Marshal.dump and Marshal.load.
+* 1 minor enhancements:
+  * Support rspec v2.
+  * some bug fixes.
+
 === 0.0.7
 * 1 bug fix:
   * Fix bugs of Methods MPFR::SquareMatrix, MPFR::ColumnVector, and MPFR::RowVector.

data/Manifest.txt

@@ -27,6 +27,7 @@ spec/mpfr/conversion_spec.rb
 spec/mpfr/exception_spec.rb
 spec/mpfr/functions_spec.rb
 spec/mpfr/generate_number_modulue.rb
+spec/mpfr/marshal_spec.rb
 spec/mpfr/precision_roundmode_spec.rb
 spec/mpfr/rounding_spec.rb
 spec/mpfr/set_value_spec.rb
@@ -35,9 +36,9 @@ spec/mpfr/string_spec.rb
 spec/mpfr_matrix/generate_matrix_arguments.rb
 spec/mpfr_matrix/mpfr_matrix_alloc_spec.rb
 spec/mpfr_matrix/mpfr_matrix_arithmetic_spec.rb
+spec/mpfr_matrix/mpfr_matrix_marshal_spec.rb
 spec/mpfr_matrix/mpfr_matrix_set_element_spec.rb
 spec/mpfr_matrix/mpfr_matrix_string_spec.rb
 spec/mpfr_matrix/mpfr_square_matrix_spec.rb
 spec/mpfr_matrix/spec_helper.rb
-spec/spec.opts
 tasks/extconf.rake

data/Rakefile

@@ -3,6 +3,7 @@ gem 'hoe', '>= 2.1.0'
 require 'hoe'
 require 'fileutils'
 require './lib/mpfr/version.rb'
+# require 'rake/clean'
 
 Hoe.plugin :newgem
 # Hoe.plugin :website

data/ext/mpfr/extconf.rb

@@ -2,8 +2,19 @@ require 'mkmf'
 
 $CFLAGS += " -Wall"
 
+REQUIRE_VERSION = 3
+
 dir_config("mpfr")
 dir_config("gmp")
 if have_header('mpfr.h') && have_library('mpfr') && have_header('gmp.h') && have_library('gmp')
-  create_makefile("mpfr")
+  if have_macro('MPFR_VERSION_MAJOR', 'mpfr.h') do |src|
+      src + <<SRC
+#if MPFR_VERSION_MAJOR < #{REQUIRE_VERSION}
+# error
+>>>>>> MPFR_VERSION_MAJOR requires larger than #{REQUIRE_VERSION} <<<<<<
+#endif
+SRC
+    end
+    create_makefile("mpfr")
+  end
 end

data/ext/mpfr/ruby_mpfr.c

@@ -1,10 +1,17 @@
 #include "ruby_mpfr.h"
 
+#define MPFR_DUMP_NUMBER 'A'
+#define MPFR_DUMP_PZERO  'B'
+#define MPFR_DUMP_MZERO  'C'
+#define MPFR_DUMP_PINF   'D'
+#define MPFR_DUMP_MINF   'E'
+#define MPFR_DUMP_NAN    'F'
+
 static ID eqq, to_s, new, class, method_defined, object_id;
 static VALUE __mpfr_class__, __sym_to_s__, __sym_to_str__;
 
 /* ------------------------------ Precision and Rounding Mode Start ------------------------------ */
-#define VALID_RND(rnd) (rnd >= 0 && rnd <= 3)
+#define VALID_RND(rnd_mode) (rnd_mode >= MPFR_RNDN && rnd_mode < ((MPFR_RNDA)+1))
 #define SPECIAL_FUNC_STATE "@@special_func_state"
 
 /* Convert VALUE rnd (rounding mode number) to C integer and */
@@ -456,6 +463,21 @@ static VALUE r_mpfr_nan(int argc, VALUE *argv, VALUE self)
   return val_ret;
 }
 
+/* Return zero. This method takes one optional argument meaning precision. */
+static VALUE r_mpfr_zero(int argc, VALUE *argv, VALUE self)
+{
+  mp_prec_t prec = r_mpfr_prec_from_optional_argument(0, 2, argc, argv);
+  MPFR *ptr_return;
+  VALUE val_ret;
+  int sign = 1;
+  if (argc >= 1) {
+    sign = NUM2INT(argv[0]);
+  }
+  r_mpfr_make_struct_init2(val_ret, ptr_return, prec);
+  mpfr_set_zero(ptr_return, sign);
+  return val_ret;
+}
+
 /* Return plus infinity. This method takes one optional argument meaning precision. */
 static VALUE r_mpfr_pinf(int argc, VALUE *argv, VALUE self)
 {
@@ -543,6 +565,21 @@ static VALUE r_mpfr_set_nan(VALUE self)
   return self;
 }
 
+/* Set the value of self to zero. */
+static VALUE r_mpfr_set_zero(int argc, VALUE *argv, VALUE self)
+{
+  int sign = 1;
+  MPFR *ptr_self;
+  r_mpfr_get_struct(ptr_self, self);
+  if (argc == 1) {
+    sign = NUM2INT(argv[0]);
+  } else if (argc > 1) {
+    rb_raise(rb_eArgError, "Invalid number of arguments.");
+  }
+  mpfr_set_zero(ptr_self, sign);
+  return self;
+}
+
 /* Swap the values self and p1 efficiently. p1 must be MPFR object. */
 static VALUE r_mpfr_swap(VALUE self, VALUE other)
 {
@@ -1197,6 +1234,14 @@ static VALUE r_mpfr_nonzero_p(VALUE self)
   return (mpfr_zero_p(ptr_self) == 0 ? Qtrue : Qfalse);
 }
 
+/* Return true if self is regular number, nil otherwise. */
+static VALUE r_mpfr_regular_p(VALUE self)
+{
+  MPFR *ptr_self;
+  r_mpfr_get_struct(ptr_self, self);
+  return (mpfr_regular_p(ptr_self) != 0 ? Qtrue : Qfalse);
+}
+
 /* mpfr_sgn(self). */
 static VALUE r_mpfr_sgn(VALUE self)
 {
@@ -1280,7 +1325,7 @@ static VALUE r_mpfr_unordered_p(VALUE self, VALUE other)
 /* ------------------------------ Integer Related Functions Start  ------------------------------ */
 
 /* mpfr_rint(ret, self, rnd) */
-static VALUE r_mpfr_m_rint(int argc, VALUE *argv, VALUE self)
+static VALUE r_mpfr_fr_rint(int argc, VALUE *argv, VALUE self)
 {
   mp_rnd_t rnd;
   mp_prec_t prec;
@@ -1294,7 +1339,7 @@ static VALUE r_mpfr_m_rint(int argc, VALUE *argv, VALUE self)
 }
 
 /* mpfr_ceil(ret, self) */
-static VALUE r_mpfr_m_ceil(int argc, VALUE *argv, VALUE self)
+static VALUE r_mpfr_fr_ceil(int argc, VALUE *argv, VALUE self)
 {
   mp_prec_t prec = r_mpfr_prec_from_optional_argument(0, 1, argc, argv);
   MPFR *ptr_self, *ptr_return;
@@ -1306,7 +1351,7 @@ static VALUE r_mpfr_m_ceil(int argc, VALUE *argv, VALUE self)
 }
 
 /* mpfr_floor(ret, self) */
-static VALUE r_mpfr_m_floor(int argc, VALUE *argv, VALUE self)
+static VALUE r_mpfr_fr_floor(int argc, VALUE *argv, VALUE self)
 {
   mp_prec_t prec = r_mpfr_prec_from_optional_argument(0, 1, argc, argv);
   MPFR *ptr_self, *ptr_return;
@@ -1318,7 +1363,7 @@ static VALUE r_mpfr_m_floor(int argc, VALUE *argv, VALUE self)
 }
 
 /* mpfr_round(ret, self) */
-static VALUE r_mpfr_m_round(int argc, VALUE *argv, VALUE self)
+static VALUE r_mpfr_fr_round(int argc, VALUE *argv, VALUE self)
 {
   mp_prec_t prec = r_mpfr_prec_from_optional_argument(0, 1, argc, argv);
   MPFR *ptr_self, *ptr_return;
@@ -1330,7 +1375,7 @@ static VALUE r_mpfr_m_round(int argc, VALUE *argv, VALUE self)
 }
 
 /* mpfr_trunc(ret, self) */
-static VALUE r_mpfr_m_trunc(int argc, VALUE *argv, VALUE self)
+static VALUE r_mpfr_fr_trunc(int argc, VALUE *argv, VALUE self)
 {
   mp_prec_t prec = r_mpfr_prec_from_optional_argument(0, 1, argc, argv);
   MPFR *ptr_self, *ptr_return;
@@ -1609,7 +1654,7 @@ static VALUE r_mpfr_prec_round(int argc, VALUE *argv, VALUE self)
   return val_ret;
 }
 
-/* mpfr_prec_round(ret, prec, rnd) */
+/* mpfr_prec_round(self, prec, rnd) */
 static VALUE r_mpfr_prec_round2(int argc, VALUE *argv, VALUE self)
 {
   mp_rnd_t rnd;
@@ -1621,6 +1666,26 @@ static VALUE r_mpfr_prec_round2(int argc, VALUE *argv, VALUE self)
   return self;
 }
 
+/* mpfr_can_round(self, err, rnd1, rnd2, prec) */
+static VALUE r_mpfr_can_round(VALUE self, VALUE err, VALUE rnd1, VALUE rnd2, VALUE prec)
+{
+  MPFR *ptr_self;
+  r_mpfr_get_struct(ptr_self, self);
+  if (mpfr_can_round(ptr_self, NUM2INT(err), NUM2INT(rnd1), NUM2INT(rnd2), NUM2INT(prec))) {
+    return Qtrue;
+  }
+  return Qnil;
+}
+
+/* mpfr_min_prec(self) */
+static VALUE r_mpfr_min_prec(VALUE self)
+{
+  MPFR *ptr_self;
+  r_mpfr_get_struct(ptr_self, self);
+  return INT2NUM(mpfr_min_prec(ptr_self));
+}
+
+
 /* ------------------------------ Rounding Mode Related Functions End ------------------------------ */
 
 /* ------------------------------ Special Functions Start ------------------------------ */
@@ -1960,6 +2025,51 @@ static VALUE r_mpfr_math_sinh_cosh(int argc, VALUE *argv, VALUE self)
   return rb_ary_new3(2, val_ret1, val_ret2);
 }
 
+/* mpfr_sech(ret, p1, rnd). */
+static VALUE r_mpfr_math_sech(int argc, VALUE *argv, VALUE self)
+{
+  mp_rnd_t rnd;
+  mp_prec_t prec;
+  r_mpfr_get_rnd_prec_from_optional_arguments(&rnd, &prec, 1, 3, argc, argv);
+  MPFR *ptr_arg1, *ptr_return;
+  VALUE val_ret;
+  volatile VALUE tmp_argv0 = r_mpfr_new_fr_obj(argv[0]);
+  r_mpfr_get_struct(ptr_arg1, tmp_argv0);
+  r_mpfr_make_struct_init2(val_ret, ptr_return, prec);
+  r_mpfr_set_special_func_state(mpfr_sech(ptr_return, ptr_arg1, rnd));
+  return val_ret;
+}
+
+/* mpfr_csch(ret, p1, rnd). */
+static VALUE r_mpfr_math_csch(int argc, VALUE *argv, VALUE self)
+{
+  mp_rnd_t rnd;
+  mp_prec_t prec;
+  r_mpfr_get_rnd_prec_from_optional_arguments(&rnd, &prec, 1, 3, argc, argv);
+  MPFR *ptr_arg1, *ptr_return;
+  VALUE val_ret;
+  volatile VALUE tmp_argv0 = r_mpfr_new_fr_obj(argv[0]);
+  r_mpfr_get_struct(ptr_arg1, tmp_argv0);
+  r_mpfr_make_struct_init2(val_ret, ptr_return, prec);
+  r_mpfr_set_special_func_state(mpfr_csch(ptr_return, ptr_arg1, rnd));
+  return val_ret;
+}
+
+/* mpfr_coth(ret, p1, rnd). */
+static VALUE r_mpfr_math_coth(int argc, VALUE *argv, VALUE self)
+{
+  mp_rnd_t rnd;
+  mp_prec_t prec;
+  r_mpfr_get_rnd_prec_from_optional_arguments(&rnd, &prec, 1, 3, argc, argv);
+  MPFR *ptr_arg1, *ptr_return;
+  VALUE val_ret;
+  volatile VALUE tmp_argv0 = r_mpfr_new_fr_obj(argv[0]);
+  r_mpfr_get_struct(ptr_arg1, tmp_argv0);
+  r_mpfr_make_struct_init2(val_ret, ptr_return, prec);
+  r_mpfr_set_special_func_state(mpfr_coth(ptr_return, ptr_arg1, rnd));
+  return val_ret;
+}
+
 /* mpfr_acosh(ret, p1, rnd). */
 static VALUE r_mpfr_math_acosh(int argc, VALUE *argv, VALUE self)
 {
@@ -2129,6 +2239,21 @@ static VALUE r_mpfr_math_lgamma(int argc, VALUE *argv, VALUE self)
   return rb_ary_new3(2, val_ret, INT2FIX(singp));
 }
 
+/* mpfr_digamma(ret, p1, rnd). */
+static VALUE r_mpfr_math_digamma(int argc, VALUE *argv, VALUE self)
+{
+  mp_rnd_t rnd;
+  mp_prec_t prec;
+  r_mpfr_get_rnd_prec_from_optional_arguments(&rnd, &prec, 1, 3, argc, argv);
+  MPFR *ptr_arg1, *ptr_return;
+  VALUE val_ret;
+  volatile VALUE tmp_argv0 = r_mpfr_new_fr_obj(argv[0]);
+  r_mpfr_get_struct(ptr_arg1, tmp_argv0);
+  r_mpfr_make_struct_init2(val_ret, ptr_return, prec);
+  r_mpfr_set_special_func_state(mpfr_digamma(ptr_return, ptr_arg1, rnd));
+  return val_ret;
+}
+
 /* mpfr_zeta(ret, p1, rnd). */
 static VALUE r_mpfr_math_zeta(int argc, VALUE *argv, VALUE self)
 {
@@ -2346,6 +2471,21 @@ static VALUE r_mpfr_math_hypot(int argc, VALUE *argv, VALUE self)
   return val_ret;
 }
 
+/* mpfr_ai(ret, p1, rnd). */
+static VALUE r_mpfr_math_ai(int argc, VALUE *argv, VALUE self)
+{
+  mp_rnd_t rnd;
+  mp_prec_t prec;
+  r_mpfr_get_rnd_prec_from_optional_arguments(&rnd, &prec, 1, 3, argc, argv);
+  MPFR *ptr_arg1, *ptr_return;
+  VALUE val_ret;
+  volatile VALUE tmp_argv0 = r_mpfr_new_fr_obj(argv[0]);
+  r_mpfr_get_struct(ptr_arg1, tmp_argv0);
+  r_mpfr_make_struct_init2(val_ret, ptr_return, prec);
+  r_mpfr_set_special_func_state(mpfr_ai(ptr_return, ptr_arg1, rnd));
+  return val_ret;
+}
+
 /* mpfr_const_log2(ret, rnd). */
 static VALUE r_mpfr_math_const_log2(int argc, VALUE *argv, VALUE self)
 {
@@ -2465,7 +2605,110 @@ static VALUE r_mpfr_math_max(int argc, VALUE *argv, VALUE self)
   return val_ret;
 }
 
+/* Version string which mpfr_get_version() returns. */
+static VALUE r_mpfr_get_version(VALUE self)
+{
+  return rb_str_new2(mpfr_get_version());
+}
+
+/* String which mpfr_get_patches() returns. */
+static VALUE r_mpfr_get_patches(VALUE self)
+{
+  return rb_str_new2(mpfr_get_patches());
+}
+
+/* Return true if MPFR was compiled as thread safe using compiler-level Thread Local Storage. Otherwise, nil. */
+static VALUE r_mpfr_buildopt_tls_p(VALUE self)
+{
+  return mpfr_buildopt_tls_p() == 0 ? Qnil : Qtrue;
+}
+
+/* Return true if MPFR was compiled with decimal float support. Otherwise, nil. */
+static VALUE r_mpfr_buildopt_decimal_p(VALUE self)
+{
+  return mpfr_buildopt_decimal_p() == 0 ? Qnil : Qtrue;
+}
+
+static VALUE r_mpfr_marshal_dump(VALUE self)
+{
+  MPFR *ptr_s;
+  r_mpfr_get_struct(ptr_s, self);
+  char *ret_str;
+  if (mpfr_regular_p(ptr_s)) {
+    mpz_t m;
+    mp_exp_t e;
+    mpz_init(m);
+    e = mpfr_get_z_2exp(m, ptr_s);
+    mpfr_asprintf(&ret_str, "%c%ld\t%ld\t%Zd", MPFR_DUMP_NUMBER, mpfr_get_prec(ptr_s), (long int)e, m);
+    mpz_clear(m);
+  } else {
+    char type;
+    if (mpfr_zero_p(ptr_s)) {
+      if (mpfr_sgn(ptr_s) >= 0) {
+	type = MPFR_DUMP_PZERO;
+      } else {
+	type = MPFR_DUMP_MZERO;
+      }
+    } else if (mpfr_nan_p(ptr_s)) {
+      type = MPFR_DUMP_NAN;
+    } else if (mpfr_sgn(ptr_s) >= 0) {
+      type = MPFR_DUMP_PINF;
+    } else {
+      type = MPFR_DUMP_MINF;
+    }
+    mpfr_asprintf(&ret_str, "%c%ld", type, mpfr_get_prec(ptr_s));
+  }
+  VALUE ret_val = rb_str_new2(ret_str);
+  mpfr_free_str(ret_str);
+  return ret_val;
+}
 
+static VALUE r_mpfr_marshal_load(VALUE self, VALUE dump_string)
+{
+  long int prec;
+  MPFR *ptr_s;
+  r_mpfr_get_struct(ptr_s, self);
+  Check_Type(dump_string, T_STRING);
+  char *dump, type;
+  dump = RSTRING_PTR(dump_string);
+  type = dump[0];
+  dump++;
+  if (type == MPFR_DUMP_NUMBER) {
+    mpz_t m;
+    long int e;
+    mpz_init(m);
+    sscanf(dump, "%ld\t%ld\t", &prec, &e);
+    int i;
+    i = 0;
+    while (i < 2) {
+      if (dump[0] == '\t') {
+	i++;
+      }
+      dump++;
+    }
+    mpz_set_str(m, dump, 10);
+    mpfr_init2(ptr_s, prec);
+    mpfr_set_z_2exp(ptr_s, m, e, MPFR_RNDN);
+    mpz_clear(m);
+  } else {
+    sscanf(dump, "%ld", &prec);
+    mpfr_init2(ptr_s, prec);
+    if (type == MPFR_DUMP_PZERO) {
+      mpfr_set_zero(ptr_s, +1);
+    } else if (type == MPFR_DUMP_MZERO){
+      mpfr_set_zero(ptr_s, -1);
+    } else if (type == MPFR_DUMP_NAN) {
+      mpfr_set_nan(ptr_s);
+    } else if (type == MPFR_DUMP_PINF) {
+      mpfr_set_inf(ptr_s, +1);
+    } else if (type == MPFR_DUMP_MINF) {
+      mpfr_set_inf(ptr_s, -1);
+    } else {
+      rb_raise(rb_eArgError, "Invalid dumped data for marshal_load.");
+    }
+  }
+  return self;
+}
 
 void Init_mpfr()
 {
@@ -2494,7 +2737,7 @@ void Init_mpfr()
     Except for some methods,
     the names of methods inherit that of MPFR functions in C language and
     each name of method is removed 'mpfr_' from that of correspnding MPFR function in C language.
-    The methods m_rint, m_ceil, m_floor, m_round and m_trunc are exceptions and
+    The methods fr_rint, fr_ceil, fr_floor, fr_round and fr_trunc are exceptions and
     respectively corresponts to mpfr_rint, mpfr_ceil, mpfr_floor, mpfr_round and mpfr_trunc.
     Note that there are also methods rint, ceil, floor, round and trunc
     which have respectively defferent actions
@@ -2535,41 +2778,44 @@ void Init_mpfr()
 
   /* ------------------------------ Constants Start ------------------------------ */
 
-  /* Version string which mpfr_get_version() returns. */
-  rb_define_const(r_mpfr_class, "MPFR_VERSION", rb_str_new2(mpfr_get_version()));
-  /* String which mpfr_get_patches() returns. */
-  rb_define_const(r_mpfr_class, "MPFR_PATCHES", rb_str_new2(mpfr_get_patches()));
-  /* Integer which is macro MPFR_VERSION. */
-  rb_define_const(r_mpfr_class, "MPFR_VERSION2", INT2NUM(MPFR_VERSION));
-  /* Integer which is macro MPFR_VERSION_MAJOR. */
+  /* Integer that is macro MPFR_VERSION. */
+  rb_define_const(r_mpfr_class, "MPFR_VERSION", INT2NUM(MPFR_VERSION));
+  /* Integer that is macro MPFR_VERSION_MAJOR. */
   rb_define_const(r_mpfr_class, "MPFR_VERSION_MAJOR", INT2NUM(MPFR_VERSION_MAJOR));
-  /* Integer which is macro MPFR_VERSION_MINOR. */
+  /* Integer that is macro MPFR_VERSION_MINOR. */
   rb_define_const(r_mpfr_class, "MPFR_VERSION_MINOR", INT2NUM(MPFR_VERSION_MINOR));
-  /* Integer which is macro MPFR_VERSION_PATCHLEVEL. */
+  /* Integer that is macro MPFR_VERSION_PATCHLEVEL. */
   rb_define_const(r_mpfr_class, "MPFR_VERSION_PATCHLEVEL", INT2NUM(MPFR_VERSION_PATCHLEVEL));
-  /* String whichi is macro MPFR_VERSION_STRING. */
+  /* String that is macro MPFR_VERSION_STRING. */
   rb_define_const(r_mpfr_class, "MPFR_VERSION_STRING", rb_str_new2(MPFR_VERSION_STRING));
 
-  /* Integer which is macro MPFR_PREC_MAX. */
+  /* Integer that is macro MPFR_PREC_MAX. */
   rb_define_const(r_mpfr_class, "PREC_MAX", INT2NUM(MPFR_PREC_MAX));
-  /* Integer which is macro MPFR_PREC_MIN. */
+  /* Integer that is macro MPFR_PREC_MIN. */
   rb_define_const(r_mpfr_class, "PREC_MIN", INT2NUM(MPFR_PREC_MIN));
-  /* Integer which is macro MPFR_EMAX_DEFAULT. */
+  /* Integer that is macro MPFR_EMAX_DEFAULT. */
   rb_define_const(r_mpfr_class, "EMAX_DEFAULT", INT2NUM(MPFR_EMAX_DEFAULT));
   /* Integer whichi is MPFR_EMIN_DEFAULT. */
   rb_define_const(r_mpfr_class, "EMIN_DEFAULT", INT2NUM(MPFR_EMIN_DEFAULT));
 
-  /* Integer which is macro GMP_RNDN. */
-  rb_define_const(r_mpfr_class, "RNDN", INT2NUM(GMP_RNDN));
-  /* Integer which is macro GMP_RNDZ. */
-  rb_define_const(r_mpfr_class, "RNDZ", INT2NUM(GMP_RNDZ));
-  /* Integer which is macro GMP_RNDU. */
-  rb_define_const(r_mpfr_class, "RNDU", INT2NUM(GMP_RNDU));
-  /* Integer which is macro GMP_RNDD. */
-  rb_define_const(r_mpfr_class, "RNDD", INT2NUM(GMP_RNDD));
+  /* Integer that is macro MPFR_RNDN. */
+  rb_define_const(r_mpfr_class, "RNDN", INT2NUM(MPFR_RNDN));
+  /* Integer that is macro MPFR_RNDZ. */
+  rb_define_const(r_mpfr_class, "RNDZ", INT2NUM(MPFR_RNDZ));
+  /* Integer that is macro MPFR_RNDU. */
+  rb_define_const(r_mpfr_class, "RNDU", INT2NUM(MPFR_RNDU));
+  /* Integer that is macro MPFR_RNDD. */
+  rb_define_const(r_mpfr_class, "RNDD", INT2NUM(MPFR_RNDD));
+  /* Integer that is macro MPFR_RNDD. */
+  rb_define_const(r_mpfr_class, "RNDA", INT2NUM(MPFR_RNDA));
 
   /* ------------------------------ Constants End ------------------------------ */
 
+  rb_define_singleton_method(r_mpfr_class, "get_version", r_mpfr_get_version, 0);
+  rb_define_singleton_method(r_mpfr_class, "get_patches", r_mpfr_get_patches, 0);
+  rb_define_singleton_method(r_mpfr_class, "buildopt_tls?", r_mpfr_buildopt_tls_p, 0);
+  rb_define_singleton_method(r_mpfr_class, "buildopt_decimal?", r_mpfr_buildopt_decimal_p, 0);
+
   /* ------------------------------ Precision and Rounding Mode Start ------------------------------ */
 
   rb_define_singleton_method(r_mpfr_class, "set_default_prec", r_mpfr_set_default_prec, 1);
@@ -2620,6 +2866,7 @@ void Init_mpfr()
   rb_define_method(r_mpfr_class, "coerce", r_mpfr_coerce, 1);
 
   rb_define_singleton_method(r_mpfr_class, "nan", r_mpfr_nan, -1);
+  rb_define_singleton_method(r_mpfr_class, "zero", r_mpfr_zero, -1);
   rb_define_singleton_method(r_mpfr_class, "pinf", r_mpfr_pinf, -1);
   rb_define_singleton_method(r_mpfr_class, "minf", r_mpfr_minf, -1);
 
@@ -2633,6 +2880,7 @@ void Init_mpfr()
   rb_define_method(r_mpfr_class, "set_fixnum_2exp", r_mpfr_set_fixnum_2exp, -1);
   rb_define_method(r_mpfr_class, "set_inf", r_mpfr_set_inf, 1);
   rb_define_method(r_mpfr_class, "set_nan", r_mpfr_set_nan, 0);
+  rb_define_method(r_mpfr_class, "set_zero", r_mpfr_set_zero, -1);
   rb_define_method(r_mpfr_class, "swap", r_mpfr_swap, 1);
 
   /* ------------------------------ Assignment Functions End ------------------------------ */
@@ -2642,6 +2890,9 @@ void Init_mpfr()
   rb_define_method(r_mpfr_class, "to_strf", r_mpfr_to_strf, 1);
   rb_define_method(r_mpfr_class, "to_s", r_mpfr_to_s, 0);
   rb_define_method(r_mpfr_class, "inspect", r_mpfr_inspect, 0);
+
+  rb_define_method(r_mpfr_class, "marshal_dump", r_mpfr_marshal_dump, 0);
+  rb_define_method(r_mpfr_class, "marshal_load", r_mpfr_marshal_load, 1);
   /* ------------------------------ Methods related to string End ------------------------------ */
 
   /* ------------------------------ Conversion functions Start ------------------------------ */
@@ -2699,16 +2950,17 @@ void Init_mpfr()
   rb_define_alias(r_mpfr_class, "finite?", "number_p");
   rb_define_alias(r_mpfr_class, "zero?", "zero_p");
   rb_define_method(r_mpfr_class, "nonzero?", r_mpfr_nonzero_p, 1);
+  rb_define_method(r_mpfr_class, "regular?", r_mpfr_regular_p, 1);
 
   /* ------------------------------ Comparison Functions Start ------------------------------ */
 
   /* ------------------------------ Integer Related Functions Start ------------------------------ */
 
-  rb_define_method(r_mpfr_class, "m_rint", r_mpfr_m_rint, -1);
-  rb_define_method(r_mpfr_class, "m_ceil", r_mpfr_m_ceil, -1);
-  rb_define_method(r_mpfr_class, "m_floor", r_mpfr_m_floor, -1);
-  rb_define_method(r_mpfr_class, "m_round", r_mpfr_m_round, -1);
-  rb_define_method(r_mpfr_class, "m_trunc", r_mpfr_m_trunc, -1);
+  rb_define_method(r_mpfr_class, "fr_rint", r_mpfr_fr_rint, -1);
+  rb_define_method(r_mpfr_class, "fr_ceil", r_mpfr_fr_ceil, -1);
+  rb_define_method(r_mpfr_class, "fr_floor", r_mpfr_fr_floor, -1);
+  rb_define_method(r_mpfr_class, "fr_round", r_mpfr_fr_round, -1);
+  rb_define_method(r_mpfr_class, "fr_trunc", r_mpfr_fr_trunc, -1);
   rb_define_method(r_mpfr_class, "rint_ceil", r_mpfr_rint_ceil, -1);
   rb_define_method(r_mpfr_class, "rint_floor", r_mpfr_rint_floor, -1);
   rb_define_method(r_mpfr_class, "rint_round", r_mpfr_rint_round, -1);
@@ -2743,6 +2995,8 @@ void Init_mpfr()
 
   rb_define_method(r_mpfr_class, "prec_round", r_mpfr_prec_round, -1);
   rb_define_method(r_mpfr_class, "prec_round!", r_mpfr_prec_round2, -1);
+  rb_define_method(r_mpfr_class, "can_round", r_mpfr_can_round, 4);
+  rb_define_method(r_mpfr_class, "min_prec", r_mpfr_min_prec, 0);
 
   /* ------------------------------ Rounding Mode Related Functions End ------------------------------ */
 
@@ -2813,6 +3067,9 @@ void Init_mpfr()
   rb_define_module_function(r_mpfr_math, "sinh", r_mpfr_math_sinh, -1);
   rb_define_module_function(r_mpfr_math, "tanh", r_mpfr_math_tanh, -1);
   rb_define_module_function(r_mpfr_math, "sinh_cosh", r_mpfr_math_sinh_cosh, -1);
+  rb_define_module_function(r_mpfr_math, "sech", r_mpfr_math_sech, -1);
+  rb_define_module_function(r_mpfr_math, "csch", r_mpfr_math_csch, -1);
+  rb_define_module_function(r_mpfr_math, "coth", r_mpfr_math_coth, -1);
   rb_define_module_function(r_mpfr_math, "acosh", r_mpfr_math_acosh, -1);
   rb_define_module_function(r_mpfr_math, "asinh", r_mpfr_math_asinh, -1);
   rb_define_module_function(r_mpfr_math, "atanh", r_mpfr_math_atanh, -1);
@@ -2824,6 +3081,7 @@ void Init_mpfr()
   rb_define_module_function(r_mpfr_math, "gamma", r_mpfr_math_gamma, -1);
   rb_define_module_function(r_mpfr_math, "lngamma", r_mpfr_math_lngamma, -1);
   rb_define_module_function(r_mpfr_math, "lgamma", r_mpfr_math_lgamma, -1);
+  rb_define_module_function(r_mpfr_math, "digamma", r_mpfr_math_digamma, -1);
   rb_define_module_function(r_mpfr_math, "zeta", r_mpfr_math_zeta, -1);
   rb_define_module_function(r_mpfr_math, "erf", r_mpfr_math_erf, -1);
   rb_define_module_function(r_mpfr_math, "erfc", r_mpfr_math_erfc, -1);
@@ -2837,6 +3095,7 @@ void Init_mpfr()
   rb_define_module_function(r_mpfr_math, "fms", r_mpfr_math_fms, -1);
   rb_define_module_function(r_mpfr_math, "agm", r_mpfr_math_agm, -1);
   rb_define_module_function(r_mpfr_math, "hypot", r_mpfr_math_hypot, -1);
+  rb_define_module_function(r_mpfr_math, "ai", r_mpfr_math_ai, -1);
   rb_define_module_function(r_mpfr_math, "const_log2", r_mpfr_math_const_log2, -1);
   rb_define_module_function(r_mpfr_math, "const_pi", r_mpfr_math_const_pi, -1);
   rb_define_module_function(r_mpfr_math, "const_euler", r_mpfr_math_const_euler, -1);
@@ -2867,3 +3126,4 @@ void Init_mpfr()
   __sym_to_str__ = rb_eval_string(":to_str");
 
 }
+

data/ext/mpfr_matrix/mpfr/extconf.rb

@@ -1,7 +1,21 @@
 require 'mkmf'
 
+$CFLAGS += " -Wall"
+
+REQUIRE_VERSION = 3
+
 dir_config("mpfr")
 dir_config("gmp")
+
 if have_header('mpfr.h') && have_library('mpfr') && have_header('gmp.h') && have_library('gmp')
-  create_makefile("mpfr/matrix")
+  if have_macro('MPFR_VERSION_MAJOR', 'mpfr.h') do |src|
+      src + <<SRC
+#if MPFR_VERSION_MAJOR < #{REQUIRE_VERSION}
+# error
+>>>>>> MPFR_VERSION_MAJOR requires larger than #{REQUIRE_VERSION} <<<<<<
+#endif
+SRC
+    end
+    create_makefile("mpfr/matrix")
+  end
 end

data/ext/mpfr_matrix/mpfr/ruby_mpfr_matrix.c

@@ -1,5 +1,12 @@
 #include "ruby_mpfr_matrix.h"
 
+#define MPFR_DUMP_NUMBER 'A'
+#define MPFR_DUMP_PZERO  'B'
+#define MPFR_DUMP_MZERO  'C'
+#define MPFR_DUMP_PINF   'D'
+#define MPFR_DUMP_MINF   'E'
+#define MPFR_DUMP_NAN    'F'
+
 static ID eqq;
 
 void r_mpfr_matrix_free(void *ptr)
@@ -185,6 +192,108 @@ static VALUE r_mpfr_matrix_initialize_copy (VALUE self, VALUE other)
   return Qtrue;
 }
 
+static VALUE r_mpfr_matrix_marshal_dump(VALUE self)
+{
+  MPFRMatrix *ptr;
+  r_mpfr_get_matrix_struct(ptr, self);
+  mpz_t m;
+  mp_exp_t e;
+  mpz_init(m);
+  int i;
+  MPFR *ptr_el;
+  char *tmp_str, type;
+  VALUE ret_ary;
+  ret_ary = rb_ary_new();
+  rb_ary_push(ret_ary, INT2FIX(ptr->row));
+  rb_ary_push(ret_ary, INT2FIX(ptr->column));
+
+  for (i = 0; i < ptr->size; i++) {
+    ptr_el = ptr->data + i;
+    if (mpfr_regular_p(ptr_el)) {
+      e = mpfr_get_z_2exp(m, ptr_el);
+      mpfr_asprintf(&tmp_str, "%c%ld\t%ld\t%Zd", MPFR_DUMP_NUMBER, mpfr_get_prec(ptr_el), (long int)e, m);
+    } else {
+      if (mpfr_zero_p(ptr_el)) {
+	if (mpfr_sgn(ptr_el) >= 0) {
+	  type = MPFR_DUMP_PZERO;
+	} else {
+	  type = MPFR_DUMP_MZERO;
+	}
+      } else if (mpfr_nan_p(ptr_el)) {
+	type = MPFR_DUMP_NAN;
+      } else if (mpfr_sgn(ptr_el) >= 0) {
+	type = MPFR_DUMP_PINF;
+      } else {
+	type = MPFR_DUMP_MINF;
+      }
+      mpfr_asprintf(&tmp_str, "%c%ld", type, mpfr_get_prec(ptr_el));
+    }
+    rb_ary_push(ret_ary, rb_str_new2(tmp_str));
+    mpfr_free_str(tmp_str);
+  }
+
+  mpz_clear(m);
+  return ret_ary;
+}
+
+static VALUE r_mpfr_matrix_marshal_load(VALUE self, VALUE dump_ary)
+{
+  MPFRMatrix *ptr;
+  r_mpfr_get_matrix_struct(ptr, self);
+
+  ptr->row = NUM2INT(rb_ary_entry(dump_ary, 0));
+  ptr->column = NUM2INT(rb_ary_entry(dump_ary, 1));
+  ptr->size = ptr->row * ptr->column;
+  ptr->data = ALLOC_N(MPFR, ptr->size);
+  int i, j;
+  long int prec, e;
+  MPFR *ptr_el;
+  char *dump, type;
+  VALUE dump_element;
+  mpz_t m;
+
+  for(i = 0; i < ptr->size; i++){
+    dump_element = rb_ary_entry(dump_ary, i + 2);
+    ptr_el = ptr->data + i;
+    Check_Type(dump_element, T_STRING);
+    dump = RSTRING_PTR(dump_element);
+    type = dump[0];
+    dump++;
+    if (type == MPFR_DUMP_NUMBER) {
+      mpz_init(m);
+      sscanf(dump, "%ld\t%ld\t", &prec, &e);
+      j = 0;
+      while (j < 2) {
+	if (dump[0] == '\t') {
+	  j++;
+	}
+	dump++;
+      }
+      mpz_set_str(m, dump, 10);
+      mpfr_init2(ptr_el, prec);
+      mpfr_set_z_2exp(ptr_el, m, e, MPFR_RNDN);
+      mpz_clear(m);
+    } else {
+      sscanf(dump, "%ld", &prec);
+      mpfr_init2(ptr_el, prec);
+      if (type == MPFR_DUMP_PZERO) {
+	mpfr_set_zero(ptr_el, +1);
+      } else if (type == MPFR_DUMP_MZERO){
+	mpfr_set_zero(ptr_el, -1);
+      } else if (type == MPFR_DUMP_NAN) {
+	mpfr_set_nan(ptr_el);
+      } else if (type == MPFR_DUMP_PINF) {
+	mpfr_set_inf(ptr_el, +1);
+      } else if (type == MPFR_DUMP_MINF) {
+	mpfr_set_inf(ptr_el, -1);
+      } else {
+	rb_raise(rb_eArgError, "Invalid dumped data for marshal_load.");
+      }
+    }
+  }
+  return self;
+}
+
 /* Allocation function for MPFR::ColumnVector. */
 static VALUE r_mpfr_col_vector_alloc (VALUE self)
 {
@@ -273,7 +382,7 @@ static VALUE r_mpfr_row_vector_global_new(int argc, VALUE arg)
 {
   MPFRMatrix *ptr;
   VALUE val;
-  r_mpfr_make_row_struct(val, ptr);
+  r_mpfr_make_row_vector_struct(val, ptr);
   r_mpfr_row_vector_set_initial_value(ptr, arg);
   return val;
 }
@@ -1031,6 +1140,9 @@ void Init_matrix()
   rb_define_private_method(r_mpfr_matrix, "initialize", r_mpfr_matrix_initialize, -1);
   rb_define_private_method(r_mpfr_matrix, "initialize_copy", r_mpfr_matrix_initialize_copy, 1);
 
+  rb_define_method(r_mpfr_matrix, "marshal_dump", r_mpfr_matrix_marshal_dump, 0);
+  rb_define_method(r_mpfr_matrix, "marshal_load", r_mpfr_matrix_marshal_load, 1);
+
   rb_define_singleton_method(tmp_r_mpfr_class, "SquareMatrix", r_mpfr_square_matrix_global_new, 1);
   rb_define_alloc_func(r_mpfr_square_matrix, r_mpfr_square_matrix_alloc);
   rb_define_private_method(r_mpfr_square_matrix, "initialize", r_mpfr_square_matrix_initialize, 1);

data/lib/mpfr/version.rb

@@ -1 +1 @@
-RUBY_MPFR_VERSION = '0.0.7'
+RUBY_MPFR_VERSION = '0.0.8'

data/spec/mpfr/allocate_spec.rb

@@ -1,4 +1,4 @@
-require File.dirname(__FILE__) + '/spec_helper.rb'
+require File.expand_path(File.dirname(__FILE__)) + '/spec_helper.rb'
 
 describe MPFR, 'when allocating objects from integer' do
   it "should equal bignum" do

data/spec/mpfr/arithmetic_spec.rb

@@ -1,4 +1,4 @@
-require File.dirname(__FILE__) + '/spec_helper.rb'
+require File.expand_path(File.dirname(__FILE__)) + '/spec_helper.rb'
 
 MPFR.set_default_prec(256)
 

data/spec/mpfr/comparison_spec.rb

@@ -1,4 +1,4 @@
-require File.dirname(__FILE__) + '/spec_helper.rb'
+require File.expand_path(File.dirname(__FILE__)) + '/spec_helper.rb'
 
 MPFR.set_default_prec(256)
 

data/spec/mpfr/constant_spec.rb

@@ -1,10 +1,8 @@
-require File.dirname(__FILE__) + '/spec_helper.rb'
+require File.expand_path(File.dirname(__FILE__)) + '/spec_helper.rb'
 
 describe MPFR, "when calling constant of MPFR" do
   it "should be defined" do
     MPFR::MPFR_VERSION.should be_true
-    MPFR::MPFR_PATCHES.should be_true
-    MPFR::MPFR_VERSION2.should be_true
     MPFR::MPFR_VERSION_MAJOR.should be_true
     MPFR::MPFR_VERSION_MINOR.should be_true
     MPFR::MPFR_VERSION_PATCHLEVEL.should be_true

data/spec/mpfr/conversion_spec.rb

@@ -1,4 +1,4 @@
-require File.dirname(__FILE__) + '/spec_helper.rb'
+require File.expand_path(File.dirname(__FILE__)) + '/spec_helper.rb'
 
 MPFR.set_default_prec(90)
 

data/spec/mpfr/exception_spec.rb

@@ -1,4 +1,4 @@
-require File.dirname(__FILE__) + '/spec_helper.rb'
+require File.expand_path(File.dirname(__FILE__)) + '/spec_helper.rb'
 
 describe MPFR, "when executing methods which are exception related functions" do
   it "should get and set emin and emax" do

data/spec/mpfr/functions_spec.rb

@@ -1,4 +1,4 @@
-require File.dirname(__FILE__) + '/spec_helper.rb'
+require File.expand_path(File.dirname(__FILE__)) + '/spec_helper.rb'
 
 MPFR.set_default_prec(256)
 

data/spec/mpfr/marshal_spec.rb

@@ -0,0 +1,57 @@
+require File.expand_path(File.dirname(__FILE__)) + '/spec_helper.rb'
+
+describe MPFR do
+  context "when marshaling" do
+    before(:all) do
+      MPFR.set_default_prec(256)
+    end
+
+    it "should restore nan" do
+      n = MPFR.nan(128)
+      n2 = Marshal.load(Marshal.dump(n))
+      n2.nan?.should be_true
+      n2.get_prec.should == 128
+    end
+
+    it "should restore plus infinity" do
+      n = MPFR.pinf(128)
+      n2 = Marshal.load(Marshal.dump(n))
+      n2.infinite?.should be_true
+      n2.get_prec.should == 128
+      n2.sgn.should > 0
+    end
+
+    it "should restore minus infinity" do
+      n = MPFR.minf(128)
+      n2 = Marshal.load(Marshal.dump(n))
+      n2.infinite?.should be_true
+      n2.get_prec.should == 128
+      n2.sgn.should < 0
+    end
+
+    it "should restore plus zero" do
+      n = MPFR.zero(1, 128)
+      n2 = Marshal.load(Marshal.dump(n))
+      n2.zero?.should be_true
+      n2.get_prec.should == 128
+      n2.sgn.should >= 0
+    end
+
+    it "should restore minus zero" do
+      n = MPFR.zero(-1, 128)
+      n2 = Marshal.load(Marshal.dump(n))
+      n2.zero?.should be_true
+      n2.get_prec.should == 128
+      n2.sgn.should <= 0
+    end
+
+    it "should restore regular numbers" do
+      GenerateNumber.string(1000).each do |s|
+        n = MPFR(s, MPFR::RNDN, 128)
+        n2 = Marshal.load(Marshal.dump(n))
+        n2.should == n
+        n2.get_prec.should == 128
+      end
+    end
+  end
+end

data/spec/mpfr/precision_roundmode_spec.rb

@@ -1,4 +1,4 @@
-require File.dirname(__FILE__) + '/spec_helper.rb'
+require File.expand_path(File.dirname(__FILE__)) + '/spec_helper.rb'
 
 def check_rnd(a)
   MPFR.new(a).should == MPFR.new(a, MPFR.get_default_rounding_mode)
@@ -6,9 +6,11 @@ def check_rnd(a)
   z = MPFR.new(a, MPFR::RNDZ)
   u = MPFR.new(a, MPFR::RNDU)
   d = MPFR.new(a, MPFR::RNDD)
+  a = MPFR.new(a, MPFR::RNDA)
   (n > 0 ? d : u).should == z
   [(a - z).abs, (a - u).abs, (a - d).abs].min.should >= (a - n).abs
   (d <= u && (n == z || n == u || n == d)).should be_true
+  z.abs.should <= a.abs
 end
 
 def check_set_rnd(rnd1, rnd2)
@@ -24,11 +26,11 @@ describe MPFR, 'when getting rounding mode' do
     args = GenerateNumber.float(100)
     args += [Math::PI, -Math::PI, Math::E, -Math::E, 0, 10, 27, -9, -293578]
     args.each{ |a| check_rnd(a) }
-    [MPFR::RNDN, MPFR::RNDZ, MPFR::RNDU, MPFR::RNDD].each{ |rnd| args.each{ |a| check_rnd(a) } }
+    [MPFR::RNDN, MPFR::RNDZ, MPFR::RNDU, MPFR::RNDD, MPFR::RNDA].each{ |rnd| args.each{ |a| check_rnd(a) } }
   end
 
   it "should be set rounding mode" do
-    rnds = [MPFR::RNDN, MPFR::RNDZ, MPFR::RNDU, MPFR::RNDD]
+    rnds = [MPFR::RNDN, MPFR::RNDZ, MPFR::RNDU, MPFR::RNDD, MPFR::RNDA]
     for i in 0...rnds.size
       for j in (i+1)...rnds.size
         check_set_rnd(rnds[i], rnds[j])

data/spec/mpfr/rounding_spec.rb

@@ -1,4 +1,4 @@
-require File.dirname(__FILE__) + '/spec_helper.rb'
+require File.expand_path(File.dirname(__FILE__)) + '/spec_helper.rb'
 
 describe MPFR, "when rounding number by prec_round" do
   before(:all) do
@@ -38,7 +38,7 @@ describe MPFR, "when rounding number by prec_round" do
 
   it "should be rounded by prec_round!" do
     @args.each do |a|
-      [MPFR::RNDU, MPFR::RNDD, MPFR::RNDN, MPFR::RNDZ].each do |rnd|
+      [MPFR::RNDU, MPFR::RNDD, MPFR::RNDN, MPFR::RNDZ, MPFR::RNDA].each do |rnd|
         b = a.dup
         old_object_id = b.object_id
         b.prec_round!(rnd, @prec)

data/spec/mpfr/set_value_spec.rb

@@ -1,4 +1,4 @@
-require File.dirname(__FILE__) + '/spec_helper.rb'
+require File.expand_path(File.dirname(__FILE__)) + '/spec_helper.rb'
 
 MPFR.set_default_prec(300)
 

data/spec/mpfr/spec_helper.rb

@@ -1,13 +1,13 @@
 begin
-  require 'spec'
+  require 'rspec'
 rescue LoadError
   require 'rubygems' unless ENV['NO_RUBYGEMS']
   gem 'rspec'
-  require 'spec'
+  require 'rspec'
 end
 
-$:.unshift(File.dirname(__FILE__) + '/../../lib')
-$:.unshift(File.dirname(__FILE__) + '/../../ext')
-$:.unshift(File.dirname(__FILE__))
-require 'mpfr'
+$:.unshift(File.expand_path(File.dirname(__FILE__)) + '/../../lib')
+$:.unshift(File.expand_path(File.dirname(__FILE__)) + '/../../ext')
+$:.unshift(File.expand_path(File.dirname(__FILE__)))
+require 'mpfr/mpfr'
 require "generate_number_modulue"

data/spec/mpfr/string_spec.rb

@@ -1,4 +1,4 @@
-require File.dirname(__FILE__) + '/spec_helper.rb'
+require File.expand_path(File.dirname(__FILE__)) + '/spec_helper.rb'
 
 describe MPFR, "when converted to string" do
   before(:all) do

data/spec/mpfr_matrix/mpfr_matrix_alloc_spec.rb

@@ -1,4 +1,4 @@
-require File.dirname(__FILE__) + '/spec_helper.rb'
+require File.expand_path(File.dirname(__FILE__)) + '/spec_helper.rb'
 
 describe "initialization of matrix" do
   before(:all) do

data/spec/mpfr_matrix/mpfr_matrix_arithmetic_spec.rb

@@ -1,4 +1,4 @@
-require File.dirname(__FILE__) + '/spec_helper.rb'
+require File.expand_path(File.dirname(__FILE__)) + '/spec_helper.rb'
 
 describe MPFR::Matrix, "when making matrix negative" do
   before(:all) do

data/spec/mpfr_matrix/mpfr_matrix_marshal_spec.rb

@@ -0,0 +1,65 @@
+require File.expand_path(File.dirname(__FILE__)) + '/spec_helper.rb'
+
+describe MPFR::Matrix, "when marshaling" do
+  before(:all) do
+    MPFR.set_default_prec(128)
+  end
+
+  it "should be restore" do
+    row = 8
+    column = 3
+    args = GenerateNumber.float_matrix_arguments(1000, row, column).map{ |a| MPFR::Matrix.new(a) }
+    args.each do |m|
+      m2 = Marshal.load(Marshal.dump(m))
+      m2.should == m
+    end
+  end
+end
+
+describe MPFR::SquareMatrix, "when marshaling" do
+  before(:all) do
+    MPFR.set_default_prec(128)
+  end
+
+  it "should be restore" do
+    row = 4
+    column = 4
+    args = GenerateNumber.float_matrix_arguments(1000, row, column).map{ |a| MPFR::SquareMatrix.new(a) }
+    args.each do |m|
+      m2 = Marshal.load(Marshal.dump(m))
+      m2.should == m
+    end
+  end
+end
+
+describe MPFR::ColumnVector, "when marshaling" do
+  before(:all) do
+    MPFR.set_default_prec(128)
+  end
+
+  it "should be restore" do
+    row = 1
+    column = 6
+    args = GenerateNumber.float_matrix_arguments(1000, row, column).map{ |a| MPFR::ColumnVector.new(a[0]) }
+    args.each do |m|
+      m2 = Marshal.load(Marshal.dump(m))
+      m2.should == m
+    end
+  end
+end
+
+describe MPFR::RowVector, "when marshaling" do
+  before(:all) do
+    MPFR.set_default_prec(128)
+  end
+
+  it "should be restore" do
+    row = 1
+    column = 6
+    args = GenerateNumber.float_matrix_arguments(1000, row, column).map{ |a| MPFR::RowVector.new(a[0]) }
+    args.each do |m|
+      m2 = Marshal.load(Marshal.dump(m))
+      m2.should == m
+    end
+  end
+end

data/spec/mpfr_matrix/mpfr_matrix_set_element_spec.rb

@@ -1,4 +1,4 @@
-require File.dirname(__FILE__) + '/spec_helper.rb'
+require File.expand_path(File.dirname(__FILE__)) + '/spec_helper.rb'
 
 describe MPFR::Matrix, "when setting number to particular element" do
   before(:all) do

data/spec/mpfr_matrix/mpfr_matrix_string_spec.rb

@@ -1,4 +1,4 @@
-require File.dirname(__FILE__) + '/spec_helper.rb'
+require File.expand_path(File.dirname(__FILE__)) + '/spec_helper.rb'
 
 describe MPFR::Matrix, "when MPFR instance is converted to string by some methods." do
   

data/spec/mpfr_matrix/mpfr_square_matrix_spec.rb

@@ -1,4 +1,4 @@
-require File.dirname(__FILE__) + '/spec_helper.rb'
+require File.expand_path(File.dirname(__FILE__)) + '/spec_helper.rb'
 
 describe MPFR::SquareMatrix, "when calculating LU decomposition" do
   before(:all) do

data/spec/mpfr_matrix/spec_helper.rb

@@ -1,9 +1,9 @@
 begin
-  require 'spec'
+  require 'rspec'
 rescue LoadError
   require 'rubygems' unless ENV['NO_RUBYGEMS']
   gem 'rspec'
-  require 'spec'
+  require 'rspec'
 end
 
 $:.unshift(File.dirname(__FILE__) + '/../../lib')

data/tasks/extconf.rake

@@ -1,18 +1,18 @@
-begin
-  require 'spec'
-rescue LoadError
-  require 'rubygems' unless ENV['NO_RUBYGEMS']
-  require 'spec'
-end
-begin
-  require 'spec/rake/spectask'
-rescue LoadError
-  puts <<-EOS
-To use rspec for testing you must install rspec gem:
-    gem install rspec
-EOS
-  exit(0)
-end
+# begin
+#   require 'rspec'
+# rescue LoadError
+#   require 'rubygems' unless ENV['NO_RUBYGEMS']
+#   require 'rspec'
+# end
+# begin
+#   require 'rspec/rake/spectask'
+# rescue LoadError
+#   puts <<-EOS
+# To use rspec for testing you must install rspec gem:
+#     gem install rspec
+# EOS
+#   exit(0)
+# end
 
 desc "Run 'make realclean' for extended libraries"
 task "ext:realclean" do

metadata

@@ -1,7 +1,12 @@
 --- !ruby/object:Gem::Specification 
 name: ruby-mpfr
 version: !ruby/object:Gem::Version 
-  version: 0.0.7
+  prerelease: false
+  segments: 
+  - 0
+  - 0
+  - 8
+  version: 0.0.8
 platform: ruby
 authors: 
 - Takayuki YAMAGUCHI
@@ -9,39 +14,39 @@ autorequire:
 bindir: bin
 cert_chain: []
 
-date: 2010-01-21 00:00:00 +09:00
+date: 2010-11-14 00:00:00 +09:00
 default_executable: 
 dependencies: 
 - !ruby/object:Gem::Dependency 
   name: rubyforge
-  type: :development
-  version_requirement: 
-  version_requirements: !ruby/object:Gem::Requirement 
+  prerelease: false
+  requirement: &id001 !ruby/object:Gem::Requirement 
+    none: false
     requirements: 
     - - ">="
       - !ruby/object:Gem::Version 
-        version: 2.0.3
-    version: 
-- !ruby/object:Gem::Dependency 
-  name: gemcutter
+        segments: 
+        - 2
+        - 0
+        - 4
+        version: 2.0.4
   type: :development
-  version_requirement: 
-  version_requirements: !ruby/object:Gem::Requirement 
-    requirements: 
-    - - ">="
-      - !ruby/object:Gem::Version 
-        version: 0.3.0
-    version: 
+  version_requirements: *id001
 - !ruby/object:Gem::Dependency 
   name: hoe
-  type: :development
-  version_requirement: 
-  version_requirements: !ruby/object:Gem::Requirement 
+  prerelease: false
+  requirement: &id002 !ruby/object:Gem::Requirement 
+    none: false
     requirements: 
     - - ">="
       - !ruby/object:Gem::Version 
-        version: 2.5.0
-    version: 
+        segments: 
+        - 2
+        - 6
+        - 2
+        version: 2.6.2
+  type: :development
+  version_requirements: *id002
 description: |-
   ruby-mpfr is a library to use MPFR[http://www.mpfr.org/] which is a C library for
   multiple-precision floating-point computations.
@@ -93,6 +98,7 @@ files:
 - spec/mpfr/exception_spec.rb
 - spec/mpfr/functions_spec.rb
 - spec/mpfr/generate_number_modulue.rb
+- spec/mpfr/marshal_spec.rb
 - spec/mpfr/precision_roundmode_spec.rb
 - spec/mpfr/rounding_spec.rb
 - spec/mpfr/set_value_spec.rb
@@ -101,11 +107,11 @@ files:
 - spec/mpfr_matrix/generate_matrix_arguments.rb
 - spec/mpfr_matrix/mpfr_matrix_alloc_spec.rb
 - spec/mpfr_matrix/mpfr_matrix_arithmetic_spec.rb
+- spec/mpfr_matrix/mpfr_matrix_marshal_spec.rb
 - spec/mpfr_matrix/mpfr_matrix_set_element_spec.rb
 - spec/mpfr_matrix/mpfr_matrix_string_spec.rb
 - spec/mpfr_matrix/mpfr_square_matrix_spec.rb
 - spec/mpfr_matrix/spec_helper.rb
-- spec/spec.opts
 - tasks/extconf.rake
 has_rdoc: true
 homepage: http://rubyforge.org/projects/ruby-mpfr/
@@ -117,23 +123,26 @@ rdoc_options:
 - README.rdoc
 require_paths: 
 - lib
-- ext
 required_ruby_version: !ruby/object:Gem::Requirement 
+  none: false
   requirements: 
   - - ">="
     - !ruby/object:Gem::Version 
+      segments: 
+      - 0
       version: "0"
-  version: 
 required_rubygems_version: !ruby/object:Gem::Requirement 
+  none: false
   requirements: 
   - - ">="
     - !ruby/object:Gem::Version 
+      segments: 
+      - 0
       version: "0"
-  version: 
 requirements: []
 
 rubyforge_project: ruby-mpfr
-rubygems_version: 1.3.5
+rubygems_version: 1.3.7
 signing_key: 
 specification_version: 3
 summary: ruby-mpfr is a library to use MPFR[http://www.mpfr.org/] which is a C library for multiple-precision floating-point computations

data/spec/spec.opts

@@ -1 +0,0 @@
---colour