crmf 0.1.1 → 0.1.2

data/ext/crmf/crmf.c

@@ -1581,7 +1581,7 @@ static VALUE crmf_add_rz(VALUE self, VALUE x, VALUE y)
 	double dx = NUM2DBL(x);
 	double dy = NUM2DBL(y);
 	fesetround(FE_TOWARDZERO);
-	double dr = dx + dy;
+	volatile double dr = dx + dy;
 	fesetround(FE_TONEAREST);
 	return DBL2NUM(dr);
 }
@@ -1591,7 +1591,7 @@ static VALUE crmf_add_rd(VALUE self, VALUE x, VALUE y)
 	double dx = NUM2DBL(x);
 	double dy = NUM2DBL(y);
 	fesetround(FE_DOWNWARD);
-	double dr = dx + dy;
+	volatile double dr = dx + dy;
 	fesetround(FE_TONEAREST);
 	return DBL2NUM(dr);
 }
@@ -1601,7 +1601,7 @@ static VALUE crmf_add_ru(VALUE self, VALUE x, VALUE y)
 	double dx = NUM2DBL(x);
 	double dy = NUM2DBL(y);
 	fesetround(FE_UPWARD);
-	double dr = dx + dy;
+	volatile double dr = dx + dy;
 	fesetround(FE_TONEAREST);
 	return DBL2NUM(dr);
 }
@@ -1611,7 +1611,7 @@ static VALUE crmf_add_rn(VALUE self, VALUE x, VALUE y)
 	double dx = NUM2DBL(x);
 	double dy = NUM2DBL(y);
 	fesetround(FE_TONEAREST);
-	double dr = dx + dy;
+	volatile double dr = dx + dy;
 	return DBL2NUM(dr);
 }
 
@@ -1621,7 +1621,7 @@ static VALUE crmf_sub_rz(VALUE self, VALUE x, VALUE y)
 	double dx = NUM2DBL(x);
 	double dy = NUM2DBL(y);
 	fesetround(FE_TOWARDZERO);
-	double dr = dx - dy;
+	volatile double dr = dx - dy;
 	fesetround(FE_TONEAREST);
 	return DBL2NUM(dr);
 }
@@ -1631,7 +1631,7 @@ static VALUE crmf_sub_rd(VALUE self, VALUE x, VALUE y)
 	double dx = NUM2DBL(x);
 	double dy = NUM2DBL(y);
 	fesetround(FE_DOWNWARD);
-	double dr = dx - dy;
+	volatile double dr = dx - dy;
 	fesetround(FE_TONEAREST);
 	return DBL2NUM(dr);
 }
@@ -1641,7 +1641,7 @@ static VALUE crmf_sub_ru(VALUE self, VALUE x, VALUE y)
 	double dx = NUM2DBL(x);
 	double dy = NUM2DBL(y);
 	fesetround(FE_UPWARD);
-	double dr = dx - dy;
+	volatile double dr = dx - dy;
 	fesetround(FE_TONEAREST);
 	return DBL2NUM(dr);
 }
@@ -1651,7 +1651,7 @@ static VALUE crmf_sub_rn(VALUE self, VALUE x, VALUE y)
 	double dx = NUM2DBL(x);
 	double dy = NUM2DBL(y);
 	fesetround(FE_TONEAREST);
-	double dr = dx - dy;
+	volatile double dr = dx - dy;
 	return DBL2NUM(dr);
 }
 
@@ -1661,7 +1661,7 @@ static VALUE crmf_mul_rz(VALUE self, VALUE x, VALUE y)
 	double dx = NUM2DBL(x);
 	double dy = NUM2DBL(y);
 	fesetround(FE_TOWARDZERO);
-	double dr = dx * dy;
+	volatile double dr = dx * dy;
 	fesetround(FE_TONEAREST);
 	return DBL2NUM(dr);
 }
@@ -1671,7 +1671,7 @@ static VALUE crmf_mul_rd(VALUE self, VALUE x, VALUE y)
 	double dx = NUM2DBL(x);
 	double dy = NUM2DBL(y);
 	fesetround(FE_DOWNWARD);
