ruby-mpc 0.0.2

data/ext/mpc/ruby_mpc.h

@@ -0,0 +1,43 @@
+#ifndef _RUBY_MPC_H_
+#define _RUBY_MPC_H_
+
+#include <stdio.h>
+#include <ruby.h>
+#include <mpfr.h>
+#include <mpc.h>
+#include "ruby_mpfr.h"
+
+typedef __mpc_struct MPC;
+
+VALUE r_mpc_class, r_mpc_math;
+
+#define r_mpc_make_struct(ruby_var, c_var) { ruby_var = Data_Make_Struct(r_mpc_class, MPC, 0, r_mpc_free, c_var); }
+#define r_mpc_make_struct_init(ruby_var, c_var) { r_mpc_make_struct(ruby_var, c_var); mpc_init2(c_var, mpfr_get_default_prec()); }
+#define r_mpc_make_struct_init2(ruby_var, c_var, prec) { r_mpc_make_struct(ruby_var, c_var); mpc_init2(c_var, prec); }
+#define r_mpc_get_struct(c_var, ruby_var) { Data_Get_Struct(ruby_var, MPC, c_var); }
+
+#define r_mpc_temp_alloc(c_var) { c_var=ALLOC_N(MPC, 1); }
+#define r_mpc_temp_alloc_init(c_var) { r_mpc_temp_alloc(c_var); mpc_init2(c_var, mpfr_get_default_prec()); }
+#define r_mpc_temp_alloc_init2(c_var, prec) { r_mpc_temp_alloc(c_var); mpc_init2(c_var, prec); }
+#define r_mpc_temp_free(c_var) { mpc_clear(c_var); free(c_var); }
+
+#define r_mpc_rnd_from_value(rnd) (mp_rnd_t)NUM2INT(rnd)
+
+void r_mpc_free(void *ptr);
+VALUE r_mpc_new_c_obj(VALUE obj);
+/* VALUE r_mpc_new_robj(MPC *ptr); */
+/* void r_mpc_set_robj(MPC *ptr, VALUE obj); */
+
+/* Arguments */
+mpc_rnd_t r_mpc_rnd_from_optional_argument(int min, int max, int argc, VALUE *argv);
+mpc_rnd_t r_mpc_prec_from_optional_argument(int min, int max, int argc, VALUE *argv);
+void r_mpc_get_rnd_prec_from_optional_arguments(mpc_rnd_t *rnd, mp_prec_t *prec, int min, int max,
+						int argc, VALUE *argv);
+
+void rb_mpc_extended_set_default_rounding_mode (mpc_rnd_t rnd_mode);
+mpc_rnd_t rb_mpc_extended_get_default_rounding_mode (void);
+mp_prec_t rb_mpc_extended_get_max_prec (mpc_t x);
+void r_mpc_set_c_function_state(int num);
+void r_mpc_set_fr_function_state(int num);
+
+#endif /* _RUBY_MPC_H_ */

data/ext/mpc/ruby_mpfr.h

