gsl 1.15.3 → 1.16.0.6

data/ext/cqp.c

@@ -1,283 +0,0 @@
-#ifdef HAVE_GSL_GSL_CQP_H
-
-#include "rb_gsl.h"
-#include "gsl/gsl_cqp.h"
-
-static const gsl_cqpminimizer_type* type_by_string(VALUE t);
-static const gsl_cqpminimizer_type* get_type(VALUE t)
-{
-
-	switch (TYPE(t)) {
-	case T_STRING:
-		return type_by_string(t);
-		break;
-	default:
-		rb_raise(rb_eTypeError, "Wrong argument type %s.", rb_class2name(CLASS_OF(t)));
-	}
-}
-
-static const gsl_cqpminimizer_type* type_by_string(VALUE t)
-{
-	char *name;	
-	name = STR2CSTR(t);
-	if (strcmp(name, "mg_pdip") == 0) {
-		return gsl_cqpminimizer_mg_pdip;
-	} else {
-		rb_raise(rb_eRuntimeError, "Unknown minimizer type %s.", name);
-	}
-	return NULL; /* never reach here */
-}
-
-static VALUE rb_cqpminimizer_alloc(VALUE klass,  VALUE t, VALUE n, VALUE me, VALUE mi)
-{
-	gsl_cqpminimizer *m;
-	m = gsl_cqpminimizer_alloc(get_type(t), (size_t) FIX2INT(n), (size_t) FIX2INT(me), (size_t) FIX2INT(mi));
-	return Data_Wrap_Struct(klass, 0, gsl_cqpminimizer_free, m);
-}
-
-static VALUE rb_cqpminimizer_set(VALUE obj, VALUE data)
-{
-	gsl_cqpminimizer *m;
-	gsl_cqp_data *d;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	Data_Get_Struct(data, gsl_cqp_data, d);
-	gsl_cqpminimizer_set(m, d);
-	return Qtrue;
-}
-
-static VALUE rb_cqpminimizer_name(VALUE obj)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return rb_str_new2(gsl_cqpminimizer_name(m));
-}
-
-static VALUE rb_cqpminimizer_iterate(VALUE obj)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return INT2FIX(gsl_cqpminimizer_iterate(m));
-}
-
-static VALUE rb_cqpminimizer_x(VALUE obj)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return Data_Wrap_Struct(cgsl_vector_view, 0, NULL, gsl_cqpminimizer_x(m));
-}
-
-static VALUE rb_cqpminimizer_lm_eq(VALUE obj)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return Data_Wrap_Struct(cgsl_vector_view, 0, NULL, gsl_cqpminimizer_lm_eq(m));
-}
-static VALUE rb_cqpminimizer_lm_ineq(VALUE obj)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return Data_Wrap_Struct(cgsl_vector_view, 0, NULL, gsl_cqpminimizer_lm_ineq(m));
-}
-static VALUE rb_cqpminimizer_f(VALUE obj)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return rb_float_new(gsl_cqpminimizer_f(m));
-}
-static VALUE rb_cqpminimizer_gap(VALUE obj)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return rb_float_new(gsl_cqpminimizer_gap(m));
-}
-static VALUE rb_cqpminimizer_residuals_norm(VALUE obj)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return rb_float_new(gsl_cqpminimizer_residuals_norm(m));
-}
-/*
-static VALUE rb_cqpminimizer_minimum(VALUE obj)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return rb_float_new(gsl_cqpminimizer_minimum(m));
-}
-*/
-static VALUE rb_cqpminimizer_test_convergence(VALUE obj, VALUE g, VALUE r)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return INT2FIX(gsl_cqpminimizer_test_convergence(m, NUM2DBL(g), NUM2DBL(r)));
-}
-static VALUE rb_cqpminimizer_test_infeasibility(VALUE obj, VALUE e)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return INT2FIX(gsl_cqp_minimizer_test_infeasibility(m, NUM2DBL(e)));
-}
-
-static VALUE rb_cqp_data_alloc(VALUE klass)
-{
-	gsl_cqp_data *d;
-	d = (gsl_cqp_data*) malloc(sizeof(gsl_cqp_data));
-	return Data_Wrap_Struct(klass, 0, free, d);
-}
-
-static VALUE rb_cqp_data_Q(VALUE obj)
-{
-	gsl_cqp_data *d;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	return Data_Wrap_Struct(cgsl_matrix_view, 0, NULL, 	d->Q);
-}
-
-static VALUE rb_cqp_data_q(VALUE obj)
-{
-	gsl_cqp_data *d;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	return Data_Wrap_Struct(cgsl_vector_view, 0, NULL, 	d->q);
-}
-
-static VALUE rb_cqp_data_A(VALUE obj)
-{
-	gsl_cqp_data *d;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	return Data_Wrap_Struct(cgsl_matrix_view, 0, NULL, 	d->A);
-}
-
-static VALUE rb_cqp_data_b(VALUE obj)
-{
-	gsl_cqp_data *d;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	return Data_Wrap_Struct(cgsl_vector_view, 0, NULL, 	d->b);
-}
-
-static VALUE rb_cqp_data_C(VALUE obj)
-{
-	gsl_cqp_data *d;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	return Data_Wrap_Struct(cgsl_matrix_view, 0, NULL, 	d->C);
-}
-
-static VALUE rb_cqp_data_d(VALUE obj)
-{
-	gsl_cqp_data *d;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	return Data_Wrap_Struct(cgsl_vector_view, 0, NULL, 	d->d);
-}
-
-static VALUE rb_cqp_data_set_Q(VALUE obj, VALUE mm)
-{
-	gsl_cqp_data *d;
-	gsl_matrix *m;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	CHECK_MATRIX(mm);
-	Data_Get_Struct(mm, gsl_matrix, m);	
-	d->Q = m;
-	return Qtrue;
-}
-
-static VALUE rb_cqp_data_set_q(VALUE obj, VALUE vv)
-{
-	gsl_cqp_data *d;
-	gsl_vector *v;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	CHECK_VECTOR(vv);
-	Data_Get_Struct(vv, gsl_vector, v);
-	d->q = v;
-	return Qtrue;
-}
-
-static VALUE rb_cqp_data_set_A(VALUE obj, VALUE mm)
-{
-	gsl_cqp_data *d;
-	gsl_matrix *m;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	CHECK_MATRIX(mm);
-	Data_Get_Struct(mm, gsl_matrix, m);
-	d->A = m;
-	return Qtrue;
-}
-
-static VALUE rb_cqp_data_set_b(VALUE obj, VALUE vv)
-{
-	gsl_cqp_data *d;
-	gsl_vector *v;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	CHECK_VECTOR(vv);
-	Data_Get_Struct(vv, gsl_vector, v);
-	d->b = v;
-	return Qtrue;
-}
-
-static VALUE rb_cqp_data_set_C(VALUE obj, VALUE mm)
-{
-	gsl_cqp_data *d;
-	gsl_matrix *m;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	CHECK_MATRIX(mm);
-	Data_Get_Struct(mm, gsl_matrix, m);
-	d->C = m;
-	return Qtrue;
-}
-
-static VALUE rb_cqp_data_set_d(VALUE obj, VALUE vv)
-{
-	gsl_cqp_data *d;
-	gsl_vector *v;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	CHECK_VECTOR(vv);
-	Data_Get_Struct(vv, gsl_vector, v);
-	d->d = v;
-	return Qtrue;
-}
-
-void Init_cqp(VALUE module)
-{
-	VALUE mCQP, cMinimizer, cData;
-	
-	mCQP = rb_define_module_under(module, "CQP");
-	cMinimizer = rb_define_class_under(mCQP, "Minimizer", cGSL_Object);
-	cData = rb_define_class_under(mCQP, "Data", cGSL_Object);	
-
-	rb_define_singleton_method(cMinimizer, "alloc", rb_cqpminimizer_alloc, 4);
-	
-	rb_define_method(cMinimizer, "set", rb_cqpminimizer_set, 1);
-	rb_define_method(cMinimizer, "name", rb_cqpminimizer_name, 0);
-	rb_define_method(cMinimizer, "iterate", rb_cqpminimizer_iterate, 0);	
-	rb_define_method(cMinimizer, "x", rb_cqpminimizer_x, 0);
-	rb_define_method(cMinimizer, "lm_eq", rb_cqpminimizer_lm_eq, 0);
-	rb_define_method(cMinimizer, "lm_ineq", rb_cqpminimizer_lm_ineq, 0);
-	rb_define_method(cMinimizer, "f", rb_cqpminimizer_f, 0);
-	rb_define_method(cMinimizer, "gap", rb_cqpminimizer_gap, 0);
-	rb_define_method(cMinimizer, "residuals_norm", rb_cqpminimizer_residuals_norm, 0);	
-/*	rb_define_method(cMinimizer, "minimum", rb_cqpminimizer_minimum, 0);	*/
-	rb_define_method(cMinimizer, "test_convergence", rb_cqpminimizer_test_convergence, 2);
-	rb_define_method(cMinimizer, "test_infeasibility", rb_cqpminimizer_test_infeasibility, 1);	
-
-	/*****/	
-	rb_define_singleton_method(cData, "alloc", rb_cqp_data_alloc, 0);	
-	
-	rb_define_method(cData, "Q", rb_cqp_data_Q, 0);
-	rb_define_method(cData, "q", rb_cqp_data_q, 0);
-	rb_define_method(cData, "A", rb_cqp_data_A, 0);
-	rb_define_method(cData, "b", rb_cqp_data_b, 0);
-	rb_define_method(cData, "C", rb_cqp_data_C, 0);
-	rb_define_method(cData, "d", rb_cqp_data_d, 0);
-	
-	rb_define_method(cData, "set_Q", rb_cqp_data_set_Q, 1);
-	rb_define_method(cData, "set_q", rb_cqp_data_set_q, 1);
-	rb_define_method(cData, "set_A", rb_cqp_data_set_A, 1);
-	rb_define_method(cData, "set_b", rb_cqp_data_set_b, 1);
-	rb_define_method(cData, "set_C", rb_cqp_data_set_C, 1);
-	rb_define_method(cData, "set_d", rb_cqp_data_set_d, 1);
-	rb_define_alias(cData, "Q=", "set_Q");
-	rb_define_alias(cData, "q=", "set_q");
-	rb_define_alias(cData, "A=", "set_A");
-	rb_define_alias(cData, "b=", "set_b");
-	rb_define_alias(cData, "C=", "set_C");
-	rb_define_alias(cData, "d=", "set_d");	
-}
-
-#endif
-