-	double dr = dx * dy;
+	volatile double dr = dx * dy;
 	fesetround(FE_TONEAREST);
 	return DBL2NUM(dr);
 }
@@ -1681,7 +1681,7 @@ static VALUE crmf_mul_ru(VALUE self, VALUE x, VALUE y)
 	double dx = NUM2DBL(x);
 	double dy = NUM2DBL(y);
 	fesetround(FE_UPWARD);
-	double dr = dx * dy;
+	volatile double dr = dx * dy;
 	fesetround(FE_TONEAREST);
 	return DBL2NUM(dr);
 }
@@ -1691,7 +1691,7 @@ static VALUE crmf_mul_rn(VALUE self, VALUE x, VALUE y)
 	double dx = NUM2DBL(x);
 	double dy = NUM2DBL(y);
 	fesetround(FE_TONEAREST);
-	double dr = dx * dy;
+	volatile double dr = dx * dy;
 	return DBL2NUM(dr);
 }
 
@@ -1701,7 +1701,7 @@ static VALUE crmf_div_rz(VALUE self, VALUE x, VALUE y)
 	double dx = NUM2DBL(x);
 	double dy = NUM2DBL(y);
 	fesetround(FE_TOWARDZERO);
-	double dr = dx / dy;
+	volatile double dr = dx / dy;
 	fesetround(FE_TONEAREST);
 	return DBL2NUM(dr);
 }
@@ -1711,7 +1711,7 @@ static VALUE crmf_div_rd(VALUE self, VALUE x, VALUE y)
 	double dx = NUM2DBL(x);
 	double dy = NUM2DBL(y);
 	fesetround(FE_DOWNWARD);
-	double dr = dx / dy;
+	volatile double dr = dx / dy;
 	fesetround(FE_TONEAREST);
 	return DBL2NUM(dr);
 }
@@ -1721,7 +1721,7 @@ static VALUE crmf_div_ru(VALUE self, VALUE x, VALUE y)
 	double dx = NUM2DBL(x);
 	double dy = NUM2DBL(y);
 	fesetround(FE_UPWARD);
-	double dr = dx / dy;
+	volatile double dr = dx / dy;
 	fesetround(FE_TONEAREST);
 	return DBL2NUM(dr);
 }
@@ -1731,7 +1731,7 @@ static VALUE crmf_div_rn(VALUE self, VALUE x, VALUE y)
 	double dx = NUM2DBL(x);
 	double dy = NUM2DBL(y);
 	fesetround(FE_TONEAREST);
-	double dr = dx / dy;
+	volatile double dr = dx / dy;
 	return DBL2NUM(dr);
 }
 

data/ext/crmf/extconf.rb

@@ -26,24 +26,28 @@ unless have_library('mpfr') then
 end
 
 ext_dir = File.absolute_path(File.join('..', '..'), __FILE__)
-crlibm_src = File.join(ext_dir, 'crlibm-1.0beta4')
-crlibm_archive = crlibm_src + '.tar.gz'
+crlibm_src_dir = File.join(ext_dir, 'crlibm-1.0beta5')
 crlibm_install_dir = File.join(ext_dir, 'crlibm')
 crlibm_lib_dir = File.join(crlibm_install_dir, 'lib')
 crlibm_inc_dir = File.join(crlibm_install_dir, 'include')
 crlibm_bin_dir = File.join(crlibm_install_dir, 'bin')
 unless File.directory?(crlibm_install_dir) then
