gsl4r 0.0.2 → 0.0.3

data/lib/t.rb

@@ -1,10 +1,19 @@
-include ::GSL4r::Complex
-a=GSL_Complex.create(1,1)
-
-a.gsl_complex_arccsc_real(5.0)
-a.arccsc_real(5.0)
-
-#a.gsl_complex_abs(a)
-#a.gsl_complex_abs()
-#a.abs()
-
+require 'rubygems'
+require 'gsl4r'
+require 'gsl4r/vector'
+include FFI
+op1 = MemoryPointer.new :pointer
+op1 = GSL4r::Vector::Methods::gsl_vector_alloc(50)
+puts op1.get_ulong(0)
+ob1=GSL4r::Vector::GSL_Vector.new(op1)
+ob2=GSL4r::Vector::GSL_Vector_Cast.new(op1)
+puts ob1.length
+puts ob2.length
+ob1.set( (1..50).to_a.collect { |i| i=i*5.0 } )
+puts ob2.values
+#op1.put_array_of_double(1, (1..50).to_a.collect { |i| i=i*5.0 } )
+#ob1[:data].put_array_of_double(0, (1..50).to_a.collect { |i| i=i*5.0 } )
+#puts ob2.set( (1..50).to_a.collect { |i| i=i*5.0 } )
+#v = ob2.values
+#v[1] = 3000.0
+#puts v

data/test/README

@@ -1,8 +1,12 @@
-This directory is empty for now, as all tests are generated automatically by running:
+Most tests are generated automatically by running at the top of the source tree:
 
 rake test
 
-This will create a C programs that make calls to GSL functions, then writing the
+This will create C programs that make calls to GSL functions, then writing the
 results as Ruby programs that cross check the FFI based wrapper answers against
 the C programs' results using Test::Unit.
 
+Other tests include:
+
+block_class_test.rb - Test basic operations dealing with gsl_block types
+vector_class_test.rb - Test basic operations dealing with gsl_vector types

data/test/block_class_test.rb

@@ -0,0 +1,104 @@
+$: << File.join('..','lib')
+
+require 'rubygems'
+require 'ffi'
+
+require 'test/unit'
+require 'test/unit/autorunner'
+
+require 'gsl4r'
+require 'gsl4r/block'
+
+include GSL4r::Block
+include FFI
+
+class BlockTests < Test::Unit::TestCase
+
+  SIZE=50
+
+  def test_gsl_block_alloc_and_free()
+    assert_nothing_raised do
+      blkptr = MemoryPointer.new :pointer
+      blkptr = GSL4r::Block::Methods::gsl_block_alloc( SIZE )
+      GSL4r::Block::Methods::gsl_block_free( blkptr )
+    end
+  end
+
+  # freed by leaving scope
+  def test_gsl_block_calloc()
+    assert_nothing_raised do
+      blkptr = MemoryPointer.new :pointer
+      blkptr = GSL4r::Block::Methods::gsl_block_calloc( SIZE )
+      blk = GSL_Block.new( blkptr )
+
+      v = blk.values
+      assert v.length == SIZE
+      v.each { |i| assert i == 0 }
+
+    end
+  end
+
+  def test_gsl_block_new()
+    assert_nothing_raised do
+      blkptr = MemoryPointer.new :pointer
+      blkptr = GSL4r::Block::Methods::gsl_block_alloc( SIZE )
+
+      blk = GSL_Block.new( blkptr )
+      assert blk.length == SIZE
+    end
+  end
+
+  class ::GSL4r::Block::GSL_Block_Monitor < ::GSL4r::Block::GSL_Block
+    def self.release(ptr)
+      super
+      $release_count = $release_count + 1
+    end
+  end
+
+  class ::GSL4r::Block::GSL_Block_Cast_Monitor < ::GSL4r::Block::GSL_Block_Cast
+    def self.release(ptr)
+      super
+      $release_count = $release_count + 1
+    end
+  end
+
+  # tests using a cast object to retain a reference
+  # to a memory pointer, without freeing it when the
+  # object is garbage collected.
+  def test_gsl_block_cast()
+    $release_count = 0
+
+    assert_nothing_raised do
+      blkptr = MemoryPointer.new :pointer
+      blkptr = GSL4r::Block::Methods::gsl_block_alloc( SIZE )
+
+      b1=GSL4r::Block::GSL_Block_Monitor.new( blkptr )
+      b2=GSL4r::Block::GSL_Block_Cast_Monitor.new( blkptr )
+
+      assert b1.length == b2.length
+    end
+
+    b1=nil
+    b2=nil
+    GC.start
+
+    assert $release_count == 1
+
+  end
+
+  def test_gsl_block_set_values()
+    assert_nothing_raised do
+      blkptr = MemoryPointer.new :pointer
+      blkptr = GSL4r::Block::Methods::gsl_block_alloc( SIZE )
+      blk = GSL_Block.new( blkptr )
+
+      blk.set( (1..SIZE).to_a.collect { |i| i=i*i } )
+      v = blk.values
+      (1..SIZE).each { |i| assert v[i-1] == i*i }
+
+    end
+  end
+
+
+
+end

data/test/complex_test.rb

@@ -5,6 +5,216 @@ require 'gsl4r/complex'
 include GSL4r::Complex
 class ComplexTests < Test::Unit::TestCase
   EPSILON = 5.0e-15
+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_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_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_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_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_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_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_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_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_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_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_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_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_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_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
+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_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_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_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_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_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_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_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_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_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_add_real()
   v12 = GSL_Complex.create
   v12.set(2.0,2.0)
@@ -277,214 +487,4 @@ def test_gsl_complex_arccoth()
   r_r2.set(0.238877861256859,-0.259573057123261)
   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_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_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_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_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_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_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_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_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_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_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_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_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_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_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
-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_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_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_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_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_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_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_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_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_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
 end