ruby-mpfi 0.0.2

data/History.txt

@@ -0,0 +1,4 @@
+=== 0.0.1 2009-12-13
+
+* 1 major enhancement:
+  * Initial release

data/Manifest.txt

@@ -0,0 +1,62 @@
+History.txt
+Manifest.txt
+PostInstall.txt
+README.rdoc
+Rakefile
+ext/mpfi/extconf.rb
+ext/mpfi/func_mpfi_extention.c
+ext/mpfi/func_mpfi_extention.h
+ext/mpfi/make_c_source.rb
+ext/mpfi/ruby_mpfi.c
+ext/mpfi/ruby_mpfi.h
+ext/mpfi/ruby_mpfr.h
+ext/mpfi/yasnippet_mpfi.el
+ext/mpfi_complex/mpfi/extconf.rb
+ext/mpfi_complex/mpfi/func_mpfi_extention.h
+ext/mpfi_complex/mpfi/func_ruby_mpfi_complex.c
+ext/mpfi_complex/mpfi/func_ruby_mpfi_complex.h
+ext/mpfi_complex/mpfi/ruby_mpfi.h
+ext/mpfi_complex/mpfi/ruby_mpfi_complex.c
+ext/mpfi_complex/mpfi/ruby_mpfi_complex.h
+ext/mpfi_complex/mpfi/ruby_mpfr.h
+ext/mpfi_matrix/mpfi/extconf.rb
+ext/mpfi_matrix/mpfi/func_mpfi_extention.h
+ext/mpfi_matrix/mpfi/func_mpfi_matrix.c
+ext/mpfi_matrix/mpfi/func_mpfi_matrix.h
+ext/mpfi_matrix/mpfi/func_mpfr_matrix.h
+ext/mpfi_matrix/mpfi/func_ruby_mpfi_complex.h
+ext/mpfi_matrix/mpfi/ruby_mpfi.h
+ext/mpfi_matrix/mpfi/ruby_mpfi_complex.h
+ext/mpfi_matrix/mpfi/ruby_mpfi_matrix.c
+ext/mpfi_matrix/mpfi/ruby_mpfi_matrix.h
+ext/mpfi_matrix/mpfi/ruby_mpfr.h
+ext/mpfi_matrix/mpfi/ruby_mpfr_matrix.h
+lib/mpfi/matrix.rb
+lib/mpfi/version.rb
+ruby-mpfi.gemspec
+script/console
+script/destroy
+script/generate
+spec/mpfi/generate_number_module.rb
+spec/mpfi/mpfi_alloc_spec.rb
+spec/mpfi/mpfi_diam_arithmetic_spec.rb
+spec/mpfi/mpfi_interval_arithmetic_spec.rb
+spec/mpfi/mpfi_interval_functions_spec.rb
+spec/mpfi/mpfi_math_functions_spec.rb
+spec/mpfi/mpfi_set_operation_spec.rb
+spec/mpfi/ruby-mpfi_spec.rb
+spec/mpfi/spec_helper.rb
+spec/mpfi_complex/spec_helper.rb
+spec/mpfi_matrix/generate_matrix_arguments.rb
+spec/mpfi_matrix/mpfi_matrix_alloc_spec.rb
+spec/mpfi_matrix/mpfi_matrix_arithmetic_spec.rb
+spec/mpfi_matrix/mpfi_matrix_interval_func_spec.rb
+spec/mpfi_matrix/mpfi_matrix_set_element_spec.rb
+spec/mpfi_matrix/mpfi_matrix_set_operation_spec.rb
+spec/mpfi_matrix/mpfi_matrix_string_spec.rb
+spec/mpfi_matrix/mpfi_matrix_subdivision_spec.rb
+spec/mpfi_matrix/mpfi_square_matrix_spec.rb
+spec/mpfi_matrix/mpfi_vector_spec.rb
+spec/mpfi_matrix/spec_helper.rb
+spec/spec.opts
+tasks/extconf.rake

data/PostInstall.txt

@@ -0,0 +1,7 @@
+
+For more information on ruby-mpfi, see http://ruby-mpfi.rubyforge.org
+
+NOTE: Change this information in PostInstall.txt 
+You can also delete it if you don't want it.
+
+

data/README.rdoc

@@ -0,0 +1,48 @@
+= ruby-mpfi
+
+* http://github.com/#{github_username}/#{project_name}
+
+== DESCRIPTION:
+
+FIX (describe your package)
+
+== FEATURES/PROBLEMS:
+
+* FIX (list of features or problems)
+
+== SYNOPSIS:
+
+  FIX (code sample of usage)
+
+== REQUIREMENTS:
+
+* FIX (list of requirements)
+
+== INSTALL:
+
+* FIX (sudo gem install, anything else)
+
+== LICENSE:
+
+(The MIT License)
+
+Copyright (c) 2009 FIXME full name
+
+Permission is hereby granted, free of charge, to any person obtaining
+a copy of this software and associated documentation files (the
+'Software'), to deal in the Software without restriction, including
+without limitation the rights to use, copy, modify, merge, publish,
+distribute, sublicense, and/or sell copies of the Software, and to
+permit persons to whom the Software is furnished to do so, subject to
+the following conditions:
+
+The above copyright notice and this permission notice shall be
+included in all copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED 'AS IS', WITHOUT WARRANTY OF ANY KIND,
+EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
+MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.
+IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY
+CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT,
+TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
+SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.