-  cmd = "tar xvf '#{crlibm_archive}' -C '#{ext_dir}'"
-  abort "-----\n#{caller[0]}\nCannot untar CRLIBM.\nCommand: \"#{cmd}\"\n-----" unless system(cmd)
-  Dir.chdir(crlibm_src) do
-    cmd = "./configure --prefix='#{crlibm_install_dir}'#{RUBY_PLATFORM =~ /x86_64/ ? ' --enable-sse2 --enable-double-extended' : ''}"
-    abort "-----\n#{caller[0]}\nCannot configure CRLIBM.\nCommand: \"#{cmd}\"\n-----" unless system(cmd)
+  Dir.chdir(crlibm_src_dir) do
+    unless File.file?(File.join(crlibm_src_dir, 'configure')) then
+      ENV['CFLAGS'] = '-g -O2 -fPIC'
+      cmd = "./prepare --prefix='#{crlibm_install_dir}'#{RUBY_PLATFORM =~ /x86_64/ ? ' --enable-sse2 --enable-double-extended' : ''}"
+      abort "-----\n#{caller[0]}\nCannot configure CRLIBM.\nCommand: \"#{cmd}\"\n-----" unless system(cmd)
+    end
+    #cmd = "../configure --prefix='#{crlibm_install_dir}'#{RUBY_PLATFORM =~ /x86_64/ ? ' --enable-sse2 --enable-double-extended' : ''}"
+    #abort "-----\n#{caller[0]}\nCannot configure CRLIBM.\nCommand: \"#{cmd}\"\n-----" unless system(cmd)
     cmd = "make"
     abort "-----\n#{caller[0]}\nCannot make CRLIBM.\nCommand: \"#{cmd}\"\n-----" unless system(cmd)
     cmd = "make -j1 install"
     abort "-----\n#{caller[0]}\nCannot locally install CRLIBM.\nCommand: \"#{cmd}\"\n-----" unless system(cmd)
+    cmd = "make clean"
+    abort "-----\n#{caller[0]}\nCannot cleanup CRLIBM.\nCommand: \"#{cmd}\"\n-----" unless system(cmd)
   end
-  cmd = "rm -rf '#{crlibm_src}' '#{crlibm_bin_dir}'"
+  cmd = "rm -rf '#{crlibm_bin_dir}'"
   abort "-----\n#{caller[0]}\nCannot cleanup CRLIBM.\nCommand: \"#{cmd}\"\n-----" unless system(cmd)
 end
 dir_config('crlibm', crlibm_inc_dir, crlibm_lib_dir)

data/lib/crmf/version.rb

@@ -1,4 +1,4 @@
 module CRMF
   # @return [String] CRMF gem version
-  VERSION = '0.1.1'
+  VERSION = '0.1.2'
 end

data/tests/perf.rb

@@ -1,5 +1,6 @@
 #!/usr/bin/env ruby
 require '../lib/crmf.rb'
+require 'colorize'
 
 def print_time(t)
   if t >= 3600 then
@@ -124,231 +125,111 @@ funcs.each do |m|
   end
   puts "libm   vs mpfr: #{nb_diff_libm_mpfr}/#{nb_tests} (#{nb_diff_libm_mpfr*100.0/nb_tests}%)#{estr}"
   nb_diff_crlibm_mpfr = nb_tests.times.count{|i| res_crlibm[i] != res_mpfr[i]}
-  puts "crlibm vs mpfr: #{nb_diff_crlibm_mpfr}/#{nb_tests} (#{nb_diff_crlibm_mpfr*100.0/nb_tests}%)"
+  str = "crlibm vs mpfr: #{nb_diff_crlibm_mpfr}/#{nb_tests} (#{nb_diff_crlibm_mpfr*100.0/nb_tests}%)"
+  if nb_diff_crlibm_mpfr > 0 then
+    puts str.light_red
+  else
+    puts str
+  end
 end
 
 puts
 
 nbtm1 = nb_tests - 1
 