@@ -0,0 +1,38 @@
+#ifndef _RUBY_MPFR_H_
+#define _RUBY_MPFR_H_
+
+#include <ruby.h>
+#include <mpfr.h>
+#include <stdio.h>
+
+typedef __mpfr_struct MPFR;
+
+VALUE r_mpfr_class, r_mpfr_math;
+
+#define r_mpfr_make_struct(ruby_var, c_var) { ruby_var = Data_Make_Struct(r_mpfr_class, MPFR, 0, r_mpfr_free, c_var); }
+#define r_mpfr_make_struct_init(ruby_var, c_var) { r_mpfr_make_struct(ruby_var, c_var); mpfr_init(c_var); }
+#define r_mpfr_make_struct_init2(ruby_var, c_var, prec) { r_mpfr_make_struct(ruby_var, c_var); mpfr_init2(c_var, prec); }
+#define r_mpfr_get_struct(c_var, ruby_var) { Data_Get_Struct(ruby_var, MPFR, c_var); }
+
+#define r_mpfr_temp_alloc(var) { var=ALLOC_N(MPFR, 1); }
+#define r_mpfr_temp_alloc_init(var) { r_mpfr_temp_alloc(var); mpfr_init(var); }
+#define r_mpfr_temp_alloc_init2(var, prec) { r_mpfr_temp_alloc(var); mpfr_init2(var, prec); }
+#define r_mpfr_temp_free(var) { mpfr_clear(var); free(var); }
+
+#define r_mpfr_check_number(c_val) { if(mpfr_number_p(c_val) == 0) rb_raise(rb_eArgError, "Not an ordinary number."); }
+#define r_mpfr_check_positive_number(c_val) { if(mpfr_number_p(c_val) == 0 && mpfr_sgn(c_val) <= 0) rb_raise(rb_eArgError, "Not a positive number."); }
+#define r_mpfr_check_non_negative_number(c_val) { if(mpfr_number_p(c_val) == 0 && mpfr_sgn(c_val) < 0) rb_raise(rb_eArgError, "Not a non negative number."); }
+#define r_mpfr_check_negative_number(c_val) { if(mpfr_number_p(c_val) == 0 && mpfr_sgn(c_val) >= 0) rb_raise(rb_eArgError, "Not a negative number."); }
+#define r_mpfr_check_non_positive_number(c_val) { if(mpfr_number_p(c_val) == 0 && mpfr_sgn(c_val) > 0) rb_raise(rb_eArgError, "Not a non positive number."); }
+
+void r_mpfr_free(void *ptr);
+VALUE r_mpfr_new_fr_obj(VALUE obj);
+void r_mpfr_set_robj(MPFR *ptr, VALUE obj, mp_rnd_t rnd);
+
+mp_rnd_t r_mpfr_rnd_from_value(VALUE rnd);
+mp_rnd_t r_mpfr_rnd_from_optional_argument(int min, int max, int argc, VALUE *argv);
+mp_rnd_t r_mpfr_prec_from_optional_argument(int min, int max, int argc, VALUE *argv);
+void r_mpfr_get_rnd_prec_from_optional_arguments(mp_rnd_t *rnd, mp_prec_t *prec, int min, int max, int argc, VALUE *argv);
+
+#endif /* _RUBY_MPFR_H_ */
+

data/lib/mpc/version.rb

@@ -0,0 +1,3 @@
+module MPC
+  VERSION = '0.0.2'
+end

data/ruby-mpc.gemspec