data/Rakefile

@@ -0,0 +1,26 @@
+require 'rubygems'
+gem 'hoe', '>= 2.1.0'
+require 'hoe'
+require 'fileutils'
+require './lib/mpfi/version.rb'
+
+Hoe.plugin :newgem
+# Hoe.plugin :website
+# Hoe.plugin :cucumberfeatures
+
+# Generate all the Rake tasks
+# Run 'rake -T' to see list of generated tasks (from gem root directory)
+$hoe = Hoe.spec 'ruby-mpfi' do
+  self.developer 'Takayuki YAMAGUCHI', 'd@ytak.info'
+  self.post_install_message = 'PostInstall.txt' # TODO remove if post-install message not required
+  self.rubyforge_name       = self.name # TODO this is default value
+  # self.extra_deps         = [['activesupport','>= 2.0.2']]
+  self.spec_extras[:extensions] = ["ext/mpfi/extconf.rb", "ext/mpfi_complex/mpfi/extconf.rb", "ext/mpfi_matrix/mpfi/extconf.rb"]
+end
+
+require 'newgem/tasks'
+Dir['tasks/**/*.rake'].each { |t| load t }
+
+# TODO - want other tests/tasks run by default? Add them to the list
+# remove_task :default
+# task :default => [:spec, :features]

data/ext/mpfi/extconf.rb

@@ -0,0 +1,10 @@
+require 'mkmf'
+
+$CFLAGS += " -Wall"
+
+dir_config("mpfr")
+dir_config("mpfi")
+dir_config("gmp")
+if have_header('mpfr.h') && have_library('mpfr') && have_header('mpfi.h') && have_library('mpfi') && have_header('gmp.h') && have_library('gmp')
+  create_makefile("mpfi")
+end

data/ext/mpfi/func_mpfi_extention.c

@@ -0,0 +1,52 @@
+#include <mpfr.h>
+#include <mpfi.h>
+#include <mpfi_io.h>
+#include "ruby_mpfi.h"
+#include "func_mpfi_extention.h"
+
+void mpfi_mid_interval(mpfi_t ret, mpfi_t x){
+  mpfr_t left, right;
+  mpfr_init2(left, mpfi_get_prec(x));
+  mpfr_init2(right, mpfi_get_prec(x));
+
+  mpfr_add(left, r_mpfi_left_ptr(x), r_mpfi_right_ptr(x), GMP_RNDD);
+  mpfr_add(right, r_mpfi_left_ptr(x), r_mpfi_right_ptr(x), GMP_RNDU);
+
+  mpfr_div_ui(r_mpfi_left_ptr(ret), left, 2, GMP_RNDD);
+  mpfr_div_ui(r_mpfi_right_ptr(ret), right, 2, GMP_RNDU);
+
+  mpfr_clear(left);
+  mpfr_clear(right);
+}
+
+/* Retrun 0 if this function puts subdivision to *ret. */
+/* Otherwise, return -1. */
+int mpfi_subdivision(int num, mpfi_t *ret, mpfi_t x){
+  int i, ret_val = -1;
+  mpfr_t l;
+  mpfr_init(l);
+  mpfr_sub(l, r_mpfi_right_ptr(x), r_mpfi_left_ptr(x), GMP_RNDD);
+  mpfr_div_si(l, l, num, GMP_RNDD);
+  
+  mpfr_t x_diam;
+  mpfr_init(x_diam);
+  mpfi_diam_abs(x_diam, x);
+  if(mpfr_cmp(x_diam, l) > 0){
+    mpfr_set(r_mpfi_left_ptr(*ret), r_mpfi_left_ptr(x), GMP_RNDN);
+    mpfr_add(r_mpfi_right_ptr(*ret), r_mpfi_left_ptr(*ret), l, GMP_RNDU);
+
+    for(i = 1; i < num - 1; i ++){
+      mpfr_set(r_mpfi_left_ptr(*ret + i), r_mpfi_right_ptr(*ret + i - 1), GMP_RNDN);
+      mpfr_add(r_mpfi_right_ptr(*ret + i), r_mpfi_left_ptr(*ret + i), l, GMP_RNDU);
+    }
+
+    mpfr_set(r_mpfi_left_ptr(*ret + i), r_mpfi_right_ptr(*ret + i - 1), GMP_RNDN);
+    mpfr_set(r_mpfi_right_ptr(*ret + i), r_mpfi_right_ptr(x), GMP_RNDN);
+    ret_val = 0;
+  }
+  
+  mpfr_clear(x_diam);
+  mpfr_clear(l);
+  return ret_val;
+}
+

