gsl4r 0.0.1

data/lib/gsl4r/f

@@ -0,0 +1,174 @@
+[Function] 
+[Function] 
+gsl_complex gsl_complex_pow (gsl complex z, gsl complex a ) 
+  The function returns the complex number z raised to the complex power a, za . This 
+  is computed as exp(log(z) 
+      ∗ a) using complex logarithms and complex exponentials. 
+  [Function] 
+  gsl_complex gsl_complex_pow_real (gsl complex z, double x ) 
+  This function returns the complex number z raised to the real power x, zx . 
+  [Function] 
+  gsl_complex gsl_complex_exp (gsl complex z ) 
+  This function returns the complex exponential of the complex number z, exp(z).
+  Chapter 5: Complex Numbers 24 
+  [Function] 
+  gsl_complex gsl_complex_log (gsl complex z ) 
+  This function returns the complex natural logarithm (base e) of the complex number 
+  z, log(z). The branch cut is the negative real axis. 
+  [Function] 
+  gsl_complex gsl_complex_log10 (gsl complex z ) 
+  This function returns the complex base-10 logarithm of the complex number z, 
+  log10 (z). 
+  [Function] 
+  gsl_complex gsl_complex_log_b (gsl complex z, gsl complex b ) 
+  This function returns the complex base-b logarithm of the complex number z, logb (z). 
+  This quantity is computed as the ratio log(z)/ log(b). 
+  5.5 Complex Trigonometric Functions 
+  [Function] 
+  gsl_complex gsl_complex_sin (gsl complex z ) 
+  This function returns the complex sine of the complex number z, sin(z) = (exp(iz) 
+      − 
+      exp( 
+	−iz))/(2i). 
+  [Function] 
+  gsl_complex gsl_complex_cos (gsl complex z ) 
+  This function returns the complex cosine of the complex number z, cos(z) = (exp(iz)+ 
+      exp( 
+	−iz))/2. 
+  [Function] 
+  gsl_complex gsl_complex_tan (gsl complex z ) 
+  This function returns the complex tangent of the complex number z, 
+  tan(z) = sin(z)/ cos(z). 
+  [Function] 
+  gsl_complex gsl_complex_sec (gsl complex z ) 
+  This function returns the complex secant of the complex number z, sec(z) = 1/ cos(z). 
+  [Function] 
+  gsl_complex gsl_complex_csc (gsl complex z ) 
+  This function returns the complex cosecant of the complex number z, csc(z) = 
+  1/ sin(z). 
+  [Function] 
+  gsl_complex gsl_complex_cot (gsl complex z ) 
+  This function returns the complex cotangent of the complex number z, cot(z) = 
+  1/ tan(z). 
+  5.6 Inverse Complex Trigonometric Functions 
+  [Function] 
+  gsl_complex gsl_complex_arcsin (gsl complex z ) 
+  This function returns the complex arcsine of the complex number z, arcsin(z). The 
+  branch cuts are on the real axis, less than 
+  −1 and greater than 1. 
+  [Function] 
+  gsl_complex gsl_complex_arcsin_real (double z ) 
+  This function returns the complex arcsine of the real number z, arcsin(z). For z 
+  between 
+  −1 and 1, the function returns a real value in the range [−π/2, π/2]. For z 
+  less than 
+  −1 the result has a real part of −π/2 and a positive imaginary part. For z 
+  greater than 1 the result has a real part of π/2 and a negative imaginary part.
+  Chapter 5: Complex Numbers 25 
+  [Function] 
+  gsl_complex gsl_complex_arccos (gsl complex z ) 
+  This function returns the complex arccosine of the complex number z, arccos(z). The 
+  branch cuts are on the real axis, less than 
+  −1 and greater than 1. 
+  [Function] 
+  gsl_complex gsl_complex_arccos_real (double z ) 
+  This function returns the complex arccosine of the real number z, arccos(z). For z 
+  between 
+  −1 and 1, the function returns a real value in the range [0, π]. For z less 
+  than 
+  −1 the result has a real part of π and a negative imaginary part. For z greater 
+  than 1 the result is purely imaginary and positive. 
+  [Function] 
+  gsl_complex gsl_complex_arctan (gsl complex z ) 
+  This function returns the complex arctangent of the complex number z, arctan(z). 
+  The branch cuts are on the imaginary axis, below 
+  −i and above i. 
+  [Function] 
+  gsl_complex gsl_complex_arcsec (gsl complex z ) 
+  This function returns the complex arcsecant of the complex number z, arcsec(z) = 
+  arccos(1/z). 
+  [Function] 
+  gsl_complex gsl_complex_arcsec_real (double z ) 
+  This function returns the complex arcsecant of the real number z, 
+  arcsec(z) = arccos(1/z). 
+  [Function] 
+  gsl_complex gsl_complex_arccsc (gsl complex z ) 
+  This function returns the complex arccosecant of the complex number z, arccsc(z) = 
+  arcsin(1/z). 
+  [Function] 
+  gsl_complex gsl_complex_arccsc_real (double z ) 
+  This function returns the complex arccosecant of the real number z, arccsc(z) = 
+  arcsin(1/z). 
+  [Function] 
+  gsl_complex gsl_complex_arccot (gsl complex z ) 
+  This function returns the complex arccotangent of the complex number z, arccot(z) = 
+  arctan(1/z). 
+  5.7 Complex Hyperbolic Functions 
+  [Function] 
+  gsl_complex gsl_complex_sinh (gsl complex z ) 
+  This function returns the complex hyperbolic sine of the complex number z, sinh(z) = 
+  (exp(z) 
+   − exp(−z))/2. 
+  [Function] 
+  gsl_complex gsl_complex_cosh (gsl complex z ) 
+  This function returns the complex hyperbolic cosine of the complex number z, 
+  cosh(z) = (exp(z) + exp( 
+	−z))/2. 
+  [Function] 
+  gsl_complex gsl_complex_tanh (gsl complex z ) 
+  This function returns the complex hyperbolic tangent of the complex number z, 
+  tanh(z) = sinh(z)/ cosh(z). 
+  [Function] 
+  gsl_complex gsl_complex_sech (gsl complex z ) 
+  This function returns the complex hyperbolic secant of the complex number z, 
+  sech(z) = 1/ cosh(z).
+  Chapter 5: Complex Numbers 26 
+  [Function] 
+  gsl_complex gsl_complex_csch (gsl complex z ) 
+  This function returns the complex hyperbolic cosecant of the complex number z, 
+  csch(z) = 1/ sinh(z). 
+  [Function] 
+  gsl_complex gsl_complex_coth (gsl complex z ) 
+  This function returns the complex hyperbolic cotangent of the complex number z, 
+  coth(z) = 1/ tanh(z). 
+  5.8 Inverse Complex Hyperbolic Functions 
+  [Function] 
+  gsl_complex gsl_complex_arcsinh (gsl complex z ) 
+  This function returns the complex hyperbolic arcsine of the complex number z, 
+  arcsinh(z). The branch cuts are on the imaginary axis, below 
+  −i and above i. 
+  [Function] 
+  gsl_complex gsl_complex_arccosh (gsl complex z ) 
+  This function returns the complex hyperbolic arccosine of the complex number z, 
+  arccosh(z). The branch cut is on the real axis, less than 1. Note that in this case 
+  we use the negative square root in formula 4.6.21 of Abramowitz & Stegun giving 
+  arccosh(z) = log(z 
+      − √z 
+      2 
+      − 1). 
+  [Function] 
+  gsl_complex gsl_complex_arccosh_real (double z ) 
+  This function returns the complex hyperbolic arccosine of the real number z, 
+  arccosh(z). 
+  [Function] 
+  gsl_complex gsl_complex_arctanh (gsl complex z ) 
+  This function returns the complex hyperbolic arctangent of the complex number z, 
+  arctanh(z). The branch cuts are on the real axis, less than 
+  −1 and greater than 1. 
+  [Function] 
+  gsl_complex gsl_complex_arctanh_real (double z ) 
+  This function returns the complex hyperbolic arctangent of the real number z, 
+  arctanh(z). 
+  [Function] 
+  gsl_complex gsl_complex_arcsech (gsl complex z ) 
+  This function returns the complex hyperbolic arcsecant of the complex number z, 
+  arcsech(z) = arccosh(1/z). 
+  [Function] 
+  gsl_complex gsl_complex_arccsch (gsl complex z ) 
+  This function returns the complex hyperbolic arccosecant of the complex number z, 
+  arccsch(z) = arcsin(1/z). 
+  [Function] 
+  gsl_complex gsl_complex_arccoth (gsl complex z ) 
+  This function returns the complex hyperbolic arccotangent of the complex number z, 
+  arccoth(z) = arctanh(1/z). 
+