-puts " Add N ".center(80, '-')
-
-tik = Time.now
-nb_tests.times do |i|
-  res_libm[i] = f[i] + f[nbtm1-i] 
-end
-tok = Time.now
-t_libm = tok - tik
-puts "libm:   #{print_time(t_libm)} (#{print_time(t_libm / nb_tests)}/t)"
-
-tik = Time.now
-nb_tests.times do |i|
-  res_mpfr[i] = Float.add(f[i], f[nbtm1-i])
-end
-tok = Time.now
-t_mpfr = tok - tik
-puts "mpfr:   #{print_time(t_mpfr)} (#{print_time(t_mpfr / nb_tests)}/t)"
-
-tik = Time.now
-nb_tests.times do |i|
-  res_crlibm[i] = CRMF.add_rn(f[i], f[nbtm1-i])
-end
-tok = Time.now
-t_crlibm = tok - tik
-puts "crlibm: #{print_time(t_crlibm)} (#{print_time(t_crlibm / nb_tests)}/t)"
-
-max_diff_pct = 0.0
-nb_tests.times do |i|
-  next if (not(res_libm[i].finite?) and not(res_mpfr[i].finite?))
-  error_pct = 100.0 * (res_libm[i] - res_mpfr[i]).abs / res_mpfr[i]
-  max_diff_pct = [max_diff_pct, error_pct].max
-end
-nb_diff_libm_mpfr = nb_tests.times.count{|i| res_libm[i] != res_mpfr[i]}
-if nb_diff_libm_mpfr > 0 then
-  estr = ", max error: #{max_diff_pct}%"
-else
-  estr = ''
-end
-puts "libm   vs mpfr: #{nb_diff_libm_mpfr}/#{nb_tests} (#{nb_diff_libm_mpfr*100.0/nb_tests}%)#{estr}"
-
-max_diff_pct = 0.0
-nb_tests.times do |i|
-  next if (not(res_crlibm[i].finite?) and not(res_mpfr[i].finite?))
-  error_pct = 100.0 * (res_crlibm[i] - res_mpfr[i]).abs / res_mpfr[i]
-  max_diff_pct = [max_diff_pct, error_pct].max
-end
-nb_diff_crlibm_mpfr = nb_tests.times.count{|i| res_crlibm[i] != res_mpfr[i]}
-if nb_diff_crlibm_mpfr > 0 then
-  estr = ", max error: #{max_diff_pct}%"
-else
-  estr = ''
-end
-puts "crlibm vs mpfr: #{nb_diff_crlibm_mpfr}/#{nb_tests} (#{nb_diff_crlibm_mpfr*100.0/nb_tests}%)#{estr}"
-
-
-puts " Add D ".center(80, '-')
-
-tik = Time.now
-nb_tests.times do |i|
-  res_libm[i] = f[i] + f[nbtm1-i] 
-end
-tok = Time.now
-t_libm = tok - tik
-puts "libm:   #{print_time(t_libm)} (#{print_time(t_libm / nb_tests)}/t)"
-
-tik = Time.now
-nb_tests.times do |i|
-  res_mpfr[i] = Float.add(f[i], f[nbtm1-i], :down)
-end
-tok = Time.now
-t_mpfr = tok - tik
-puts "mpfr:   #{print_time(t_mpfr)} (#{print_time(t_mpfr / nb_tests)}/t)"
-
-tik = Time.now
-nb_tests.times do |i|
-  res_crlibm[i] = CRMF.add_rd(f[i], f[nbtm1-i])
-end
-tok = Time.now
-t_crlibm = tok - tik
-puts "crlibm: #{print_time(t_crlibm)} (#{print_time(t_crlibm / nb_tests)}/t)"
-
-max_diff_pct = 0.0
-nb_tests.times do |i|
-  next if (not(res_libm[i].finite?) and not(res_mpfr[i].finite?))
-  error_pct = 100.0 * (res_libm[i] - res_mpfr[i]).abs / res_mpfr[i]
-  max_diff_pct = [max_diff_pct, error_pct].max
-end
-nb_diff_libm_mpfr = nb_tests.times.count{|i| res_libm[i] != res_mpfr[i]}
-if nb_diff_libm_mpfr > 0 then
-  estr = ", max error: #{max_diff_pct}%"
-else
-  estr = ''
-end
-puts "libm   vs mpfr: #{nb_diff_libm_mpfr}/#{nb_tests} (#{nb_diff_libm_mpfr*100.0/nb_tests}%)#{estr}"
-
-max_diff_pct = 0.0
-nb_tests.times do |i|
-  next if (not(res_crlibm[i].finite?) and not(res_mpfr[i].finite?))
-  error_pct = 100.0 * (res_crlibm[i] - res_mpfr[i]).abs / res_mpfr[i]
-  max_diff_pct = [max_diff_pct, error_pct].max
-end
-nb_diff_crlibm_mpfr = nb_tests.times.count{|i| res_crlibm[i] != res_mpfr[i]}
-if nb_diff_crlibm_mpfr > 0 then
-  estr = ", max error: #{max_diff_pct}%"
-else
-  estr = ''
-end
-puts "crlibm vs mpfr: #{nb_diff_crlibm_mpfr}/#{nb_tests} (#{nb_diff_crlibm_mpfr*100.0/nb_tests}%)#{estr}"
-
-
-
-puts " Add U ".center(80, '-')
-
-tik = Time.now
-nb_tests.times do |i|
-  res_libm[i] = f[i] + f[nbtm1-i] 
-end
-tok = Time.now
-t_libm = tok - tik
-puts "libm:   #{print_time(t_libm)} (#{print_time(t_libm / nb_tests)}/t)"
-
-tik = Time.now
-nb_tests.times do |i|
-  res_mpfr[i] = Float.add(f[i], f[nbtm1-i], :up)
-end
-tok = Time.now
-t_mpfr = tok - tik
-puts "mpfr:   #{print_time(t_mpfr)} (#{print_time(t_mpfr / nb_tests)}/t)"
-
-tik = Time.now
-nb_tests.times do |i|
-  res_crlibm[i] = CRMF.add_ru(f[i], f[nbtm1-i])
-end
-tok = Time.now
-t_crlibm = tok - tik
-puts "crlibm: #{print_time(t_crlibm)} (#{print_time(t_crlibm / nb_tests)}/t)"
-
-max_diff_pct = 0.0
-nb_tests.times do |i|
-  next if (not(res_libm[i].finite?) and not(res_mpfr[i].finite?))
-  error_pct = 100.0 * (res_libm[i] - res_mpfr[i]).abs / res_mpfr[i]
-  max_diff_pct = [max_diff_pct, error_pct].max
-end
-nb_diff_libm_mpfr = nb_tests.times.count{|i| res_libm[i] != res_mpfr[i]}
-if nb_diff_libm_mpfr > 0 then
-  estr = ", max error: #{max_diff_pct}%"
-else
-  estr = ''
-end
-puts "libm   vs mpfr: #{nb_diff_libm_mpfr}/#{nb_tests} (#{nb_diff_libm_mpfr*100.0/nb_tests}%)#{estr}"
-
-max_diff_pct = 0.0
-nb_tests.times do |i|
-  next if (not(res_crlibm[i].finite?) and not(res_mpfr[i].finite?))
-  error_pct = 100.0 * (res_crlibm[i] - res_mpfr[i]).abs / res_mpfr[i]
-  max_diff_pct = [max_diff_pct, error_pct].max
-end
-nb_diff_crlibm_mpfr = nb_tests.times.count{|i| res_crlibm[i] != res_mpfr[i]}
-if nb_diff_crlibm_mpfr > 0 then
-  estr = ", max error: #{max_diff_pct}%"
-else
-  estr = ''
-end
-puts "crlibm vs mpfr: #{nb_diff_crlibm_mpfr}/#{nb_tests} (#{nb_diff_crlibm_mpfr*100.0/nb_tests}%)#{estr}"
-
-
-puts " Add Z ".center(80, '-')
-
-tik = Time.now
-nb_tests.times do |i|
-  res_libm[i] = f[i] + f[nbtm1-i] 
-end
-tok = Time.now
-t_libm = tok - tik
-puts "libm:   #{print_time(t_libm)} (#{print_time(t_libm / nb_tests)}/t)"
-
-tik = Time.now
-nb_tests.times do |i|
-  res_mpfr[i] = Float.add(f[i], f[nbtm1-i], :zero)
-end
-tok = Time.now
-t_mpfr = tok - tik
-puts "mpfr:   #{print_time(t_mpfr)} (#{print_time(t_mpfr / nb_tests)}/t)"
-
-tik = Time.now
-nb_tests.times do |i|
-  res_crlibm[i] = CRMF.add_rz(f[i], f[nbtm1-i])
-end
-tok = Time.now
-t_crlibm = tok - tik
-puts "crlibm: #{print_time(t_crlibm)} (#{print_time(t_crlibm / nb_tests)}/t)"
-
-max_diff_pct = 0.0
-nb_tests.times do |i|
-  next if (not(res_libm[i].finite?) and not(res_mpfr[i].finite?))
-  error_pct = 100.0 * (res_libm[i] - res_mpfr[i]).abs / res_mpfr[i]
-  max_diff_pct = [max_diff_pct, error_pct].max
-end
-nb_diff_libm_mpfr = nb_tests.times.count{|i| res_libm[i] != res_mpfr[i]}
-if nb_diff_libm_mpfr > 0 then
-  estr = ", max error: #{max_diff_pct}%"
-else
-  estr = ''
-end
-puts "libm   vs mpfr: #{nb_diff_libm_mpfr}/#{nb_tests} (#{nb_diff_libm_mpfr*100.0/nb_tests}%)#{estr}"
-
-max_diff_pct = 0.0
-nb_tests.times do |i|
-  next if (not(res_crlibm[i].finite?) and not(res_mpfr[i].finite?))
-  error_pct = 100.0 * (res_crlibm[i] - res_mpfr[i]).abs / res_mpfr[i]
-  max_diff_pct = [max_diff_pct, error_pct].max
-end
-nb_diff_crlibm_mpfr = nb_tests.times.count{|i| res_crlibm[i] != res_mpfr[i]}
-if nb_diff_crlibm_mpfr > 0 then
-  estr = ", max error: #{max_diff_pct}%"
-else
-  estr = ''
+[
+  [:add, :+], 
+  [:sub, :-], 
+  [:mul, :*], 
+  [:div, :/] 
+].each do |m_m|
+  [:z, :d, :u, :n].each do |r|
+    puts " #{m_m[0].to_s.capitalize} #{r.to_s.upcase} ".center(80, '-')
+    rnd = case r
+          when :z
+            :zero
+          when :d
+            :down
+          when :u
+            :up
+          when :n
+            :nearest
+          end
+
+    tik = Time.now
+    case m_m[0]
+    when :add
+      nb_tests.times do |i|
+        res_libm[i] = f[i] + f[nbtm1-i] 
+      end
+    when :sub
+      nb_tests.times do |i|
+        res_libm[i] = f[i] - f[nbtm1-i] 
+      end
+    when :mul
+      nb_tests.times do |i|
+        res_libm[i] = f[i] * f[nbtm1-i] 
+      end
+    when :div
+      nb_tests.times do |i|
+        res_libm[i] = f[i] / f[nbtm1-i] 
+      end
+    end
+    tok = Time.now
+    t_libm = tok - tik
+    puts "libm: #{print_time(t_libm)} (#{print_time(t_libm / nb_tests)}/t)"
+
+    method = Float.method(m_m[0])
+    tik = Time.now
+    nb_tests.times do |i|
+      res_mpfr[i] = method.call(f[i], f[nbtm1-i], rnd)
+    end
+    tok = Time.now
+    t_mpfr = tok - tik
+    puts "mpfr: #{print_time(t_mpfr)} (#{print_time(t_mpfr / nb_tests)}/t)"
+
+    method = CRMF.method("#{m_m[0]}_r#{r}")
+    tik = Time.now
+    nb_tests.times do |i|
+      res_crlibm[i] = method.call(f[i], f[nbtm1-i])
+    end
+    tok = Time.now
+    t_crlibm = tok - tik
+    puts "FPU:  #{print_time(t_crlibm)} (#{print_time(t_crlibm / nb_tests)}/t)"
+
+    max_diff_pct = 0.0
+    nb_tests.times do |i|
+      next if (not(res_libm[i].finite?) and not(res_mpfr[i].finite?))
+      error_pct = 100.0 * (res_libm[i] - res_mpfr[i]).abs / res_mpfr[i]
+      max_diff_pct = [max_diff_pct, error_pct].max
+    end
+    nb_diff_libm_mpfr = nb_tests.times.count{|i| res_libm[i] != res_mpfr[i]}
+    if nb_diff_libm_mpfr > 0 then
+      estr = ", max error: #{max_diff_pct}%"
+    else
+      estr = ''
+    end
+    puts "libm vs mpfr: #{nb_diff_libm_mpfr}/#{nb_tests} (#{nb_diff_libm_mpfr*100.0/nb_tests}%)#{estr}"
+
+    max_diff_pct = 0.0
+    nb_tests.times do |i|
+      next if (not(res_crlibm[i].finite?) and not(res_mpfr[i].finite?))
+      error_pct = 100.0 * (res_crlibm[i] - res_mpfr[i]).abs / res_mpfr[i]
+      max_diff_pct = [max_diff_pct, error_pct].max
+    end
+    nb_diff_crlibm_mpfr = nb_tests.times.count{|i| res_crlibm[i] != res_mpfr[i]}
+    if nb_diff_crlibm_mpfr > 0 then
+      estr = ", max error: #{max_diff_pct}%"
+    else
+      estr = ''
+    end
+    str = "FPU  vs mpfr: #{nb_diff_crlibm_mpfr}/#{nb_tests} (#{nb_diff_crlibm_mpfr*100.0/nb_tests}%)#{estr}"
+    if nb_diff_crlibm_mpfr > 0 then
+      puts str.light_red
+    else
+      puts str
+    end
+  end
 end