data/ext/extconf.rb

@@ -1,295 +0,0 @@
-require 'mkmf'
-
-
-module GSL
-  class Version
-    def initialize(str)
-      @str = str
-      @ary = str.split(".").collect { |elm| elm.to_i }
-    end
-    def to_s; @str; end
-    def inspect; @str; end
-    def >=(ver)
-      ary2 = ver.split(".").collect { |elm| elm.to_i }
-      if @ary[0] > ary2[0]; return true; end			
-      if @ary[0] < ary2[0]; return false; end
-      if @ary[1] > ary2[1]; return true; end
-      if @ary[1] < ary2[1]; return false; end
-      if @ary.size < ary2.size; return false; end
-      if @ary.size == 3 and ary2.size == 3
-        if @ary[2] < ary2[2]; return false; end
-      end		
-      return true
-    end
-    def <(ver)
-      ary2 = ver.split(".").collect { |elm| elm.to_i }
-      if @ary[0] >= ary2[0]; return false; end
-      if @ary[0] >= ary2[0]; return false; end
-      return true
-    end
-  end
-end
-
-if /mingw/ =~ RUBY_PLATFORM
-  GSL_CONFIG = "sh gsl-config"
-else
-  GSL_CONFIG = "gsl-config"
-end
-
-def gsl_config()
-  print("checking gsl cflags... ")
-  IO.popen("#{GSL_CONFIG} --cflags") do |f|
-    cflags = f.gets.chomp
-    puts(cflags)
-    $CFLAGS += " " + cflags
-  end
-  
-  IO.popen("#{GSL_CONFIG} --libs") do |f|
-    libs = f.gets.chomp
-    dir_config("cblas")
-    dir_config("atlas")
-    if have_library("cblas") and have_library("atlas")
-      libs.gsub!("-lgslcblas", "-lcblas -latlas")
-      $LOCAL_LIBS += " " + libs.gsub(" -lgslcblas", "")
-      print("checking gsl libs... ")
-      puts(libs)
-    else
-      print("checking gsl libs... ")
-      puts(libs)
-      $LOCAL_LIBS += " " + libs
-    end
-  end
-
-end
-
-def check_version(configfile)
-  
-  print("checking gsl version... ")
-  IO.popen("#{GSL_CONFIG} --version") do |f|
-    ver = GSL::Version.new(f.gets.chomp)
-    puts(ver)
-    configfile.printf("#ifndef GSL_VERSION\n#define GSL_VERSION \"#{ver}\"\n#endif\n")
-
-    if ver >= "0.9.4"
-      configfile.printf("#ifndef GSL_0_9_4_LATER\n#define GSL_0_9_4_LATER\n#endif\n")
-    else
-      configfile.close
-      raise("Ruby/GSL requires gsl-0.9.4 or later.")
-    end
-    if ver >= "1.0"
-      configfile.printf("#ifndef GSL_1_0_LATER\n#define GSL_1_0_LATER\n#endif\n")
-    end
-    if ver >= "1.1"
-      configfile.printf("#ifndef GSL_1_1_LATER\n#define GSL_1_1_LATER\n#endif\n")
-    end
-    if ver >= "1.1.1"
-      configfile.printf("#ifndef GSL_1_1_1_LATER\n#define GSL_1_1_1_LATER\n#endif\n")
-    end
-    if ver >= "1.2"
-      configfile.printf("#ifndef GSL_1_2_LATER\n#define GSL_1_2_LATER\n#endif\n")
-    end
-    if ver >= "1.3"
-      configfile.printf("#ifndef GSL_1_3_LATER\n#define GSL_1_3_LATER\n#endif\n")
-    end
-    if ver >= "1.4"
-      configfile.printf("#ifndef GSL_1_4_LATER\n#define GSL_1_4_LATER\n#endif\n")
-    end
-    if ver >= "1.4.90"
-      configfile.printf("#ifndef GSL_1_4_9_LATER\n#define GSL_1_4_9_LATER\n#endif\n")
-    end
-    
-    if ver >= "1.5.90"
-      configfile.printf("#ifndef GSL_1_6_LATER\n#define GSL_1_6_LATER\n#endif\n")
-    end
-
-    if ver >= "1.7.90"
-      configfile.printf("#ifndef GSL_1_8_LATER\n#define GSL_1_8_LATER\n#endif\n")
-    end
-    if ver >= "1.8.90"
-      configfile.printf("#ifndef GSL_1_9_LATER\n#define GSL_1_9_LATER\n#endif\n")
-    end
-
-    if ver >= "1.9.90"
-      configfile.printf("#ifndef GSL_1_10_LATER\n#define GSL_1_10_LATER\n#endif\n")
-    end    
-    if ver < "1.4"
-      configfile.printf("#ifndef GSL_CONST_OLD\n#define GSL_CONST_OLD\n#endif\n")
-    end
-
-    if ver >= "1.11"
-      configfile.printf("#ifndef GSL_1_11_LATER\n#define GSL_1_11_LATER\n#endif\n")
-    end    
-
-   if ver >= "1.12.90"
-      configfile.printf("#ifndef GSL_1_13_LATER\n#define GSL_1_13_LATER\n#endif\n")
-    end
-
-    if ver >= "1.14"
-      configfile.printf("#ifndef GSL_1_14_LATER\n#define GSL_1_14_LATER\n#endif\n")
-    end
-
-    if ver >= "1.15"
-      configfile.printf("#ifndef GSL_1_15_LATER\n#define GSL_1_15_LATER\n#endif\n")
-    end
-
-  end
-end
-
-#####
-
-$CFLAGS ||= ''
-$CFLAGS += " -Wall -I../include "
-
-begin
-  RB_GSL_CONFIG = File.open("../include/rb_gsl_config.h", "w")
-  RB_GSL_CONFIG.printf("#ifndef ___RB_GSL_CONFIG_H___\n")
-  RB_GSL_CONFIG.printf("#define ___RB_GSL_CONFIG_H___\n\n")
-
-  check_version(RB_GSL_CONFIG)
-
-  gsl_config()
-
-  have_func("round")
-
-# Check GSL extensions
-
-  if have_header("rngextra/rngextra.h")
-    have_library("rngextra")
-  end
-
-  if have_header("qrngextra/qrngextra.h")
-    have_library("qrngextra")
-  end
-
-  if have_header("ool/ool_version.h")
-    have_library("ool")
-  end
-  
-  if have_header("tensor/tensor.h")
-    have_library("tensor")
-  end
-
-  if have_header("jacobi.h")
-    have_library("jacobi")
-  end
-  if have_header("gsl/gsl_cqp.h")
-    have_library("cqp")
-  end
-  if have_header("gsl/gsl_multimin_fsdf.h")
-    have_library("bundle_method")
-  end
-     
-  if have_library("gsl", "gsl_poly_solve_quartic")
-    RB_GSL_CONFIG.printf("#ifndef HAVE_POLY_SOLVE_QUARTIC\n#define HAVE_POLY_SOLVE_QUARTIC\n#endif\n")
-  end
-  if have_library("gsl", "gsl_eigen_francis")
-    RB_GSL_CONFIG.printf("#ifndef HAVE_EIGEN_FRANCIS\n#define HAVE_EIGEN_FRANCIS\n#endif\n")
-  end
-
-  if have_header("ndlinear/gsl_multifit_ndlinear.h")
-    have_library("ndlinear")
-  end
-
-# Added 2009/Apr/20
-  if have_header("alf/alf.h")
-    have_library("alf")
-  end
-  
-  begin
-    print("checking rb-gsl version...")
-    IO.popen("cat ../VERSION") do |f|
-      ver = GSL::Version.new(f.gets.chomp)
-      puts(ver)
-      RB_GSL_CONFIG.printf("#ifndef RUBY_GSL_VERSION\n#define RUBY_GSL_VERSION \"#{ver}\"\n#endif\n")
-    end
-  end
-
-  RUBY_VERSION2 = GSL::Version.new(RUBY_VERSION)
-  
-  puts("checking ruby version... #{RUBY_VERSION2}")
-  if RUBY_VERSION2 >= "1.8"
-    RB_GSL_CONFIG.printf("#ifndef RUBY_1_8_LATER\n#define RUBY_1_8_LATER\n#endif\n")
-
-    if find_executable("graph")
-      RB_GSL_CONFIG.printf("#ifndef HAVE_GNU_GRAPH\n#define HAVE_GNU_GRAPH\n#endif\n")
-    end
-  else
-    path = (path || ENV['PATH']).split(File::PATH_SEPARATOR)  
-    flag = 0
-    print("checking for GNU graph... ")
-    path.each do |dir|    
-      if File.executable?(file = File.join(dir, "graph")) 
-        puts("yes")
-        RB_GSL_CONFIG.printf("#ifndef HAVE_GNU_GRAPH\n#define HAVE_GNU_GRAPH\n#endif\n")
-        flag = 1
-        break
-      end
-    end
-    puts("no") if flag == 0
-  end
-  if RUBY_VERSION2 >= "1.9"
-    RB_GSL_CONFIG.printf("#ifndef RUBY_1_9_LATER\n#define RUBY_1_9_LATER\n#endif\n")
-# Added 2010/Sep/29
-    if RUBY_VERSION2 >= "1.9.2"
-      RB_GSL_CONFIG.printf("#ifndef RUBY_1_9_2_LATER\n#define RUBY_1_9_2_LATER\n#endif\n")
-    end
-  end
-
-  RB_GSL_CONFIG.printf("\n#endif\n")
-  RB_GSL_CONFIG.close
-  
-rescue
-  raise("Check GSL>=0.9.4 is installed, and the command \"gsl-config\" is in search path.")
-end
-
-#narray_config = dir_config("narray")
-narray_config = dir_config('narray',$sitearchdir,$sitearchdir)
-# Try to find narray with RubyGems
-begin
-  require 'rubygems'
-  na_gemspec=Gem::Specification.find_by_path('narray.h')
-  if na_gemspec
-    narray_config = File.join(na_gemspec.full_gem_path, na_gemspec.require_path)
-    $CPPFLAGS = " -I#{narray_config} "+$CPPFLAGS
-  end
-rescue LoadError
-end
-have_narray_h = have_header("narray.h")
-if narray_config
-  if RUBY_PLATFORM =~ /cygwin|mingw/
-#    have_library("narray") || raise("ERROR: narray import library is not found") 
-  have_library("narray")
-  end
-end
-
-tamu_anova_config = dir_config('tamu_anova',$sitearchdir,$sitearchdir)
-have_tamu_anova_h = have_header("tamu_anova/tamu_anova.h")
-if tamu_anova_config
-  have_library("tamuanova")
-#  if RUBY_PLATFORM =~ /cygwin|mingw/
-#    have_library("tamuanova") || raise("ERROR: tamu_anova import library is not found")
-#  end
-end
-
-File.open("../lib/gsl.rb", "w") do |file|
-  if have_narray_h
-    file.print("require('narray')\n")
-  end
-#  file.print("require('rb_gsl')\ninclude GSL\n")
-  file.print("require('rb_gsl')\n")  
-  file.print("require('gsl/oper.rb')\n")
-end
-
-File.open("../lib/rbgsl.rb", "w") do |file|
-  if have_narray_h
-    file.print("require('narray')\n")
-  end
-  file.print("require('rb_gsl')\n")
-  file.print("require('gsl/oper.rb')\n")
-end
-
-srcs = Dir.glob("*.c") - ["vector_source.c", "matrix_source.c", "tensor_source.c", "poly_source.c", "block_source.c"]
-
-$objs = srcs.collect { |f| f.sub(".c", ".o") }
-
-create_makefile("rb_gsl")