data/lib/gsl4r/harness.rb

@@ -0,0 +1,58 @@
+
+#
+# == Other Info
+#
+# Author::      Colby Gutierrez-Kraybill
+# Version::     $Id$
+#
+
+require 'rubygems'
+
+module GSL4r
+  module Harness
+    attr_accessor :c_includes, :c_libs, :c_tests, :c_src_name, :c_binary
+    attr_accessor :r_header, :r_footer
+
+    TEST_DIR = "test"
+
+    def write_c_tests
+      f = File.new("#{TEST_DIR}/#{@c_src_name}", "w")
+      f.puts "/* Auto generated by #{self.class.name} */"
+      f.puts "#include <stdio.h>"
+      @c_includes.each { |i|
+	f.puts "#include \"#{i}\""
+      }
+      f.puts "int main( int argc, char **argv )\n{\n"
+      
+      f.puts "  puts(\"#{@r_header}\");"
+
+      @c_tests.each { |t|
+	t_fqmn_a = self.class.name.split("::")
+	t_fqmn_a.pop
+	src = ""
+	eval <<-end_eval
+	   src = ::#{t_fqmn_a.join("::")}::Methods::#{t}
+	end_eval
+	f.puts " /* #{t} */"
+	f.puts src
+      }
+
+      f.puts "  puts(\"#{@r_footer}\");"
+
+      f.puts "  return(0);\n}\n"
+      f.close
+    end
+
+    def compile_c_tests
+      compile_s = "#{@c_compiler} #{@c_flags.join(" ")} " +
+	"-o #{TEST_DIR}/#{@c_binary} #{TEST_DIR}/#{@c_src_name}"
+      p compile_s
+      `#{compile_s}`
+    end
+
+    def run_c_tests( filename )
+      `#{TEST_DIR}/#{@c_binary} > test/#{filename}`
+    end
+
+  end
+end