-puts "crlibm vs mpfr: #{nb_diff_crlibm_mpfr}/#{nb_tests} (#{nb_diff_crlibm_mpfr*100.0/nb_tests}%)#{estr}"
 
 

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: crmf
 version: !ruby/object:Gem::Version
-  version: 0.1.1
+  version: 0.1.2
 platform: ruby
 authors:
 - Théotime Bollengier
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2022-07-21 00:00:00.000000000 Z
+date: 2022-07-22 00:00:00.000000000 Z
 dependencies: []
 description: CRMF is a Ruby C extension which provides correctly rounded math functions
   for Ruby floats, using MPFR and CRlibm. Provided rounding modes are toward zero,
@@ -23,7 +23,108 @@ files:
 - LICENSE
 - README.md
 - crmf.gemspec
-- ext/crlibm-1.0beta4.tar.gz
+- ext/crlibm-1.0beta5/AUTHORS
+- ext/crlibm-1.0beta5/CMakeLists.txt
+- ext/crlibm-1.0beta5/COPYING
+- ext/crlibm-1.0beta5/COPYING.LIB
+- ext/crlibm-1.0beta5/ChangeLog
+- ext/crlibm-1.0beta5/Makefile.am
+- ext/crlibm-1.0beta5/NEWS
+- ext/crlibm-1.0beta5/README
+- ext/crlibm-1.0beta5/README.DEV
+- ext/crlibm-1.0beta5/README.md
+- ext/crlibm-1.0beta5/TODO
+- ext/crlibm-1.0beta5/VERSION
+- ext/crlibm-1.0beta5/acos-td.c
+- ext/crlibm-1.0beta5/acos-td.h
+- ext/crlibm-1.0beta5/asin-td.c
+- ext/crlibm-1.0beta5/asin-td.h
+- ext/crlibm-1.0beta5/asincos.c
+- ext/crlibm-1.0beta5/asincos.h
+- ext/crlibm-1.0beta5/atan-itanium.c
+- ext/crlibm-1.0beta5/atan-pentium.c
+- ext/crlibm-1.0beta5/atan-pentium.h
+- ext/crlibm-1.0beta5/atan_accurate.c
+- ext/crlibm-1.0beta5/atan_accurate.h
+- ext/crlibm-1.0beta5/atan_fast.c
+- ext/crlibm-1.0beta5/atan_fast.h
+- ext/crlibm-1.0beta5/configure.ac
+- ext/crlibm-1.0beta5/crlibm.h
+- ext/crlibm-1.0beta5/crlibm.spec
+- ext/crlibm-1.0beta5/crlibm_private.c
+- ext/crlibm-1.0beta5/crlibm_private.h
+- ext/crlibm-1.0beta5/csh_fast.c
+- ext/crlibm-1.0beta5/csh_fast.h
+- ext/crlibm-1.0beta5/double-extended.h
+- ext/crlibm-1.0beta5/exp-itanium.c
+- ext/crlibm-1.0beta5/exp-td-standalone.c
+- ext/crlibm-1.0beta5/exp-td.c
+- ext/crlibm-1.0beta5/exp-td.h
+- ext/crlibm-1.0beta5/exp_build_coeffs/exp_fast_table.c
+- ext/crlibm-1.0beta5/expm1-standalone.c
+- ext/crlibm-1.0beta5/expm1.c
+- ext/crlibm-1.0beta5/expm1.h
+- ext/crlibm-1.0beta5/interval.h
+- ext/crlibm-1.0beta5/log-de.c
+- ext/crlibm-1.0beta5/log-de.h
+- ext/crlibm-1.0beta5/log-de2.c
+- ext/crlibm-1.0beta5/log-de2.h
+- ext/crlibm-1.0beta5/log-td.c
+- ext/crlibm-1.0beta5/log-td.h
+- ext/crlibm-1.0beta5/log.c
+- ext/crlibm-1.0beta5/log.h
+- ext/crlibm-1.0beta5/log10-td.c
+- ext/crlibm-1.0beta5/log10-td.h
+- ext/crlibm-1.0beta5/log1p.c
+- ext/crlibm-1.0beta5/log2-td.c
+- ext/crlibm-1.0beta5/log2-td.h
+- ext/crlibm-1.0beta5/log2_accurate.c
+- ext/crlibm-1.0beta5/log2_accurate.h
+- ext/crlibm-1.0beta5/log_accurate.c
+- ext/crlibm-1.0beta5/log_accurate.h
+- ext/crlibm-1.0beta5/log_fast.c
+- ext/crlibm-1.0beta5/log_fast.h
+- ext/crlibm-1.0beta5/pow.c
+- ext/crlibm-1.0beta5/pow.h
+- ext/crlibm-1.0beta5/prepare
+- ext/crlibm-1.0beta5/rem_pio2_accurate.c
+- ext/crlibm-1.0beta5/rem_pio2_accurate.h
+- ext/crlibm-1.0beta5/scs_lib/AUTHORS
+- ext/crlibm-1.0beta5/scs_lib/COPYING
+- ext/crlibm-1.0beta5/scs_lib/ChangeLog
+- ext/crlibm-1.0beta5/scs_lib/Doxyfile.dev
+- ext/crlibm-1.0beta5/scs_lib/Doxyfile.user
+- ext/crlibm-1.0beta5/scs_lib/INSTALL
+- ext/crlibm-1.0beta5/scs_lib/Makefile.am
+- ext/crlibm-1.0beta5/scs_lib/NEWS
+- ext/crlibm-1.0beta5/scs_lib/README
+- ext/crlibm-1.0beta5/scs_lib/README.DEV
+- ext/crlibm-1.0beta5/scs_lib/TODO
+- ext/crlibm-1.0beta5/scs_lib/VERSION
+- ext/crlibm-1.0beta5/scs_lib/addition_scs.c
+- ext/crlibm-1.0beta5/scs_lib/division_scs.c
+- ext/crlibm-1.0beta5/scs_lib/double2scs.c
+- ext/crlibm-1.0beta5/scs_lib/main.dox
+- ext/crlibm-1.0beta5/scs_lib/multiplication_scs.c
+- ext/crlibm-1.0beta5/scs_lib/poly_fct.c
+- ext/crlibm-1.0beta5/scs_lib/print_scs.c
+- ext/crlibm-1.0beta5/scs_lib/rand_scs.c
+- ext/crlibm-1.0beta5/scs_lib/scs.h
+- ext/crlibm-1.0beta5/scs_lib/scs2double.c
+- ext/crlibm-1.0beta5/scs_lib/scs2mpf.c
+- ext/crlibm-1.0beta5/scs_lib/scs2mpfr.c
+- ext/crlibm-1.0beta5/scs_lib/scs_private.c
+- ext/crlibm-1.0beta5/scs_lib/scs_private.h
+- ext/crlibm-1.0beta5/scs_lib/wrapper_scs.h
+- ext/crlibm-1.0beta5/scs_lib/zero_scs.c
+- ext/crlibm-1.0beta5/trigo_accurate.c
+- ext/crlibm-1.0beta5/trigo_accurate.h
+- ext/crlibm-1.0beta5/trigo_fast.c
+- ext/crlibm-1.0beta5/trigo_fast.h
+- ext/crlibm-1.0beta5/trigpi.c
+- ext/crlibm-1.0beta5/trigpi.h
+- ext/crlibm-1.0beta5/triple-double.c
+- ext/crlibm-1.0beta5/triple-double.h
 - ext/crmf/crmf.c
 - ext/crmf/crmf.h
 - ext/crmf/crmf.map