data/ext/mpfi/func_mpfi_extention.h

@@ -0,0 +1,2 @@
+void mpfi_mid_interval(mpfi_t ret, mpfi_t x);
+int mpfi_subdivision(int num, mpfi_t *ret, mpfi_t x);

data/ext/mpfi/make_c_source.rb

@@ -0,0 +1,115 @@
+
+def make_sepcial_functions(one_arguments, no_argument)
+  one_arguments.each do |name|
+    puts <<SRC
+static VALUE r_mpfi_math_#{name} (int argc, VALUE *argv, VALUE self){
+  R_MPFI *ptr_a0, *ptr_ret;
+  r_mpfi_get_convert_struct(argv[0], ptr_a0);
+  VALUE val_ret;
+  r_mpfi_make_struct_init2(val_ret, ptr_ret, r_mpfr_prec_from_optional_argument(1, 2, argc, argv));
+  r_mpfi_set_function_state(mpfi_#{name}(ptr_ret, ptr_a0));
+  return val_ret;
+}
+
+SRC
+  end
+
+  no_argument.each do |name|
+    puts <<SRC
+static VALUE r_mpfi_math_#{name} (int argc, VALUE *argv, VALUE self){
+  R_MPFI *ptr_ret;
+  VALUE val_ret;
+  r_mpfi_make_struct_init2(val_ret, ptr_ret, r_mpfr_prec_from_optional_argument(0, 1, argc, argv));
+  r_mpfi_set_function_state(mpfi_#{name}(ptr_ret));
+  return val_ret;
+}
+
+SRC
+  end
+
+  (one_arguments + no_argument).each do |name|
+    puts %!rb_define_module_function(r_mpfi_math, "#{name}", r_mpfi_math_#{name}, -1);!
+  end
+  
+end
+
+def make_comparison(positive_returned_value, non_zero_returned_value)
+  positive_returned_value.each do |name|
+    puts <<SRC
+static VALUE r_mpfi_#{name} (VALUE self){
+  R_MPFI *ptr_self;
+  r_mpfi_get_struct(ptr_self, self);
+  if (mpfi_#{name}(ptr_self) > 0) {
+    return Qtrue;
+  }else{
+    return Qnil;
+  }
+}
+
+SRC
+  end
+
+  non_zero_returned_value.each do |name|
+    puts <<SRC
+static VALUE r_mpfi_#{name} (VALUE self){
+  R_MPFI *ptr_self;
+  r_mpfi_get_struct(ptr_self, self);
+  if (mpfi_#{name}(ptr_self) != 0) {
+    return Qtrue;
+  }else{
+    return Qnil;
+  }
+}
+
+SRC
+  end
+
+  (positive_returned_value + non_zero_returned_value).each do |name|
+    puts %!rb_define_module_function(r_mpfi_class, "#{name}", r_mpfi_#{name}, 0);!
+  end
+end
+
+# Parse options.
+require 'optparse'
+
+help_message =<<HELP
+Make C source for extended library mpfr.
+HELP
+
+category = nil
+
+begin
+  OptionParser.new(help_message) do |opt|
+    # on(short [, klass [, pat [, desc]]]) {|v| ...}
+    # on(long [, klass [, pat [, desc]]) {|v| ...}
+    # on(short, long [, klass [, pat [, desc]]]) {|v| ...}
+    opt.on('--special') { |v| category = :special }
+    opt.on('--comparison') { |v| category = :comparison }
+    opt.parse!(ARGV)
+  end
+rescue OptionParser::InvalidOption
+  $stderr.print <<MES
+error: Invalid Option
+#{help_message}
+MES
+  exit(2)
+rescue OptionParser::InvalidArgument
+  $stderr.print <<MES
+error: Invalid Argument
+#{help_message}
+MES
+  exit(2)
+end
+
+case category
+when :special
+  one_arguments = ['log', 'exp', 'exp2', 'cos', 'sin', 'tan', 'acos', 'asin', 'atan', 'cosh', 'sinh', 'tanh',
+                   'acosh', 'asinh', 'atanh', 'log1p', 'expm1', 'log2', 'log10']
+  no_argument = ['const_log2', 'const_pi', 'const_euler']
+  make_sepcial_functions(one_arguments, no_argument)
+when :comparison
+  positive_returned_value = ['is_pos', 'is_strictly_pos', 'is_nonneg', 'is_neg', 'is_strictly_neg', 'is_nonpos', 'is_zero']
+  non_zero_returned_value = ['nan_p', 'inf_p', 'bounded_p']
+  make_comparison(positive_returned_value, non_zero_returned_value)
+end
+