data/lib/gsl4r/platform.rb

@@ -0,0 +1,8 @@
+#
+# This file is automatically modified at the
+# time of gem installation
+#
+
+module GSL4r
+  GSL_SHARED_LIB_DIR = "GSL_CONFIG_LOCATION"
+end

data/lib/gsl4r/util.rb

@@ -0,0 +1,85 @@
+
+#
+# == Other Info
+#
+# Author::      Colby Gutierrez-Kraybill
+# Version::     $Id$
+#
+
+require 'rubygems'
+
+module GSL4r
+  module Util
+
+    $c_var_num = 0
+
+    def attach_gsl_function( method_name, args, return_var, args_type=nil, return_type=nil )
+
+      # This function is attached to the extended ::FFI::Library
+      # module from the calling namespace, e.g. ::GSL4r::Complex::Methods
+      attach_function method_name, args, return_var
+
+      if ( args_type != nil )
+        # prepare c and ruby args code
+	c_src = ""
+	c_call_vars = []
+	c_return_name = "c_r#{$c_var_num}"
+	r_src = []
+	if ( ! args_type.is_a?(Array) )
+	  args_type = Array.new([args_type])
+	end
+        args_type.each { |a_t|
+	  c_var_name = "v#{$c_var_num += 1}"
+	  c_src << (a_t.respond_to?("c_type") ?
+		    "  #{a_t.c_type} #{c_var_name};\n" : "#{a_t.to_s} #{c_var_name} ")
+	  c_src << (a_t.respond_to?("c_assignment") ?
+		    "  #{a_t.c_assignment("#{c_var_name}")}\n" : "= (#{a_t.to_s})2.0;\n")
+	  c_call_vars << "#{c_var_name}"
+
+	  r_src << (a_t.respond_to?("r_type") ?
+		    "  #{c_var_name} = #{a_t.r_type}.create" : "")
+	  r_src << (a_t.respond_to?("r_assignment") ?
+		    "  #{a_t.r_assignment("#{c_var_name}")}" : "  #{c_var_name} = 2.0")
+	} # args_type.each
+
+	# prepare c return type
+	c_src << (return_type.respond_to?("c_type") ?
+		  "  #{return_type.c_type} #{c_return_name};\n" :
+		  "  #{return_type.to_s} #{c_return_name};\n")
+
+	# prepare c call
+	c_src << "  #{c_return_name} = #{method_name}(#{c_call_vars.join(",")});\n"
+	
+	# now generate the ruby code for the unit test
+	c_src << "  puts(" << %Q{\\"def test_#{method_name}()\\"} << ");\n"
+
+	# TODO, Need to insert ruby object instantiation code here!
+	#
+	r_src.each { |v|
+	  c_src << "  puts(" << %Q{\\"#{v}\\"} << ");\n"
+	}
+
+	r_r1 = "r_r1" # ruby result
+	c_src << "  puts(" << %Q{\\"  #{r_r1} = ::#{self.to_s}::#{method_name}(#{c_call_vars.join(",")})\\"} << ");\n"
+	if ( return_type.respond_to?("c_to_r_assignment") )
+	  r_r2 = "r_r2" # ruby result comparitor
+	  c_src << "  puts(" << %Q{\\"  #{r_r2} = #{return_type.r_type}.new\\"} << ");\n"
+	  c_src << "  #{return_type.c_to_r_assignment(r_r2,c_return_name)}"
+	  c_src << "  printf(" << %Q{\\"  assert r_r1.equals(r_r2)\\\\n\\"} << ");\n"
+	else
+	  c_src << "  printf(" << %Q{\\"  assert_in_delta r_r1, %.15g, EPSILON\\\\n\\"} << ", #{c_return_name});\n"
+	end
+
+	c_src << "  puts(" << %Q{\\"end\\"} << ");"
+
+
+	eval <<-end_eval
+        def c_test_#{method_name}
+  	  # Build list of arguments and their values
+	  "#{c_src}"
+        end
+	end_eval
+      end
+    end
+  end
+end

data/lib/gsl4r/vector.rb

@@ -0,0 +1,29 @@
+
+#
+# == Other Info
+#
+# Author::	Colby Gutierrez-Kraybill
+# Version::	$Id$
+#
+
+require 'rubygems'
+require 'ffi'
+
+module GSL4r
+
+  extend ::FFI::Library
+
+  ffi_lib 'gsl'
+
+  class Vector
+    attach_function :gsl_vector_alloc, [ :size_t ], :pointer
+    attach_function :gsl_vector_calloc, [ :size_t ], :pointer
+    attach_function :gsl_vector_free, [ :pointer ], :void
+
+    # initializing
+    attach_function :gsl_vector_set_all, [ :pointer, :double ], :void
+    attach_function :gsl_vector_set_zero, [ :pointer ], :void
+    attach_function :gsl_vector_set_basis, [ :pointer, :size_t ], :int
+  end
+
+end

data/test/complex_test.rb

@@ -0,0 +1,490 @@
+$: << File.join('..','lib')
+require 'test/unit'
+require 'test/unit/autorunner'
+require 'gsl4r/complex'
+include GSL4r::Complex
+class ComplexTests < Test::Unit::TestCase
+  EPSILON = 5.0e-15
+def test_gsl_complex_tanh()
+  v64 = GSL_Complex.create
+  v64.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_tanh(v64)
+  r_r2 = GSL_Complex.new
+  r_r2.set(1.02383559457047,-0.0283929528682323)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_pow_real()
+  v37 = GSL_Complex.create
+  v37.set(2.0,2.0)
+
+  v38 = 2.0
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_pow_real(v37,v38)
+  r_r2 = GSL_Complex.new
+  r_r2.set(4.89858719658941e-16,8)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_sinh()
+  v60 = GSL_Complex.create
+  v60.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_sinh(v60)
+  r_r2 = GSL_Complex.new
+  r_r2.set(-1.50930648532362,3.42095486111701)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_sqrt()
+  v31 = GSL_Complex.create
+  v31.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_sqrt(v31)
+  r_r2 = GSL_Complex.new
+  r_r2.set(1.55377397403004,0.643594252905582)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_arccsc()
+  v56 = GSL_Complex.create
+  v56.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arccsc(v56)
+  r_r2 = GSL_Complex.new
+  r_r2.set(0.24452216513554,-0.254895573340551)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_div_imag()
+  v26 = GSL_Complex.create
+  v26.set(2.0,2.0)
+
+  v27 = 2.0
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_div_imag(v26,v27)
+  r_r2 = GSL_Complex.new
+  r_r2.set(1,-1)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_abs2()
+  v3 = GSL_Complex.create
+  v3.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_abs2(v3)
+  assert_in_delta r_r1, 8, EPSILON
+end
+def test_gsl_complex_arccos()
+  v52 = GSL_Complex.create
+  v52.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arccos(v52)
+  r_r2 = GSL_Complex.new
+  r_r2.set(0.816547182096851,-1.73432452148797)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_arccoth()
+  v73 = GSL_Complex.create
+  v73.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arccoth(v73)
+  r_r2 = GSL_Complex.new
+  r_r2.set(0.238877861256859,-0.259573057123261)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_div_real()
+  v18 = GSL_Complex.create
+  v18.set(2.0,2.0)
+
+  v19 = 2.0
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_div_real(v18,v19)
+  r_r2 = GSL_Complex.new
+  r_r2.set(1,1)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_tan()
+  v48 = GSL_Complex.create
+  v48.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_tan(v48)
+  r_r2 = GSL_Complex.new
+  r_r2.set(-0.0283929528682323,1.02383559457047)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_arcsech()
+  v69 = GSL_Complex.create
+  v69.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arcsech(v69)
+  r_r2 = GSL_Complex.new
+  r_r2.set(0.254895573340551,-1.32627416165936)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_div()
+  v10 = GSL_Complex.create
+  v10.set(2.0,2.0)
+  v11 = GSL_Complex.create
+  v11.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_div(v10,v11)
+  r_r2 = GSL_Complex.new
+  r_r2.set(1,0)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_sin()
+  v44 = GSL_Complex.create
+  v44.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_sin(v44)
+  r_r2 = GSL_Complex.new
+  r_r2.set(3.42095486111701,-1.50930648532362)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_coth()
+  v65 = GSL_Complex.create
+  v65.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_coth(v65)
+  r_r2 = GSL_Complex.new
+  r_r2.set(0.975968735117309,0.0270655117325545)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_exp()
+  v39 = GSL_Complex.create
+  v39.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_exp(v39)
+  r_r2 = GSL_Complex.new
+  r_r2.set(-3.07493232063936,6.71884969742825)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_cosh()
+  v61 = GSL_Complex.create
+  v61.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_cosh(v61)
+  r_r2 = GSL_Complex.new
+  r_r2.set(-1.56562583531574,3.29789483631124)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_sqrt_real()
+
+  v32 = 2.0
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_sqrt_real(v32)
+  r_r2 = GSL_Complex.new
+  r_r2.set(1.4142135623731,0)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_arccsc_real()
+
+  v57 = 2.0
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arccsc_real(v57)
+  r_r2 = GSL_Complex.new
+  r_r2.set(0.523598775598299,0)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_conjugate()
+  v28 = GSL_Complex.create
+  v28.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_conjugate(v28)
+  r_r2 = GSL_Complex.new
+  r_r2.set(2,-2)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_arccos_real()
+
+  v53 = 2.0
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arccos_real(v53)
+  r_r2 = GSL_Complex.new
+  r_r2.set(0,1.31695789692482)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_add_imag()
+  v20 = GSL_Complex.create
+  v20.set(2.0,2.0)
+
+  v21 = 2.0
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_add_imag(v20,v21)
+  r_r2 = GSL_Complex.new
+  r_r2.set(2,4)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_cot()
+  v49 = GSL_Complex.create
+  v49.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_cot(v49)
+  r_r2 = GSL_Complex.new
+  r_r2.set(-0.0270655117325545,-0.975968735117309)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_arccsch()
+  v70 = GSL_Complex.create
+  v70.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arccsch(v70)
+  r_r2 = GSL_Complex.new
+  r_r2.set(0.254895573340551,-0.24452216513554)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_add_real()
+  v12 = GSL_Complex.create
+  v12.set(2.0,2.0)
+
+  v13 = 2.0
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_add_real(v12,v13)
+  r_r2 = GSL_Complex.new
+  r_r2.set(4,2)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_cos()
+  v45 = GSL_Complex.create
+  v45.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_cos(v45)
+  r_r2 = GSL_Complex.new
+  r_r2.set(-1.56562583531574,-3.29789483631124)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_arcsinh()
+  v66 = GSL_Complex.create
+  v66.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arcsinh(v66)
+  r_r2 = GSL_Complex.new
+  r_r2.set(1.73432452148797,0.754249144698046)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_add()
+  v4 = GSL_Complex.create
+  v4.set(2.0,2.0)
+  v5 = GSL_Complex.create
+  v5.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_add(v4,v5)
+  r_r2 = GSL_Complex.new
+  r_r2.set(4,4)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_log()
+  v40 = GSL_Complex.create
+  v40.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_log(v40)
+  r_r2 = GSL_Complex.new
+  r_r2.set(1.03972077083992,0.785398163397448)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_sech()
+  v62 = GSL_Complex.create
+  v62.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_sech(v62)
+  r_r2 = GSL_Complex.new
+  r_r2.set(-0.117475142661415,-0.24745418582076)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_polar()
+
+  v33 = 2.0
+
+  v34 = 2.0
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_polar(v33,v34)
+  r_r2 = GSL_Complex.new
+  r_r2.set(-0.832293673094285,1.81859485365136)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_arctan()
+  v58 = GSL_Complex.create
+  v58.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arctan(v58)
+  r_r2 = GSL_Complex.new
+  r_r2.set(1.31122326967164,0.238877861256859)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_inverse()
+  v29 = GSL_Complex.create
+  v29.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_inverse(v29)
+  r_r2 = GSL_Complex.new
+  r_r2.set(0.25,-0.25)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_arcsec()
+  v54 = GSL_Complex.create
+  v54.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arcsec(v54)
+  r_r2 = GSL_Complex.new
+  r_r2.set(1.32627416165936,0.254895573340551)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_sub_imag()
+  v22 = GSL_Complex.create
+  v22.set(2.0,2.0)
+
+  v23 = 2.0
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_sub_imag(v22,v23)
+  r_r2 = GSL_Complex.new
+  r_r2.set(2,0)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_arcsin()
+  v50 = GSL_Complex.create
+  v50.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arcsin(v50)
+  r_r2 = GSL_Complex.new
+  r_r2.set(0.754249144698046,1.73432452148797)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_arg()
+  v1 = GSL_Complex.create
+  v1.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arg(v1)
+  assert_in_delta r_r1, 0.785398163397448, EPSILON
+end
+def test_gsl_complex_arctanh()
+  v71 = GSL_Complex.create
+  v71.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arctanh(v71)
+  r_r2 = GSL_Complex.new
+  r_r2.set(0.238877861256859,1.31122326967164)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_sub_real()
+  v14 = GSL_Complex.create
+  v14.set(2.0,2.0)
+
+  v15 = 2.0
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_sub_real(v14,v15)
+  r_r2 = GSL_Complex.new
+  r_r2.set(0,2)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_sec()
+  v46 = GSL_Complex.create
+  v46.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_sec(v46)
+  r_r2 = GSL_Complex.new
+  r_r2.set(-0.117475142661415,0.24745418582076)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_arccosh()
+  v67 = GSL_Complex.create
+  v67.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arccosh(v67)
+  r_r2 = GSL_Complex.new
+  r_r2.set(1.73432452148797,0.816547182096851)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_sub()
+  v6 = GSL_Complex.create
+  v6.set(2.0,2.0)
+  v7 = GSL_Complex.create
+  v7.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_sub(v6,v7)
+  r_r2 = GSL_Complex.new
+  r_r2.set(0,0)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_log10()
+  v41 = GSL_Complex.create
+  v41.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_log10(v41)
+  r_r2 = GSL_Complex.new
+  r_r2.set(0.451544993495972,0.34109408846046)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_csch()
+  v63 = GSL_Complex.create
+  v63.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_csch(v63)
+  r_r2 = GSL_Complex.new
+  r_r2.set(-0.107954592221385,-0.244687073586957)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_pow()
+  v35 = GSL_Complex.create
+  v35.set(2.0,2.0)
+  v36 = GSL_Complex.create
+  v36.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_pow(v35,v36)
+  r_r2 = GSL_Complex.new
+  r_r2.set(-1.4525046270557,-0.8098895463353)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_arccot()
+  v59 = GSL_Complex.create
+  v59.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arccot(v59)
+  r_r2 = GSL_Complex.new
+  r_r2.set(0.259573057123261,-0.238877861256859)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_negative()
+  v30 = GSL_Complex.create
+  v30.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_negative(v30)
+  r_r2 = GSL_Complex.new
+  r_r2.set(-2,-2)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_arcsec_real()
+
+  v55 = 2.0
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arcsec_real(v55)
+  r_r2 = GSL_Complex.new
+  r_r2.set(1.0471975511966,0)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_mul_imag()
+  v24 = GSL_Complex.create
+  v24.set(2.0,2.0)
+
+  v25 = 2.0
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_mul_imag(v24,v25)
+  r_r2 = GSL_Complex.new
+  r_r2.set(-4,4)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_arcsin_real()
+
+  v51 = 2.0
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arcsin_real(v51)
+  r_r2 = GSL_Complex.new
+  r_r2.set(1.5707963267949,-1.31695789692482)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_abs()
+  v2 = GSL_Complex.create
+  v2.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_abs(v2)
+  assert_in_delta r_r1, 2.82842712474619, EPSILON
+end
+def test_gsl_complex_arctanh_real()
+
+  v72 = 2.0
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arctanh_real(v72)
+  r_r2 = GSL_Complex.new
+  r_r2.set(0.549306144334055,-1.5707963267949)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_mul_real()
+  v16 = GSL_Complex.create
+  v16.set(2.0,2.0)
+
+  v17 = 2.0
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_mul_real(v16,v17)
+  r_r2 = GSL_Complex.new
+  r_r2.set(4,4)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_csc()
+  v47 = GSL_Complex.create
+  v47.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_csc(v47)
+  r_r2 = GSL_Complex.new
+  r_r2.set(0.244687073586957,0.107954592221385)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_arccosh_real()
+
+  v68 = 2.0
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arccosh_real(v68)
+  r_r2 = GSL_Complex.new
+  r_r2.set(1.31695789692482,0)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_mul()
+  v8 = GSL_Complex.create
+  v8.set(2.0,2.0)
+  v9 = GSL_Complex.create
+  v9.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_mul(v8,v9)
+  r_r2 = GSL_Complex.new
+  r_r2.set(0,8)
+  assert r_r1.equals(r_r2)
+end
+def test_gsl_complex_log_b()
+  v42 = GSL_Complex.create
+  v42.set(2.0,2.0)
+  v43 = GSL_Complex.create
+  v43.set(2.0,2.0)
+  r_r1 = ::GSL4r::Complex::Methods::gsl_complex_log_b(v42,v43)
+  r_r2 = GSL_Complex.new
+  r_r2.set(1,0)
+  assert r_r1.equals(r_r2)
+end
+end