@@ -0,0 +1,35 @@
+# -*- encoding: utf-8 -*-
+
+Gem::Specification.new do |s|
+  s.name = %q{ruby-mpc}
+  s.version = "0.0.2"
+
+  s.required_rubygems_version = Gem::Requirement.new(">= 0") if s.respond_to? :required_rubygems_version=
+  s.authors = ["Takayuki YAMAGUCHI"]
+  s.date = %q{2009-12-17}
+  s.description = %q{FIX (describe your package)}
+  s.email = ["d@ytak.info"]
+  s.extensions = ["ext/mpc/extconf.rb"]
+  s.extra_rdoc_files = ["History.txt", "Manifest.txt", "PostInstall.txt"]
+  s.files = ["History.txt", "Manifest.txt", "PostInstall.txt", "README.rdoc", "Rakefile", "ext/mpc/extconf.rb", "ext/mpc/ruby_mpc.c", "ext/mpc/ruby_mpc.h", "ext/mpc/ruby_mpfr.h", "lib/mpc/version.rb", "ruby-mpc.gemspec", "script/console", "script/destroy", "script/generate", "spec/mpc/arithmetic_op_spec.rb", "spec/mpc/const_spec.rb", "spec/mpc/functions_spec.rb", "spec/mpc/generate_complex_number.rb", "spec/mpc/spec_helper.rb", "spec/spec.opts", "tasks/extconf.rake"]
+  s.homepage = %q{http://github.com/#{github_username}/#{project_name}}
+  s.post_install_message = %q{PostInstall.txt}
+  s.rdoc_options = ["--main", "README.rdoc"]
+  s.require_paths = ["lib", "ext"]
+  s.rubyforge_project = %q{ruby-mpc}
+  s.rubygems_version = %q{1.3.5}
+  s.summary = %q{FIX (describe your package)}
+
+  if s.respond_to? :specification_version then
+    current_version = Gem::Specification::CURRENT_SPECIFICATION_VERSION
+    s.specification_version = 3
+
+    if Gem::Version.new(Gem::RubyGemsVersion) >= Gem::Version.new('1.2.0') then
+      s.add_development_dependency(%q<hoe>, [">= 2.4.0"])
+    else
+      s.add_dependency(%q<hoe>, [">= 2.4.0"])
+    end
+  else
+    s.add_dependency(%q<hoe>, [">= 2.4.0"])
+  end
+end

data/script/console

@@ -0,0 +1,10 @@
+#!/usr/bin/env ruby
+# File: script/console
+irb = RUBY_PLATFORM =~ /(:?mswin|mingw)/ ? 'irb.bat' : 'irb'
+
+libs =  " -r irb/completion"
+# Perhaps use a console_lib to store any extra methods I may want available in the cosole
+# libs << " -r #{File.dirname(__FILE__) + '/../lib/console_lib/console_logger.rb'}"
+libs <<  " -r #{File.dirname(__FILE__) + '/../lib/ruby-mpc.rb'}"
+puts "Loading ruby-mpc gem"
+exec "#{irb} #{libs} --simple-prompt"

data/script/destroy

@@ -0,0 +1,14 @@
+#!/usr/bin/env ruby
+APP_ROOT = File.expand_path(File.join(File.dirname(__FILE__), '..'))
+
+begin
+  require 'rubigen'
+rescue LoadError
+  require 'rubygems'
+  require 'rubigen'
+end
+require 'rubigen/scripts/destroy'
+
+ARGV.shift if ['--help', '-h'].include?(ARGV[0])
+RubiGen::Base.use_component_sources! [:rubygems, :newgem, :newgem_theme, :test_unit]
+RubiGen::Scripts::Destroy.new.run(ARGV)

data/script/generate

@@ -0,0 +1,14 @@
+#!/usr/bin/env ruby
+APP_ROOT = File.expand_path(File.join(File.dirname(__FILE__), '..'))
+
+begin
+  require 'rubigen'
+rescue LoadError
+  require 'rubygems'
+  require 'rubigen'
+end
+require 'rubigen/scripts/generate'
+
+ARGV.shift if ['--help', '-h'].include?(ARGV[0])
+RubiGen::Base.use_component_sources! [:rubygems, :newgem, :newgem_theme, :test_unit]
+RubiGen::Scripts::Generate.new.run(ARGV)

data/spec/mpc/arithmetic_op_spec.rb

@@ -0,0 +1,30 @@
+require File.dirname(__FILE__) + '/spec_helper.rb'
+
+require "generate_complex_number.rb"
+
+MPFR.set_default_prec(256)
+
+ARITHMETIC_NUM_OF_CALCULATE = 1000
+
+def check_mpc_arithmetic_op(arg1, arg2, error = 10e-10, &block)
+  ruby_res = yield(Complex(*(arg1.map{ |a| a.to_f })), Complex(*(arg2.map{ |a| a.to_f })))
+  mpc_res = yield(MPC.new(*arg1), MPC.new(*arg2))
+  (MPFR.new(ruby_res.real) - mpc_res.real).abs.should < error
+  (MPFR.new(ruby_res.imag) - mpc_res.imag).abs.should < error
+end
+
+describe MPC, "when calculating" do
+  it "should be complex number" do
+    args = GenerateComplex.float_arguments(ARITHMETIC_NUM_OF_CALCULATE)
+    args.each_index do |i|
+      if i > 0
+        check_mpc_arithmetic_op(args[i-1], args[i]){ |a, b| a + b }
+        check_mpc_arithmetic_op(args[i-1], args[i]){ |a, b| a - b }
+        check_mpc_arithmetic_op(args[i-1], args[i]){ |a, b| a * b }
+        check_mpc_arithmetic_op(args[i-1], args[i]){ |a, b| a / b }
+      end
+    end
+  end
+
+end
+

data/spec/mpc/const_spec.rb

@@ -0,0 +1,23 @@
+require File.dirname(__FILE__) + '/spec_helper.rb'
+
+describe MPC, "when returning constants" do
+  it "should have constants" do
+    constants = MPC.constants
+    [:RNDNN, :RNDNZ, :RNDNU, :RNDND, :RNDZN, :RNDZZ, :RNDZU, :RNDZD,
+     :RNDUN, :RNDUZ, :RNDUU, :RNDUD, :RNDDN, :RNDDZ, :RNDDU, :RNDDD].each do |cst|
+      constants.include?(cst).should be_true
+      eval("MPC::#{cst}").should be_instance_of Fixnum
+    end
+
+    [:VERSION_MAJOR, :VERSION_MINOR, :VERSION_PATCHLEVEL].each do |cst|
+      constants.include?(cst).should be_true
+      eval("MPC::#{cst}").should be_instance_of Fixnum
+    end
+
+    constants.include?(:VERSION_STRING).should be_true
+    MPC::VERSION_STRING.should be_instance_of String
+  end
+
+end
+
+

data/spec/mpc/functions_spec.rb

@@ -0,0 +1,30 @@
+require File.dirname(__FILE__) + '/spec_helper.rb'
+
+require "generate_complex_number.rb"
+
+MPFR.set_default_prec(512)
+
+FUNCTIONS_NUM_OF_CALCULATE = 100
+
+def check_mpc_function(arg1, error = 10e-10, &block)
+  ruby_res = yield(Complex(*(arg1.map{ |a| a.to_f })))
+  mpc_res = yield(MPC.new(*arg1))
+  if MPC === mpc_res
+    (MPFR.new(ruby_res.real) - mpc_res.real).abs.should < error
+    (MPFR.new(ruby_res.imag) - mpc_res.imag).abs.should < error
+  elsif MPFR === mpc_res
+    (MPFR.new(ruby_res) - mpc_res).abs.should < error
+  end
+end
+
+describe MPC, "when calculating some functions" do
+  it "should be the same result as Complex class" do
+    args = GenerateComplex.float_arguments(FUNCTIONS_NUM_OF_CALCULATE)
+    args.each do |a|
+      check_mpc_function(a) { |b| b.abs }
+      check_mpc_function(a) { |b| b.arg }
+      check_mpc_function(a) { |b| b.conj }
+    end
+  end
+end
+  

data/spec/mpc/generate_complex_number.rb

@@ -0,0 +1,29 @@
+module GenerateComplex
+
+  def self.float_arguments(number, max = 100)
+    ret = Array.new(number){ |j| Array.new(2){ |i| rand(max) - rand } }
+    for i in 0...(ret.size)
+      case i % 4
+      when 1
+        ret[i].map!{ |a| -a }
+      when 2, 3
+        ret[i].map!{ |a| (rand > 0.5 ? -a : a) }
+      end
+    end
+    ret
+  end
+
+  MAX_NUMBER = 1000
+  
+  def self.string_arguments(number)
+    prec = MPFR.get_default_prec
+    max = 2 ** prec
+    Array.new(number) do |i|
+      sign1 = ((rand > 0.5 ? '-' : ''))
+      sign2 = ((rand > 0.5 ? '-' : ''))
+      ["#{sign1}#{rand(MAX_NUMBER)}.#{rand(max)}", "#{sign2}#{rand(MAX_NUMBER)}.#{rand(max)}"]
+    end
+  end
+  
+
+end

data/spec/mpc/spec_helper.rb

@@ -0,0 +1,16 @@
+begin
+  require 'spec'
+rescue LoadError
+  require 'rubygems' unless ENV['NO_RUBYGEMS']
+  gem 'rspec'
+  require 'spec'
+end
+
+$:.unshift(File.dirname(__FILE__) + '/../../lib')
+$:.unshift(*Dir.glob(File.dirname(__FILE__) + '/../../ext/*'))
+$:.unshift(File.dirname(__FILE__))
+
+require 'mpfr'
+require 'mpc'
+require "generate_complex_number"
+

data/spec/spec.opts

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

data/tasks/extconf.rake

@@ -0,0 +1,36 @@
+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
+
+desc "Run 'make realclean' for extended libraries"
+task "ext:realclean" do
+  Dir.glob("ext/**/Makefile").each do |path|
+    system("cd #{File.dirname(path)}; make realclean")
+  end
+end
+
+desc "Run 'make clean' for extended libraries"
+task "ext:clean" do
+  Dir.glob("ext/**/Makefile").each do |path|
+    system("cd #{File.dirname(path)}; make clean")
+  end
+end
+
+desc "Run 'make realclean' for extended libraries"
+task 'ext:make' do
+  Dir.glob("ext/**/extconf.rb").each do |path|
+    system("cd #{File.dirname(path)}; ruby extconf.rb && make")
+  end
+end

metadata

@@ -0,0 +1,89 @@
+--- !ruby/object:Gem::Specification 
+name: ruby-mpc
+version: !ruby/object:Gem::Version 
+  version: 0.0.2
+platform: ruby
+authors: 
+- Takayuki YAMAGUCHI
+autorequire: 
+bindir: bin
+cert_chain: []
+
+date: 2009-12-17 00:00:00 +09:00
+default_executable: 
+dependencies: 
+- !ruby/object:Gem::Dependency 
+  name: hoe
+  type: :development
+  version_requirement: 
+  version_requirements: !ruby/object:Gem::Requirement 
+    requirements: 
+    - - ">="
+      - !ruby/object:Gem::Version 
+        version: 2.4.0
+    version: 
+description: FIX (describe your package)
+email: 
+- d@ytak.info
+executables: []
+
+extensions: 
+- ext/mpc/extconf.rb
+extra_rdoc_files: 
+- History.txt
+- Manifest.txt
+- PostInstall.txt
+files: 
+- History.txt
+- Manifest.txt
+- PostInstall.txt
+- README.rdoc
+- Rakefile
+- ext/mpc/extconf.rb
+- ext/mpc/ruby_mpc.c
+- ext/mpc/ruby_mpc.h
+- ext/mpc/ruby_mpfr.h
+- lib/mpc/version.rb
+- ruby-mpc.gemspec
+- script/console
+- script/destroy
+- script/generate
+- spec/mpc/arithmetic_op_spec.rb
+- spec/mpc/const_spec.rb
+- spec/mpc/functions_spec.rb
+- spec/mpc/generate_complex_number.rb
+- spec/mpc/spec_helper.rb
+- spec/spec.opts
+- tasks/extconf.rake
+has_rdoc: true
+homepage: http://github.com/#{github_username}/#{project_name}
+licenses: []
+
+post_install_message: PostInstall.txt
+rdoc_options: 
+- --main
+- README.rdoc
+require_paths: 
+- lib
+- ext
+required_ruby_version: !ruby/object:Gem::Requirement 
+  requirements: 
+  - - ">="
+    - !ruby/object:Gem::Version 
+      version: "0"
+  version: 
+required_rubygems_version: !ruby/object:Gem::Requirement 
+  requirements: 
+  - - ">="
+    - !ruby/object:Gem::Version 
+      version: "0"
+  version: 
+requirements: []
+
+rubyforge_project: ruby-mpc
+rubygems_version: 1.3.5
+signing_key: 
+specification_version: 3
+summary: FIX (describe your package)
+test_files: []
+